Changeset 6404 for branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/SHARED/namelist_ref
- Timestamp:
- 2016-03-29T11:24:48+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/SHARED/namelist_ref
r6401 r6404 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 = .true. ! 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 !----------------------------------------------------------------------- … … 595 598 / 596 599 !----------------------------------------------------------------------- 597 &nam_tide ! tide parameters (#ifdef key_tide)600 &nam_tide ! tide parameters ("key_tide") 598 601 !----------------------------------------------------------------------- 599 602 ln_tide_pot = .true. ! use tidal potential forcing … … 637 640 / 638 641 !----------------------------------------------------------------------- 639 &nambdy_dta ! open boundaries - external data("key_bdy")640 !----------------------------------------------------------------------- 641 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !642 ! ! ! (if <0 months) ! name ! (logical)! (T/F ) ! 'monthly' ! filename ! pairing ! filename !643 bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig' , .true., .false. , 'daily' , '' , '' , ''644 bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx' , .true., .false. , 'daily' , '' , '' , ''645 bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty' , .true., .false. , 'daily' , '' , '' , ''646 bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx' , .true., .false. , 'daily' , '' , '' , ''647 bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty' , .true., .false. , 'daily' , '' , '' , ''648 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true., .false. , 'daily' , '' , '' , ''649 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' , '' , '' , '' 650 653 ! for lim2 651 ! bn_frld = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true., .false. , 'daily' , '' , '' , ''652 ! bn_hicif = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true., .false. , 'daily' , '' , '' , ''653 ! 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' , '' , '' , '' 654 657 ! for lim3 655 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 656 ! bn_ht_i = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 657 ! bn_ht_s = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 658 cn_dir = 'bdydta/' 659 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. ! 660 664 / 661 665 !----------------------------------------------------------------------- 662 666 &nambdy_tide ! tidal forcing at open boundaries 663 667 !----------------------------------------------------------------------- 664 filtide = 'bdydta/amm12_bdytide_' 665 ln_bdytide_2ddta = .false. 666 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. ! 667 671 / 668 672 !!====================================================================== … … 675 679 ! 676 680 !----------------------------------------------------------------------- 677 &nambfr ! bottom friction 681 &nambfr ! bottom friction (default: linear) 678 682 !----------------------------------------------------------------------- 679 683 nn_bfr = 1 ! type of bottom friction : = 0 : free slip, = 1 : linear friction … … 698 702 / 699 703 !----------------------------------------------------------------------- 700 &nambbc ! bottom temperature boundary condition 701 !----------------------------------------------------------------------- 702 ! ! ! (if <0 months) ! 703 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 704 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 705 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' , '' , '' , '' 706 709 ! 707 cn_dir = './' ! root directory for the location of the runoff files 708 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 710 ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom 709 711 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 710 712 ! = 1 constant flux 711 713 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 712 714 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 713 714 / 715 !----------------------------------------------------------------------- 716 &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") 717 719 !----------------------------------------------------------------------- 718 720 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) … … 725 727 !! Tracer (T & S ) namelists 726 728 !!====================================================================== 727 !! nameos equation of state728 !! namtra_adv advection scheme729 !! nameos equation of state 730 !! namtra_adv advection scheme 729 731 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 730 !! namtra_ldf lateral diffusion scheme 731 !! 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 732 735 !!====================================================================== 733 736 ! … … 740 743 ! = 1, S-EOS (simplified eos) 741 744 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 742 !!745 ! 743 746 ! ! S-EOS coefficients : 744 !! 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 745 748 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 746 749 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 752 755 / 753 756 !----------------------------------------------------------------------- 754 &namtra_adv ! advection scheme for tracer 755 !----------------------------------------------------------------------- 756 ln_traadv_cen2 = .false. ! 2nd order centered scheme 757 ln_traadv_tvd = .true. ! TVD scheme 758 ln_traadv_muscl = .false. ! MUSCL scheme 759 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 760 ln_traadv_ubs = .false. ! UBS scheme 761 ln_traadv_qck = .false. ! QUICKEST scheme 762 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 763 ln_traadv_tvd_zts= .false. ! TVD scheme with sub-timestepping of vertical tracer advection 764 / 765 !----------------------------------------------------------------------- 766 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 767 !----------------------------------------------------------------------- 768 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 769 777 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 770 778 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 776 784 rn_rho_c_mle = 0.01 ! delta rho criterion used to calculate MLD for FK 777 785 / 778 !----------------------------------------------------------------------- -----------779 &namtra_ldf ! lateral diffusion scheme for tracers 780 !----------------------------------------------------------------------- -----------786 !----------------------------------------------------------------------- 787 &namtra_ldf ! lateral diffusion scheme for tracers (default: NO diffusion) 788 !----------------------------------------------------------------------- 781 789 ! ! Operator type: 782 ln_traldf_lap = .true. ! laplacian operator 783 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 784 793 ! ! Direction of action: 785 ln_traldf_level = .false. ! iso-level 786 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 787 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 788 ! ! Griffies parameters (all need "key_ldfslp") 789 ln_traldf_grif = .false. ! use griffies triads 790 ln_traldf_gdia = .false. ! output griffies eddy velocities 791 ln_triad_iso = .false. ! pure lateral mixing in ML 792 ln_botmix_grif = .false. ! lateral mixing on bottom 793 ! ! Coefficients 794 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 795 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 796 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 797 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 798 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 799 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 800 ! (normally=0; not used with Griffies) 801 rn_slpmax = 0.01 ! slope limit 802 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 803 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 804 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 805 / 806 !----------------------------------------------------------------------- 807 &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) 808 834 !----------------------------------------------------------------------- 809 835 ln_tradmp = .true. ! add a damping termn (T) or not (F) … … 811 837 ! =1 no damping in the mixing layer (kz criteria) 812 838 ! =2 no damping in the mixed layer (rho crieria) 813 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) 814 840 / 815 841 … … 820 846 !! namdyn_vor advection scheme 821 847 !! namdyn_hpg hydrostatic pressure gradient 822 !! namdyn_spg surface pressure gradient (CPP key only)848 !! namdyn_spg surface pressure gradient 823 849 !! namdyn_ldf lateral diffusion scheme 824 850 !!====================================================================== 825 851 ! 826 852 !----------------------------------------------------------------------- 827 &namdyn_adv ! formulation of the momentum advection 853 &namdyn_adv ! formulation of the momentum advection (default: vector form) 828 854 !----------------------------------------------------------------------- 829 855 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) … … 834 860 / 835 861 !----------------------------------------------------------------------- 836 &nam_vvl ! vertical coordinate options 862 &nam_vvl ! vertical coordinate options (default: zstar) 837 863 !----------------------------------------------------------------------- 838 864 ln_vvl_zstar = .true. ! zstar vertical coordinate … … 848 874 / 849 875 !----------------------------------------------------------------------- 850 &namdyn_vor ! option of physics/algorithm (not control by CPP keys)876 &namdyn_vor ! option of physics/algorithm (default: NO) 851 877 !----------------------------------------------------------------------- 852 878 ln_dynvor_ene = .false. ! enstrophy conserving scheme 853 879 ln_dynvor_ens = .false. ! energy conserving scheme 854 880 ln_dynvor_mix = .false. ! mixed scheme 855 ln_dynvor_een = .true. ! energy & enstrophy scheme 856 ln_dynvor_een_old = .false. ! energy & enstrophy scheme - original formulation 857 / 858 !----------------------------------------------------------------------- 859 &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) 860 887 !----------------------------------------------------------------------- 861 888 ln_hpg_zco = .false. ! z-coordinate - full steps … … 865 892 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 866 893 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 867 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 868 ! centered time scheme (F) 869 / 870 !----------------------------------------------------------------------- 871 !namdyn_spg ! surface pressure gradient (CPP key only) 872 !----------------------------------------------------------------------- 873 ! ! explicit free surface ("key_dynspg_exp") 874 ! ! filtered free surface ("key_dynspg_flt") 875 ! ! split-explicit free surface ("key_dynspg_ts") 876 877 !----------------------------------------------------------------------- 878 &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) 879 911 !----------------------------------------------------------------------- 880 912 ! ! Type of the operator : 881 ln_dynldf_lap = .true. ! laplacian operator 882 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 883 916 ! ! Direction of action : 884 ln_dynldf_lev el = .false.! iso-level885 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)886 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 887 920 ! ! Coefficient 888 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 889 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 890 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 891 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 892 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 893 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 894 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 895 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) 896 934 / 897 935 … … 902 940 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 903 941 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 904 !! namzdf_ kpp KPP dependent vertical mixing ("key_zdfkpp")942 !! namzdf_gls GLS vertical mixing ("key_zdfgls") 905 943 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 906 944 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 934 972 rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer 935 973 rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer 936 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. 937 975 / 938 976 !----------------------------------------------------------------------- … … 964 1002 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 965 1003 / 966 !------------------------------------------------------------------------967 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:968 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")969 ln_kpprimix = .true. ! shear instability mixing970 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]971 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]972 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability973 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]974 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]975 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]976 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv977 nn_ave = 1 ! constant (=0) or profile (=1) background on avt978 /979 1004 !----------------------------------------------------------------------- 980 1005 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") … … 1015 1040 !! *** Miscellaneous namelists *** 1016 1041 !!====================================================================== 1017 !! namsol elliptic solver / island / free surface1018 1042 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 1019 1043 !! namctl Control prints & Benchmark 1020 !! namc1d 1D configuration options ("key_c1d")1021 !! namc1d_uvd data: U & V currents ("key_c1d")1022 !! namc1d_dyndmp U & V newtonian damping ("key_c1d")1023 1044 !! namsto Stochastic parametrization of EOS 1024 1045 !!====================================================================== 1025 1046 ! 1026 !-----------------------------------------------------------------------1027 &namsol ! elliptic solver / island / free surface1028 !-----------------------------------------------------------------------1029 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)1030 ! =2 successive-over-relaxation (sor)1031 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test1032 rn_eps = 1.e-6 ! absolute precision of the solver1033 nn_nmin = 300 ! minimum of iterations for the SOR solver1034 nn_nmax = 800 ! maximum of iterations for the SOR solver1035 nn_nmod = 10 ! frequency of test for the SOR solver1036 rn_resmax = 1.e-10 ! absolute precision for the SOR solver1037 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain)1038 /1039 1047 !----------------------------------------------------------------------- 1040 1048 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) … … 1062 1070 ! (no physical validity of the results) 1063 1071 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 1064 / 1065 !----------------------------------------------------------------------- 1066 &namc1d_uvd ! data: U & V currents ("key_c1d") 1067 !----------------------------------------------------------------------- 1068 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 1069 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 1070 sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' 1071 sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' 1072 ! 1073 cn_dir = './' ! root directory for the location of the files 1074 ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F) 1075 ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F) 1076 / 1077 !----------------------------------------------------------------------- 1078 &namc1d_dyndmp ! U & V newtonian damping ("key_c1d") 1079 !----------------------------------------------------------------------- 1080 ln_dyndmp = .false. ! add a damping term (T) or not (F) 1081 / 1082 !----------------------------------------------------------------------- 1083 &namsto ! Stochastic parametrization of EOS 1084 !----------------------------------------------------------------------- 1085 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) 1086 1086 ln_rstseed = .true. ! read seed of RNG from restart file 1087 1087 cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input) 1088 1088 cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output) 1089 1090 ln_sto_eos = .false. ! stochastic equation of state1091 nn_sto_eos = 1 ! number of independent random walks1092 rn_eos_stdxy = 1.4 ! random walk horz. standard deviation (in grid points)1093 rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points)1094 rn_eos_tcor = 1440.0 ! random walk time correlation (in timesteps)1095 nn_eos_ord = 1 ! order of autoregressive processes1096 nn_eos_flt = 0 ! passes of Laplacian filter1097 rn_eos_lim = 2.0 ! limitation factor (default = 3.0)1098 1089 / 1099 1090 … … 1101 1092 !! *** Diagnostics namelists *** 1102 1093 !!====================================================================== 1103 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")1104 1094 !! namtrd dynamics and/or tracer trends 1105 1095 !! namptr Poleward Transport Diagnostics 1096 !! namhsb Heat and salt budgets 1106 1097 !! namflo float parameters ("key_float") 1107 !! namhsb Heat and salt budgets 1108 !!====================================================================== 1109 ! 1110 !----------------------------------------------------------------------- 1111 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 1112 !----------------------------------------------------------------------- 1113 nn_nchunks_i= 4 ! number of chunks in i-dimension 1114 nn_nchunks_j= 4 ! number of chunks in j-dimension 1115 nn_nchunks_k= 31 ! number of chunks in k-dimension 1116 ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 1117 ! is optimal for postprocessing which works exclusively with horizontal slabs 1118 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 1119 ! (F) ignore chunking information and produce netcdf3-compatible files 1120 / 1121 !----------------------------------------------------------------------- 1122 &namtrd ! diagnostics on dynamics and/or tracer trends 1123 ! ! 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 1124 1106 !----------------------------------------------------------------------- 1125 1107 ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE 1126 1108 ln_dyn_trd = .false. ! (T) 3D momentum trend output 1127 ln_dyn_mxl = . FALSE. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)1128 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) 1129 1111 ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends 1130 1112 ln_PE_trd = .false. ! (T) 3D Potential Energy trends 1131 ln_tra_trd = . FALSE. ! (T) 3D tracer trend output1113 ln_tra_trd = .false. ! (T) 3D tracer trend output 1132 1114 ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 1133 1115 nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step) … … 1140 1122 !!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 1141 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 / 1142 1135 !----------------------------------------------------------------------- 1143 1136 &namflo ! float parameters ("key_float") … … 1155 1148 / 1156 1149 !----------------------------------------------------------------------- 1157 &namptr ! Poleward Transport Diagnostic 1158 !----------------------------------------------------------------------- 1159 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1160 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 1161 / 1162 !----------------------------------------------------------------------- 1163 &namhsb ! Heat and salt budgets 1164 !----------------------------------------------------------------------- 1165 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 1166 / 1167 !----------------------------------------------------------------------- 1168 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 1150 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 1169 1151 !----------------------------------------------------------------------- 1170 1152 nit000_han = 1 ! First time step used for harmonic analysis … … 1175 1157 / 1176 1158 !----------------------------------------------------------------------- 1177 &namdct ! transports through s ections1159 &namdct ! transports through some sections 1178 1160 !----------------------------------------------------------------------- 1179 1161 nn_dct = 15 ! time step frequency for transports computing … … 1183 1165 ! 0 < n : debug section number n 1184 1166 / 1185 1186 !!====================================================================== 1187 !! *** Observation & Assimilation namelists *** 1188 !!====================================================================== 1189 !! 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 1190 1183 !! nam_asminc assimilation increments ('key_asminc') 1191 1184 !!====================================================================== 1192 1185 ! 1193 1186 !----------------------------------------------------------------------- 1194 &namobs ! observation usage switch ('key_diaobs') 1195 !----------------------------------------------------------------------- 1196 ln_t3d = .false. ! Logical switch for T profile observations 1197 ln_s3d = .false. ! Logical switch for S profile observations 1198 ln_ena = .false. ! Logical switch for ENACT insitu data set 1199 ln_cor = .false. ! Logical switch for Coriolis insitu data set 1200 ln_profb = .false. ! Logical switch for feedback insitu data set 1201 ln_sla = .false. ! Logical switch for SLA observations 1202 ln_sladt = .false. ! Logical switch for AVISO SLA data 1203 ln_slafb = .false. ! Logical switch for feedback SLA data 1204 ln_ssh = .false. ! Logical switch for SSH observations 1205 ln_sst = .false. ! Logical switch for SST observations 1206 ln_reysst = .false. ! Logical switch for Reynolds observations 1207 ln_ghrsst = .false. ! Logical switch for GHRSST observations 1208 ln_sstfb = .false. ! Logical switch for feedback SST data 1209 ln_sss = .false. ! Logical switch for SSS observations 1210 ln_seaice = .false. ! Logical switch for Sea Ice observations 1211 ln_vel3d = .false. ! Logical switch for velocity observations 1212 ln_velavcur= .false ! Logical switch for velocity daily av. cur. 1213 ln_velhrcur= .false ! Logical switch for velocity high freq. cur. 1214 ln_velavadcp = .false. ! Logical switch for velocity daily av. ADCP 1215 ln_velhradcp = .false. ! Logical switch for velocity high freq. ADCP 1216 ln_velfb = .false. ! Logical switch for feedback velocity data 1217 ln_grid_global = .false. ! Global distribtion of observations 1218 ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table 1219 grid_search_file = 'grid_search' ! Grid search lookup file header 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 1220 1203 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files 1221 enactfiles = 'enact.nc' ! ENACT input observation file names (specify full array in namelist_cfg) 1222 coriofiles = 'corio.nc' ! Coriolis input observation file name 1223 profbfiles = 'profiles_01.nc' ! Profile feedback input observation file name 1224 ln_profb_enatim = .false ! Enact feedback input time setting switch 1225 slafilesact = 'sla_act.nc' ! Active SLA input observation file names 1226 slafilespas = 'sla_pass.nc' ! Passive SLA input observation file names 1227 slafbfiles = 'sla_01.nc' ! slafbfiles: Feedback SLA input observation file names 1228 sstfiles = 'ghrsst.nc' ! GHRSST input observation file names 1229 sstfbfiles = 'sst_01.nc' ! Feedback SST input observation file names 1230 seaicefiles = 'seaice_01.nc' ! Sea Ice input observation file names 1231 velavcurfiles = 'velavcurfile.nc' ! Vel. cur. daily av. input file name 1232 velhrcurfiles = 'velhrcurfile.nc' ! Vel. cur. high freq. input file name 1233 velavadcpfiles = 'velavadcpfile.nc' ! Vel. ADCP daily av. input file name 1234 velhradcpfiles = 'velhradcpfile.nc' ! Vel. ADCP high freq. input file name 1235 velfbfiles = 'velfbfile.nc' ! Vel. feedback input observation file name 1236 dobsini = 20000101.000000 ! Initial date in window YYYYMMDD.HHMMSS 1237 dobsend = 20010101.000000 ! Final date in window YYYYMMDD.HHMMSS 1238 n1dint = 0 ! Type of vertical interpolation method 1239 n2dint = 0 ! Type of horizontal interpolation method 1240 ln_nea = .false. ! Rejection of observations near land switch 1241 nmsshc = 0 ! MSSH correction scheme 1242 mdtcorr = 1.61 ! MDT correction 1243 mdtcutoff = 65.0 ! MDT cutoff for computed correction 1244 ln_altbias = .false. ! Logical switch for alt bias 1245 ln_ignmis = .true. ! Logical switch for ignoring missing files 1246 endailyavtypes = 820 ! ENACT daily average types - array (use namelist_cfg to set more values) 1247 ln_grid_global = .true. 1248 ln_grid_search_lookup = .false. 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' 1249 1222 / 1250 1223 !----------------------------------------------------------------------- … … 1267 1240 / 1268 1241 !----------------------------------------------------------------------- 1269 &namsbc_wave ! External fields from wave model 1270 !----------------------------------------------------------------------- 1271 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 1272 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 1273 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' , '' 1274 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1275 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1276 sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' , '' , '' , '' 1277 ! 1278 cn_dir_cdg = './' ! root directory for the location of drag coefficient files 1279 / 1280 !----------------------------------------------------------------------- 1281 &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) 1282 !----------------------------------------------------------------------- 1283 ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 1284 ln_neptsimp = .false. ! yes/no use simplified neptune 1285 1286 ln_smooth_neptvel = .false. ! yes/no smooth zunep, zvnep 1287 rn_tslse = 1.2e4 ! value of lengthscale L at the equator 1288 rn_tslsp = 3.0e3 ! value of lengthscale L at the pole 1289 ! Specify whether to ramp down the Neptune velocity in shallow 1290 ! water, and if so the depth range controlling such ramping down 1291 ln_neptramp = .true. ! ramp down Neptune velocity in shallow water 1292 rn_htrmin = 100.0 ! min. depth of transition range 1293 rn_htrmax = 200.0 ! max. depth of transition range 1294 / 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.