[MITgcm-devel] depressed Eta and sIceLoad

[Dimitris Menemenlis]

Martin, my vote is to add all the fixes, as you describe in your e-mail, without 
caring for backward compatibility.  Although I have not yet tried it, I agree 
that the evap is potentially the most serious bug, especially for those 
integrations that include Southern Ocean.  Regarding too much snow accumulation, 
if the problem is with the NCEP or CORE forcing fields, then we should be fixing 
the NCEP or CORE forcing fields, not the model.

The link:
http://ecco2.jpl.nasa.gov/data2/cube/cube38/pickup/HSNOW.jpg
is live once again, it was down when you last tried to access it.
This particular cubed-sphere integration is driven by NCEP reanalysis.
As Jinlun pointed out the "too-much snow" problem is mostly around Antarctica.
In the Arctic it is less of a problem, although there are some small regions 
with unrealistic 25 m snow around Greenland and in the Canadian Archipelago.

The companion figures to above, that show complete contents of seaice_pickup 
from same integration and time in Arctic and Antarctic are
http://ecco2.jpl.nasa.gov/data2/cube/cube38/pickup/pickup_seaiceArctic.jpg
http://ecco2.jpl.nasa.gov/data2/cube/cube38/pickup/pickup_seaiceAntarctic.jpg
A yet different problem is ice that becomes too thick, up to 25 m in some 
isolated spots.  Shall we use this same occasion to uncomment

#ifdef DO_WE_NEED_THIS
c          HEFF(I,J,1,bi,bj)=MIN(MAX_HEFF,HEFF(I,J,1,bi,bj))
#endif

in growth.F.  We can set a large default MAX_HEFF in seaice_readparms.F

Dimitris