Mesh Corner Elements coordinate error

[chlnddev]

Mesh Corner Elements coordinate error

[WPringle]

Hi Drak,

I'm not sure what you mean by coordinate error and why the triangles should look like that in Blue Kenue...

I ran your script with the gshhs shoreline and the grb bathymetry and the coordinates were fine.

What I'm seeing here is actually a problem with the interpolation of the bathymetry. When you use the gdat to interpolate sometimes it can miss the corners because the point is just outside of the extracted bathy data. A simple way to avoid this is to use the 'nan','fill' option, i.e.,
m = interp(m,gdat,'nan','fill');

or you can just interpolate directly from the dem:
dem = 'gbr100_16jan.nc';
m = interp(m,dem);

When the dem is used to interpolate it will extract the region of interest from the dem with a buffer region so it will not have issues at the corners.

I noticed an issue with using the grb bathymetry to interpolate directly due to the y-coordinates being in the upside down direction - I will need to make sure it works for this case too.

[DrakonianMight]

Thanks for getting back to me so fast!

Yes interpolating directly from the DEM solved the issue.

FYI, with Blue Kenue it appears that when it ingests the .14 format it doesn't pull in the coordinates properly if the elevation isn't there/ is not what it's expecting.

Great job with the module though, it is taking me a fraction of the time to make some nice mesh's.

[WPringle]

Yeah that's interesting about Blue Kenue, I figured something like that but the default bathymetry is usually just zero so I didn't think that a bathymetry value of zero would cause such a strange visualization.

Anyway, great to hear it worked and looking forward to seeing more meshes.