Create warning for interpolation from geodata

CHLNDDEV / OceanMesh2D

A two-dimensional triangular mesh generator with pre- and post-processing utilities written in pure MATLAB (no toolboxes required) designed specifically to build models that solve shallow-water equations or wave equations in a coastal environment (ADCIRC, FVCOM, WaveWatch3, SWAN, SCHISM, Telemac, etc.).

https://github.com/sponsors/krober10nd

GNU General Public License v3.0

182 stars 65 forks source link

Create warning for interpolation from geodata #148

Open krober10nd opened 3 years ago

krober10nd commented 3 years ago

Altering user geodata may have been downsampled to fit in RAM and it's best to probably just interpolate directly from the DEM.

m = interp(m,'DEM.nc');

HamishB commented 3 years ago

just to note, I tried changing the hardcoded AVAILABLE_MEMORY=4 GB size in geodata.m to 100 GB on our big-memory workstation but was still getting downsampling stride messages. Maybe that was just because the DEM is much finer resolution than the minimum element size; is it supposed to be like that?

krober10nd commented 3 years ago

yes. If the DEM is much finer than the minimum element size, then there's little point in reading in everything. I guess we could put some more advanced up and down sampling methods.

krober10nd commented 3 years ago

Just to note: I would strongly recommend not interpolating from the geodata object and instead using the original NetCDF file with msh.interp.

Building a mesh from a geodata object is fine, but interpolation is more critical to get right and sensitive to downsampling.