| 1 | |
|---|
| 2 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 3 | !! NEMO/OPA Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref |
|---|
| 4 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 5 | ! |
|---|
| 6 | !----------------------------------------------------------------------- |
|---|
| 7 | &namrun ! parameters of the run |
|---|
| 8 | !----------------------------------------------------------------------- |
|---|
| 9 | cn_exp = _AUTOBLOCKER_ ! Experience name |
|---|
| 10 | nn_it000 = _AUTOBLOCKER_ ! First time step |
|---|
| 11 | nn_itend = _AUTOBLOCKER_ ! Last time step |
|---|
| 12 | nn_date0 = _AUTOBLOCKER_ ! Date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) |
|---|
| 13 | nn_leapy = _AUTOBLOCKER_ ! Leap year calendar (1) or not (0) |
|---|
| 14 | ln_rstart = _AUTOBLOCKER_ ! start from rest (F) or from a restart file (T) |
|---|
| 15 | ln_rstart_ts = _AUTOBLOCKER_ ! start from rest for current only (F) or from a restart file (T) |
|---|
| 16 | nn_rstctl = _AUTOBLOCKER_ ! Restart control => activated only if ln_rstart = T |
|---|
| 17 | ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist |
|---|
| 18 | ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart |
|---|
| 19 | ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart |
|---|
| 20 | cn_ocerst_in = "restartopa" ! Suffix of ocean restart name (input) |
|---|
| 21 | cn_ocerst_indir = "." ! directory from which to read input ocean restarts |
|---|
| 22 | cn_ocerst_out = "restart" ! Suffix of ocean restart name (output) |
|---|
| 23 | cn_ocerst_outdir = "." ! directory in which to write output ocean restarts |
|---|
| 24 | nn_istate = 0 ! Output the initial state (1) or not (0) |
|---|
| 25 | nn_stock = _AUTOBLOCKER_ ! Frequency of creation of a restart file (modulo referenced to 1) |
|---|
| 26 | nn_write = 5475 ! Requency of write in the output file (modulo referenced to nn_it000) |
|---|
| 27 | ln_mskland = .true. ! Masks land points in NetCDF outputs |
|---|
| 28 | ln_mskutil = .true. ! Outputs without halos |
|---|
| 29 | ln_cfmeta = .true. ! output additional data to netCDF files required for compliance with the CF metadata standard |
|---|
| 30 | / |
|---|
| 31 | !----------------------------------------------------------------------- |
|---|
| 32 | &namcfg ! parameters of the configuration |
|---|
| 33 | !----------------------------------------------------------------------- |
|---|
| 34 | cp_cfg = "orca" ! name of the configuration |
|---|
| 35 | jp_cfg = 1 ! resolution of the configuration |
|---|
| 36 | jpidta = 362 ! 1st lateral dimension ( >= jpi ) |
|---|
| 37 | jpjdta = 332 ! 2nd " " ( >= jpj ) |
|---|
| 38 | jpkdta = 75 ! number of levels ( >= jpk ) |
|---|
| 39 | jpiglo = 362 ! 1st dimension of global domain --> i =jpidta |
|---|
| 40 | jpjglo = 332 ! 2nd - - --> j =jpjdta |
|---|
| 41 | jperio = 6 ! lateral cond. type (between 0 and 6) |
|---|
| 42 | / |
|---|
| 43 | !----------------------------------------------------------------------- |
|---|
| 44 | &namzgr ! vertical coordinate |
|---|
| 45 | !----------------------------------------------------------------------- |
|---|
| 46 | / |
|---|
| 47 | !----------------------------------------------------------------------- |
|---|
| 48 | &namzgr_sco ! s-coordinate or hybrid z-s-coordinate |
|---|
| 49 | !----------------------------------------------------------------------- |
|---|
| 50 | / |
|---|
| 51 | !----------------------------------------------------------------------- |
|---|
| 52 | &namdom ! space and time domain (bathymetry, mesh, timestep) |
|---|
| 53 | !----------------------------------------------------------------------- |
|---|
| 54 | nn_closea = 1 ! remove (=0) or keep (=1) closed seas and lakes (ORCA) |
|---|
| 55 | ln_clodyn = _AUTOBLOCKER_ ! |
|---|
| 56 | ! |
|---|
| 57 | jphgr_msh = 0 ! type of horizontal mesh |
|---|
| 58 | ppglam0 = 999999.0 ! longitude of first raw and column T-point (jphgr_msh = 1) |
|---|
| 59 | ppgphi0 = 999999.0 ! latitude of first raw and column T-point (jphgr_msh = 1) |
|---|
| 60 | ppe1_deg = 999999.0 ! zonal grid-spacing (degrees) |
|---|
| 61 | ppe2_deg = 999999.0 ! meridional grid-spacing (degrees) |
|---|
| 62 | ppe1_m = 999999.0 ! zonal grid-spacing (degrees) |
|---|
| 63 | ppe2_m = 999999.0 ! meridional grid-spacing (degrees) |
|---|
| 64 | ppsur = -3958.951371276829 ! ORCA r4, r2 and r05 coefficients |
|---|
| 65 | ppa0 = 103.9530096000000 ! (default coefficients) |
|---|
| 66 | ppa1 = 2.415951269000000 ! |
|---|
| 67 | ppkth = 15.35101370000000 ! |
|---|
| 68 | ppacr = 7.0 ! |
|---|
| 69 | ppdzmin = 999999.0 ! Minimum vertical spacing |
|---|
| 70 | pphmax = 999999.0 ! Maximum depth |
|---|
| 71 | ppa2 = 100.7609285000000 ! Double tanh function parameters |
|---|
| 72 | ppkth2 = 48.02989372000000 ! |
|---|
| 73 | ppacr2 = 13.00000000000 ! |
|---|
| 74 | rn_rdt = 2700. ! time step for the dynamics (and tracer if nn_acc=0) |
|---|
| 75 | rn_hmin = 20. |
|---|
| 76 | nn_msh = _AUTO_ ! AUTO - Create (=1) a mesh file or not (=0) |
|---|
| 77 | / |
|---|
| 78 | !----------------------------------------------------------------------- |
|---|
| 79 | &namsplit |
|---|
| 80 | !----------------------------------------------------------------------- |
|---|
| 81 | ln_bt_fw = .FALSE. ! leap-frog integration of barotropic equations |
|---|
| 82 | ln_bt_av = .TRUE. ! Time filtering of barotropic variables |
|---|
| 83 | ln_bt_nn_auto = .TRUE. ! Set nn_baro automatically to be just below |
|---|
| 84 | ! a user defined maximum courant number (rn_bt_cmax) |
|---|
| 85 | nn_baro = 30 ! Number of iterations of barotropic mode |
|---|
| 86 | ! during rn_rdt seconds. Only used if ln_bt_nn_auto=F |
|---|
| 87 | rn_bt_cmax = 0.8 ! Maximum courant number allowed if ln_bt_nn_auto=T |
|---|
| 88 | nn_bt_flt = 1 ! Time filter choice |
|---|
| 89 | ! = 0 None |
|---|
| 90 | ! = 1 Boxcar over nn_baro barotropic steps |
|---|
| 91 | ! = 2 Boxcar over 2*nn_baro " |
|---|
| 92 | / |
|---|
| 93 | !----------------------------------------------------------------------- |
|---|
| 94 | &namcrs ! Grid coarsening for dynamics output and/or |
|---|
| 95 | ! passive tracer coarsened online simulations |
|---|
| 96 | !----------------------------------------------------------------------- |
|---|
| 97 | / |
|---|
| 98 | !----------------------------------------------------------------------- |
|---|
| 99 | &namtsd ! data : Temperature & Salinity |
|---|
| 100 | !----------------------------------------------------------------------- |
|---|
| 101 | ln_tsd_tradmp = .false. ! damping of ocean T & S toward T &S input data (T) or not (F) |
|---|
| 102 | sn_tem = 'conservative_temperature_WOA13_decav_Reg1L75_clim', -1 ,'votemper' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc' , '' , '' |
|---|
| 103 | sn_sal = 'absolute_salinity_WOA13_decav_Reg1L75_clim' , -1 ,'vosaline' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc' , '' , '' |
|---|
| 104 | / |
|---|
| 105 | !----------------------------------------------------------------------- |
|---|
| 106 | &namsbc ! Surface Boundary Condition (surface module) |
|---|
| 107 | !----------------------------------------------------------------------- |
|---|
| 108 | nn_fsbc = 2 ! frequency of surface boundary condition computation |
|---|
| 109 | ! (also = the frequency of sea-ice model call) |
|---|
| 110 | ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) |
|---|
| 111 | ln_cpl = .true. ! atmosphere coupled formulation ( requires key_oasis3 ) |
|---|
| 112 | nn_limflx = 2 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) |
|---|
| 113 | ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled |
|---|
| 114 | ! = 0 Average per-category fluxes (forced and coupled mode) |
|---|
| 115 | ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled |
|---|
| 116 | ! = 2 Redistribute a single flux over categories (coupled mode only) |
|---|
| 117 | nn_ice_embd = 1 ! AUTO - |
|---|
| 118 | ! =0 levitating ice (no mass exchange, concentration/dilution effect) |
|---|
| 119 | ! =1 levitating ice with mass and salt exchange but no presure effect |
|---|
| 120 | ! =2 embedded sea-ice (full salt and mass exchanges and pressure) |
|---|
| 121 | ln_rnf = .false. ! runoffs (T => fill namsbc_rnf) |
|---|
| 122 | ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) |
|---|
| 123 | nn_fwb = 0 ! FreshWater Budget: =0 unchecked |
|---|
| 124 | nn_isf = 3 ! ice shelf melting/freezing (/=0 => fill namsbc_isf) |
|---|
| 125 | ! 3 = rnf file for isf |
|---|
| 126 | !----------------------------------------------------------------------- |
|---|
| 127 | &namsbc_core ! namsbc_core CORE bulk formulae |
|---|
| 128 | !----------------------------------------------------------------------- |
|---|
| 129 | / |
|---|
| 130 | !----------------------------------------------------------------------- |
|---|
| 131 | &namtra_qsr ! penetrative solar radiation |
|---|
| 132 | !----------------------------------------------------------------------- |
|---|
| 133 | sn_chl ='merged_ESACCI_BIOMER4V1R1_CHL_REG05', -1 , 'CHLA' , .true. , .true. , 'yearly' , 'weights_reg05_2_eorca1_bilinear.nc' , '' , '' |
|---|
| 134 | ln_traqsr = .true. ! Light penetration (T) or not (F) |
|---|
| 135 | ln_qsr_rgb = .false. ! RGB (Red-Green-Blue) light penetration |
|---|
| 136 | ln_qsr_2bd = .false. ! 2 bands light penetration |
|---|
| 137 | ln_qsr_bio = .true. ! bio-model light penetration |
|---|
| 138 | / |
|---|
| 139 | !----------------------------------------------------------------------- |
|---|
| 140 | &namsbc_rnf ! runoffs namelist surface boundary condition |
|---|
| 141 | !----------------------------------------------------------------------- |
|---|
| 142 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 143 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 144 | sn_rnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 145 | sn_i_rnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', -1 , 'Icb_flux', .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 146 | sn_cnf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc', 0 , 'socoeff' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 147 | sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 148 | sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 149 | sn_dep_rnf = 'runoffs_eORCA1.0_depths.nc' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 150 | |
|---|
| 151 | ln_rnf_icb = .false. ! read in iceberg flux |
|---|
| 152 | ln_rnf_mouth = .false. ! specific treatment at rivers mouths |
|---|
| 153 | ln_rnf_depth = .true. ! read in depth information for runoff |
|---|
| 154 | ln_rnf_tem = .false. ! read in temperature information for runoff |
|---|
| 155 | ln_rnf_sal = .false. ! read in salinity information for runoff |
|---|
| 156 | ln_rnf_depth_ini = .false.! compute depth at initialisation from runoff file |
|---|
| 157 | rn_rnf_max = 0.05 ! max value of the runoff climatology over global domain ( if ln_rnf_depth_ini = .true ) |
|---|
| 158 | rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) |
|---|
| 159 | nn_rnf_depth_file = _AUTO_ ! create (=1) a runoff depth file or not (=0) |
|---|
| 160 | / |
|---|
| 161 | !----------------------------------------------------------------------- |
|---|
| 162 | &namsbc_isf ! Top boundary layer (ISF) |
|---|
| 163 | !----------------------------------------------------------------------- |
|---|
| 164 | ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! |
|---|
| 165 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
|---|
| 166 | ! ! |
|---|
| 167 | sn_rnfisf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sornfisf', .false. , .true. , 'yearly' , '' , '' |
|---|
| 168 | sn_depmax_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sodepmax_isf' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 169 | sn_depmin_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' , -12 ,'sodepmin_isf' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 170 | / |
|---|
| 171 | !----------------------------------------------------------------------- |
|---|
| 172 | &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk |
|---|
| 173 | !----------------------------------------------------------------------- |
|---|
| 174 | / |
|---|
| 175 | !----------------------------------------------------------------------- |
|---|
| 176 | &namsbc_ssr ! surface boundary condition : sea surface restoring |
|---|
| 177 | !----------------------------------------------------------------------- |
|---|
| 178 | ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! |
|---|
| 179 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
|---|
| 180 | sn_sst = 'sst_conservative_temperature_WOA13_decav_Reg1L75_clim', -1 ,'sosstsst' , .true. , .true. , 'yearly' , 'weights_WOA13d1_2_eorca1_bilinear.nc' , '' , '' |
|---|
| 181 | sn_sss = 'sss_absolute_salinity_WOA13_decav_Reg1L75_clim', -1., 'sosaline', .true. , .true. , 'yearly' , 'weights_WOA13d1_2_eorca1_bilinear.nc' , '' |
|---|
| 182 | |
|---|
| 183 | cn_dir = './' ! root directory for the location of the runoff files |
|---|
| 184 | nn_sstr = 0 ! add a retroaction term in the surface heat not (=0) or flux (=1) or (conditional nudging =2) |
|---|
| 185 | nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2) or (conditional nudging =3) |
|---|
| 186 | ! or to SSS only (=1) or no damping term (=0) |
|---|
| 187 | rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] |
|---|
| 188 | rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] |
|---|
| 189 | ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) |
|---|
| 190 | rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] |
|---|
| 191 | nn_icedmp = 0 ! Cntrl of surface restoration under ice nn_icedmp |
|---|
| 192 | ! ( 0 = no restoration under ice ) |
|---|
| 193 | ! ( 1 = restoration everywhere ) |
|---|
| 194 | ! ( > 1 = reinforced damping (x nn_icedmp) under ice ) |
|---|
| 195 | ln_sssd_bnd = .false. ! flag to bound S-S* |
|---|
| 196 | rn_sssd_bnd = 0.01 ! ABS(Max./Min.) S-S* threshold [psu] |
|---|
| 197 | ! |
|---|
| 198 | ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 199 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 200 | sn_msk_sst = 'eORCA1.2_mask_SST_Merged_good' , -1 , 'mask_part' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 201 | sn_msk_sss = 'eORCA1.2_mask_SSS_FDS_yearly' , -12 , 'mask_part' , .false. , .true. , 'yearly' , '' , '' , '' |
|---|
| 202 | ln_vargamma = .true. ! false = fix gamma (value=rn_dqdt), true= variable gamma (proportional to mld) |
|---|
| 203 | nn_mask_sst = 2 ! (associated with nn_sstr=2) : 0 = global |
|---|
| 204 | ! : 1 = fix (regional mask) in time: one yearly value, read in sn_msk_sst |
|---|
| 205 | ! : 2 = variable (regional mask) in time: read different value every month, read in sn_msk_sst |
|---|
| 206 | nn_mask_sss = 1 ! (associated ONLY with nn_sssr=3) : 0 = global |
|---|
| 207 | / |
|---|
| 208 | !----------------------------------------------------------------------- |
|---|
| 209 | &namsbc_alb ! albedo parameters |
|---|
| 210 | !----------------------------------------------------------------------- |
|---|
| 211 | nn_ice_alb = 1 ! parameterization of ice/snow albedo |
|---|
| 212 | ! 0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo |
|---|
| 213 | ! 1: "home made" based on Brandt et al. (JClim 2005) and Grenfell & Perovich (JGR 2004), |
|---|
| 214 | ! giving cloud-sky albedo |
|---|
| 215 | rn_alb_sdry = 0.87 ! dry snow albedo : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky) |
|---|
| 216 | rn_alb_smlt = 0.82 ! melting snow albedo : 0.65 ( '' ) ; 0.75 ( '' ) ; obs 0.72-0.82 ( '' ) |
|---|
| 217 | rn_alb_idry = 0.65 ! dry ice albedo : 0.72 ( '' ) ; 0.60 ( '' ) ; obs 0.54-0.65 ( '' ) |
|---|
| 218 | rn_alb_imlt = 0.58 ! bare puddled ice albedo : 0.53 ( '' ) ; 0.50 ( '' ) ; obs 0.49-0.58 ( '' ) |
|---|
| 219 | / |
|---|
| 220 | !----------------------------------------------------------------------- |
|---|
| 221 | &namsbc_cpl ! coupling parameters |
|---|
| 222 | !----------------------------------------------------------------------- |
|---|
| 223 | ! ! description ! multiple ! vector ! vector ! vector ! |
|---|
| 224 | ! ! ! categories ! reference ! orientation ! grids ! |
|---|
| 225 | ! send |
|---|
| 226 | sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' |
|---|
| 227 | sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' |
|---|
| 228 | sn_snd_thick = 'none' , 'no' , '' , '' , '' |
|---|
| 229 | sn_snd_crt = 'mixed oce-ice' , 'no' , 'cartesian' , 'eastward-northward' , 'T' |
|---|
| 230 | sn_snd_co2 = 'none' , 'no' , '' , '' , '' |
|---|
| 231 | ! receive |
|---|
| 232 | sn_rcv_w10m = 'coupled' , 'no' , '' , '' , '' |
|---|
| 233 | sn_rcv_taumod = 'none' , 'no' , '' , '' , '' |
|---|
| 234 | sn_rcv_tau = 'mixed oce-ice' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' |
|---|
| 235 | sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' |
|---|
| 236 | sn_rcv_qsr = 'conservative' , 'no' , '' , '' , '' |
|---|
| 237 | sn_rcv_qns = 'conservative' , 'no' , '' , '' , '' |
|---|
| 238 | sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' |
|---|
| 239 | sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' |
|---|
| 240 | sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' |
|---|
| 241 | sn_rcv_co2 = 'none' , 'no' , '' , '' , '' |
|---|
| 242 | sn_rcv_icb = 'coupled' , 'no' , '' , '' , '' |
|---|
| 243 | sn_rcv_isf = 'coupled' , 'no' , '' , '' , '' |
|---|
| 244 | |
|---|
| 245 | / |
|---|
| 246 | !----------------------------------------------------------------------- |
|---|
| 247 | &namberg ! iceberg parameters |
|---|
| 248 | !----------------------------------------------------------------------- |
|---|
| 249 | ln_icebergs = .false. |
|---|
| 250 | ln_bergdia = .false. ! Calculate budgets |
|---|
| 251 | nn_verbose_level = 0 ! Turn on more verbose output if level > 0 |
|---|
| 252 | nn_verbose_write = 120 ! Timesteps between verbose messages |
|---|
| 253 | nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage |
|---|
| 254 | ! Initial mass required for an iceberg of each class |
|---|
| 255 | rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 |
|---|
| 256 | ! Proportion of calving mass to apportion to each class |
|---|
| 257 | rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 |
|---|
| 258 | ! Ratio between effective and real iceberg mass (non-dim) |
|---|
| 259 | ! i.e. number of icebergs represented at a point |
|---|
| 260 | rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 |
|---|
| 261 | ! thickness of newly calved bergs (m) |
|---|
| 262 | rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. |
|---|
| 263 | rn_rho_bergs = 850. ! Density of icebergs |
|---|
| 264 | rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs |
|---|
| 265 | ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics |
|---|
| 266 | rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits |
|---|
| 267 | rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) |
|---|
| 268 | ln_passive_mode = .false. ! iceberg - ocean decoupling |
|---|
| 269 | nn_test_icebergs = 8 ! Create test icebergs of this class (-1 = no) |
|---|
| 270 | ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) |
|---|
| 271 | !rn_test_box = 108.0, 116.0, -66.0, -58.0 |
|---|
| 272 | rn_test_box = -180.0, 180.0, 70.0, 90.0 ! |
|---|
| 273 | rn_speed_limit = 0. ! CFL speed limit for a berg |
|---|
| 274 | |
|---|
| 275 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
|---|
| 276 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! |
|---|
| 277 | sn_icb = 'calving' , -1 , 'calvingmask', .true. , .true. , 'yearly' , '' , '' , '' |
|---|
| 278 | |
|---|
| 279 | cn_dir = './' |
|---|
| 280 | / |
|---|
| 281 | !----------------------------------------------------------------------- |
|---|
| 282 | &namlbc ! lateral momentum boundary condition |
|---|
| 283 | !----------------------------------------------------------------------- |
|---|
| 284 | rn_shlat = 0.0 ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat |
|---|
| 285 | / |
|---|
| 286 | !----------------------------------------------------------------------- |
|---|
| 287 | &namcla ! cross land advection |
|---|
| 288 | !----------------------------------------------------------------------- |
|---|
| 289 | / |
|---|
| 290 | !----------------------------------------------------------------------- |
|---|
| 291 | &nambfr ! bottom friction |
|---|
| 292 | !----------------------------------------------------------------------- |
|---|
| 293 | nn_bfr = 2 ! type of bottom friction : = 0 : free slip, = 1 : linear friction |
|---|
| 294 | / |
|---|
| 295 | !----------------------------------------------------------------------- |
|---|
| 296 | &nambbc ! bottom temperature boundary condition |
|---|
| 297 | !----------------------------------------------------------------------- |
|---|
| 298 | sn_qgh ='Goutorbe_ghflux.nc', -12. , 'gh_flux' , .false. , .true. , 'yearly' , 'weights_Goutorbe1_2_eorca1_bilinear.nc' , '' , '' |
|---|
| 299 | ! |
|---|
| 300 | cn_dir = './' ! root directory for the location of the runoff files |
|---|
| 301 | nn_geoflx = 2 ! geothermal heat flux: = 0 no flux |
|---|
| 302 | / |
|---|
| 303 | !----------------------------------------------------------------------- |
|---|
| 304 | &nambbl ! bottom boundary layer scheme |
|---|
| 305 | !----------------------------------------------------------------------- |
|---|
| 306 | / |
|---|
| 307 | !----------------------------------------------------------------------- |
|---|
| 308 | &nameos ! ocean physical parameters |
|---|
| 309 | !----------------------------------------------------------------------- |
|---|
| 310 | / |
|---|
| 311 | !----------------------------------------------------------------------- |
|---|
| 312 | &namtra_adv ! advection scheme for tracer |
|---|
| 313 | !----------------------------------------------------------------------- |
|---|
| 314 | ln_traadv_tvd = .true. ! TVD scheme |
|---|
| 315 | ln_traadv_ubs = .false. ! UBS scheme |
|---|
| 316 | / |
|---|
| 317 | !----------------------------------------------------------------------- |
|---|
| 318 | &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) |
|---|
| 319 | !----------------------------------------------------------------------- |
|---|
| 320 | / |
|---|
| 321 | !---------------------------------------------------------------------------------- |
|---|
| 322 | &namtra_ldf ! lateral diffusion scheme for tracers |
|---|
| 323 | !---------------------------------------------------------------------------------- |
|---|
| 324 | ln_traldf_grif = .false. ! griffies skew flux formulation (require "key_ldfslp") |
|---|
| 325 | ln_traldf_gdia = .false. ! griffies operator strfn diagnostics (require "key_ldfslp") |
|---|
| 326 | ln_botmix_grif = .false. ! griffies operator with lateral mixing on bottom (require "key_ldfslp") |
|---|
| 327 | rn_aht_0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] |
|---|
| 328 | rn_aeiv_0 = 1000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") |
|---|
| 329 | / |
|---|
| 330 | !----------------------------------------------------------------------- |
|---|
| 331 | &namtra_dmp ! tracer: T & S newtonian damping |
|---|
| 332 | !----------------------------------------------------------------------- |
|---|
| 333 | ln_tradmp = .false. ! add a damping termn (T) or not (F) |
|---|
| 334 | / |
|---|
| 335 | !----------------------------------------------------------------------- |
|---|
| 336 | &namdyn_adv ! formulation of the momentum advection |
|---|
| 337 | !----------------------------------------------------------------------- |
|---|
| 338 | ln_dynadv_vec = .true. ! vector form (T) or flux form (F) |
|---|
| 339 | ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme |
|---|
| 340 | ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme |
|---|
| 341 | nn_dynkeg = 1 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction |
|---|
| 342 | / |
|---|
| 343 | !----------------------------------------------------------------------- |
|---|
| 344 | &nam_vvl ! vertical coordinate options |
|---|
| 345 | !----------------------------------------------------------------------- |
|---|
| 346 | / |
|---|
| 347 | !----------------------------------------------------------------------- |
|---|
| 348 | &namdyn_vor ! option of physics/algorithm (not control by CPP keys) |
|---|
| 349 | !----------------------------------------------------------------------- |
|---|
| 350 | / |
|---|
| 351 | !----------------------------------------------------------------------- |
|---|
| 352 | &namdyn_hpg ! Hydrostatic pressure gradient option |
|---|
| 353 | !----------------------------------------------------------------------- |
|---|
| 354 | ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) |
|---|
| 355 | ln_hpg_sco = .true. ! s-coordinate (standard jacobian formulation) |
|---|
| 356 | !ln_hpg_isf = .true. ! s-coordinate (sco ) adapted to isf |
|---|
| 357 | ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) |
|---|
| 358 | ! centered time scheme (F) |
|---|
| 359 | / |
|---|
| 360 | !----------------------------------------------------------------------- |
|---|
| 361 | &namdyn_ldf ! lateral diffusion on momentum |
|---|
| 362 | !----------------------------------------------------------------------- |
|---|
| 363 | rn_ahm_0_lap = 20000. ! horizontal laplacian eddy viscosity [m2/s] |
|---|
| 364 | / |
|---|
| 365 | !----------------------------------------------------------------------- |
|---|
| 366 | &namzdf ! vertical physics |
|---|
| 367 | !----------------------------------------------------------------------- |
|---|
| 368 | / |
|---|
| 369 | !----------------------------------------------------------------------- |
|---|
| 370 | &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") |
|---|
| 371 | !----------------------------------------------------------------------- |
|---|
| 372 | nn_etau = 0 ! penetration of tke below the mixed layer (ML) due to internal & intertial waves |
|---|
| 373 | ! = 0 no penetration |
|---|
| 374 | ! = 1 add a tke source below the ML |
|---|
| 375 | ! = 2 add a tke source just at the base of the ML |
|---|
| 376 | ! = 3 as = 1 applied on HF part of the stress ("key_oasis3") |
|---|
| 377 | nn_mxl0 = 2 ! type of scaling under sea-ice |
|---|
| 378 | ! = 0 no scaling under sea-ice |
|---|
| 379 | ! = 1 scaling with constant sea-ice thickness |
|---|
| 380 | ! = 2 scaling with mean sea-ice thickness |
|---|
| 381 | ! = 3 scaling with maximum sea-ice thickness |
|---|
| 382 | rn_hice = 10. ! max constant ice thickness value when scaling under sea-ice ( nn_mxl0=1) |
|---|
| 383 | ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002) |
|---|
| 384 | rn_lc = 0.20 ! coef. associated to Langmuir cells |
|---|
| 385 | / |
|---|
| 386 | !----------------------------------------------------------------------- |
|---|
| 387 | &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") |
|---|
| 388 | !----------------------------------------------------------------------- |
|---|
| 389 | / |
|---|
| 390 | !----------------------------------------------------------------------- |
|---|
| 391 | &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") |
|---|
| 392 | !----------------------------------------------------------------------- |
|---|
| 393 | / |
|---|
| 394 | !----------------------------------------------------------------------- |
|---|
| 395 | &namzdf_tmx_new ! new tidal mixing parameterization ("key_zdftmx_new") |
|---|
| 396 | !----------------------------------------------------------------------- |
|---|
| 397 | nn_zpyc = 2 ! pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) |
|---|
| 398 | ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency |
|---|
| 399 | ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F) |
|---|
| 400 | / |
|---|
| 401 | !----------------------------------------------------------------------- |
|---|
| 402 | &namsol ! elliptic solver / island / free surface |
|---|
| 403 | !----------------------------------------------------------------------- |
|---|
| 404 | / |
|---|
| 405 | !----------------------------------------------------------------------- |
|---|
| 406 | &nammpp ! Massively Parallel Processing ("key_mpp_mpi) |
|---|
| 407 | !----------------------------------------------------------------------- |
|---|
| 408 | ln_nnogather= .true. ! |
|---|
| 409 | jpni = 22 ! jpni number of processors following i (set automatically if < 1) |
|---|
| 410 | jpnj = 22 ! jpnj number of processors following j (set automatically if < 1) |
|---|
| 411 | jpnij = 360 ! jpnij number of local domains (set automatically if < 1), 360 for eORCA1/IPSLCM6-LR |
|---|
| 412 | / |
|---|
| 413 | !----------------------------------------------------------------------- |
|---|
| 414 | &namctl ! Control prints & Benchmark |
|---|
| 415 | !----------------------------------------------------------------------- |
|---|
| 416 | / |
|---|
| 417 | !----------------------------------------------------------------------- |
|---|
| 418 | &namptr ! Poleward Transport Diagnostic |
|---|
| 419 | !----------------------------------------------------------------------- |
|---|
| 420 | ln_diaptr = .true. ! Poleward heat and salt transport (T) or not (F) |
|---|
| 421 | ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not |
|---|
| 422 | ! (orca configuration only, need input basins mask file named "subbasins.nc" |
|---|
| 423 | / |
|---|
| 424 | !----------------------------------------------------------------------- |
|---|
| 425 | &namhsb ! Heat and salt budgets |
|---|
| 426 | !----------------------------------------------------------------------- |
|---|
| 427 | ln_diahsb = .true. |
|---|
| 428 | / |
|---|
| 429 | !----------------------------------------------------------------------- |
|---|
| 430 | &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) |
|---|
| 431 | !----------------------------------------------------------------------- |
|---|
| 432 | / |
|---|
| 433 | !----------------------------------------------------------------------- |
|---|
| 434 | &nam_vvl ! vertical coordinate options |
|---|
| 435 | !----------------------------------------------------------------------- |
|---|
| 436 | / |
|---|
| 437 | !----------------------------------------------------------------------- |
|---|
| 438 | &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") |
|---|
| 439 | !----------------------------------------------------------------------- |
|---|
| 440 | / |
|---|
| 441 | !----------------------------------------------------------------------- |
|---|
| 442 | &namtrd ! diagnostics on dynamics and/or tracer trends |
|---|
| 443 | ! ! and/or mixed-layer trends and/or barotropic vorticity |
|---|
| 444 | !----------------------------------------------------------------------- |
|---|
| 445 | ln_tra_trd = .true. ! (T) 3D tracer trend output |
|---|
| 446 | / |
|---|
| 447 | !----------------------------------------------------------------------- |
|---|
| 448 | &namsto ! Stochastic parametrization of EOS |
|---|
| 449 | !----------------------------------------------------------------------- |
|---|
| 450 | / |
|---|
