- Timestamp:
- 2017-09-08T17:09:25+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref
r8512 r8514 40 40 !------------------------------------------------------------------------------ 41 41 ! -- ice_rdgrft_strength -- ! 42 ln_str_H ib= .true. ! ice strength param.: Hibler_79 => P = pstar*<h>*exp(-c_rhg*A)42 ln_str_H79 = .true. ! ice strength param.: Hibler_79 => P = pstar*<h>*exp(-c_rhg*A) 43 43 rn_pstar = 2.0e+04 ! ice strength thickness parameter (N/m2) 44 44 rn_crhg = 20.0 ! ice strength conc. parameter (-) 45 ln_str_R ot= .false. ! ice strength param.: Rothrock_75 => P = Cf*coeff*integral(wr.h^2)45 ln_str_R75 = .false. ! ice strength param.: Rothrock_75 => P = Cf*coeff*integral(wr.h^2) 46 46 rn_perdg = 17.0 ! ridging work divided by pot. energy change in ridging 47 47 ! -- ice_rdgrft -- ! … … 94 94 ln_adv_Pra = .false. ! Advection scheme (Prather) 95 95 ln_adv_UMx = .true. ! Advection scheme (Ultimate-Macho) 96 nn_UMx = 5 ! order of the scheme for UMx (1-5 )96 nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order) 97 97 / 98 98 !------------------------------------------------------------------------------ 99 99 &namice_thd ! Ice thermodynamics 100 100 !------------------------------------------------------------------------------ 101 ln_limthd = .true. ! ice thermo (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 102 ! -- limthd_dif -- ! 103 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice (m-1) 104 nn_ice_thcon = 1 ! sea ice thermal conductivity 105 ! 0: k = k0 + beta.S/T (Untersteiner, 1964) 106 ! 1: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 107 ln_dqnsice = .true. ! change the surface non-solar flux with surface temperature (T) or not (F) 108 rn_cdsn = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 109 ! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 110 ! -- limthd_dh -- ! 111 ln_limdH = .true. ! activate ice thickness change from growing/melting (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 112 rn_betas = 0.66 ! exponent in lead-ice repratition of snow precipitation 113 ! betas = 1 -> equipartition, betas < 1 -> more on leads 114 ! -- limthd_da -- ! 115 ln_limdA = .true. ! activate lateral melting param. (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 116 rn_beta = 1.0 ! (ln_latmelt=T) coef. beta for lateral melting param. Recommended range=[0.8-1.2] 117 ! => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2) 118 ! 0.3 = best fit for western Fram Strait and Antarctica 119 ! 1.4 = best fit for eastern Fram Strait 120 rn_dmin = 8. ! (ln_latmelt=T) minimum floe diameter for lateral melting param. Recommended range=[6-10] 121 ! => 6 vs 8m = +40% melting at the peak (A~0.5) 122 ! 10 vs 8m = -20% melting 123 ! -- limthd_lac -- ! 124 ln_limdO = .true. ! activate ice growth in open-water (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 125 rn_hnewice = 0.1 ! thickness for new ice formation in open water (m) 126 ln_frazil = .false. ! Frazil ice parameterization (ice collection as a function of wind) 127 rn_maxfrazb = 1.0 ! (ln_frazil=T) maximum fraction of frazil ice collecting at the ice base 128 rn_vfrazb = 0.417 ! (ln_frazil=T) thresold drift speed for frazil ice collecting at the ice bottom (m/s) 129 rn_Cfrazb = 5.0 ! (ln_frazil=T) squeezing coefficient for frazil ice collecting at the ice bottom 101 ln_icethd = .true. ! ice thermo (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 102 ! -- limthd_dif -- ! 103 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice (m-1) 104 ln_cndi_U64 = .false. ! sea ice thermal conductivity: k = k0 + beta.S/T (Untersteiner, 1964) 105 ln_cndi_P07 = .true. ! sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 106 ln_dqns_i = .true. ! change the surface non-solar flux with surface temperature (T) or not (F) 107 rn_cnd_s = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 108 ! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 109 ! -- limthd_dh -- ! 110 ln_icedH = .true. ! activate ice thickness change from growing/melting (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 111 rn_blow_s = 0.66 ! mesure of snow blowing into the leads 112 ! = 1 => no snow blowing, < 1 => some snow blowing 113 ! -- limthd_da -- ! 114 ln_icedA = .true. ! activate lateral melting param. (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 115 rn_beta = 1.0 ! coef. beta for lateral melting param. Recommended range=[0.8-1.2] 116 ! => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2) 117 ! 0.3 = best fit for western Fram Strait and Antarctica 118 ! 1.4 = best fit for eastern Fram Strait 119 rn_dmin = 8. ! minimum floe diameter for lateral melting param. Recommended range=[6-10] 120 ! => 6 vs 8m = +40% melting at the peak (A~0.5) 121 ! 10 vs 8m = -20% melting 122 ! -- limthd_lac -- ! 123 ln_icedO = .true. ! activate ice growth in open-water (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 124 rn_hinew = 0.1 ! thickness for new ice formation in open water (m) 125 ln_frazil = .false. ! Frazil ice parameterization (ice collection as a function of wind) 126 rn_maxfraz = 1.0 ! maximum fraction of frazil ice collecting at the ice base 127 rn_vfraz = 0.417 ! thresold drift speed for frazil ice collecting at the ice bottom (m/s) 128 rn_Cfraz = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom 130 129 ! -- limitd_th -- ! 131 rn_himin = 0.1 ! minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 132 ! -- icestp -- ! 133 nn_limflx = -1 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 134 ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 135 ! = 0 Average per-category fluxes (forced and coupled mode) 136 ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 137 ! = 2 Redistribute a single flux over categories (coupled mode only) 130 rn_himin = 0.1 ! minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 131 ! -- icestp -- ! 132 nn_iceflx = -1 ! Redistribute heat flux over ice categories 133 ! =-1 Do nothing (needs N(cat) fluxes) 134 ! ==> forced mode only 135 ! = 0 Average N(cat) fluxes then apply the average over the N(cat) ice 136 ! ==> forced and coupled modes 137 ! = 1 Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity 138 ! ==> forced mode only 139 ! = 2 Redistribute a single flux over categories 140 ! ==> coupled mode only 138 141 / 139 142 !------------------------------------------------------------------------------ 140 143 &namice_sal ! Ice salinity 141 144 !------------------------------------------------------------------------------ 142 ! -- limthd_sal -- !143 ln_ limdS = .true.! activate gravity drainage and flushing (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO144 nn_icesal = 2 ! ice salinity option145 ! 1: constant ice salinity (S=rn_icesal)146 ! 2: varying salinity parameterization S(z,t)147 ! 3: prescribed salinity profile S(z), Schwarzacher, 1959148 rn_icesal = 4. ! (nn_icesal=1) ice salinity (g/kg)149 rn_sal_gd = 5. ! restoring ice salinity, gravity drainage (g/kg)150 rn_time_gd = 1.73e+6 ! restoring time scale, gravity drainage (s)151 rn_sal_fl = 2. ! restoring ice salinity, flushing (g/kg)152 rn_time_fl = 8.64e+5 ! restoring time scale, flushing (s)153 rn_simax = 20. ! maximum tolerated ice salinity (g/kg)154 rn_simin = 0.1 ! minimum tolerated ice salinity (g/kg)145 ! -- limthd_sal -- ! 146 ln_icedS = .true. ! activate gravity drainage and flushing (T) or not (F) => DO NOT TOUCH UNLESS U KNOW WHAT U DO 147 nn_icesal = 2 ! ice salinity option 148 ! 1: constant ice salinity (S=rn_icesal) 149 ! 2: varying salinity parameterization S(z,t) 150 ! 3: prescribed salinity profile S(z), Schwarzacher, 1959 151 rn_icesal = 4. ! (nn_icesal=1) ice salinity (g/kg) 152 rn_sal_gd = 5. ! restoring ice salinity, gravity drainage (g/kg) 153 rn_time_gd = 1.73e+6 ! restoring time scale, gravity drainage (s) 154 rn_sal_fl = 2. ! restoring ice salinity, flushing (g/kg) 155 rn_time_fl = 8.64e+5 ! restoring time scale, flushing (s) 156 rn_simax = 20. ! maximum tolerated ice salinity (g/kg) 157 rn_simin = 0.1 ! minimum tolerated ice salinity (g/kg) 155 158 / 156 159 !------------------------------------------------------------------------------ 157 160 &namicemp ! Melt ponds 158 161 !------------------------------------------------------------------------------ 159 ! -- limmp -- !160 ln_pnd = .false. ! active melt ponds161 ln_pnd_rad = .false. ! active melt ponds radiative coupling162 ln_pnd_fw = .false. ! active melt ponds freshwater coupling163 nn_pnd_scheme = 0! type of melt pond scheme : =0 prescribed ( Tsu=0 ), =1 empirical, =2 topographic164 rn_apnd = 0.2! prescribed pond fraction, at Tsu=0 : (0<rn_apnd<1, nn_pnd_scheme = 0)165 rn_hpnd = 0.05! prescribed pond depth, at Tsu=0 : (0<rn_apnd<1, nn_pnd_scheme = 0)162 ! -- limmp -- ! 163 ln_pnd = .false. ! active melt ponds 164 ln_pnd_rad = .false. ! active melt ponds radiative coupling 165 ln_pnd_fw = .false. ! active melt ponds freshwater coupling 166 nn_pnd_scheme = 0 ! type of melt pond scheme : =0 prescribed ( Tsu=0 ), =1 empirical, =2 topographic 167 rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 : (0<rn_apnd<1, nn_pnd_scheme = 0) 168 rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 : (0<rn_apnd<1, nn_pnd_scheme = 0) 166 169 / 167 170 !------------------------------------------------------------------------------ 168 171 &namice_ini ! Ice initialization 169 172 !------------------------------------------------------------------------------ 170 ! -- limistate -- !171 ln_iceini = .true. ! activate ice initialization (T) or not (F)172 ln_iceini_file = .false. ! netcdf file provided for initialization (T) or not (F)173 rn_thres_sst = 2.0! maximum water temperature with initial ice (degC)174 rn_hts_ini_n = 0.3! initial real snow thickness (m), North175 rn_hts_ini_s = 0.3! " " South176 rn_hti_ini_n = 3.0! initial real ice thickness (m), North177 rn_hti_ini_s = 1.0! " " South178 rn_ati_ini_n = 0.9! initial ice concentration (-), North179 rn_ati_ini_s = 0.9! " " South180 rn_smi_ini_n = 6.3! initial ice salinity (g/kg), North181 rn_smi_ini_s = 6.3! " " South182 rn_tmi_ini_n = 270.! initial ice/snw temperature (K), North183 rn_tmi_ini_s = 270.! " " South173 ! -- limistate -- ! 174 ln_iceini = .true. ! activate ice initialization (T) or not (F) 175 ln_iceini_file = .false. ! netcdf file provided for initialization (T) or not (F) 176 rn_thres_sst = 2.0 ! maximum water temperature with initial ice (degC) 177 rn_hts_ini_n = 0.3 ! initial real snow thickness (m), North 178 rn_hts_ini_s = 0.3 ! " " South 179 rn_hti_ini_n = 3.0 ! initial real ice thickness (m), North 180 rn_hti_ini_s = 1.0 ! " " South 181 rn_ati_ini_n = 0.9 ! initial ice concentration (-), North 182 rn_ati_ini_s = 0.9 ! " " South 183 rn_smi_ini_n = 6.3 ! initial ice salinity (g/kg), North 184 rn_smi_ini_s = 6.3 ! " " South 185 rn_tmi_ini_n = 270. ! initial ice/snw temperature (K), North 186 rn_tmi_ini_s = 270. ! " " South 184 187 185 188 sn_hti = 'Ice_initialization' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' … … 194 197 &namice_alb ! albedo parameters 195 198 !------------------------------------------------------------------------------ 196 nn_ice_alb = 1 ! parameterization of ice/snow albedo 197 ! 0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo 198 ! 1: "home made" based on Brandt et al. (JClim 2005) and Grenfell & Perovich (JGR 2004), 199 ! giving cloud-sky albedo 200 rn_alb_sdry = 0.85 ! dry snow albedo : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky) 201 rn_alb_smlt = 0.75 ! melting snow albedo : 0.65 ( '' ) ; 0.75 ( '' ) ; obs 0.72-0.82 ( '' ) 202 rn_alb_idry = 0.60 ! dry ice albedo : 0.72 ( '' ) ; 0.60 ( '' ) ; obs 0.54-0.65 ( '' ) 203 rn_alb_imlt = 0.50 ! bare puddled ice albedo : 0.53 ( '' ) ; 0.50 ( '' ) ; obs 0.49-0.58 ( '' ) 204 rn_alb_dpnd = 0.27 ! ponded ice albedo : 0.25 ( '' ) ; 0.27 ( '' ) ; obs 0.10-0.30 ( '' ) 199 nn_ice_alb = 1 ! parameterization of ice/snow albedo 200 ! 0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo 201 ! 1: "home made" based on Brandt et al. (JClim 2005) 202 ! and Grenfell & Perovich (JGR 2004), giving cloud-sky albedo 203 ! 2: as 1 with melt ponds 204 rn_alb_sdry = 0.85 ! dry snow albedo : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky) 205 rn_alb_smlt = 0.75 ! melting snow albedo : 0.65 ( '' ) ; 0.75 ( '' ) ; obs 0.72-0.82 ( '' ) 206 rn_alb_idry = 0.60 ! dry ice albedo : 0.72 ( '' ) ; 0.60 ( '' ) ; obs 0.54-0.65 ( '' ) 207 rn_alb_imlt = 0.50 ! bare puddled ice albedo : 0.53 ( '' ) ; 0.50 ( '' ) ; obs 0.49-0.58 ( '' ) 208 rn_alb_dpnd = 0.27 ! ponded ice albedo : 0.25 ( '' ) ; 0.27 ( '' ) ; obs 0.10-0.30 ( '' ) 205 209 / 206 210 !------------------------------------------------------------------------------ 207 211 &namice_dia ! Diagnostics 208 212 !------------------------------------------------------------------------------ 209 ln_ limdiachk = .false.! check online the heat, mass & salt budgets (T) or not (F)210 ln_ limdiahsb = .false.! output the heat, mass & salt budgets (T) or not (F)211 ln_ limctl = .false.! ice points output for debug (T or F)212 iiceprt = 10! i-index for debug213 jiceprt = 10! j-index for debug214 / 213 ln_icediachk = .false. ! check online the heat, mass & salt budgets (T) or not (F) 214 ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F) 215 ln_icectl = .false. ! ice points output for debug (T or F) 216 iiceprt = 10 ! i-index for debug 217 jiceprt = 10 ! j-index for debug 218 /
Note: See TracChangeset
for help on using the changeset viewer.