[MITgcm-devel] another stop in seaice_check: mcPheePiston

[Martin Losch]

Hi there,

I just ran across another check that made the model stop in seaice_check:
when the default
      SEAICE_mcPheePiston  =
     &      MCPHEE_TAPER_FAC * STANTON_NUMBER * USTAR_BASE
= 8.4000e-04 m/s
is larger than SEAICE_mcphee_max=dRf(kSurface)/SEAICE_deltaTtherm
the model stops. This happens to me for my very coarse cs32 experiment (basically global_ocean.cs32 with seaice), when I use the default parameters, because in this configuration dz/dt =  50m/86400s = 5.7870e-04 m/s.

Is this a useful stop? Why shouldn't mcPheePiston be larger than dz/dt? Please advise and possibly correct seaice_check.F

Martin