[MITgcm-support] Icefront returnin NaN based on topology.

[Menemenlis, Dimitris (3248)]

One trick, if you want to avoid computing shelfice potential, is to set
the reference density "rhoConst" to a value that is approximately equal
to the mean density of sea water in the ice shelf cavity.
Then a reasonable value of shelficeLoadAnomaly will be
computed automatically in shelfice_init_fixed.F

Note that this reference pressure is only important for reducing
model initialization shock.  After that there is no difference
to within a time-independent constant.

On Feb 19, 2014, at 7:45 AM, Dustin Carroll <dcarroll at uoregon.edu> wrote:

> Hi Harry,
> 
> If you haven’t already, you can also check out the verification/icefront/2D_example, there are m files that will show you how to compute the shelfice potential for a 2-D setup.
> 
> Cheers,
> Dustin