[MITgcm-devel] sea ice clean up this week

[Gael Forget]

Dear all,

here is a summary of Friday's commit and implication for known users (of Patrick's table).

I unified the freezing point treatment and the old SEAICE_VARIABLE_FREEZING_POINT is now the default.
Ian, Martin & Matt : you can recover the old default constant freezing point by setting
SEAICE_tempFrz0    = -1.96,
SEAICE_dTempFrz_dS = 0.,

I removed the ALLOW_SEAICE_FLOODING brackets. The run time switch was already there 
and sufficient, plus everybody seems to be using it. In case you had runs 
with "#undef ALLOW_SEAICE_FLOODING" you need to set SEAICEuseFlooding=.FALSE.

I activated the sublimation term that was missing and Martin added last summer.
Matt : if you want to reproduce old results you want to remove it, by defining 
SEAICE_DISABLE_SUBLIM (as I did in all verif exp but cs32x15).

Cheers,
Gael

On Feb 9, 2012, at 2:16 PM, Gael Forget wrote:

> 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.
> 
> 
> <seaice_whodoeswhat_3.pptx>
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel