[MITgcm-devel] question regarding seaice_solve4temp.F

[Jean-Michel Campin]

Hi Ian,

I have a question regarding this part of the code, formerly in
seaice_budget_ice_if.F and now in seaice_solve4temp.F (non-legacy):
> c    If the search falls below 50 Kelvin then kick the search back up to
> c    TMELT.  Note that a solution to the equation is for a large negative
> c    value of ice surface temperature since the longwave outgoing radiation
> c    goes as the fourth power of temperature.
>           IF (tsurfLoc(I,J) .LT. 50.0 ) THEN
>                tsurfLoc(I,J) = TMELT
>           ENDIF
Does this fix a problem that can happen with the new saturation vapor 
pressure (SVP) formulation ? 
It seems to me that with the new (SVP), the search for tsurf looks
quite safe, with always dFiDTs1 > 0. May be I am missing something ?
If you could provide an example where this problem show-up, it would be nice.

Or was it something you added when trying to add the latent-heat derivative 
relative to tsurf with the old SVP (Maykut polynomia) ?
In this case, I can understand that low temp can be a problem.

Cheers,
Jean-Michel