--- MITgcm/verification/offline_exf_seaice/input_ad/data.seaice 2012/01/17 15:29:27 1.10 +++ MITgcm/verification/offline_exf_seaice/input_ad/data.seaice 2013/03/25 02:37:36 1.16 @@ -1,27 +1,48 @@ # SEAICE parameters &SEAICE_PARM01 -#- old DIFF1 default: - DIFF1 = 4.E-3, - SEAICE_initialHEFF = 1.0, - SEAICE_deltaTtherm = 3600., - SEAICE_deltaTdyn = 3600., - SEAICEuseDYNAMICS =.FALSE., -#-- According to Martin, SEAICE_clipVelocities is not recommended -# SEAICE_clipVelocities=.TRUE., - SEAICEadvSnow =.FALSE., -#-- above: to reproduce old results - LSR_ERROR = 1.E-12, - SEAICEwriteState = .TRUE., + SEAICErestoreUnderIce=.TRUE., +#same freezing point as pkg/thsice + SEAICE_tempFrz0 = 0.0, + SEAICE_dTempFrz_dS = -5.4e-02, +#seaice growth parameters + SEAICE_frazilFrac = 1., +#- same mcPheePiston value as in thsice with surf-current=0.2 m/s + SEAICE_mcPheePiston= 8.7854425e-5, + SEAICEuseFlooding = .TRUE., + SEAICE_salt0 = 4., +#SEAICE_saltFrac = 0.3, +#simplest (constant) albedo scheme + SEAICE_dryIceAlb = 0.6, + SEAICE_wetIceAlb = 0.6, + SEAICE_drySnowAlb = 0.6, + SEAICE_wetSnowAlb = 0.6, +#SEAICE_waterAlbedo = 0.1, +# + SEAICE_strength = 2.6780e+04, + OCEAN_drag = 8.1541e-04, + SEAICE_waterDrag = 5.3508, + SEAICEuseDYNAMICS =.FALSE., + LSR_ERROR = 1.E-12, + SOLV_MAX_ITERS = 1500, + LSR_mixIniGuess = 1, + SEAICEadvScheme = 77, + SEAICEadvSnow = .TRUE., +#SEAICEadvSalt = .TRUE., + SEAICE_no_Slip = .FALSE., +# + HsnowFile = 'const_00.bin', + HsaltFile = 'const_00.bin', + AreaFile = 'ice0_area.bin', + HeffFile = 'ice0_heff.bin', +#SEAICE_initialHEFF = 0.2, +#uIceFile = 'const_00.bin', + SEAICEwriteState = .TRUE., + SEAICE_monFreq = 43200., & &SEAICE_PARM02 - mult_ice = 1., -# choose which seaice cost term you want + mult_ice = 1.E-15, cost_ice_flag = 1, - costIceStart1 = 20000101, - costIceStart2 = 00000, - costIceEnd1 = 20000201, - costIceEnd2 = 00000, & &SEAICE_PARM03