[MITgcm-devel] sea ice clean up this week

[Gael Forget]

Dear all,

below is (A) a summary of the last two day clean-up steps (cosmetic changes) 
and (B) what it means for known users (based on Patrick's table; attached).
In case you follow testreport, you already know that the modifs have not 
changed anything to results. Notice : there is more clean-up coming in the 
next couple days, so you may want to wait until doing an update, up to you.

Cheers,
Gael

(A) 

retired CPPs : 
	MCPHEE_OCEAN_ICE_HEAT_FLUX
	GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR
	SEAICE_OCN_MELT_ACT_ON_AREA, SEAICE_DO_OPEN_WATER_GROWTH
	SEAICE_DO_OPEN_WATER_MELT, FENTY_AREA_EXPANSION_CONTRACTION

added run time params. I include a copy of the doc in SEAICE_PARAMS.h and 
tried to trace back where the various options came from (text in parentheses).
	SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt
		T/F 	just like the old CPP versions
		T	use open water heat flux directly to grow ice	(Ian way)
		F 	cool ocean, and grow later if needed		(Zhang way)
	SEAICEturbFluxFormula :: selects formula for ocean-ice turbulent flux
		1=direct spec. of SEAICE_availHeatFrac/Frz		(Dimitri/Ian way)
		2=spec. via SEAICE_gamma_t/_frz				(Martin/Ian way)
		3=McPhee with linear tapering					(Ian latest way)
		4=McPhee with step function					(Ian other latest way)
	SEAICE_areaLossFormula :: selects formula for ice cover loss from melt
		1=from all but only melt conributions by ATM and OCN	(Ian way)
		2=from net melt-growth>0 by ATM and OCN		(Gael/Ian way)
		3=from predicted melt by ATM					(Zhang way)
	SEAICE_areaGainFormula :: selects formula for ice cover gain from open water growth
		1=from growth by ATM						(Ian way)
		2=from predicted growth by ATM				(Zhang way)

as a side note, I put stops in seaice_check.F	 for combinaisons
that are obviously problematic. Tyipically, if you do NOT use  
SEAICE_doOpenWaterGrowth you need to use the Zhang way
for SEAICE_areaLossFormula (3) and SEAICE_areaGainFormula (2),
and this is the only case allowed with SEAICE_GROWTH_LEGACY.

(B) 

You should double check, but based on Patrick's table, 
here is what you will need to add to your namelists : 

Martin : 	in SEAICE_GROWTH_LEGACY case you should be all set. Otherwise
		if TRUE SEAICE_doOpenWaterGrowth then SEAICE_areaLossFormula = 2
		if FALSE SEAICE_doOpenWaterGrowth then SEAICE_areaLossFormula = 3, SEAICE_areaGainFormula=2

Ian et al:	SEAICE_doOpenWaterMelt=.TRUE.
		SEAICEturbFluxFormula=3
		ps: def/undef SEAICE_OCN_MELT_ACT_ON_AREA was doing nothing 
		      since FENTY_AREA_EXPANSION_CONTRACTION was def.

Me : TBD.


-------------- next part --------------
A non-text attachment was scrubbed...
Name: seaice_whodoeswhat_3.pptx
Type: application/vnd.openxmlformats-officedocument.presentationml.presentation
Size: 69786 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20120209/0c5a38f1/attachment-0001.pptx>
-------------- next part --------------