Doubts about Seagrid

[chaos4u2]

I have a bathymetry file way too big, about 2 million data, and when I use the appliaction seagrid in matlab it almost blows up and goes very very slow.

My question is how would you do to have less data. I would take 1 of each 4 data value for example; but there must me a better way to do it. Isn't it?

And 2nd question: when I have my grid, how do you tell ROMS that in your domain everything is closed except 2 zones. I know there's the lateral boundary conditions, but they are north, south, etc... How do you concretise the coordenates of your domain where currents exists?


Thank you in advanced!!!

[lanerolle]

What we do is to generate a grid via Seagrid with a default/dummy bathymetric value.

Thereafter, we take the grid NetCDF file and write the bathymetry into it via say Matlab as a separate procedure. I have dealt with way more than 2 million points (may be x 10 or x 50 more than that) and loop over the wet-points of the ROMS grid and interpolate the bathymetry via an inverse distance squared algorithm which only looks at the soundings within a given radius around the grid point in question - you can do this Matlab via the inpolygon function.

Hope this helps.

[chaos4u2]

Thank you for responding!!

Another question: I have 2 files of coastline. Seagrid doesn't allow me to load both at the same time. When I charge one, the other disappears. So I put it together both in one file. But when I charge it, there's a line between the last point of the first file and the first point of the second file.

How can I charge both files without joining this 2 different figures?

[miguel.solano]

I have never tried to load two coastlines at the same time in Seagrid, but I do know that coastline files are just a longitude, latitude pair where coastline segments are separated by Nans. So I would just try to separate file 1 and 2 by a pair of Nans so that Seagrid knows these are separated coastline segments.

Hope this helps

[chaos4u2]

Thank you very much. It worked as you said with the NaNs.

[kate]

And 2nd question: when I have my grid, how do you tell ROMS that in your domain everything is closed except 2 zones. I know there's the lateral boundary conditions, but they are north, south, etc... How do you concretise the coordenates of your domain where currents exists?

If you have an open segment on a boundary that is otherwise closed, you tell it that it's open, then mask out some of the boundary points with the land mask.