[MITgcm-devel] Unrealistic low SST with SEAICE_GROWTH_LEGACY undef

[Martin Losch]

Hi Jean-Michel,

here's my data.seaice (with very recent code). If have started another run with this set of options:
C-    Ian suggests this:
#if (!defined SEAICE_SOLVE4TEMP_LEGACY && !defined SEAICE_GROWTH_LEGACY)
# undef SEAICE_HEAT_CONSERV_FIX
C# define MCPHEE_OCEAN_ICE_HEAT_FLUX
# undef MCPHEE_OCEAN_ICE_HEAT_FLUX
# define GRADIENT_MIXED_LAYER_TURBULENCE_FACTOR
# define SEAICE_DO_OPEN_WATER_GROWTH
# define SEAICE_DO_OPEN_WATER_MELT
# define FENTY_AREA_EXPANSION_CONTRACTION
# define SEAICE_ADD_SUBLIMATION_TO_FWBUDGET
#endif
(so MCPHEE_OCEAN_ICE_HEAT_FLUX undefined, so that the GRADIENT_* flag is not used either)
and after 4months, the minimum temperatures (based on daily means) are -1.98456degC which is what I expect. So the preliminar conclusion is that this MCPHEE-flag is a major player in this.

M.

cat data.seaice 
# SEAICE parameters
 &SEAICE_PARM01
 SEAICEwriteState   = .TRUE.,
 SEAICE_initialHEFF = 0.0,
#SEAICE_waterTurnAngle=25.,
# lead closing parameter, default is 0.5
 HO = 1.,
 SEAICEuseFlooding = .true.,
 SEAICEadvSnow     = .true.,
 SEAICEadvAge      = .true.,
 SEAICEadvScheme   = 33,
 DIFF1             = 0.,
 SEAICE_clipVelocities = .false.,
 SEAICE_no_slip    = .true.
# I think we need this to have some salinity in the ice
 SIsalFRAC          = 0.3,
#LSR_error         = 1.e-7,
 SEAICE_gamma_t   = 259200.,
# for AOMIP stupid convention
 SEAICE_strength    = 15000.,
#SEAICE_strength    = 27500.,
 SEAICErestoreUnderIce = .true.,
 SEAICE_drag = 0.001,
 ocean_drag  = 0.001,
 SEAICE_dryIceAlb   = 0.70,
 SEAICE_wetIceAlb   = 0.68,
 SEAICE_drySnowAlb  = 0.81,
 SEAICE_wetSnowAlb  = 0.77,
 SEAICE_waterAlbedo = 0.1,
#SEAICE_waterAlbedo = 0.066,
 SEAICE_rhoIce      = 900.,
# CORE
# SEAICE_rhoAir  = 1.22,
 SEAICE_dalton = 1.63e-03,
#SEAICE_lhSublim = 2.839e06,
 SEAICE_rhoIce  = 900.,
 SEAICE_rhoSnow = 300.,
 &
# stupid dummy namelist
 &SEAICE_PARM03
 &

On Jan 31, 2012, at 4:19 PM, Jean-Michel Campin wrote:

> Hi,
> 
> I am started a new thread (since I thought the other one deviate on 
> different subjets).
> 
> Martin and Gael, can you send the data.seaice that leads to 
> unrealistic low SST (~ -20.C) ?
> I know you alreday detailed some of SEAICE_OPTIONS.h, but may be 
> could be useful to have it together with data.seaice.
> 
> Gael, would be nice also to have the modified data.seaice which
> fix the Pb, with some idea about how low the SST gets.
> 
> Thanks,
> Jean-Michel
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel