Changeset 13472 for NEMO/trunk/cfgs/SHARED
- Timestamp:
- 2020-09-16T15:05:19+02:00 (4 years ago)
- Location:
- NEMO/trunk/cfgs/SHARED
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/cfgs/SHARED/field_def_nemo-ice.xml
r12377 r13472 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" /> … … 173 176 <field id="frq_m" unit="-" /> 174 177 178 <!-- rheology convergence tests --> 179 <field id="uice_cvg" long_name="sea ice velocity convergence" standard_name="sea_ice_velocity_convergence" unit="m/s" /> 180 175 181 <!-- ================= --> 176 182 <!-- Add-ons for SIMIP --> … … 211 217 <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" /> 212 218 <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" /> 219 <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" /> 213 220 <field id="dmsspr" long_name="snow mass change through snow fall" standard_name="snowfall_flux" unit="kg/m2/s" /> 214 221 <field id="dmsmel" long_name="snow mass change through melt" standard_name="surface_snow_melt_flux" unit="kg/m2/s" /> … … 289 296 <field id="iceapnd_cat" long_name="Ice melt pond concentration per category" unit="" /> 290 297 <field id="icehpnd_cat" long_name="Ice melt pond thickness per category" unit="m" detect_missing_value="true" /> 298 <field id="icehlid_cat" long_name="Ice melt pond lid thickness per category" unit="m" detect_missing_value="true" /> 291 299 <field id="iceafpnd_cat" long_name="Ice melt pond fraction per category" unit="" /> 300 <field id="iceaepnd_cat" long_name="Ice melt pond effective fraction per category" unit="" /> 292 301 <field id="icemask_cat" long_name="Fraction of time step with sea ice (per category)" unit="" /> 293 302 <field id="iceage_cat" long_name="Ice age per category" unit="days" detect_missing_value="true" /> … … 300 309 <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> 301 310 <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> 311 312 <!-- heat diffusion convergence tests --> 313 <field id="tice_cvgerr" long_name="sea ice temperature convergence error" standard_name="sea_ice_temperature_convergence_err" unit="K" /> 314 <field id="tice_cvgstp" long_name="sea ice temperature convergence iterations" standard_name="sea_ice_temperature_convergence_stp" unit="" /> 302 315 303 316 </field_group> <!-- SBC_3D --> … … 560 573 <field field_ref="dmisum" name="sidmassmelttop" /> 561 574 <field field_ref="dmibom" name="sidmassmeltbot" /> 575 <field field_ref="dmilam" name="sidmassmeltlat" /> 562 576 <field field_ref="dmsspr" name="sndmasssnf" /> 563 577 <field field_ref="dmsmel" name="sndmassmelt" /> -
NEMO/trunk/cfgs/SHARED/field_def_nemo-oce.xml
r13214 r13472 129 129 <!-- AGRIF sponge --> 130 130 <field id="agrif_spt" long_name=" AGRIF t-sponge coefficient" unit=" " /> 131 132 <!-- additions to diawri.F90 --> 133 <field id="socegrad" long_name="module of salinity gradient" unit="psu/m" grid_ref="grid_T_3D"/> 134 <field id="socegrad2" long_name="square of module of salinity gradient" unit="psu2/m2" grid_ref="grid_T_3D"/> 135 <field id="ke" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D" /> 136 <field id="ke_int" long_name="vertical integration of kinetic energy" unit="m3/s2" /> 137 <field id="relvor" long_name="relative vorticity" unit="s-1" grid_ref="grid_T_3D"/> 138 <field id="absvor" long_name="absolute vorticity" unit="s-1" grid_ref="grid_T_3D"/> 139 <field id="potvor" long_name="potential vorticity" unit="s-1" grid_ref="grid_T_3D"/> 140 <field id="salt2c" long_name="Salt content vertically integrated" unit="1e-3*kg/m2" /> 131 141 132 142 <!-- t-eddy viscosity coefficients (ldfdyn) --> … … 177 187 <field id="alpha" long_name="thermal expansion" unit="degC-1" grid_ref="grid_T_3D" /> 178 188 <field id="beta" long_name="haline contraction" unit="1e3" grid_ref="grid_T_3D" /> 179 <field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="s-1" grid_ref="grid_T_3D" />180 189 <field id="rhop" long_name="potential density (sigma0)" standard_name="sea_water_sigma_theta" unit="kg/m3" grid_ref="grid_T_3D" /> 181 190 182 191 <!-- Energy - horizontal divergence --> 183 <field id="eken" long_name="kinetic energy" standard_name="specific_kinetic_energy_of_sea_water" unit="m2/s2" grid_ref="grid_T_3D" />184 192 <field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D" /> 185 193 … … 655 663 <field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" /> 656 664 665 <!-- EOS --> 666 <field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="s-2" /> 667 657 668 </field_group> 658 669 … … 710 721 <field id="masstr_strait" long_name="Sea water transport across line" grid_ref="grid_4strait" > u_masstr_strait + v_masstr_strait </field> 711 722 </field_group> 712 713 723 714 724 <!-- variables available with ln_floats --> -
NEMO/trunk/cfgs/SHARED/namelist_ice_ref
r12377 r13472 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 !------------------------------------------------------------------------------ … … 94 95 rn_creepl = 2.0e-9 ! creep limit [1/s] 95 96 rn_ecc = 2.0 ! eccentricity of the elliptical yield curve 96 nn_nevp = 1 20 ! number of EVP subcycles97 nn_nevp = 100 ! number of EVP subcycles 97 98 rn_relast = 0.333 ! ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast 98 ! advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 99 ! advised value: 1/3 (nn_nevp=100) or 1/9 (nn_nevp=300) 100 nn_rhg_chkcvg = 0 ! check convergence of rheology 101 ! = 0 no check 102 ! = 1 check at the main time step (output xml: uice_cvg) 103 ! = 2 check at both main and rheology time steps (additional output: ice_cvg.nc) 104 ! this option 2 asks a lot of communications between cpu 99 105 / 100 106 !------------------------------------------------------------------------------ 101 107 &namdyn_adv ! Ice advection 102 108 !------------------------------------------------------------------------------ 103 ln_adv_Pra = .true. ! Advection scheme (Prather)104 ln_adv_UMx = .false. 109 ln_adv_Pra = .true. ! Advection scheme (Prather) 110 ln_adv_UMx = .false. ! Advection scheme (Ultimate-Macho) 105 111 nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order) 106 112 / … … 109 115 !------------------------------------------------------------------------------ 110 116 rn_cio = 5.0e-03 ! ice-ocean drag coefficient (-) 111 rn_blow_s = 0.66 ! mesure of snow blowing into the leads 117 nn_snwfra = 2 ! calculate the fraction of ice covered by snow (for zdf and albedo) 118 ! = 0 fraction = 1 (if snow) or 0 (if no snow) 119 ! = 1 fraction = 1-exp(-0.2*rhos*hsnw) [MetO formulation] 120 ! = 2 fraction = hsnw / (hsnw+0.02) [CICE formulation] 121 rn_snwblow = 0.66 ! mesure of snow blowing into the leads 112 122 ! = 1 => no snow blowing, < 1 => some snow blowing 113 123 nn_flxdist = -1 ! Redistribute heat flux over ice categories … … 118 128 ln_cndflx = .false. ! Use conduction flux as surface boundary conditions (i.e. for Jules coupling) 119 129 ln_cndemulate = .false. ! emulate conduction flux (if not provided in the inputs) 130 nn_qtrice = 1 ! Solar flux transmitted thru the surface scattering layer: 131 ! = 0 Grenfell and Maykut 1977 (depends on cloudiness and is 0 when there is snow) 132 ! = 1 Lebrun 2019 (equals 0.3 anytime with different melting/dry snw conductivities) 120 133 / 121 134 !------------------------------------------------------------------------------ … … 126 139 ln_icedO = .true. ! activate ice growth in open-water (T) or not (F) 127 140 ln_icedS = .true. ! activate brine drainage (T) or not (F) 141 ! 142 ln_leadhfx = .true. ! heat in the leads is used to melt sea-ice before warming the ocean 128 143 / 129 144 !------------------------------------------------------------------------------ … … 135 150 rn_cnd_s = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971) 136 151 ! Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 137 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice [1/m] 152 rn_kappa_i = 1.0 ! radiation attenuation coefficient in sea ice [1/m] 153 rn_kappa_s = 10.0 ! nn_qtrice = 0: radiation attenuation coefficient in snow [1/m] 154 rn_kappa_smlt = 7.0 ! nn_qtrice = 1: radiation attenuation coefficient in melting snow [1/m] 155 rn_kappa_sdry = 10.0 ! radiation attenuation coefficient in dry snow [1/m] 156 ln_zdf_chkcvg = .false. ! check convergence of heat diffusion scheme (outputs: tice_cvgerr, tice_cvgstp) 138 157 / 139 158 !------------------------------------------------------------------------------ … … 175 194 &namthd_pnd ! Melt ponds 176 195 !------------------------------------------------------------------------------ 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 196 ln_pnd = .true. ! activate melt ponds or not 197 ln_pnd_LEV = .true. ! level ice melt ponds (from Flocco et al 2007,2010 & Holland et al 2012) 198 rn_apnd_min = 0.15 ! minimum ice fraction that contributes to melt pond. range: 0.0 -- 0.15 ?? 199 rn_apnd_max = 0.85 ! maximum ice fraction that contributes to melt pond. range: 0.7 -- 0.85 ?? 200 ln_pnd_CST = .false. ! constant melt ponds 201 rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC 202 rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC 203 ln_pnd_lids = .true. ! frozen lids on top of the ponds (only for ln_pnd_LEV) 204 ln_pnd_alb = .true. ! effect of melt ponds on ice albedo 183 205 / 184 206 !------------------------------------------------------------------------------ … … 186 208 !------------------------------------------------------------------------------ 187 209 ln_iceini = .true. ! activate ice initialization (T) or not (F) 188 ln_iceini_file = .false. ! netcdf file provided for initialization (T) or not (F) 210 nn_iceini_file = 0 ! 0 = Initialise sea ice based on SSTs 211 ! 1 = Initialise sea ice from single category netcdf file 212 ! 2 = Initialise sea ice from multi category restart file 189 213 rn_thres_sst = 2.0 ! max temp. above Tfreeze with initial ice = (sst - tfreeze) 190 214 rn_hti_ini_n = 3.0 ! initial ice thickness (m), North … … 206 230 rn_hpd_ini_n = 0.05 ! initial pond depth (m), North 207 231 rn_hpd_ini_s = 0.05 ! " " South 208 ! -- for ln_iceini_file = T 232 rn_hld_ini_n = 0.0 ! initial pond lid depth (m), North 233 rn_hld_ini_s = 0.0 ! " " South 234 ! -- for nn_iceini_file = 1 209 235 sn_hti = 'Ice_initialization' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' 210 236 sn_hts = 'Ice_initialization' , -12 ,'hts' , .false. , .true., 'yearly' , '' , '', '' … … 217 243 sn_apd = 'NOT USED' , -12 ,'apd' , .false. , .true., 'yearly' , '' , '', '' 218 244 sn_hpd = 'NOT USED' , -12 ,'hpd' , .false. , .true., 'yearly' , '' , '', '' 245 sn_hld = 'NOT USED' , -12 ,'hld' , .false. , .true., 'yearly' , '' , '', '' 219 246 cn_dir='./' 220 247 / … … 238 265 ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F) 239 266 ln_icectl = .false. ! ice points output for debug (T or F) 240 iiceprt = 10 !i-index for debug241 jiceprt = 10 !j-index for debug242 / 267 iiceprt = 10 ! i-index for debug 268 jiceprt = 10 ! j-index for debug 269 / -
NEMO/trunk/cfgs/SHARED/namelist_ref
r13461 r13472 296 296 sn_uoatm = 'NOT USED' , 6. , 'UOATM' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , 'Uoceatm', '' 297 297 sn_voatm = 'NOT USED' , 6. , 'VOATM' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , 'Voceatm', '' 298 sn_cc = 'NOT USED' , 24. , 'CC' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 298 299 sn_hpgi = 'NOT USED' , 24. , 'uhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG' , '' 299 300 sn_hpgj = 'NOT USED' , 24. , 'vhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG' , '' … … 335 336 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 336 337 !----------------------------------------------------------------------- 337 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data 338 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 339 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 340 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 338 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data 339 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 340 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 341 ln_scale_ice_flux = .false. ! use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 342 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 341 343 !_____________!__________________________!____________!_____________!______________________!________! 342 344 ! ! description ! multiple ! vector ! vector ! vector ! … … 732 734 bn_aip = 'NOT USED' , 24. , 'siapnd' , .true. , .false., 'daily' , '' , '' , '' 733 735 bn_hip = 'NOT USED' , 24. , 'sihpnd' , .true. , .false., 'daily' , '' , '' , '' 736 bn_hil = 'NOT USED' , 24. , 'sihlid' , .true. , .false., 'daily' , '' , '' , '' 734 737 ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 735 738 rn_ice_tem = 270. ! arbitrary temperature of incoming sea ice … … 738 741 rn_ice_apnd = 0.2 ! -- pond fraction = a_ip/a_i -- 739 742 rn_ice_hpnd = 0.05 ! -- pond depth -- 743 rn_ice_hlid = 0.0 ! -- pond lid depth -- 740 744 / 741 745 !----------------------------------------------------------------------- … … 765 769 ! 766 770 ln_drgimp = .true. ! implicit top/bottom friction flag 771 ln_drgice_imp = .false. ! implicit ice-ocean drag 767 772 / 768 773 !----------------------------------------------------------------------- … … 1135 1140 rn_bshear = 1.e-20 ! background shear (>0) currently a numerical threshold (do not change it) 1136 1141 nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 1137 nn_mxl = 2! mixing length: = 0 bounded by the distance to surface and bottom1142 nn_mxl = 3 ! mixing length: = 0 bounded by the distance to surface and bottom 1138 1143 ! ! = 1 bounded by the local vertical scale factor 1139 1144 ! ! = 2 first vertical derivative of mixing length bounded by 1 1140 1145 ! ! = 3 as =2 with distinct dissipative an mixing length scale 1141 1146 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) 1142 nn_mxlice = 0! type of scaling under sea-ice1147 nn_mxlice = 2 ! type of scaling under sea-ice 1143 1148 ! = 0 no scaling under sea-ice 1144 1149 ! = 1 scaling with constant sea-ice thickness 1145 ! = 2 1146 ! = 3 1150 ! = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 1151 ! = 3 scaling with maximum sea-ice thickness 1147 1152 rn_mxlice = 10. ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 1148 1153 rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value … … 1157 1162 ! = 0 constant 10 m length scale 1158 1163 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 1159 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1164 nn_eice = 1 ! attenutaion of langmuir & surface wave breaking under ice 1165 ! ! = 0 no impact of ice cover on langmuir & surface wave breaking 1166 ! ! = 1 weigthed by 1-TANH(10*fr_i) 1167 ! ! = 2 weighted by 1-fr_i 1168 ! ! = 3 weighted by 1-MIN(1,4*fr_i) 1160 1169 / 1161 1170 !----------------------------------------------------------------------- … … 1170 1179 rn_charn = 70000. ! Charnock constant for wb induced roughness length 1171 1180 rn_hsro = 0.02 ! Minimum surface roughness 1181 rn_hsri = 0.03 ! Ice-ocean roughness 1172 1182 rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met>1) 1173 1183 nn_z0_met = 2 ! Method for surface roughness computation (0/1/2/3) 1174 ! ! =3 requires ln_wave=T 1184 ! ! = 3 requires ln_wave=T 1185 nn_z0_ice = 1 ! attenutaion of surface wave breaking under ice 1186 ! ! = 0 no impact of ice cover 1187 ! ! = 1 roughness uses rn_hsri and is weigthed by 1-TANH(10*fr_i) 1188 ! ! = 2 roughness uses rn_hsri and is weighted by 1-fr_i 1189 ! ! = 3 roughness uses rn_hsri and is weighted by 1-MIN(1,4*fr_i) 1175 1190 nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum) 1176 1191 nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum)
Note: See TracChangeset
for help on using the changeset viewer.