[MITgcm-devel] sea ice clean up this week

Jean-Michel Campin jmc at ocean.mit.edu
Fri Feb 10 11:56:51 EST 2012


Hi, 

Just a comment about this:
> PS. I am surprised to see a MS-PPT document in this list, check this out:
> http://www.gnu.org/philosophy/no-word-attachments.html
I tried to read it with my OpenOffice, and the 1rst raw was just black,
which makes it less useful. I had to ask Gael to sent a pdf which I am
attaching here.

Cheers,
Jean-Michel

On Fri, Feb 10, 2012 at 08:39:17AM +0100, Martin Losch wrote:
> Hi Gael,
> 
> thanks for this information. I'd appreciate it if this could end up somewhere in the documentation and/or in the header section of one of the seaice-files, e.g. SEAICE_OPTIONS.h or seaice_growth.F. You've almost done is already, but what would be useful is list of corresponding parameter values, e.g.
> 
> MCPHEE_OCEAN_ICE_HEAT_FLUX + 	GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR defined => SEAICEturbFluxFormula = 3
> MCPHEE_OCEAN_ICE_HEAT_FLUX defined  => SEAICEturbFluxFormula = 4
> 
> SEAICE_DO_OPEN_WATER_GROWTH + SEAICE_DO_OPEN_WATER_MELT  defined => SEAICE_areaLossFormula = 2, SEAICE_areaGainFormula=2
> 
> etc.
> 
> The configuration I sent to Patrick is the one that I am/we are routinely using, but I am also playing with all sorts of other combinations of flags (which I thought would be too much to send, my mistake). As a consquence I am currently a lttle confused, how I can convert my different CPP-flag combination into the appropriate run time flags. Would that be possible to do?
> 
> Martin
> 
> PS. I am surprised to see a MS-PPT document in this list, check this out:
> http://www.gnu.org/philosophy/no-word-attachments.html
> 
> 
> On Feb 9, 2012, at 8: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
> 
> 
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
-------------- next part --------------
A non-text attachment was scrubbed...
Name: seaice_whodoeswhat_3-1.pdf
Type: application/pdf
Size: 39874 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20120210/ffdae856/attachment-0001.pdf>


More information about the MITgcm-devel mailing list