- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/CONFIG/SHARED/namelist_ref
- Property svn:keywords set to Id
r5501 r6808 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! namelist_ref 1 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 4 !! NEMO/OPA : 1 - run manager (namrun) … … 4 6 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 5 7 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 6 !! namsbc_apr, namsbc_ssr, namsbc_alb )7 !! 4 - lateral boundary (namlbc, nam cla, namobc, namagrif, nambdy, nambdy_tide)8 !! namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave) 9 !! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 8 10 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 9 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ dmp)11 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 10 12 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 11 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ kpp, namzdf_ddm, namzdf_tmx)13 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 12 14 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto) 13 !! 10 - miscellaneous (nam sol, nammpp, namctl)15 !! 10 - miscellaneous (nammpp, namctl) 14 16 !! 11 - Obs & Assim (namobs, nam_asminc) 15 17 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 29 31 nn_itend = 5475 ! last time step (std 5475) 30 32 nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 33 nn_time0 = 0 ! initial time of day in hhmm 31 34 nn_leapy = 0 ! Leap year calendar (1) or not (0) 32 35 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 33 nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T 34 nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T 35 ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 36 ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 37 ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart 38 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 39 cn_ocerst_indir = "." ! directory from which to read input ocean restarts 40 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 41 cn_ocerst_outdir = "." ! directory in which to write output ocean restarts 36 nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T 37 nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T 38 ! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 39 ! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 40 ! ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart 41 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 42 cn_ocerst_indir = "." ! directory from which to read input ocean restarts 43 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 44 cn_ocerst_outdir= "." ! directory in which to write output ocean restarts 45 ln_iscpl = .false. ! cavity evolution forcing or coupling to ice sheet model 42 46 nn_istate = 0 ! output the initial state (1) or not (0) 43 47 ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) … … 45 49 nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 46 50 nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000) 47 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T)48 51 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 49 52 ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard 50 ln_clobber = . false.! clobber (overwrite) an existing file53 ln_clobber = .true. ! clobber (overwrite) an existing file 51 54 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 52 55 / … … 59 62 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 60 63 !! namdom space and time domain (bathymetry, mesh, timestep) 64 !! namcrs coarsened grid (for outputs and/or TOP) ("key_crs") 65 !! namc1d 1D configuration options ("key_c1d") 66 !! namc1d_uvd 1D data (currents) ("key_c1d") 67 !! namc1d_dyndmp 1D newtonian damping applied on currents ("key_c1d") 61 68 !! namtsd data: temperature & salinity 62 69 !!====================================================================== 63 70 ! 64 71 !----------------------------------------------------------------------- 65 &namcfg ! parameters of the configuration 66 !----------------------------------------------------------------------- 67 cp_cfg = "default" ! name of the configuration 68 cp_cfz = "no zoom" ! name of the zoom of configuration 69 jp_cfg = 0 ! resolution of the configuration 70 jpidta = 10 ! 1st lateral dimension ( >= jpi ) 71 jpjdta = 12 ! 2nd " " ( >= jpj ) 72 jpkdta = 31 ! number of levels ( >= jpk ) 73 jpiglo = 10 ! 1st dimension of global domain --> i =jpidta 74 jpjglo = 12 ! 2nd - - --> j =jpjdta 75 jpizoom = 1 ! left bottom (i,j) indices of the zoom 76 jpjzoom = 1 ! in data domain indices 77 jperio = 0 ! lateral cond. type (between 0 and 6) 78 ! = 0 closed ; = 1 cyclic East-West 79 ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 80 ! = 4 cyclic East-West AND North fold T-point pivot 81 ! = 5 North fold F-point pivot 82 ! = 6 cyclic East-West AND North fold F-point pivot 83 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 84 ! in netcdf input files, as the start j-row for reading 85 / 86 !----------------------------------------------------------------------- 87 &namzgr ! vertical coordinate 88 !----------------------------------------------------------------------- 89 ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) 90 ln_zps = .true. ! z-coordinate - partial steps (T/F) 91 ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F) 92 ln_isfcav = .false. ! ice shelf cavity (T/F) 72 &namcfg ! parameters of the configuration 73 !----------------------------------------------------------------------- 74 cp_cfg = "default" ! name of the configuration 75 cp_cfz = "no zoom" ! name of the zoom of configuration 76 jp_cfg = 0 ! resolution of the configuration 77 jpidta = 10 ! 1st lateral dimension ( >= jpi ) 78 jpjdta = 12 ! 2nd " " ( >= jpj ) 79 jpkdta = 31 ! number of levels ( >= jpk ) 80 jpiglo = 10 ! 1st dimension of global domain --> i =jpidta 81 jpjglo = 12 ! 2nd - - --> j =jpjdta 82 jpizoom = 1 ! left bottom (i,j) indices of the zoom 83 jpjzoom = 1 ! in data domain indices 84 jperio = 0 ! lateral cond. type (between 0 and 6) 85 ! = 0 closed ; = 1 cyclic East-West 86 ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 87 ! = 4 cyclic East-West AND North fold T-point pivot 88 ! = 5 North fold F-point pivot 89 ! = 6 cyclic East-West AND North fold F-point pivot 90 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 91 ! in netcdf input files, as the start j-row for reading 92 / 93 !----------------------------------------------------------------------- 94 &namzgr ! vertical coordinate (default: NO selection) 95 !----------------------------------------------------------------------- 96 ln_zco = .false. ! z-coordinate - full steps 97 ln_zps = .false. ! z-coordinate - partial steps 98 ln_sco = .false. ! s- or hybrid z-s-coordinate 99 ln_isfcav = .false. ! ice shelf cavity 100 ln_linssh = .false. ! linear free surface 93 101 / 94 102 !----------------------------------------------------------------------- 95 103 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 96 104 !----------------------------------------------------------------------- 97 ln_s_sh94 = . true. ! Song & Haidvogel 1994 hybrid S-sigma (T)|105 ln_s_sh94 = .false. ! Song & Haidvogel 1994 hybrid S-sigma (T)| 98 106 ln_s_sf12 = .false. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 99 107 ln_sigcrit = .false. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch … … 125 133 nn_msh = 1 ! create (=1) a mesh file or not (=0) 126 134 rn_hmin = -3. ! min depth of the ocean (>0) or min number of ocean level (<0) 135 rn_isfhmin = 1.00 ! treshold (m) to discriminate grounding ice to floating ice 127 136 rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of 128 137 rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 … … 130 139 rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0) 131 140 rn_atfp = 0.1 ! asselin time filter parameter 132 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k)133 ! =0, not used, rdt = rdttra134 rn_rdtmin = 28800. ! minimum time step on tracers (used if nn_acc=1)135 rn_rdtmax = 28800. ! maximum time step on tracers (used if nn_acc=1)136 rn_rdth = 800. ! depth variation of tracer time step (used if nn_acc=1)137 141 ln_crs = .false. ! Logical switch for coarsening module 138 142 jphgr_msh = 0 ! type of horizontal mesh … … 161 165 / 162 166 !----------------------------------------------------------------------- 163 &namsplit ! time splitting parameters ("key_dynspg_ts") 164 !----------------------------------------------------------------------- 165 ln_bt_fw = .TRUE. ! Forward integration of barotropic equations 166 ln_bt_av = .TRUE. ! Time filtering of barotropic variables 167 ln_bt_nn_auto = .TRUE. ! Set nn_baro automatically to be just below 168 ! a user defined maximum courant number (rn_bt_cmax) 169 nn_baro = 30 ! Number of iterations of barotropic mode 170 ! during rn_rdt seconds. Only used if ln_bt_nn_auto=F 171 rn_bt_cmax = 0.8 ! Maximum courant number allowed if ln_bt_nn_auto=T 172 nn_bt_flt = 1 ! Time filter choice 173 ! = 0 None 174 ! = 1 Boxcar over nn_baro barotropic steps 175 ! = 2 Boxcar over 2*nn_baro " " 176 / 177 !----------------------------------------------------------------------- 178 &namcrs ! Grid coarsening for dynamics output and/or 179 ! passive tracer coarsened online simulations 167 &namcrs ! coarsened grid (for outputs and/or TOP) ("key_crs") 180 168 !----------------------------------------------------------------------- 181 169 nn_factx = 3 ! Reduction factor of x-direction … … 199 187 / 200 188 !----------------------------------------------------------------------- 189 &namc1d_dyndmp ! U & V newtonian damping ("key_c1d") 190 !----------------------------------------------------------------------- 191 ln_dyndmp = .false. ! add a damping term (T) or not (F) 192 / 193 !----------------------------------------------------------------------- 194 &namc1d_uvd ! data: U & V currents ("key_c1d") 195 !----------------------------------------------------------------------- 196 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 197 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 198 sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' 199 sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' 200 ! 201 cn_dir = './' ! root directory for the location of the files 202 ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F) 203 ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F) 204 / 205 !----------------------------------------------------------------------- 201 206 &namtsd ! data : Temperature & Salinity 202 !-----------------------------------------------------------------------203 207 !----------------------------------------------------------------------- 204 208 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 211 215 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T &S input data (T) or not (F) 212 216 / 217 213 218 !!====================================================================== 214 219 !! *** Surface Boundary Condition namelists *** 215 220 !!====================================================================== 216 221 !! namsbc surface boundary condition 217 !! namsbc_ana analytical formulation 218 !! namsbc_flx flux formulation 219 !! namsbc_clio CLIO bulk formulae formulation 220 !! namsbc_core CORE bulk formulae formulation 221 !! namsbc_mfs MFS bulk formulae formulation 222 !! namsbc_cpl CouPLed formulation ("key_oasis3" )222 !! namsbc_ana analytical formulation (ln_ana =T) 223 !! namsbc_flx flux formulation (ln_flx =T) 224 !! namsbc_clio CLIO bulk formulae formulation (ln_blk_clio=T) 225 !! namsbc_core CORE bulk formulae formulation (ln_blk_core=T) 226 !! namsbc_mfs MFS bulk formulae formulation (ln_blk_mfs =T) 227 !! namsbc_cpl CouPLed formulation ("key_oasis3" ) 223 228 !! namsbc_sas StAndalone Surface module 224 !! namtra_qsr penetrative solar radiation 225 !! namsbc_rnf river runoffs 226 !! namsbc_isf ice shelf melting/freezing 227 !! namsbc_apr Atmospheric Pressure 228 !! namsbc_ssr sea surface restoring term (for T and/or S) 229 !! namtra_qsr penetrative solar radiation (ln_traqsr =T) 230 !! namsbc_rnf river runoffs (ln_rnf =T) 231 !! namsbc_isf ice shelf melting/freezing (nn_isf >0) 232 !! namsbc_iscpl coupling option between land ice model and ocean 233 !! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T) 234 !! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T) 229 235 !! namsbc_alb albedo parameters 236 !! namsbc_wave external fields from wave model (ln_wave =T) 237 !! namberg iceberg floats (ln_icebergs=T) 230 238 !!====================================================================== 231 239 ! … … 234 242 !----------------------------------------------------------------------- 235 243 nn_fsbc = 5 ! frequency of surface boundary condition computation 236 ! (also = the frequency of sea-ice model call) 244 ! (also = the frequency of sea-ice & iceberg model call) 245 ! Type of air-sea fluxes 237 246 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 238 247 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) … … 240 249 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core) 241 250 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 251 ! Type of coupling (Ocean/Ice/Atmosphere) : 242 252 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) 243 253 ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 ) … … 246 256 ! =1 opa-sas OASIS coupling: multi executable configuration, OPA component 247 257 ! =2 opa-sas OASIS coupling: multi executable configuration, SAS component 248 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr )249 nn_ice = 2 ! =0 no ice boundary condition ,250 ! =1 use observed ice-cover ,251 ! =2 ice-model used ("key_lim3" or "key_lim2")252 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect)253 ! =1 levitating ice with mass and salt exchange but no presure effect254 ! =2 embedded sea-ice (full salt and mass exchanges and pressure)255 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave256 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf)257 nn_isf = 0 ! ice shelf melting/freezing (/=0 => fill namsbc_isf)258 ! 0 =no isf 1 = presence of ISF259 ! 2 = bg03 parametrisation 3 = rnf file for isf260 ! 4 = ISF fwf specified261 ! option 1 and 4 need ln_isfcav = .true. (domzgr)262 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)263 nn_fwb = 2 ! FreshWater Budget: =0 unchecked264 ! =1 global mean of e-p-r set to zero at each time step265 ! =2 annual global mean of e-p-r set to zero266 ln_wave = .false. ! Activate coupling with wave (either Stokes Drift or Drag coefficient, or both) (T => fill namsbc_wave)267 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave)268 ln_sdw = .false. ! Computation of 3D stokes drift (T => fill namsbc_wave)269 nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,270 ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)271 258 nn_limflx = -1 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 272 259 ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled … … 274 261 ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 275 262 ! = 2 Redistribute a single flux over categories (coupled mode only) 263 ! Sea-ice : 264 nn_ice = 2 ! =0 no ice boundary condition , 265 ! =1 use observed ice-cover , 266 ! =2 ice-model used ("key_lim3", "key_lim2", "key_cice") 267 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 268 ! =1 levitating ice with mass and salt exchange but no presure effect 269 ! =2 embedded sea-ice (full salt and mass exchanges and pressure) 270 ! Misc. options of sbc : 271 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr ) 272 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 273 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 274 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 275 nn_fwb = 2 ! FreshWater Budget: =0 unchecked 276 ! =1 global mean of e-p-r set to zero at each time step 277 ! =2 annual global mean of e-p-r set to zero 278 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 279 ln_isf = .false. ! ice shelf (T => fill namsbc_isf) 280 ln_wave = .false. ! coupling with surface wave (T => fill namsbc_wave) 281 nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 282 ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 276 283 / 277 284 !----------------------------------------------------------------------- … … 316 323 &namsbc_core ! namsbc_core CORE bulk formulae 317 324 !----------------------------------------------------------------------- 318 ! ! file name 319 ! ! 320 sn_wndi = 'u_10.15JUNE2009_fill' 321 sn_wndj = 'v_10.15JUNE2009_fill' 322 sn_qsr = 'ncar_rad.15JUNE2009_fill' 323 sn_qlw = 'ncar_rad.15JUNE2009_fill' 324 sn_tair = 't_10.15JUNE2009_fill' 325 sn_humi = 'q_10.15JUNE2009_fill' 326 sn_prec = 'ncar_precip.15JUNE2009_fill' 327 sn_snow = 'ncar_precip.15JUNE2009_fill' 328 sn_tdif = 'taudif_core' 325 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 326 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 327 sn_wndi = 'u_10.15JUNE2009_fill' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Uwnd' , '' 328 sn_wndj = 'v_10.15JUNE2009_fill' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Vwnd' , '' 329 sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 330 sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 331 sn_tair = 't_10.15JUNE2009_fill' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 332 sn_humi = 'q_10.15JUNE2009_fill' , 6 , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 333 sn_prec = 'ncar_precip.15JUNE2009_fill', -1 , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 334 sn_snow = 'ncar_precip.15JUNE2009_fill', -1 , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 335 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 329 336 330 337 cn_dir = './' ! root directory for the location of the bulk files 331 338 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 332 rn_zqt = 10. 333 rn_zu = 10. 339 rn_zqt = 10. ! Air temperature and humidity reference height (m) 340 rn_zu = 10. ! Wind vector reference height (m) 334 341 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 335 342 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) … … 340 347 &namsbc_mfs ! namsbc_mfs MFS bulk formulae 341 348 !----------------------------------------------------------------------- 342 ! ! file name ! frequency (hours) ! variable 343 ! ! ! (if <0 months) ! name 344 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''345 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''346 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false. , 'daily' ,'bilinear.nc', '' ,''347 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''348 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''349 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false. , 'daily' ,'bilinear.nc', '' ,''350 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' ,''349 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 350 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 351 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 352 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 353 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false., 'daily' ,'bilinear.nc', '' , '' 354 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 355 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 356 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false., 'daily' ,'bilinear.nc', '' , '' 357 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' , '' 351 358 352 359 cn_dir = './ECMWF/' ! root directory for the location of the bulk files … … 355 362 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 356 363 !----------------------------------------------------------------------- 357 ! ! description 358 ! ! 364 ! ! description ! multiple ! vector ! vector ! vector ! 365 ! ! ! categories ! reference ! orientation ! grids ! 359 366 ! send 360 sn_snd_temp = 361 sn_snd_alb = 362 sn_snd_thick = 'none' , 'no', '' , '' , ''363 sn_snd_crt = 364 sn_snd_co2 = 367 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' 368 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' 369 sn_snd_thick = 'none' , 'no' , '' , '' , '' 370 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 371 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 365 372 ! receive 366 sn_rcv_w10m = 367 sn_rcv_taumod = 368 sn_rcv_tau = 369 sn_rcv_dqnsdt = 370 sn_rcv_qsr = 371 sn_rcv_qns = 372 sn_rcv_emp = 373 sn_rcv_rnf = 374 sn_rcv_cal = 375 sn_rcv_co2 = 373 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 374 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 375 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 376 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 377 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 378 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 379 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 380 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 381 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 382 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 376 383 ! 377 384 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentialy sending/receiving data … … 397 404 / 398 405 !----------------------------------------------------------------------- 399 &namtra_qsr ! penetrative solar radiation 406 &namtra_qsr ! penetrative solar radiation (ln_traqsr=T) 400 407 !----------------------------------------------------------------------- 401 408 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 404 411 405 412 cn_dir = './' ! root directory for the location of the runoff files 406 ln_traqsr = .true. ! Light penetration (T) or not (F)407 413 ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration 408 414 ln_qsr_2bd = .false. ! 2 bands light penetration … … 415 421 / 416 422 !----------------------------------------------------------------------- 417 &namsbc_rnf ! runoffs namelist surface boundary condition 423 &namsbc_rnf ! runoffs namelist surface boundary condition (ln_rnf=T) 418 424 !----------------------------------------------------------------------- 419 425 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 439 445 / 440 446 !----------------------------------------------------------------------- 441 &namsbc_isf ! Top boundary layer (ISF) 442 !----------------------------------------------------------------------- 443 ! ! file name ! frequency (hours) ! variable ! time interp ol. ! clim ! 'yearly'/ ! weights ! rotation!444 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing!447 &namsbc_isf ! Top boundary layer (ISF) (nn_isf >0) 448 !----------------------------------------------------------------------- 449 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 450 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 445 451 ! nn_isf == 4 446 sn_qisf = 'rnfisf' , -12 ,'sohflisf', .false. , .true. , 'yearly' , '' , '' 447 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' 452 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' , '' 448 453 ! nn_isf == 3 449 sn_rnfisf = 'runoffs' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' ,''454 sn_rnfisf = 'rnfisf' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' , '' 450 455 ! nn_isf == 2 and 3 451 sn_depmax_isf = 'runoffs' , -12 ,'sozisfmax' , .false. , .true. , 'yearly' , '' ,''452 sn_depmin_isf = 'runoffs' , -12 ,'sozisfmin' , .false. , .true. , 'yearly' , '' ,''456 sn_depmax_isf='rnfisf' , -12 ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 457 sn_depmin_isf='rnfisf' , -12 ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' 453 458 ! nn_isf == 2 454 sn_Leff_isf = 'rnfisf' , 0 ,'Leff' , .false. , .true. , 'yearly' , '' , '' 459 sn_Leff_isf = 'rnfisf' , -12 ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 460 ! 455 461 ! for all case 456 ln_divisf = .true. ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) 462 nn_isf = 1 ! ice shelf melting/freezing 463 ! 1 = presence of ISF 2 = bg03 parametrisation 464 ! 3 = rnf file for isf 4 = ISF fwf specified 465 ! option 1 and 4 need ln_isfcav = .true. (domzgr) 457 466 ! only for nn_isf = 1 or 2 458 rn_gammat0 = 1.0e-4 ! gammat coefficient used in blk formula 459 rn_gammas0 = 1.0e-4 ! gammas coefficient used in blk formula 467 rn_gammat0 = 1.e-4 ! gammat coefficient used in blk formula 468 rn_gammas0 = 1.e-4 ! gammas coefficient used in blk formula 469 ! only for nn_isf = 1 or 4 470 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 471 ! 0 => thickness of the tbl = thickness of the first wet cell 460 472 ! only for nn_isf = 1 461 nn_isfblk = 1 ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 462 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 463 ! 0 => thickness of the tbl = thickness of the first wet cell 464 ln_conserve = .true. ! conservative case (take into account meltwater advection) 473 nn_isfblk = 1 ! 1 ISOMIP like: 2 equations formulation (Hunter et al., 2006) 474 ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 465 475 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) 466 476 ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 467 ! if you want to keep the cd as in global config, adjust rn_gammat0 to compensate 468 ! 2 = velocity and stability dependent Gamma Holland et al. 1999 477 ! 2 = velocity and stability dependent Gamma (Holland et al. 1999) 478 / 479 !----------------------------------------------------------------------- 480 &namsbc_iscpl ! land ice / ocean coupling option 481 !----------------------------------------------------------------------- 482 nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells) 483 ln_hsb = .false. ! activate conservation module (conservation exact after a time of rn_fiscpl) 484 nn_fiscpl = 43800 ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 469 485 / 470 486 !----------------------------------------------------------------------- 471 487 &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk 472 488 !----------------------------------------------------------------------- 473 ! ! file name 474 ! ! ! (if <0 months) ! name ! (logical)! (T/F) ! 'monthly' ! filename ! pairing ! filename !475 sn_apr = 'patm' 489 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 490 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 491 sn_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' , '' 476 492 477 493 cn_dir = './' ! root directory for the location of the bulk files … … 481 497 / 482 498 !----------------------------------------------------------------------- 483 &namsbc_ssr ! surface boundary condition : sea surface restoring 499 &namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr=T) 484 500 !----------------------------------------------------------------------- 485 501 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 507 523 / 508 524 !----------------------------------------------------------------------- 509 &namberg ! iceberg parameters 510 !----------------------------------------------------------------------- 511 ln_icebergs = .false. 525 &namsbc_wave ! External fields from wave model (ln_wave=T) 526 !----------------------------------------------------------------------- 527 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 528 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 529 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff', .true. , .false., 'daily' , '' , '' , '' 530 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false., 'daily' , '' , '' , '' 531 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false., 'daily' , '' , '' , '' 532 sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false., 'daily' , '' , '' , '' 533 ! 534 cn_dir_cdg = './' ! root directory for the location of drag coefficient files 535 ln_cdgw = .false. ! Neutral drag coefficient read from wave model 536 ln_sdw = .false. ! Computation of 3D stokes drift 537 / 538 !----------------------------------------------------------------------- 539 &namberg ! iceberg parameters (default: No iceberg) 540 !----------------------------------------------------------------------- 541 ln_icebergs = .false. ! iceberg floats or not 512 542 ln_bergdia = .true. ! Calculate budgets 513 543 nn_verbose_level = 1 ! Turn on more verbose output if level > 0 … … 534 564 rn_speed_limit = 0. ! CFL speed limit for a berg 535 565 536 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !537 ! ! ! (if <0 months) ! name ! (logical)! (T/F ) ! 'monthly' ! filename ! pairing ! filename !538 sn_icb = 'calving' , -1 , 'calvingmask', .true., .true. , 'yearly' , '' , '' , ''566 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 567 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 568 sn_icb = 'calving', -1 , 'calvingmask', .true. , .true. , 'yearly' , '' , '' , '' 539 569 540 570 cn_dir = './' … … 545 575 !!====================================================================== 546 576 !! namlbc lateral momentum boundary condition 547 !! namcla cross land advection548 577 !! namobc open boundaries parameters ("key_obc") 549 578 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") … … 555 584 &namlbc ! lateral momentum boundary condition 556 585 !----------------------------------------------------------------------- 586 ! ! free slip ! partial slip ! no slip ! strong slip 557 587 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 558 ! free slip ! partial slip ! no slip ! strong slip 559 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical eqs. 560 / 561 !----------------------------------------------------------------------- 562 &namcla ! cross land advection 563 !----------------------------------------------------------------------- 564 nn_cla = 0 ! advection between 2 ocean pts separates by land 565 / 566 !----------------------------------------------------------------------- 567 &namobc ! open boundaries parameters ("key_obc") 568 !----------------------------------------------------------------------- 569 ln_obc_clim = .false. ! climatological obc data files (T) or not (F) 570 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 571 ln_obc_fla = .false. ! Flather open boundary condition 572 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 573 ! = 1 the obc data are read in 'obc.dta' files 574 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data 575 ! set to monthly if obc datafile hold 1 month of data 576 rn_dpein = 1. ! damping time scale for inflow at east open boundary 577 rn_dpwin = 1. ! - - - west - - 578 rn_dpnin = 1. ! - - - north - - 579 rn_dpsin = 1. ! - - - south - - 580 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary 581 rn_dpwob = 15. ! - - - west - - 582 rn_dpnob = 3000. ! - - - north - - 583 rn_dpsob = 15. ! - - - south - - 584 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R) 585 ! = 1 the total volume remains constant 588 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs. 586 589 / 587 590 !----------------------------------------------------------------------- … … 592 595 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 593 596 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 594 / 595 !----------------------------------------------------------------------- 596 &nam_tide ! tide parameters (#ifdef key_tide) 597 ln_chk_bathy = .FALSE. ! 598 / 599 !----------------------------------------------------------------------- 600 &nam_tide ! tide parameters ("key_tide") 597 601 !----------------------------------------------------------------------- 598 602 ln_tide_pot = .true. ! use tidal potential forcing … … 636 640 / 637 641 !----------------------------------------------------------------------- 638 &nambdy_dta ! open boundaries - external data("key_bdy")639 !----------------------------------------------------------------------- 640 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !641 ! ! ! (if <0 months) ! name ! (logical)! (T/F ) ! 'monthly' ! filename ! pairing ! filename !642 bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig' , .true., .false. , 'daily' , '' , '' , ''643 bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx' , .true., .false. , 'daily' , '' , '' , ''644 bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty' , .true., .false. , 'daily' , '' , '' , ''645 bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx' , .true., .false. , 'daily' , '' , '' , ''646 bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty' , .true., .false. , 'daily' , '' , '' , ''647 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true., .false. , 'daily' , '' , '' , ''648 bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline' , .true., .false. , 'daily' , '' , '' , ''642 &nambdy_dta ! open boundaries - external data ("key_bdy") 643 !----------------------------------------------------------------------- 644 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 645 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 646 bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig', .true. , .false. , 'daily' , '' , '' , '' 647 bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx', .true. , .false. , 'daily' , '' , '' , '' 648 bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty', .true. , .false. , 'daily' , '' , '' , '' 649 bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx', .true. , .false. , 'daily' , '' , '' , '' 650 bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty', .true. , .false. , 'daily' , '' , '' , '' 651 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper', .true. , .false. , 'daily' , '' , '' , '' 652 bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline', .true. , .false. , 'daily' , '' , '' , '' 649 653 ! for lim2 650 ! bn_frld = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true., .false. , 'daily' , '' , '' , ''651 ! bn_hicif = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true., .false. , 'daily' , '' , '' , ''652 ! bn_hsnif = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true., .false. , 'daily' , '' , '' , ''654 ! bn_frld = 'amm12_bdyT_ice' , 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , '' 655 ! bn_hicif = 'amm12_bdyT_ice' , 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , '' 656 ! bn_hsnif = 'amm12_bdyT_ice' , 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , '' 653 657 ! for lim3 654 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 655 ! bn_ht_i = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 656 ! bn_ht_s = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 657 cn_dir = 'bdydta/' 658 ln_full_vel = .false. 658 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , '' 659 ! bn_ht_i = 'amm12_bdyT_ice' , 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , '' 660 ! bn_ht_s = 'amm12_bdyT_ice' , 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , '' 661 662 cn_dir = 'bdydta/' ! root directory for the location of the bulk files 663 ln_full_vel = .false. ! 659 664 / 660 665 !----------------------------------------------------------------------- 661 666 &nambdy_tide ! tidal forcing at open boundaries 662 667 !----------------------------------------------------------------------- 663 filtide = 'bdydta/amm12_bdytide_' 664 ln_bdytide_2ddta = .false. 665 ln_bdytide_conj = .false. 668 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 669 ln_bdytide_2ddta = .false. ! 670 ln_bdytide_conj = .false. ! 666 671 / 667 672 !!====================================================================== … … 674 679 ! 675 680 !----------------------------------------------------------------------- 676 &nambfr ! bottom friction 681 &nambfr ! bottom friction (default: linear) 677 682 !----------------------------------------------------------------------- 678 683 nn_bfr = 1 ! type of bottom friction : = 0 : free slip, = 1 : linear friction … … 697 702 / 698 703 !----------------------------------------------------------------------- 699 &nambbc ! bottom temperature boundary condition 700 !----------------------------------------------------------------------- 701 ! ! ! (if <0 months) ! 702 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 703 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 704 sn_qgh ='geothermal_heating.nc', -12. , 'heatflow' , .false. , .true. , 'yearly' , '' , '' , '' 704 &nambbc ! bottom temperature boundary condition (default: NO) 705 !----------------------------------------------------------------------- 706 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 707 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 708 sn_qgh ='geothermal_heating.nc', -12. , 'heatflow', .false. , .true. , 'yearly' , '' , '' , '' 705 709 ! 706 cn_dir = './' ! root directory for the location of the runoff files 707 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 710 ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom 708 711 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 709 712 ! = 1 constant flux 710 713 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 711 714 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 712 713 / 714 !----------------------------------------------------------------------- 715 &nambbl ! bottom boundary layer scheme 715 cn_dir = './' ! root directory for the location of the runoff files 716 / 717 !----------------------------------------------------------------------- 718 &nambbl ! bottom boundary layer scheme ("key_trabbl") 716 719 !----------------------------------------------------------------------- 717 720 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) … … 724 727 !! Tracer (T & S ) namelists 725 728 !!====================================================================== 726 !! nameos equation of state727 !! namtra_adv advection scheme729 !! nameos equation of state 730 !! namtra_adv advection scheme 728 731 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 729 !! namtra_ldf lateral diffusion scheme 730 !! namtra_dmp T & S newtonian damping 732 !! namtra_ldf lateral diffusion scheme 733 !! namtra_ldfeiv eddy induced velocity param. 734 !! namtra_dmp T & S newtonian damping 731 735 !!====================================================================== 732 736 ! … … 739 743 ! = 1, S-EOS (simplified eos) 740 744 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 741 !!745 ! 742 746 ! ! S-EOS coefficients : 743 !! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS747 ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 744 748 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 745 749 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 751 755 / 752 756 !----------------------------------------------------------------------- 753 &namtra_adv ! advection scheme for tracer 754 !----------------------------------------------------------------------- 755 ln_traadv_cen2 = .false. ! 2nd order centered scheme 756 ln_traadv_tvd = .true. ! TVD scheme 757 ln_traadv_muscl = .false. ! MUSCL scheme 758 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 759 ln_traadv_ubs = .false. ! UBS scheme 760 ln_traadv_qck = .false. ! QUICKEST scheme 761 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 762 ln_traadv_tvd_zts= .false. ! TVD scheme with sub-timestepping of vertical tracer advection 763 / 764 !----------------------------------------------------------------------- 765 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 766 !----------------------------------------------------------------------- 767 ln_mle = .true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 757 &namtra_adv ! advection scheme for tracer (default: NO advection) 758 !----------------------------------------------------------------------- 759 ln_traadv_cen = .false. ! 2nd order centered scheme 760 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 761 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 762 ln_traadv_fct = .false. ! FCT scheme 763 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 764 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 765 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 766 ! ! (number of sub-timestep = nn_fct_zts) 767 ln_traadv_mus = .false. ! MUSCL scheme 768 ln_mus_ups = .false. ! use upstream scheme near river mouths 769 ln_traadv_ubs = .false. ! UBS scheme 770 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT / COMPACT 4th order 771 ln_traadv_qck = .false. ! QUICKEST scheme 772 / 773 !----------------------------------------------------------------------- 774 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) 775 !----------------------------------------------------------------------- 776 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 768 777 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 769 778 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 775 784 rn_rho_c_mle = 0.01 ! delta rho criterion used to calculate MLD for FK 776 785 / 777 !----------------------------------------------------------------------- -----------778 &namtra_ldf ! lateral diffusion scheme for tracers 779 !----------------------------------------------------------------------- -----------786 !----------------------------------------------------------------------- 787 &namtra_ldf ! lateral diffusion scheme for tracers (default: NO diffusion) 788 !----------------------------------------------------------------------- 780 789 ! ! Operator type: 781 ln_traldf_lap = .true. ! laplacian operator 782 ln_traldf_bilap = .false. ! bilaplacian operator 790 ! ! no diffusion: set ln_traldf_lap=..._blp=F 791 ln_traldf_lap = .false. ! laplacian operator 792 ln_traldf_blp = .false. ! bilaplacian operator 783 793 ! ! Direction of action: 784 ln_traldf_level = .false. ! iso-level 785 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 786 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 787 ! ! Griffies parameters (all need "key_ldfslp") 788 ln_traldf_grif = .false. ! use griffies triads 789 ln_traldf_gdia = .false. ! output griffies eddy velocities 790 ln_triad_iso = .false. ! pure lateral mixing in ML 791 ln_botmix_grif = .false. ! lateral mixing on bottom 792 ! ! Coefficients 793 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 794 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 795 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 796 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 797 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 798 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 799 ! (normally=0; not used with Griffies) 800 rn_slpmax = 0.01 ! slope limit 801 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 802 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 803 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 804 / 805 !----------------------------------------------------------------------- 806 &namtra_dmp ! tracer: T & S newtonian damping 794 ln_traldf_lev = .false. ! iso-level 795 ln_traldf_hor = .false. ! horizontal (geopotential) 796 ln_traldf_iso = .false. ! iso-neutral (standard operator) 797 ln_traldf_triad = .false. ! iso-neutral (triad operator) 798 ! 799 ! ! iso-neutral options: 800 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 801 rn_slpmax = 0.01 ! slope limit (both operators) 802 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 803 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 804 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 805 ! 806 ! ! Coefficients: 807 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 808 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 809 ! ! = 0 constant 810 ! ! = 10 F(k) =ldf_c1d 811 ! ! = 20 F(i,j) =ldf_c2d 812 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 813 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 814 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 815 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 816 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 817 / 818 !----------------------------------------------------------------------- 819 &namtra_ldfeiv ! eddy induced velocity param. (default: NO) 820 !----------------------------------------------------------------------- 821 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 822 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 823 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 824 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 825 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 826 ! ! = 0 constant 827 ! ! = 10 F(k) =ldf_c1d 828 ! ! = 20 F(i,j) =ldf_c2d 829 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 830 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 831 / 832 !----------------------------------------------------------------------- 833 &namtra_dmp ! tracer: T & S newtonian damping (default: NO) 807 834 !----------------------------------------------------------------------- 808 835 ln_tradmp = .true. ! add a damping termn (T) or not (F) … … 810 837 ! =1 no damping in the mixing layer (kz criteria) 811 838 ! =2 no damping in the mixed layer (rho crieria) 812 cn_resto = 'resto.nc' ! Name of file containing restoration coefficientfield (use dmp_tools to create this)839 cn_resto ='resto.nc' ! Name of file containing restoration coeff. field (use dmp_tools to create this) 813 840 / 814 841 … … 819 846 !! namdyn_vor advection scheme 820 847 !! namdyn_hpg hydrostatic pressure gradient 821 !! namdyn_spg surface pressure gradient (CPP key only)848 !! namdyn_spg surface pressure gradient 822 849 !! namdyn_ldf lateral diffusion scheme 823 850 !!====================================================================== 824 851 ! 825 852 !----------------------------------------------------------------------- 826 &namdyn_adv ! formulation of the momentum advection 853 &namdyn_adv ! formulation of the momentum advection (default: vector form) 827 854 !----------------------------------------------------------------------- 828 855 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) … … 833 860 / 834 861 !----------------------------------------------------------------------- 835 &nam_vvl ! vertical coordinate options 862 &nam_vvl ! vertical coordinate options (default: zstar) 836 863 !----------------------------------------------------------------------- 837 864 ln_vvl_zstar = .true. ! zstar vertical coordinate … … 847 874 / 848 875 !----------------------------------------------------------------------- 849 &namdyn_vor ! option of physics/algorithm (not control by CPP keys)876 &namdyn_vor ! option of physics/algorithm (default: NO) 850 877 !----------------------------------------------------------------------- 851 878 ln_dynvor_ene = .false. ! enstrophy conserving scheme 852 879 ln_dynvor_ens = .false. ! energy conserving scheme 853 880 ln_dynvor_mix = .false. ! mixed scheme 854 ln_dynvor_een = .true. ! energy & enstrophy scheme 855 ln_dynvor_een_old = .false. ! energy & enstrophy scheme - original formulation 856 / 857 !----------------------------------------------------------------------- 858 &namdyn_hpg ! Hydrostatic pressure gradient option 881 ln_dynvor_een = .false. ! energy & enstrophy scheme 882 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 883 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT ACTIVATE 884 / 885 !----------------------------------------------------------------------- 886 &namdyn_hpg ! Hydrostatic pressure gradient option (default: zps) 859 887 !----------------------------------------------------------------------- 860 888 ln_hpg_zco = .false. ! z-coordinate - full steps … … 864 892 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 865 893 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 866 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 867 ! centered time scheme (F) 868 / 869 !----------------------------------------------------------------------- 870 !namdyn_spg ! surface pressure gradient (CPP key only) 871 !----------------------------------------------------------------------- 872 ! ! explicit free surface ("key_dynspg_exp") 873 ! ! filtered free surface ("key_dynspg_flt") 874 ! ! split-explicit free surface ("key_dynspg_ts") 875 876 !----------------------------------------------------------------------- 877 &namdyn_ldf ! lateral diffusion on momentum 894 / 895 !----------------------------------------------------------------------- 896 &namdyn_spg ! surface pressure gradient (default: NO) 897 !----------------------------------------------------------------------- 898 ln_dynspg_exp = .false. ! explicit free surface 899 ln_dynspg_ts = .false. ! split-explicit free surface 900 ln_bt_fw = .true. ! Forward integration of barotropic Eqs. 901 ln_bt_av = .true. ! Time filtering of barotropic variables 902 nn_bt_flt = 1 ! Time filter choice = 0 None 903 ! ! = 1 Boxcar over nn_baro sub-steps 904 ! ! = 2 Boxcar over 2*nn_baro " " 905 ln_bt_auto = .true. ! Number of sub-step defined from: 906 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed 907 nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds 908 / 909 !----------------------------------------------------------------------- 910 &namdyn_ldf ! lateral diffusion on momentum (default: NO) 878 911 !----------------------------------------------------------------------- 879 912 ! ! Type of the operator : 880 ln_dynldf_lap = .true. ! laplacian operator 881 ln_dynldf_bilap = .false. ! bilaplacian operator 913 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 914 ln_dynldf_lap = .false. ! laplacian operator 915 ln_dynldf_blp = .false. ! bilaplacian operator 882 916 ! ! Direction of action : 883 ln_dynldf_lev el = .false.! iso-level884 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)885 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp")917 ln_dynldf_lev = .false. ! iso-level 918 ln_dynldf_hor = .false. ! horizontal (geopotential) 919 ln_dynldf_iso = .false. ! iso-neutral 886 920 ! ! Coefficient 887 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 888 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 889 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 890 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 891 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 892 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 893 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 894 rn_ahm_m_lap = 40000. ! upper limit for lap ahm < min(dx^2/16rdt, rn_ahm_m_lap) 921 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 922 ! ! =-30 read in eddy_viscosity_3D.nc file 923 ! ! =-20 read in eddy_viscosity_2D.nc file 924 ! ! = 0 constant 925 ! ! = 10 F(k)=c1d 926 ! ! = 20 F(i,j)=F(grid spacing)=c2d 927 ! ! = 30 F(i,j,k)=c2d*c1d 928 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 929 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 930 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 931 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 932 ! 933 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 895 934 / 896 935 … … 901 940 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 902 941 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 903 !! namzdf_ kpp KPP dependent vertical mixing ("key_zdfkpp")942 !! namzdf_gls GLS vertical mixing ("key_zdfgls") 904 943 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 905 944 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 933 972 rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer 934 973 rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer 935 ln_mldw = .true. ! Flag to use or not the mi zed layer depth param.974 ln_mldw = .true. ! Flag to use or not the mixed layer depth param. 936 975 / 937 976 !----------------------------------------------------------------------- … … 963 1002 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 964 1003 / 965 !------------------------------------------------------------------------966 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:967 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")968 ln_kpprimix = .true. ! shear instability mixing969 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]970 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]971 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability972 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]973 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]974 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]975 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv976 nn_ave = 1 ! constant (=0) or profile (=1) background on avt977 /978 1004 !----------------------------------------------------------------------- 979 1005 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") … … 1014 1040 !! *** Miscellaneous namelists *** 1015 1041 !!====================================================================== 1016 !! namsol elliptic solver / island / free surface1017 1042 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 1018 1043 !! namctl Control prints & Benchmark 1019 !! namc1d 1D configuration options ("key_c1d")1020 !! namc1d_uvd data: U & V currents ("key_c1d")1021 !! namc1d_dyndmp U & V newtonian damping ("key_c1d")1022 1044 !! namsto Stochastic parametrization of EOS 1023 1045 !!====================================================================== 1024 1046 ! 1025 !-----------------------------------------------------------------------1026 &namsol ! elliptic solver / island / free surface1027 !-----------------------------------------------------------------------1028 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)1029 ! =2 successive-over-relaxation (sor)1030 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test1031 rn_eps = 1.e-6 ! absolute precision of the solver1032 nn_nmin = 300 ! minimum of iterations for the SOR solver1033 nn_nmax = 800 ! maximum of iterations for the SOR solver1034 nn_nmod = 10 ! frequency of test for the SOR solver1035 rn_resmax = 1.e-10 ! absolute precision for the SOR solver1036 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain)1037 /1038 1047 !----------------------------------------------------------------------- 1039 1048 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) … … 1061 1070 ! (no physical validity of the results) 1062 1071 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 1063 / 1064 !----------------------------------------------------------------------- 1065 &namc1d_uvd ! data: U & V currents ("key_c1d") 1066 !----------------------------------------------------------------------- 1067 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 1068 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 1069 sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' 1070 sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' 1071 ! 1072 cn_dir = './' ! root directory for the location of the files 1073 ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F) 1074 ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F) 1075 / 1076 !----------------------------------------------------------------------- 1077 &namc1d_dyndmp ! U & V newtonian damping ("key_c1d") 1078 !----------------------------------------------------------------------- 1079 ln_dyndmp = .false. ! add a damping term (T) or not (F) 1080 / 1081 !----------------------------------------------------------------------- 1082 &namsto ! Stochastic parametrization of EOS 1083 !----------------------------------------------------------------------- 1084 ln_rststo = .false. ! start from mean parameter (F) or from restart file (T) 1072 nn_diacfl = 0 ! Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 1073 / 1074 !----------------------------------------------------------------------- 1075 &namsto ! Stochastic parametrization of EOS (default: NO) 1076 !----------------------------------------------------------------------- 1077 ln_sto_eos = .false. ! stochastic equation of state 1078 nn_sto_eos = 1 ! number of independent random walks 1079 rn_eos_stdxy = 1.4 ! random walk horz. standard deviation (in grid points) 1080 rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points) 1081 rn_eos_tcor = 1440. ! random walk time correlation (in timesteps) 1082 nn_eos_ord = 1 ! order of autoregressive processes 1083 nn_eos_flt = 0 ! passes of Laplacian filter 1084 rn_eos_lim = 2.0 ! limitation factor (default = 3.0) 1085 ln_rststo = .false. ! start from mean parameter (F) or from restart file (T) 1085 1086 ln_rstseed = .true. ! read seed of RNG from restart file 1086 1087 cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input) 1087 1088 cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output) 1088 1089 ln_sto_eos = .false. ! stochastic equation of state1090 nn_sto_eos = 1 ! number of independent random walks1091 rn_eos_stdxy = 1.4 ! random walk horz. standard deviation (in grid points)1092 rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points)1093 rn_eos_tcor = 1440.0 ! random walk time correlation (in timesteps)1094 nn_eos_ord = 1 ! order of autoregressive processes1095 nn_eos_flt = 0 ! passes of Laplacian filter1096 rn_eos_lim = 2.0 ! limitation factor (default = 3.0)1097 1089 / 1098 1090 … … 1100 1092 !! *** Diagnostics namelists *** 1101 1093 !!====================================================================== 1102 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")1103 1094 !! namtrd dynamics and/or tracer trends 1104 1095 !! namptr Poleward Transport Diagnostics 1096 !! namhsb Heat and salt budgets 1105 1097 !! namflo float parameters ("key_float") 1106 !! namhsb Heat and salt budgets 1107 !!====================================================================== 1108 ! 1109 !----------------------------------------------------------------------- 1110 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 1111 !----------------------------------------------------------------------- 1112 nn_nchunks_i= 4 ! number of chunks in i-dimension 1113 nn_nchunks_j= 4 ! number of chunks in j-dimension 1114 nn_nchunks_k= 31 ! number of chunks in k-dimension 1115 ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 1116 ! is optimal for postprocessing which works exclusively with horizontal slabs 1117 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 1118 ! (F) ignore chunking information and produce netcdf3-compatible files 1119 / 1120 !----------------------------------------------------------------------- 1121 &namtrd ! diagnostics on dynamics and/or tracer trends 1122 ! ! and/or mixed-layer trends and/or barotropic vorticity 1098 !! nam_diaharm Harmonic analysis of tidal constituents ('key_diaharm') 1099 !! namdct transports through some sections 1100 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 1101 !!====================================================================== 1102 ! 1103 !----------------------------------------------------------------------- 1104 &namtrd ! diagnostics on dynamics and/or tracer trends (default F) 1105 ! ! and/or mixed-layer trends and/or barotropic vorticity 1123 1106 !----------------------------------------------------------------------- 1124 1107 ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE 1125 1108 ln_dyn_trd = .false. ! (T) 3D momentum trend output 1126 ln_dyn_mxl = . FALSE. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)1127 ln_vor_trd = . FALSE. ! (T) 2D barotropic vorticity trends (not coded yet)1109 ln_dyn_mxl = .false. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 1110 ln_vor_trd = .false. ! (T) 2D barotropic vorticity trends (not coded yet) 1128 1111 ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends 1129 1112 ln_PE_trd = .false. ! (T) 3D Potential Energy trends 1130 ln_tra_trd = . FALSE. ! (T) 3D tracer trend output1113 ln_tra_trd = .false. ! (T) 3D tracer trend output 1131 1114 ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 1132 1115 nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step) … … 1139 1122 !!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 1140 1123 !!gm 1124 !----------------------------------------------------------------------- 1125 &namptr ! Poleward Transport Diagnostic (default F) 1126 !----------------------------------------------------------------------- 1127 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1128 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 1129 / 1130 !----------------------------------------------------------------------- 1131 &namhsb ! Heat and salt budgets (default F) 1132 !----------------------------------------------------------------------- 1133 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 1134 / 1141 1135 !----------------------------------------------------------------------- 1142 1136 &namflo ! float parameters ("key_float") … … 1154 1148 / 1155 1149 !----------------------------------------------------------------------- 1156 &namptr ! Poleward Transport Diagnostic 1157 !----------------------------------------------------------------------- 1158 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1159 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 1160 / 1161 !----------------------------------------------------------------------- 1162 &namhsb ! Heat and salt budgets 1163 !----------------------------------------------------------------------- 1164 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 1165 / 1166 !----------------------------------------------------------------------- 1167 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 1150 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 1168 1151 !----------------------------------------------------------------------- 1169 1152 nit000_han = 1 ! First time step used for harmonic analysis … … 1174 1157 / 1175 1158 !----------------------------------------------------------------------- 1176 &namdct ! transports through s ections1159 &namdct ! transports through some sections 1177 1160 !----------------------------------------------------------------------- 1178 1161 nn_dct = 15 ! time step frequency for transports computing … … 1182 1165 ! 0 < n : debug section number n 1183 1166 / 1184 1185 !!====================================================================== 1186 !! *** Observation & Assimilation namelists *** 1187 !!====================================================================== 1188 !! namobs observation and model comparison ('key_diaobs') 1167 !----------------------------------------------------------------------- 1168 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 1169 !----------------------------------------------------------------------- 1170 nn_nchunks_i= 4 ! number of chunks in i-dimension 1171 nn_nchunks_j= 4 ! number of chunks in j-dimension 1172 nn_nchunks_k= 31 ! number of chunks in k-dimension 1173 ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 1174 ! is optimal for postprocessing which works exclusively with horizontal slabs 1175 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 1176 ! (F) ignore chunking information and produce netcdf3-compatible files 1177 / 1178 1179 !!====================================================================== 1180 !! *** Observation & Assimilation *** 1181 !!====================================================================== 1182 !! namobs observation and model comparison 1189 1183 !! nam_asminc assimilation increments ('key_asminc') 1190 1184 !!====================================================================== 1191 1185 ! 1192 1186 !----------------------------------------------------------------------- 1193 &namobs ! observation usage switch ('key_diaobs') 1194 !----------------------------------------------------------------------- 1195 ln_t3d = .false. ! Logical switch for T profile observations 1196 ln_s3d = .false. ! Logical switch for S profile observations 1197 ln_ena = .false. ! Logical switch for ENACT insitu data set 1198 ! ! ln_cor Logical switch for Coriolis insitu data set 1199 ln_profb = .false. ! Logical switch for feedback insitu data set 1200 ln_sla = .false. ! Logical switch for SLA observations 1201 1202 ln_sladt = .false. ! Logical switch for AVISO SLA data 1203 1204 ln_slafb = .false. ! Logical switch for feedback SLA data 1205 ! ln_ssh Logical switch for SSH observations 1206 1207 ln_sst = .false. ! Logical switch for SST observations 1208 ln_reysst = .false. ! ln_reysst Logical switch for Reynolds observations 1209 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 1210 1211 ln_sstfb = .false. ! Logical switch for feedback SST data 1212 ! ln_sss Logical switch for SSS observations 1213 ln_seaice = .false. ! Logical switch for Sea Ice observations 1214 ! ln_vel3d Logical switch for velocity observations 1215 ! ln_velavcur Logical switch for velocity daily av. cur. 1216 ! ln_velhrcur Logical switch for velocity high freq. cur. 1217 ! ln_velavadcp Logical switch for velocity daily av. ADCP 1218 ! ln_velhradcp Logical switch for velocity high freq. ADCP 1219 ! ln_velfb Logical switch for feedback velocity data 1220 ! ln_grid_global Global distribtion of observations 1221 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 1222 ! grid_search_file Grid search lookup file header 1223 ! enactfiles ENACT input observation file names 1224 ! coriofiles Coriolis input observation file name 1225 ! ! profbfiles: Profile feedback input observation file name 1226 profbfiles = 'profiles_01.nc' 1227 ! ln_profb_enatim Enact feedback input time setting switch 1228 ! slafilesact Active SLA input observation file name 1229 ! slafilespas Passive SLA input observation file name 1230 ! ! slafbfiles: Feedback SLA input observation file name 1231 slafbfiles = 'sla_01.nc' 1232 ! sstfiles GHRSST input observation file name 1233 ! ! sstfbfiles: Feedback SST input observation file name 1234 sstfbfiles = 'sst_01.nc' 1235 ! seaicefiles Sea Ice input observation file names 1236 seaicefiles = 'seaice_01.nc' 1237 ! velavcurfiles Vel. cur. daily av. input file name 1238 ! velhvcurfiles Vel. cur. high freq. input file name 1239 ! velavadcpfiles Vel. ADCP daily av. input file name 1240 ! velhvadcpfiles Vel. ADCP high freq. input file name 1241 ! velfbfiles Vel. feedback input observation file name 1242 ! dobsini Initial date in window YYYYMMDD.HHMMSS 1243 ! dobsend Final date in window YYYYMMDD.HHMMSS 1244 ! n1dint Type of vertical interpolation method 1245 ! n2dint Type of horizontal interpolation method 1246 ! ln_nea Rejection of observations near land switch 1247 nmsshc = 0 ! MSSH correction scheme 1248 ! mdtcorr MDT correction 1249 ! mdtcutoff MDT cutoff for computed correction 1250 ln_altbias = .false. ! Logical switch for alt bias 1251 ln_ignmis = .true. ! Logical switch for ignoring missing files 1252 ! endailyavtypes ENACT daily average types 1253 ln_grid_global = .true. 1254 ln_grid_search_lookup = .false. 1187 &namobs ! observation usage switch 1188 !----------------------------------------------------------------------- 1189 ln_diaobs = .false. ! Logical switch for the observation operator 1190 ln_t3d = .false. ! Logical switch for T profile observations 1191 ln_s3d = .false. ! Logical switch for S profile observations 1192 ln_sla = .false. ! Logical switch for SLA observations 1193 ln_sst = .false. ! Logical switch for SST observations 1194 ln_sic = .false. ! Logical switch for Sea Ice observations 1195 ln_vel3d = .false. ! Logical switch for velocity observations 1196 ln_altbias = .false. ! Logical switch for altimeter bias correction 1197 ln_nea = .false. ! Logical switch for rejection of observations near land 1198 ln_grid_global = .true. ! Logical switch for global distribution of observations 1199 ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table 1200 ln_ignmis = .true. ! Logical switch for ignoring missing files 1201 ln_s_at_t = .false. ! Logical switch for computing model S at T obs if not there 1202 ln_sstnight = .false. ! Logical switch for calculating night-time average for SST obs 1203 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files 1204 cn_profbfiles = 'profiles_01.nc' ! Profile feedback input observation file names 1205 cn_slafbfiles = 'sla_01.nc' ! SLA feedback input observation file names 1206 cn_sstfbfiles = 'sst_01.nc' ! SST feedback input observation file names 1207 cn_sicfbfiles = 'sic_01.nc' ! SIC feedback input observation file names 1208 cn_velfbfiles = 'vel_01.nc' ! Velocity feedback input observation file names 1209 cn_altbiasfile = 'altbias.nc' ! Altimeter bias input file name 1210 cn_gridsearchfile = 'gridsearch.nc' ! Grid search file name 1211 rn_gridsearchres = 0.5 ! Grid search resolution 1212 rn_dobsini = 00010101.000000 ! Initial date in window YYYYMMDD.HHMMSS 1213 rn_dobsend = 00010102.000000 ! Final date in window YYYYMMDD.HHMMSS 1214 nn_1dint = 0 ! Type of vertical interpolation method 1215 nn_2dint = 0 ! Type of horizontal interpolation method 1216 nn_msshc = 0 ! MSSH correction scheme 1217 rn_mdtcorr = 1.61 ! MDT correction 1218 rn_mdtcutoff = 65.0 ! MDT cutoff for computed correction 1219 nn_profdavtypes = -1 ! Profile daily average types - array 1220 ln_sstbias = .false. 1221 cn_sstbias_files = 'sstbias.nc' 1255 1222 / 1256 1223 !----------------------------------------------------------------------- … … 1273 1240 / 1274 1241 !----------------------------------------------------------------------- 1275 &namsbc_wave ! External fields from wave model 1276 !----------------------------------------------------------------------- 1277 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 1278 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 1279 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' , '' 1280 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1281 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1282 sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' , '' , '' , '' 1283 ! 1284 cn_dir_cdg = './' ! root directory for the location of drag coefficient files 1285 / 1286 !----------------------------------------------------------------------- 1287 &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) 1288 !----------------------------------------------------------------------- 1289 ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 1290 ln_neptsimp = .false. ! yes/no use simplified neptune 1291 1292 ln_smooth_neptvel = .false. ! yes/no smooth zunep, zvnep 1293 rn_tslse = 1.2e4 ! value of lengthscale L at the equator 1294 rn_tslsp = 3.0e3 ! value of lengthscale L at the pole 1295 ! Specify whether to ramp down the Neptune velocity in shallow 1296 ! water, and if so the depth range controlling such ramping down 1297 ln_neptramp = .true. ! ramp down Neptune velocity in shallow water 1298 rn_htrmin = 100.0 ! min. depth of transition range 1299 rn_htrmax = 200.0 ! max. depth of transition range 1300 / 1242 &namdiu ! Cool skin and warm layer models 1243 !----------------------------------------------------------------------- 1244 ln_diurnal = .false. ! 1245 ln_diurnal_only = .false. ! 1246 / 1247 !----------------------------------------------------------------------- 1248 &nam_diatmb ! Top Middle Bottom Output 1249 !----------------------------------------------------------------------- 1250 ln_diatmb = .false. ! Choose Top Middle and Bottom output or not 1251 / 1252 !----------------------------------------------------------------------- 1253 &namwad ! Wetting and drying 1254 !----------------------------------------------------------------------- 1255 ln_wd = .false. ! T/F activation of wetting and drying 1256 rn_wdmin1 = 0.1 ! Minimum wet depth on dried cells 1257 rn_wdmin2 = 0.01 ! Tolerance of min wet depth on dried cells 1258 rn_wdld = 20.0 ! Land elevation below which wetting/drying is allowed 1259 nn_wdit = 10 ! Max iterations for W/D limiter 1260 / 1261 !----------------------------------------------------------------------- 1262 &nam_dia25h ! 25h Mean Output 1263 !----------------------------------------------------------------------- 1264 ln_dia25h = .false. ! Choose 25h mean output or not 1265 /
Note: See TracChangeset
for help on using the changeset viewer.