Should the wet/dry masks (e.g. rmask_wet, umask_wet) be used in step3d_t, step2d_lf_am3, etc? They are defined and read into those routines, but never actually used. I am having problems with anomalously high salinities and temperatures on the edges of wet/dry masks, and I wonder if fill values are being used to compute salinity and temperature.
I also have a separate/related problem implementing point sources at cells that go dry periodically, I wonder if applying the masks correctly would solve that problem also. Thanks...
-Neil