[MITgcm-support] sea-ice, freeze, etc.

[Dimitris Menemenlis]

Jean-Michel and I have come up with following list of suggested changes
for the sea-ice models.

1. The two sea-ice models and subroutine freeze will be moved to
forward_step, following loading of external fields and preceding
thermodynamnics.

2. By default subroutine freeze will be 2-D and act at the surface
only.  That, and its new location, cause some small changes to
output.txt for three verification experiments.

3. The sea-ice packages and subroutine freeze will modify temperature
directly but generate an array called surfaceTendencyTice that keeps
track of equivalent heat flux.

4. Net surface heat and wind stress diagnostics will be added to
timeave_statvars.  These will include relaxation and sea-ice terms.

5. Subroutine external_forcing_surf will remain.  It will be called at
very beginning (Jean-Michel's preference) or right before
thermodynamics.

6. Qnet to include Qsw by default and signs of Qnet and Qsw to change
back to "positive for increase in theta", with runtime options that
maintain backward compatibility for the verification experiments.

7. Make sure that shortwave forcing does work without KPP.

8. Surface salt flux changed to MKS units.  Sign convention will depend
on whether 6. is adopted.


-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Insitute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720