Changeset 13284 for NEMO/releases/r4.0/r4.0-HEAD/cfgs
- Timestamp:
- 2020-07-09T17:12:23+02:00 (4 years ago)
- Location:
- NEMO/releases/r4.0/r4.0-HEAD/cfgs
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/cfgs/AGRIF_DEMO/EXPREF/1_namelist_cfg
r13278 r13284 353 353 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T) 354 354 !----------------------------------------------------------------------- 355 rn_eice = 0 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4356 355 / 357 356 !!====================================================================== -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg
r13278 r13284 353 353 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T) 354 354 !----------------------------------------------------------------------- 355 rn_eice = 0 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4356 355 / 357 356 !!====================================================================== -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg
r13278 r13284 374 374 ! = 2 add a tke source just at the base of the ML 375 375 ! = 3 as = 1 applied on HF part of the stress (ln_cpl=T) 376 rn_eice = 0 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4377 376 / 378 377 !----------------------------------------------------------------------- -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/SHARED/field_def_nemo-ice.xml
r12337 r13284 49 49 <field id="icehpnd" long_name="melt pond depth" standard_name="sea_ice_meltpond_depth" unit="m" /> 50 50 <field id="icevpnd" long_name="melt pond volume" standard_name="sea_ice_meltpond_volume" unit="m" /> 51 <field id="icehlid" long_name="melt pond lid depth" standard_name="sea_ice_meltpondlid_depth" unit="m" /> 52 <field id="icevlid" long_name="melt pond lid volume" standard_name="sea_ice_meltpondlid_volume" unit="m" /> 51 53 52 54 <!-- heat --> … … 81 83 <field id="icediv" long_name="Divergence of the sea-ice velocity field" standard_name="divergence_of_sea_ice_velocity" unit="s-1" /> 82 84 <field id="iceshe" long_name="Maximum shear of sea-ice velocity field" standard_name="maximum_shear_of_sea_ice_velocity" unit="s-1" /> 83 85 <field id="beta_evp" long_name="Relaxation parameter of ice rheology (beta)" standard_name="relaxation_parameter_of_ice_rheology" unit="" /> 86 84 87 <!-- surface heat fluxes --> 85 88 <field id="qt_ice" long_name="total heat flux at ice surface" standard_name="surface_downward_heat_flux_in_air" unit="W/m2" /> … … 171 174 <field id="frq_m" unit="-" /> 172 175 176 <!-- rheology convergence tests --> 177 <field id="uice_cvg" long_name="sea ice velocity convergence" standard_name="sea_ice_velocity_convergence" unit="m/s" /> 178 173 179 <!-- ================= --> 174 180 <!-- Add-ons for SIMIP --> … … 209 215 <field id="dmisum" long_name="sea-ice mass change through surface melting" standard_name="tendency_of_sea_ice_amount_due_to_surface_melting" unit="kg/m2/s" /> 210 216 <field id="dmibom" long_name="sea-ice mass change through bottom melting" standard_name="tendency_of_sea_ice_amount_due_to_basal_melting" unit="kg/m2/s" /> 217 <field id="dmilam" long_name="sea-ice mass change through lateral melting" standard_name="tendency_of_sea_ice_amount_due_to_lateral_melting" unit="kg/m2/s" /> 211 218 <field id="dmsspr" long_name="snow mass change through snow fall" standard_name="snowfall_flux" unit="kg/m2/s" /> 212 219 <field id="dmsmel" long_name="snow mass change through melt" standard_name="surface_snow_melt_flux" unit="kg/m2/s" /> … … 287 294 <field id="iceapnd_cat" long_name="Ice melt pond concentration per category" unit="" /> 288 295 <field id="icehpnd_cat" long_name="Ice melt pond thickness per category" unit="m" detect_missing_value="true" /> 296 <field id="icehlid_cat" long_name="Ice melt pond lid thickness per category" unit="m" detect_missing_value="true" /> 289 297 <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category" unit="" /> 298 <field id="iceaepnd_cat" long_name="Ice melt pond effective fraction per category" unit="" /> 290 299 <field id="icemask_cat" long_name="Fraction of time step with sea ice (per category)" unit="" /> 291 300 <field id="iceage_cat" long_name="Ice age per category" unit="days" detect_missing_value="true" /> … … 298 307 <field id="snwthic_cat_cmip" long_name="Snow thickness in thickness categories" standard_name="snow_thickness_over_categories" detect_missing_value="true" unit="m" > snwthic_cat * icemask_cat + $missval * (1.-icemask_cat) </field> 299 308 <field id="iceconc_cat_pct_cmip" long_name="Sea-ice area fractions in thickness categories" standard_name="sea_ice_area_fraction_over_categories" detect_missing_value="true" unit="%" > iceconc_cat*100. * icemask_cat + $missval * (1.-icemask_cat) </field> 309 310 <!-- heat diffusion convergence tests --> 311 <field id="tice_cvgerr" long_name="sea ice temperature convergence error" standard_name="sea_ice_temperature_convergence_err" unit="K" /> 312 <field id="tice_cvgstp" long_name="sea ice temperature convergence iterations" standard_name="sea_ice_temperature_convergence_stp" unit="" /> 300 313 301 314 </field_group> <!-- SBC_3D --> … … 558 571 <field field_ref="dmisum" name="sidmassmelttop" /> 559 572 <field field_ref="dmibom" name="sidmassmeltbot" /> 573 <field field_ref="dmilam" name="sidmassmeltlat" /> 560 574 <field field_ref="dmsspr" name="sndmasssnf" /> 561 575 <field field_ref="dmsmel" name="sndmassmelt" /> -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/SHARED/field_def_nemo-oce.xml
r12288 r13284 36 36 <field id="toce_vmean300" field_ref="toce_e3t_vsum300" unit="degree_C" grid_ref="grid_T_vsum" detect_missing_value="true" > toce_e3t_vsum300/e3t_vsum300 </field> 37 37 38 <!--- additions to diawri.F90 ---> 39 <field id="socegrad" long_name="module of salinity gradient" unit="psu/m" grid_ref="grid_T_3D"/> 40 <field id="socegrad2" long_name="square of module of salinity gradient" unit="psu2/m2" grid_ref="grid_T_3D"/> 41 <field id="eken_int" long_name="vertical integration of kinetic energy" unit="m3/s2" /> 42 <field id="relvor" long_name="relative vorticity" unit="s-1" grid_ref="grid_T_3D"/> 43 <field id="absvor" long_name="absolute vorticity" unit="s-1" grid_ref="grid_T_3D"/> 44 <field id="potvor" long_name="potential vorticity" unit="s-1" grid_ref="grid_T_3D"/> 45 <field id="salt2c" long_name="Salt content vertically integrated" unit="1e-3*kg/m2" /> 38 46 39 47 <!-- t-eddy viscosity coefficients (ldfdyn) --> -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/SHARED/namelist_ice_ref
r12121 r13284 43 43 ln_cat_usr = .false. ! ice categories are defined by rn_catbnd below (m) 44 44 rn_catbnd = 0.,0.45,1.1,2.1,3.7,6.0 45 rn_himin = 0.1 ! minimum ice thickness (m) used in remapping 45 rn_himin = 0.1 ! minimum ice thickness (m) allowed 46 rn_himax = 99.0 ! maximum ice thickness (m) allowed 46 47 / 47 48 !------------------------------------------------------------------------------ … … 56 57 rn_ishlat = 2. ! lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2) 57 58 ln_landfast_L16 = .false. ! landfast: parameterization from Lemieux 2016 58 rn_ depfra= 0.125 ! fraction of ocean depth that ice must reach to initiate landfast59 rn_lf_depfra = 0.125 ! fraction of ocean depth that ice must reach to initiate landfast 59 60 ! recommended range: [0.1 ; 0.25] 60 rn_ icebfr = 15. ! maximum bottom stress per unit volume [N/m3]61 rn_lf relax= 1.e-5 ! relaxation time scale to reach static friction [s-1]62 rn_ tensile= 0.05 ! isotropic tensile strength [0-0.5??]61 rn_lf_bfr = 15. ! maximum bottom stress per unit volume [N/m3] 62 rn_lf_relax = 1.e-5 ! relaxation time scale to reach static friction [s-1] 63 rn_lf_tensile = 0.05 ! isotropic tensile strength [0-0.5??] 63 64 / 64 65 !------------------------------------------------------------------------------ … … 97 98 rn_relast = 0.333 ! ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast 98 99 ! advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 100 ln_rhg_chkcvg = .false. ! check convergence of rheology (outputs: file ice_cvg.nc & variable uice_cvg) 99 101 / 100 102 !------------------------------------------------------------------------------ 101 103 &namdyn_adv ! Ice advection 102 104 !------------------------------------------------------------------------------ 103 ln_adv_Pra = .true. ! Advection scheme (Prather)104 ln_adv_UMx = .false. 105 ln_adv_Pra = .true. ! Advection scheme (Prather) 106 ln_adv_UMx = .false. ! Advection scheme (Ultimate-Macho) 105 107 nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order) 106 108 / … … 109 111 !------------------------------------------------------------------------------ 110 112 rn_cio = 5.0e-03 ! ice-ocean drag coefficient (-) 111 rn_blow_s = 0.66 ! mesure of snow blowing into the leads 113 nn_snwfra = 0 ! calculate the fraction of ice covered by snow (for zdf and albedo) 114 ! = 0 fraction = 1 (if snow) or 0 (if no snow) 115 ! = 1 fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation] 116 ! = 2 fraction = hsnw / (hsnw+0.02) [CICE formulation] 117 rn_snwblow = 0.66 ! mesure of snow blowing into the leads 112 118 ! = 1 => no snow blowing, < 1 => some snow blowing 113 119 nn_flxdist = -1 ! Redistribute heat flux over ice categories … … 118 124 ln_cndflx = .false. ! Use conduction flux as surface boundary conditions (i.e. for Jules coupling) 119 125 ln_cndemulate = .false. ! emulate conduction flux (if not provided in the inputs) 126 nn_qtrice = 0 ! Solar flux transmitted thru the surface scattering layer: 127 ! = 0 Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow) 128 ! = 1 Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities) 120 129 / 121 130 !------------------------------------------------------------------------------ … … 126 135 ln_icedO = .true. ! activate ice growth in open-water (T) or not (F) 127 136 ln_icedS = .true. ! activate brine drainage (T) or not (F) 137 ! 138 ln_leadhfx = .true. ! heat in the leads is used to melt sea-ice before warming the ocean 128 139 / 129 140 !------------------------------------------------------------------------------ … … 135 146 rn_cnd_s = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 136 147 ! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 137 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice [1/m] 148 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice [1/m] 149 rn_kappa_s = 10.0 ! nn_qtrice = 0: radiation attenuation coefficient in snow [1/m] 150 rn_kappa_smlt = 7.0 ! nn_qtrice = 1: radiation attenuation coefficient in melting snow [1/m] 151 rn_kappa_sdry = 10.0 ! radiation attenuation coefficient in dry snow [1/m] 152 ln_zdf_chkcvg = .false. ! check convergence of heat diffusion scheme (outputs: tice_cvgerr, tice_cvgstp) 138 153 / 139 154 !------------------------------------------------------------------------------ … … 175 190 &namthd_pnd ! Melt ponds 176 191 !------------------------------------------------------------------------------ 177 ln_pnd = .false. ! activate melt ponds or not 178 ln_pnd_H12 = .false. ! activate evolutive melt ponds (from Holland et al 2012) 179 ln_pnd_CST = .false. ! activate constant melt ponds 180 rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC 181 rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC 182 ln_pnd_alb = .false. ! melt ponds affect albedo or not 192 ln_pnd = .false. ! activate melt ponds or not 193 ln_pnd_LEV = .false. ! level ice melt ponds (from Flocco et al 2007,2010 & Holland et al 2012) 194 rn_apnd_min = 0.15 ! minimum ice fraction that contributes to melt pond. range: 0.0 -- 0.15 ?? 195 rn_apnd_max = 0.85 ! maximum ice fraction that contributes to melt pond. range: 0.7 -- 0.85 ?? 196 ln_pnd_CST = .false. ! constant melt ponds 197 rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC 198 rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC 199 ln_pnd_lids = .true. ! frozen lids on top of the ponds (only for ln_pnd_LEV) 200 ln_pnd_alb = .true. ! effect of melt ponds on ice albedo 183 201 / 184 202 !------------------------------------------------------------------------------ … … 186 204 !------------------------------------------------------------------------------ 187 205 ln_iceini = .true. ! activate ice initialization (T) or not (F) 188 ln_iceini_file = .false. ! netcdf file provided for initialization (T) or not (F) 206 nn_iceini_file = 0 ! 0 = Initialise sea ice based on SSTs 207 ! 1 = Initialise sea ice from single category netcdf file 208 ! 2 = Initialise sea ice from multi category restart file 189 209 rn_thres_sst = 2.0 ! max temp. above Tfreeze with initial ice = (sst - tfreeze) 190 210 rn_hti_ini_n = 3.0 ! initial ice thickness (m), North … … 206 226 rn_hpd_ini_n = 0.05 ! initial pond depth (m), North 207 227 rn_hpd_ini_s = 0.05 ! " " South 208 ! -- for ln_iceini_file = T 228 rn_hld_ini_n = 0.0 ! initial pond lid depth (m), North 229 rn_hld_ini_s = 0.0 ! " " South 230 ! -- for nn_iceini_file = 1 209 231 sn_hti = 'Ice_initialization' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' 210 232 sn_hts = 'Ice_initialization' , -12 ,'hts' , .false. , .true., 'yearly' , '' , '', '' … … 217 239 sn_apd = 'NOT USED' , -12 ,'apd' , .false. , .true., 'yearly' , '' , '', '' 218 240 sn_hpd = 'NOT USED' , -12 ,'hpd' , .false. , .true., 'yearly' , '' , '', '' 241 sn_hld = 'NOT USED' , -12 ,'hld' , .false. , .true., 'yearly' , '' , '', '' 219 242 cn_dir='./' 220 243 / … … 238 261 ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F) 239 262 ln_icectl = .false. ! ice points output for debug (T or F) 240 iiceprt = 10 !i-index for debug241 jiceprt = 10 !j-index for debug242 / 263 iiceprt = 10 ! i-index for debug 264 jiceprt = 10 ! j-index for debug 265 / -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/SHARED/namelist_ref
r13278 r13284 281 281 sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 282 282 sn_slp = 'slp.15JUNE2009_fill' , 6. , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 283 sn_cc = 'NOT USED' , 24 , 'CC' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 283 284 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 284 285 / … … 286 287 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 287 288 !----------------------------------------------------------------------- 288 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data 289 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 290 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 291 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 289 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data 290 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 291 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 292 ln_scale_ice_flux = .false. ! use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 293 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 292 294 !_____________!__________________________!____________!_____________!______________________!________! 293 295 ! ! description ! multiple ! vector ! vector ! vector ! … … 645 647 bn_aip = 'NOT USED' , 24. , 'siapnd' , .true. , .false., 'daily' , '' , '' , '' 646 648 bn_hip = 'NOT USED' , 24. , 'sihpnd' , .true. , .false., 'daily' , '' , '' , '' 649 bn_hil = 'NOT USED' , 24. , 'sihlid' , .true. , .false., 'daily' , '' , '' , '' 647 650 ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 648 651 rn_ice_tem = 270. ! arbitrary temperature of incoming sea ice … … 651 654 rn_ice_apnd = 0.2 ! -- pond fraction = a_ip/a_i -- 652 655 rn_ice_hpnd = 0.05 ! -- pond depth -- 656 rn_ice_hlid = 0.0 ! -- pond lid depth -- 653 657 / 654 658 !----------------------------------------------------------------------- … … 679 683 ! 680 684 ln_drgimp = .true. ! implicit top/bottom friction flag 685 ln_drgice_imp = .false. ! implicit ice-ocean drag 681 686 / 682 687 !----------------------------------------------------------------------- … … 1054 1059 ! ! = 3 as =2 with distinct dissipative an mixing length scale 1055 1060 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) 1061 nn_mxlice = 0 ! type of scaling under sea-ice 1062 ! = 0 no scaling under sea-ice 1063 ! = 1 scaling with constant sea-ice thickness 1064 ! = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 1065 ! = 3 scaling with maximum sea-ice thickness 1066 rn_mxlice = 10. ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 1056 1067 rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value 1057 1068 ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002) … … 1065 1076 ! = 0 constant 10 m length scale 1066 1077 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 1067 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1078 nn_eice = 1 ! attenutaion of langmuir & surface wave breaking under ice 1079 ! ! = 0 no impact of ice cover on langmuir & surface wave breaking 1080 ! ! = 1 weigthed by 1-TANH(10*fr_i) 1081 ! ! = 2 weighted by 1-fr_i 1082 ! ! = 3 weighted by 1-MIN(1,4*fr_i) 1068 1083 / 1069 1084 !----------------------------------------------------------------------- … … 1078 1093 rn_charn = 70000. ! Charnock constant for wb induced roughness length 1079 1094 rn_hsro = 0.02 ! Minimum surface roughness 1095 rn_hsri = 0.03 ! Ice-ocean roughness 1080 1096 rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met>1) 1081 1097 nn_z0_met = 2 ! Method for surface roughness computation (0/1/2/3) 1082 ! ! =3 requires ln_wave=T 1098 ! ! = 3 requires ln_wave=T 1099 nn_z0_ice = 1 ! attenutaion of surface wave breaking under ice 1100 ! ! = 0 no impact of ice cover 1101 ! ! = 1 roughness uses rn_hsri and is weigthed by 1-TANH(10*fr_i) 1102 ! ! = 2 roughness uses rn_hsri and is weighted by 1-fr_i 1103 ! ! = 3 roughness uses rn_hsri and is weighted by 1-MIN(1,4*fr_i) 1083 1104 nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum) 1084 1105 nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum) -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/SPITZ12/EXPREF/namelist_cfg
r13278 r13284 217 217 ln_loglayer = .true. ! logarithmic drag: Cd = vkarmn/log(z/z0) |U| 218 218 ln_drgimp = .true. ! implicit top/bottom friction flag 219 ln_drgice_imp = .true. ! implicit ice-ocean drag 219 220 / 220 221 !----------------------------------------------------------------------- … … 340 341 nn_havtb = 1 ! horizontal shape for avtb (=1) or not (=0) 341 342 / 343 !----------------------------------------------------------------------- 344 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T) 345 !----------------------------------------------------------------------- 346 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) 347 nn_mxlice = 2 ! type of scaling under sea-ice 348 ! = 0 no scaling under sea-ice 349 ! = 1 scaling with constant sea-ice thickness 350 ! = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 351 ! = 3 scaling with maximum sea-ice thickness 352 nn_eice = 1 ! attenutaion of langmuir & surface wave breaking under ice 353 ! ! = 0 no impact of ice cover on langmuir & surface wave breaking 354 ! ! = 1 weigthed by 1-TANH(10*fr_i) 355 ! ! = 2 weighted by 1-fr_i 356 ! ! = 3 weighted by 1-MIN(1,4*fr_i) 357 / 342 358 !!====================================================================== 343 359 !! *** Diagnostics namelists *** !! -
NEMO/releases/r4.0/r4.0-HEAD/cfgs/SPITZ12/EXPREF/namelist_ice_cfg
r11731 r13284 82 82 !------------------------------------------------------------------------------ 83 83 ln_pnd = .true. ! activate melt ponds or not 84 ln_pnd_H12 = .true. ! activate evolutive melt ponds (from Holland et al 2012) 85 ln_pnd_alb = .true. ! melt ponds affect albedo or not 84 ln_pnd_LEV = .true. ! activate level ice melt ponds 86 85 / 87 86
Note: See TracChangeset
for help on using the changeset viewer.