- Timestamp:
- 2016-01-08T10:35:19+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/CONFIG/SHARED/namelist_ref
- Property svn:keywords set to Id
r4699 r6225 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : 1 - run manager (namrun, namcfg) 3 !! namelists 2 - Domain (namzgr, namzgr_sco, namdom, namtsd) 2 !! namelist_ref 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 4 !! NEMO/OPA : 1 - run manager (namrun) 5 !! namelists 2 - Domain (namcfg, namzgr, namzgr_sco, namdom, namtsd) 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)12 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, nam ptr, namhsb)13 !! 10 - miscellaneous (nam sol, nammpp, namctl)13 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 14 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto) 15 !! 10 - miscellaneous (nammpp, namctl) 14 16 !! 11 - Obs & Assim (namobs, nam_asminc) 15 17 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 18 20 !! *** Run management namelists *** 19 21 !!====================================================================== 20 !! namrun 22 !! namrun parameters of the run 21 23 !!====================================================================== 22 24 ! … … 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=.true. 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_out = "restart" ! suffix of ocean restart name (output) 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 40 46 nn_istate = 0 ! output the initial state (1) or not (0) 47 ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 41 48 nn_stock = 5475 ! frequency of creation of a restart file (modulo referenced to 1) 49 nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 42 50 nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000) 43 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T)44 51 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 45 ln_clobber = .false. ! clobber (overwrite) an existing file 52 ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard 53 ln_clobber = .true. ! clobber (overwrite) an existing file 46 54 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 47 55 / 48 56 ! 49 !-----------------------------------------------------------------------50 &namcfg ! default parameters of the configuration51 !-----------------------------------------------------------------------52 cp_cfg = "default" ! name of the configuration53 cp_cfz = '' ! name of the zoom of configuration54 jp_cfg = 0 ! resolution of the configuration55 jpidta = 10 ! 1st lateral dimension ( >= jpi )56 jpjdta = 12 ! 2nd " " ( >= jpj )57 jpkdta = 31 ! number of levels ( >= jpk )58 jpiglo = 10 ! 1st dimension of global domain --> i =jpidta59 jpjglo = 12 ! 2nd - - --> j =jpjdta60 jpizoom = 1 ! left bottom (i,j) indices of the zoom61 jpjzoom = 1 ! in data domain indices62 jperio = 0 ! lateral cond. type (between 0 and 6)63 ! = 0 closed ; = 1 cyclic East-West64 ! = 2 equatorial symmetric ; = 3 North fold T-point pivot65 ! = 4 cyclic East-West AND North fold T-point pivot66 ! = 5 North fold F-point pivot67 ! = 6 cyclic East-West AND North fold F-point pivot68 /69 57 !!====================================================================== 70 58 !! *** Domain namelists *** 71 59 !!====================================================================== 60 !! namcfg parameters of the configuration 72 61 !! namzgr vertical coordinate 73 62 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 74 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") 75 68 !! namtsd data: temperature & salinity 76 69 !!====================================================================== 77 70 ! 78 71 !----------------------------------------------------------------------- 79 &namzgr ! vertical coordinate 80 !----------------------------------------------------------------------- 81 ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) 82 ln_zps = .true. ! z-coordinate - partial steps (T/F) 83 ln_sco = .false. ! s- or hybrid z-s-coordinate (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 84 101 / 85 102 !----------------------------------------------------------------------- 86 103 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 87 104 !----------------------------------------------------------------------- 88 ln_s_sh94 = . true. ! Song & Haidvogel 1994 hybrid S-sigma (T)|105 ln_s_sh94 = .false. ! Song & Haidvogel 1994 hybrid S-sigma (T)| 89 106 ln_s_sf12 = .false. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 90 107 ln_sigcrit = .false. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch … … 97 114 !!!!!!! SH94 stretching coefficients (ln_s_sh94 = .true.) 98 115 rn_theta = 6.0 ! surface control parameter (0<=theta<=20) 99 rn_bb = 0.8 ! stretching with SH94 s-sigma 116 rn_bb = 0.8 ! stretching with SH94 s-sigma 100 117 !!!!!!! SF12 stretching coefficient (ln_s_sf12 = .true.) 101 118 rn_alpha = 4.4 ! stretching with SF12 s-sigma … … 106 123 rn_zb_b = -0.2 ! offset for calculating Zb 107 124 !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 108 rn_thetb = 1.0 ! bottom control parameter (0<=thetb<= 1) 125 rn_thetb = 1.0 ! bottom control parameter (0<=thetb<= 1) 109 126 / 110 127 !----------------------------------------------------------------------- … … 114 131 rn_bathy = 0. ! value of the bathymetry. if (=0) bottom flat at jpkm1 115 132 nn_closea = 0 ! remove (=0) or keep (=1) closed seas and lakes (ORCA) 116 nn_msh = 0! create (=1) a mesh file or not (=0)133 nn_msh = 1 ! create (=1) a mesh file or not (=0) 117 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 118 136 rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of 119 137 rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 … … 121 139 rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0) 122 140 rn_atfp = 0.1 ! asselin time filter parameter 123 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k)124 ! =0, not used, rdt = rdttra125 rn_rdtmin = 28800. ! minimum time step on tracers (used if nn_acc=1)126 rn_rdtmax = 28800. ! maximum time step on tracers (used if nn_acc=1)127 rn_rdth = 800. ! depth variation of tracer time step (used if nn_acc=1)128 141 ln_crs = .false. ! Logical switch for coarsening module 129 142 jphgr_msh = 0 ! type of horizontal mesh … … 152 165 / 153 166 !----------------------------------------------------------------------- 154 &namsplit ! time splitting parameters ("key_dynspg_ts") 155 !----------------------------------------------------------------------- 156 ln_bt_fw = .TRUE. ! Forward integration of barotropic equations 157 ln_bt_av = .TRUE. ! Time filtering of barotropic variables 158 ln_bt_nn_auto = .TRUE. ! Set nn_baro automatically to be just below 159 ! a user defined maximum courant number (rn_bt_cmax) 160 nn_baro = 30 ! Number of iterations of barotropic mode 161 ! during rn_rdt seconds. Only used if ln_bt_nn_auto=F 162 rn_bt_cmax = 0.8 ! Maximum courant number allowed if ln_bt_nn_auto=T 163 nn_bt_flt = 1 ! Time filter choice 164 ! = 0 None 165 ! = 1 Boxcar over nn_baro barotropic steps 166 ! = 2 Boxcar over 2*nn_baro " " 167 / 168 !----------------------------------------------------------------------- 169 &namcrs ! Grid coarsening for dynamics output and/or 170 ! passive tracer coarsened online simulations 167 &namcrs ! coarsened grid (for outputs and/or TOP) ("key_crs") 171 168 !----------------------------------------------------------------------- 172 169 nn_factx = 3 ! Reduction factor of x-direction … … 190 187 / 191 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 !----------------------------------------------------------------------- 192 206 &namtsd ! data : Temperature & Salinity 193 !-----------------------------------------------------------------------194 207 !----------------------------------------------------------------------- 195 208 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 202 215 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T &S input data (T) or not (F) 203 216 / 217 204 218 !!====================================================================== 205 219 !! *** Surface Boundary Condition namelists *** 206 220 !!====================================================================== 207 221 !! namsbc surface boundary condition 208 !! namsbc_ana analytical formulation 209 !! namsbc_flx flux formulation 210 !! namsbc_clio CLIO bulk formulae formulation 211 !! namsbc_core CORE bulk formulae formulation 212 !! namsbc_mfs MFS bulk formulae formulation 213 !! namsbc_cpl CouPLed formulation ("key_ coupled")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" ) 214 228 !! namsbc_sas StAndalone Surface module 215 !! namtra_qsr penetrative solar radiation 216 !! namsbc_rnf river runoffs 217 !! namsbc_apr Atmospheric Pressure 218 !! 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) 219 235 !! namsbc_alb albedo parameters 236 !! namsbc_wave external fields from wave model (ln_wave =T) 237 !! namberg iceberg floats (ln_icebergs=T) 220 238 !!====================================================================== 221 239 ! … … 224 242 !----------------------------------------------------------------------- 225 243 nn_fsbc = 5 ! frequency of surface boundary condition computation 226 ! (also = the frequency of sea-ice model call) 244 ! (also = the frequency of sea-ice & iceberg model call) 245 ! Type of air-sea fluxes 227 246 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 228 247 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) … … 230 249 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core) 231 250 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 232 ln_cpl = .false. ! Coupled formulation (T => fill namsbc_cpl ) 233 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 251 ! Type of coupling (Ocean/Ice/Atmosphere) : 252 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) 253 ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 ) 254 nn_components = 0 ! configuration of the opa-sas OASIS coupling 255 ! =0 no opa-sas OASIS coupling: default single executable configuration 256 ! =1 opa-sas OASIS coupling: multi executable configuration, OPA component 257 ! =2 opa-sas OASIS coupling: multi executable configuration, SAS component 258 nn_limflx = -1 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 259 ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 260 ! = 0 Average per-category fluxes (forced and coupled mode) 261 ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 262 ! = 2 Redistribute a single flux over categories (coupled mode only) 263 ! Sea-ice : 234 264 nn_ice = 2 ! =0 no ice boundary condition , 235 265 ! =1 use observed ice-cover , 236 ! =2 ice-model used ("key_lim3" or "key_lim2)266 ! =2 ice-model used ("key_lim3", "key_lim2", "key_cice") 237 267 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 238 268 ! =1 levitating ice with mass and salt exchange but no presure effect 239 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 ) 240 272 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 241 273 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) … … 244 276 ! =1 global mean of e-p-r set to zero at each time step 245 277 ! =2 annual global mean of e-p-r set to zero 246 ! =3 global emp set to zero and spread out over erp area 247 ln_wave = .false. ! Activate coupling with wave (either Stokes Drift or Drag coefficient, or both) (T => fill namsbc_wave) 248 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave) 249 ln_sdw = .false. ! Computation of 3D stokes drift (T => fill namsbc_wave) 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) 250 281 nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 251 282 ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 252 cn_iceflx = 'linear' ! redistribution of solar input into ice categories during coupling ice/atm.253 283 / 254 284 !----------------------------------------------------------------------- … … 265 295 &namsbc_flx ! surface boundary condition : flux formulation 266 296 !----------------------------------------------------------------------- 267 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 297 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 268 298 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 269 299 sn_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , '' , '' … … 293 323 &namsbc_core ! namsbc_core CORE bulk formulae 294 324 !----------------------------------------------------------------------- 295 ! ! file name 296 ! ! 297 sn_wndi = 'u_10.15JUNE2009_fill' 298 sn_wndj = 'v_10.15JUNE2009_fill' 299 sn_qsr = 'ncar_rad.15JUNE2009_fill' 300 sn_qlw = 'ncar_rad.15JUNE2009_fill' 301 sn_tair = 't_10.15JUNE2009_fill' 302 sn_humi = 'q_10.15JUNE2009_fill' 303 sn_prec = 'ncar_precip.15JUNE2009_fill' 304 sn_snow = 'ncar_precip.15JUNE2009_fill' 305 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' , '' , '' 306 336 307 337 cn_dir = './' ! root directory for the location of the bulk files 308 ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F)309 338 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 310 ln_bulk2z = .false. ! Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 311 rn_zqt = 3. ! Air temperature and humidity reference height (m) (ln_bulk2z) 312 rn_zu = 4. ! Wind vector reference height (m) (ln_bulk2z) 339 rn_zqt = 10. ! Air temperature and humidity reference height (m) 340 rn_zu = 10. ! Wind vector reference height (m) 313 341 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 314 342 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) 315 rn_vfac = 0. ! multiplicative factor for ocean/ice velocity 343 rn_vfac = 0. ! multiplicative factor for ocean/ice velocity 316 344 ! in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 317 345 / … … 319 347 &namsbc_mfs ! namsbc_mfs MFS bulk formulae 320 348 !----------------------------------------------------------------------- 321 ! ! file name ! frequency (hours) ! variable 322 ! ! ! (if <0 months) ! name 323 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''324 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''325 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false. , 'daily' ,'bilinear.nc', '' ,''326 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''327 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''328 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false. , 'daily' ,'bilinear.nc', '' ,''329 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' , '' , '' 330 358 331 359 cn_dir = './ECMWF/' ! root directory for the location of the bulk files 332 360 / 333 361 !----------------------------------------------------------------------- 334 &namsbc_cpl ! coupled ocean/atmosphere model ("key_ coupled")335 !----------------------------------------------------------------------- 336 ! ! description 337 ! ! 362 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 363 !----------------------------------------------------------------------- 364 ! ! description ! multiple ! vector ! vector ! vector ! 365 ! ! ! categories ! reference ! orientation ! grids ! 338 366 ! send 339 sn_snd_temp ='weighted oce and ice' , 'no' , '' , '' , ''340 sn_snd_alb ='weighted ice' , 'no' , '' , '' , ''341 sn_snd_thick = 'none' , 'no', '' , '' , ''342 sn_snd_crt ='none' , 'no' , 'spherical' , 'eastward-northward' , 'T'343 sn_snd_co2 ='coupled' , 'no' , '' , '' , ''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' , '' , '' , '' 344 372 ! receive 345 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 346 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 347 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 348 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 349 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 350 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 351 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 352 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 353 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 354 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 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' , '' , '' , '' 383 ! 384 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentialy sending/receiving data 385 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 386 ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 355 387 / 356 388 !----------------------------------------------------------------------- … … 359 391 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 360 392 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 361 sn_usp = 'sas_grid_U' , 120 , 'vozocrtx' , .true. , .true. , 'yearly' , '' , '' , '' 393 sn_usp = 'sas_grid_U' , 120 , 'vozocrtx' , .true. , .true. , 'yearly' , '' , '' , '' 362 394 sn_vsp = 'sas_grid_V' , 120 , 'vomecrty' , .true. , .true. , 'yearly' , '' , '' , '' 363 395 sn_tem = 'sas_grid_T' , 120 , 'sosstsst' , .true. , .true. , 'yearly' , '' , '' , '' 364 396 sn_sal = 'sas_grid_T' , 120 , 'sosaline' , .true. , .true. , 'yearly' , '' , '' , '' 365 397 sn_ssh = 'sas_grid_T' , 120 , 'sossheig' , .true. , .true. , 'yearly' , '' , '' , '' 366 367 ln_3d_uv = .true. ! specify whether we are supplying a 3D u,v field 398 sn_e3t = 'sas_grid_T' , 120 , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , '' 399 sn_frq = 'sas_grid_T' , 120 , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , '' 400 401 ln_3d_uve = .true. ! specify whether we are supplying a 3D u,v and e3 field 402 ln_read_frq = .false. ! specify whether we must read frq or not 368 403 cn_dir = './' ! root directory for the location of the bulk files are 369 404 / 370 405 !----------------------------------------------------------------------- 371 &namtra_qsr ! penetrative solar radiation 406 &namtra_qsr ! penetrative solar radiation (ln_traqsr=T) 372 407 !----------------------------------------------------------------------- 373 408 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 376 411 377 412 cn_dir = './' ! root directory for the location of the runoff files 378 ln_traqsr = .true. ! Light penetration (T) or not (F)379 413 ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration 380 414 ln_qsr_2bd = .false. ! 2 bands light penetration … … 387 421 / 388 422 !----------------------------------------------------------------------- 389 &namsbc_rnf ! runoffs namelist surface boundary condition 423 &namsbc_rnf ! runoffs namelist surface boundary condition (ln_rnf=T) 390 424 !----------------------------------------------------------------------- 391 425 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 398 432 399 433 cn_dir = './' ! root directory for the location of the runoff files 400 ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F)401 434 ln_rnf_mouth = .true. ! specific treatment at rivers mouths 402 435 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used … … 406 439 ln_rnf_tem = .false. ! read in temperature information for runoff 407 440 ln_rnf_sal = .false. ! read in salinity information for runoff 441 ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file 442 rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 443 rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 444 nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0) 445 / 446 !----------------------------------------------------------------------- 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 ! 451 ! nn_isf == 4 452 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' , '' 453 ! nn_isf == 3 454 sn_rnfisf = 'rnfisf' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' , '' 455 ! nn_isf == 2 and 3 456 sn_depmax_isf='rnfisf' , -12 ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 457 sn_depmin_isf='rnfisf' , -12 ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' 458 ! nn_isf == 2 459 sn_Leff_isf = 'rnfisf' , -12 ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 460 ! 461 ! for all case 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) 466 ! only for nn_isf = 1 or 2 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 472 ! only for nn_isf = 1 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) 475 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) 476 ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 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) 408 485 / 409 486 !----------------------------------------------------------------------- 410 487 &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk 411 488 !----------------------------------------------------------------------- 412 ! ! file name 413 ! ! ! (if <0 months) ! name ! (logical)! (T/F) ! 'monthly' ! filename ! pairing ! filename !414 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' , '' , '' , '' 415 492 416 493 cn_dir = './' ! root directory for the location of the bulk files 417 rn_pref = 101000. _wp! reference atmospheric pressure [N/m2]/494 rn_pref = 101000. ! reference atmospheric pressure [N/m2]/ 418 495 ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F) 419 496 ln_apr_obc = .false. ! inverse barometer added to OBC ssh data 420 497 / 421 498 !----------------------------------------------------------------------- 422 &namsbc_ssr ! surface boundary condition : sea surface restoring 499 &namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr=T) 423 500 !----------------------------------------------------------------------- 424 501 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 446 523 / 447 524 !----------------------------------------------------------------------- 448 &namberg ! iceberg parameters 449 !----------------------------------------------------------------------- 450 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 451 542 ln_bergdia = .true. ! Calculate budgets 452 543 nn_verbose_level = 1 ! Turn on more verbose output if level > 0 … … 455 546 ! Initial mass required for an iceberg of each class 456 547 rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 457 ! Proportion of calving mass to apportion to each class 548 ! Proportion of calving mass to apportion to each class 458 549 rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 459 550 ! Ratio between effective and real iceberg mass (non-dim) 460 ! i.e. number of icebergs represented at a point 551 ! i.e. number of icebergs represented at a point 461 552 rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 462 553 ! thickness of newly calved bergs (m) … … 467 558 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits 468 559 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 469 ln_passive_mode = .false. ! iceberg - ocean decoupling 560 ln_passive_mode = .false. ! iceberg - ocean decoupling 470 561 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no) 471 562 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 472 563 rn_test_box = 108.0, 116.0, -66.0, -58.0 473 rn_speed_limit = 0. ! CFL speed limit for a berg 474 475 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !476 ! ! ! (if <0 months) ! name ! (logical)! (T/F ) ! 'monthly' ! filename ! pairing ! filename !477 sn_icb = 'calving' , -1 , 'calvingmask', .true., .true. , 'yearly' , '' , '' , ''478 479 cn_dir = './' 564 rn_speed_limit = 0. ! CFL speed limit for a berg 565 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' , '' , '' , '' 569 570 cn_dir = './' 480 571 / 481 572 … … 484 575 !!====================================================================== 485 576 !! namlbc lateral momentum boundary condition 486 !! namcla cross land advection487 577 !! namobc open boundaries parameters ("key_obc") 488 578 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") … … 494 584 &namlbc ! lateral momentum boundary condition 495 585 !----------------------------------------------------------------------- 586 ! ! free slip ! partial slip ! no slip ! strong slip 496 587 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 497 ! free slip ! partial slip ! no slip ! strong slip 498 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical eqs. 499 / 500 !----------------------------------------------------------------------- 501 &namcla ! cross land advection 502 !----------------------------------------------------------------------- 503 nn_cla = 0 ! advection between 2 ocean pts separates by land 504 / 505 !----------------------------------------------------------------------- 506 &namobc ! open boundaries parameters ("key_obc") 507 !----------------------------------------------------------------------- 508 ln_obc_clim = .false. ! climatological obc data files (T) or not (F) 509 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) 510 ln_obc_fla = .false. ! Flather open boundary condition 511 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state 512 ! = 1 the obc data are read in 'obc.dta' files 513 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data 514 ! set to monthly if obc datafile hold 1 month of data 515 rn_dpein = 1. ! damping time scale for inflow at east open boundary 516 rn_dpwin = 1. ! - - - west - - 517 rn_dpnin = 1. ! - - - north - - 518 rn_dpsin = 1. ! - - - south - - 519 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary 520 rn_dpwob = 15. ! - - - west - - 521 rn_dpnob = 3000. ! - - - north - - 522 rn_dpsob = 15. ! - - - south - - 523 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R) 524 ! = 1 the total volume remains constant 588 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs. 525 589 / 526 590 !----------------------------------------------------------------------- … … 531 595 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 532 596 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 533 / 534 !----------------------------------------------------------------------- 535 &nam_tide ! tide parameters (#ifdef key_tide) 597 ln_chk_bathy = .FALSE. ! 598 / 599 !----------------------------------------------------------------------- 600 &nam_tide ! tide parameters ("key_tide") 536 601 !----------------------------------------------------------------------- 537 602 ln_tide_pot = .true. ! use tidal potential forcing 538 603 ln_tide_ramp = .false. ! 539 604 rdttideramp = 0. ! 540 clname(1) = 'M2' ! name of constituent 541 clname(2) = 'S2' 542 clname(3) = 'N2' 543 clname(4) = 'K1' 544 clname(5) = 'O1' 545 clname(6) = 'Q1' 546 clname(7) = 'M4' 547 clname(8) = 'K2' 548 clname(9) = 'P1' 549 clname(10) = 'Mf' 550 clname(11) = 'Mm' 605 clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg 551 606 / 552 607 !----------------------------------------------------------------------- … … 563 618 ! = 2, use tidal harmonic forcing data from files 564 619 ! = 3, use external data AND tidal harmonic forcing 565 cn_dyn3d = 'none' ! 620 cn_dyn3d = 'none' ! 566 621 nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state 567 622 ! = 1, bdy data are read in 'bdydata .nc' files 568 cn_tra = 'none' ! 623 cn_tra = 'none' ! 569 624 nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state 570 625 ! = 1, bdy data are read in 'bdydata .nc' files 571 cn_ice_lim = 'none' ! 626 cn_ice_lim = 'none' ! 572 627 nn_ice_lim_dta = 0 ! = 0, bdy data are equal to the initial state 573 628 ! = 1, bdy data are read in 'bdydata .nc' files … … 578 633 ln_tra_dmp =.false. ! open boudaries conditions for tracers 579 634 ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities 580 rn_time_dmp = 1. ! Damping time scale in days 635 rn_time_dmp = 1. ! Damping time scale in days 581 636 rn_time_dmp_out = 1. ! Outflow damping time scale 582 637 nn_rimwidth = 10 ! width of the relaxation zone … … 585 640 / 586 641 !----------------------------------------------------------------------- 587 &nambdy_dta ! open boundaries - external data("key_bdy")588 !----------------------------------------------------------------------- 589 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !590 ! ! ! (if <0 months) ! name ! (logical)! (T/F ) ! 'monthly' ! filename ! pairing ! filename !591 bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig' , .true., .false. , 'daily' , '' , '' , ''592 bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx' , .true., .false. , 'daily' , '' , '' , ''593 bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty' , .true., .false. , 'daily' , '' , '' , ''594 bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx' , .true., .false. , 'daily' , '' , '' , ''595 bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty' , .true., .false. , 'daily' , '' , '' , ''596 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true., .false. , 'daily' , '' , '' , ''597 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' , '' , '' , '' 598 653 ! for lim2 599 ! bn_frld = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true., .false. , 'daily' , '' , '' , ''600 ! bn_hicif = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true., .false. , 'daily' , '' , '' , ''601 ! 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' , '' , '' , '' 602 657 ! for lim3 603 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 604 ! bn_ht_i = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 605 ! bn_ht_s = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 606 cn_dir = 'bdydta/' 607 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. ! 608 664 / 609 665 !----------------------------------------------------------------------- 610 666 &nambdy_tide ! tidal forcing at open boundaries 611 667 !----------------------------------------------------------------------- 612 filtide = 'bdydta/amm12_bdytide_' 613 ln_bdytide_2ddta = .false. 614 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. ! 615 671 / 616 672 !!====================================================================== … … 623 679 ! 624 680 !----------------------------------------------------------------------- 625 &nambfr ! bottom friction 681 &nambfr ! bottom friction (default: linear) 626 682 !----------------------------------------------------------------------- 627 683 nn_bfr = 1 ! type of bottom friction : = 0 : free slip, = 1 : linear friction … … 631 687 rn_bfri2_max = 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T) 632 688 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 633 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T 634 ln_loglayer = .false. ! logarithmic formulation (non linear case) 689 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T 635 690 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 636 691 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) 692 rn_tfri1 = 4.e-4 ! top drag coefficient (linear case) 693 rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 694 rn_tfri2_max = 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T) 695 rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2) 696 rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T 697 ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file ) 698 rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T) 699 637 700 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) 638 / 639 !----------------------------------------------------------------------- 640 &nambbc ! bottom temperature boundary condition 641 !----------------------------------------------------------------------- 642 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 701 ln_loglayer = .false. ! logarithmic formulation (non linear case) 702 / 703 !----------------------------------------------------------------------- 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' , '' , '' , '' 709 ! 710 ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom 643 711 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 644 712 ! = 1 constant flux 645 713 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 646 714 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 647 / 648 !----------------------------------------------------------------------- 649 &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") 650 719 !----------------------------------------------------------------------- 651 720 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) … … 658 727 !! Tracer (T & S ) namelists 659 728 !!====================================================================== 660 !! nameos equation of state 661 !! namtra_adv advection scheme 662 !! namtra_ldf lateral diffusion scheme 663 !! namtra_dmp T & S newtonian damping 729 !! nameos equation of state 730 !! namtra_adv advection scheme 731 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 732 !! namtra_ldf lateral diffusion scheme 733 !! namtra_ldfeiv eddy induced velocity param. 734 !! namtra_dmp T & S newtonian damping 664 735 !!====================================================================== 665 736 ! … … 667 738 &nameos ! ocean physical parameters 668 739 !----------------------------------------------------------------------- 669 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 670 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 671 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 672 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 673 rn_alpha = 2.0e-4 ! thermal expension coefficient (nn_eos= 1 or 2) 674 rn_beta = 7.7e-4 ! saline expension coefficient (nn_eos= 2) 675 / 676 !----------------------------------------------------------------------- 677 &namtra_adv ! advection scheme for tracer 678 !----------------------------------------------------------------------- 679 ln_traadv_cen2 = .false. ! 2nd order centered scheme 680 ln_traadv_tvd = .true. ! TVD scheme 681 ln_traadv_muscl = .false. ! MUSCL scheme 682 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 683 ln_traadv_ubs = .false. ! UBS scheme 684 ln_traadv_qck = .false. ! QUICKEST scheme 685 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 686 / 687 !----------------------------------------------------------------------- 688 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 689 !----------------------------------------------------------------------- 690 ln_mle = .true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 740 nn_eos = -1 ! type of equation of state and Brunt-Vaisala frequency 741 ! =-1, TEOS-10 742 ! = 0, EOS-80 743 ! = 1, S-EOS (simplified eos) 744 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 745 ! 746 ! ! S-EOS coefficients : 747 ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 748 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 749 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) 750 rn_lambda1 = 5.9520e-2 ! cabbeling coeff in T^2 (=0 for linear eos) 751 rn_lambda2 = 7.4914e-4 ! cabbeling coeff in S^2 (=0 for linear eos) 752 rn_mu1 = 1.4970e-4 ! thermobaric coeff. in T (=0 for linear eos) 753 rn_mu2 = 1.1090e-5 ! thermobaric coeff. in S (=0 for linear eos) 754 rn_nu = 2.4341e-3 ! cabbeling coeff in T*S (=0 for linear eos) 755 / 756 !----------------------------------------------------------------------- 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 691 777 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 692 778 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 698 784 rn_rho_c_mle = 0.01 ! delta rho criterion used to calculate MLD for FK 699 785 / 700 !----------------------------------------------------------------------- -----------701 &namtra_ldf ! lateral diffusion scheme for tracers 702 !----------------------------------------------------------------------- -----------786 !----------------------------------------------------------------------- 787 &namtra_ldf ! lateral diffusion scheme for tracers (default: NO diffusion) 788 !----------------------------------------------------------------------- 703 789 ! ! Operator type: 704 ln_traldf_lap = .true. ! laplacian operator 705 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 706 793 ! ! Direction of action: 707 ln_traldf_level = .false. ! iso-level 708 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 709 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 710 ! ! Griffies parameters (all need "key_ldfslp") 711 ln_traldf_grif = .false. ! use griffies triads 712 ln_traldf_gdia = .false. ! output griffies eddy velocities 713 ln_triad_iso = .false. ! pure lateral mixing in ML 714 ln_botmix_grif = .false. ! lateral mixing on bottom 715 ! ! Coefficients 716 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 717 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 718 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 719 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 720 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 721 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 722 ! (normally=0; not used with Griffies) 723 rn_slpmax = 0.01 ! slope limit 724 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 725 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 726 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 727 / 728 !----------------------------------------------------------------------- 729 &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) 730 834 !----------------------------------------------------------------------- 731 835 ln_tradmp = .true. ! add a damping termn (T) or not (F) 732 nn_hdmp = -1 ! horizontal shape =-1, damping in Med and Red Seas only733 ! =XX, damping poleward of XX degrees (XX>0)734 ! + F(distance-to-coast) + Red and Med Seas735 836 nn_zdmp = 0 ! vertical shape =0 damping throughout the water column 736 837 ! =1 no damping in the mixing layer (kz criteria) 737 838 ! =2 no damping in the mixed layer (rho crieria) 738 rn_surf = 50. ! surface time scale of damping [days] 739 rn_bot = 360. ! bottom time scale of damping [days] 740 rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters] 741 nn_file = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) 839 cn_resto ='resto.nc' ! Name of file containing restoration coeff. field (use dmp_tools to create this) 742 840 / 743 841 … … 748 846 !! namdyn_vor advection scheme 749 847 !! namdyn_hpg hydrostatic pressure gradient 750 !! namdyn_spg surface pressure gradient (CPP key only)848 !! namdyn_spg surface pressure gradient 751 849 !! namdyn_ldf lateral diffusion scheme 752 850 !!====================================================================== 753 851 ! 754 852 !----------------------------------------------------------------------- 755 &namdyn_adv ! formulation of the momentum advection 853 &namdyn_adv ! formulation of the momentum advection (default: vector form) 756 854 !----------------------------------------------------------------------- 757 855 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 856 nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction 758 857 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme 759 858 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 760 / 761 !----------------------------------------------------------------------- 762 &nam_vvl ! vertical coordinate options 763 !----------------------------------------------------------------------- 764 ln_vvl_zstar = .true. ! zstar vertical coordinate 859 ln_dynzad_zts = .false. ! Use (T) sub timestepping for vertical momentum advection 860 / 861 !----------------------------------------------------------------------- 862 &nam_vvl ! vertical coordinate options (default: zstar) 863 !----------------------------------------------------------------------- 864 ln_vvl_zstar = .true. ! zstar vertical coordinate 765 865 ln_vvl_ztilde = .false. ! ztilde vertical coordinate: only high frequency variations 766 866 ln_vvl_layer = .false. ! full layer vertical coordinate … … 774 874 / 775 875 !----------------------------------------------------------------------- 776 &namdyn_vor ! option of physics/algorithm (not control by CPP keys)876 &namdyn_vor ! option of physics/algorithm (default: NO) 777 877 !----------------------------------------------------------------------- 778 878 ln_dynvor_ene = .false. ! enstrophy conserving scheme 779 879 ln_dynvor_ens = .false. ! energy conserving scheme 780 880 ln_dynvor_mix = .false. ! mixed scheme 781 ln_dynvor_een = .true. ! energy & enstrophy scheme 782 / 783 !----------------------------------------------------------------------- 784 &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) 785 887 !----------------------------------------------------------------------- 786 888 ln_hpg_zco = .false. ! z-coordinate - full steps 787 889 ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation) 788 890 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 891 ln_hpg_isf = .false. ! s-coordinate (sco ) adapted to isf 789 892 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 790 893 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 791 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 792 ! centered time scheme (F) 793 / 794 !----------------------------------------------------------------------- 795 !namdyn_spg ! surface pressure gradient (CPP key only) 796 !----------------------------------------------------------------------- 797 ! ! explicit free surface ("key_dynspg_exp") 798 ! ! filtered free surface ("key_dynspg_flt") 799 ! ! split-explicit free surface ("key_dynspg_ts") 800 801 !----------------------------------------------------------------------- 802 &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) 803 911 !----------------------------------------------------------------------- 804 912 ! ! Type of the operator : 805 ln_dynldf_lap = .true. ! laplacian operator 806 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 807 916 ! ! Direction of action : 808 ln_dynldf_lev el = .false.! iso-level809 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)810 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 811 920 ! ! Coefficient 812 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 813 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 814 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 815 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 816 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 817 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 818 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 819 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) 820 934 / 821 935 … … 826 940 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 827 941 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 828 !! namzdf_ kpp KPP dependent vertical mixing ("key_zdfkpp")942 !! namzdf_gls GLS vertical mixing ("key_zdfgls") 829 943 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 830 944 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 858 972 rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer 859 973 rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer 860 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. 861 975 / 862 976 !----------------------------------------------------------------------- … … 882 996 ! = 1 add a tke source below the ML 883 997 ! = 2 add a tke source just at the base of the ML 884 ! = 3 as = 1 applied on HF part of the stress ("key_ coupled")998 ! = 3 as = 1 applied on HF part of the stress ("key_oasis3") 885 999 rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 886 1000 nn_htau = 1 ! type of exponential decrease of tke penetration below the ML … … 888 1002 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 889 1003 / 890 !------------------------------------------------------------------------891 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:892 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")893 ln_kpprimix = .true. ! shear instability mixing894 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]895 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]896 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability897 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]898 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]899 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]900 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv901 nn_ave = 1 ! constant (=0) or profile (=1) background on avt902 /903 1004 !----------------------------------------------------------------------- 904 1005 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") 905 1006 !----------------------------------------------------------------------- 906 rn_emin = 1.e- 6! minimum value of e [m2/s2]1007 rn_emin = 1.e-7 ! minimum value of e [m2/s2] 907 1008 rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] 908 1009 ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) 909 rn_clim_galp = 0.53 ! galperin limit 910 ln_crban = .true. ! Use Craig & Banner (1994) surface wave mixing parametrisation 1010 rn_clim_galp = 0.267 ! galperin limit 911 1011 ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case 912 1012 rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux 913 1013 rn_charn = 70000. ! Charnock constant for wb induced roughness length 914 nn_tkebc_surf = 1 ! surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 915 nn_tkebc_bot = 1 ! bottom tke condition (0/1=Dir/Neum) 916 nn_psibc_surf = 1 ! surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 917 nn_psibc_bot = 1 ! bottom psi condition (0/1=Dir/Neum) 918 nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 919 nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 1014 rn_hsro = 0.02 ! Minimum surface roughness 1015 rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met=2) 1016 nn_z0_met = 2 ! Method for surface roughness computation (0/1/2) 1017 nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum) 1018 nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum) 1019 nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 1020 nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 920 1021 / 921 1022 !----------------------------------------------------------------------- … … 941 1042 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 942 1043 !! namctl Control prints & Benchmark 943 !! namsol elliptic solver / island / free surface 944 !!====================================================================== 945 ! 946 !----------------------------------------------------------------------- 947 &namsol ! elliptic solver / island / free surface 948 !----------------------------------------------------------------------- 949 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg) 950 ! =2 successive-over-relaxation (sor) 951 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test 952 rn_eps = 1.e-6 ! absolute precision of the solver 953 nn_nmin = 300 ! minimum of iterations for the SOR solver 954 nn_nmax = 800 ! maximum of iterations for the SOR solver 955 nn_nmod = 10 ! frequency of test for the SOR solver 956 rn_resmax = 1.e-10 ! absolute precision for the SOR solver 957 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain) 958 / 1044 !! namsto Stochastic parametrization of EOS 1045 !!====================================================================== 1046 ! 959 1047 !----------------------------------------------------------------------- 960 1048 &nammpp ! Massively Parallel Processing ("key_mpp_mpi) … … 982 1070 ! (no physical validity of the results) 983 1071 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 984 / 985 !----------------------------------------------------------------------- 986 &namc1d_uvd ! data: U & V currents ("key_c1d") 987 !----------------------------------------------------------------------- 988 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 989 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 990 sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' 991 sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' 992 ! 993 cn_dir = './' ! root directory for the location of the files 994 ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F) 995 ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F) 996 / 997 !----------------------------------------------------------------------- 998 &namc1d_dyndmp ! U & V newtonian damping ("key_c1d") 999 !----------------------------------------------------------------------- 1000 ln_dyndmp = .false. ! add a damping term (T) or not (F) 1001 / 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 ln_rstseed = .true. ! read seed of RNG from restart file 1087 cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input) 1088 cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output) 1089 / 1090 1002 1091 !!====================================================================== 1003 1092 !! *** Diagnostics namelists *** 1004 1093 !!====================================================================== 1005 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 1006 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 1007 !! namflo float parameters ("key_float") 1094 !! namtrd dynamics and/or tracer trends 1008 1095 !! namptr Poleward Transport Diagnostics 1009 1096 !! namhsb Heat and salt budgets 1010 !!====================================================================== 1011 ! 1012 !----------------------------------------------------------------------- 1013 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 1014 !----------------------------------------------------------------------- 1015 nn_nchunks_i= 4 ! number of chunks in i-dimension 1016 nn_nchunks_j= 4 ! number of chunks in j-dimension 1017 nn_nchunks_k= 31 ! number of chunks in k-dimension 1018 ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 1019 ! is optimal for postprocessing which works exclusively with horizontal slabs 1020 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 1021 ! (F) ignore chunking information and produce netcdf3-compatible files 1022 / 1023 !----------------------------------------------------------------------- 1024 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 1025 ! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor") 1026 !----------------------------------------------------------------------- 1027 nn_trd = 365 ! time step frequency dynamics and tracers trends 1028 nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 1029 rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 1030 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 1031 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 1032 ln_trdmld_restart = .false. ! restart for ML diagnostics 1033 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 1097 !! namflo float parameters ("key_float") 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 1106 !----------------------------------------------------------------------- 1107 ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE 1108 ln_dyn_trd = .false. ! (T) 3D momentum trend output 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) 1111 ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends 1112 ln_PE_trd = .false. ! (T) 3D Potential Energy trends 1113 ln_tra_trd = .false. ! (T) 3D tracer trend output 1114 ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 1115 nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step) 1116 / 1117 !!gm nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 1118 !!gm rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 1119 !!gm cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 1120 !!gm cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 1121 !!gm ln_trdmld_restart = .false. ! restart for ML diagnostics 1122 !!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 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) 1034 1134 / 1035 1135 !----------------------------------------------------------------------- … … 1048 1148 / 1049 1149 !----------------------------------------------------------------------- 1050 &namptr ! Poleward Transport Diagnostic 1051 !----------------------------------------------------------------------- 1052 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1053 ln_diaznl = .true. ! Add zonal means and meridional stream functions 1054 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 1055 ! (orca configuration only, need input basins mask file named "subbasins.nc" 1056 ln_ptrcomp = .true. ! Add decomposition : overturning 1057 nn_fptr = 1 ! Frequency of ptr computation [time step] 1058 nn_fwri = 15 ! Frequency of ptr outputs [time step] 1059 / 1060 !----------------------------------------------------------------------- 1061 &namhsb ! Heat and salt budgets 1062 !----------------------------------------------------------------------- 1063 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 1064 / 1065 !----------------------------------------------------------------------- 1066 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 1150 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 1067 1151 !----------------------------------------------------------------------- 1068 1152 nit000_han = 1 ! First time step used for harmonic analysis … … 1073 1157 / 1074 1158 !----------------------------------------------------------------------- 1075 &namdct ! transports through s ections1159 &namdct ! transports through some sections 1076 1160 !----------------------------------------------------------------------- 1077 1161 nn_dct = 15 ! time step frequency for transports computing … … 1081 1165 ! 0 < n : debug section number n 1082 1166 / 1083 1084 !!====================================================================== 1085 !! *** Observation & Assimilation namelists *** 1086 !!====================================================================== 1087 !! 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 1088 1183 !! nam_asminc assimilation increments ('key_asminc') 1089 1184 !!====================================================================== 1090 1185 ! 1091 1186 !----------------------------------------------------------------------- 1092 &namobs ! observation usage switch ('key_diaobs') 1093 !----------------------------------------------------------------------- 1094 ln_t3d = .false. ! Logical switch for T profile observations 1095 ln_s3d = .false. ! Logical switch for S profile observations 1096 ln_ena = .false. ! Logical switch for ENACT insitu data set 1097 ! ! ln_cor Logical switch for Coriolis insitu data set 1098 ln_profb = .false. ! Logical switch for feedback insitu data set 1099 ln_sla = .false. ! Logical switch for SLA observations 1100 1101 ln_sladt = .false. ! Logical switch for AVISO SLA data 1102 1103 ln_slafb = .false. ! Logical switch for feedback SLA data 1104 ! ln_ssh Logical switch for SSH observations 1105 1106 ln_sst = .true. ! Logical switch for SST observations 1107 ln_reysst = .true. ! ln_reysst Logical switch for Reynolds observations 1108 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 1109 1110 ln_sstfb = .false. ! Logical switch for feedback SST data 1111 ! ln_sss Logical switch for SSS observations 1112 ! ln_seaice Logical switch for Sea Ice observations 1113 ! ln_vel3d Logical switch for velocity observations 1114 ! ln_velavcur Logical switch for velocity daily av. cur. 1115 ! ln_velhrcur Logical switch for velocity high freq. cur. 1116 ! ln_velavadcp Logical switch for velocity daily av. ADCP 1117 ! ln_velhradcp Logical switch for velocity high freq. ADCP 1118 ! ln_velfb Logical switch for feedback velocity data 1119 ! ln_grid_global Global distribtion of observations 1120 ! ln_grid_search_lookup Logical switch for obs grid search w/lookup table 1121 ! grid_search_file Grid search lookup file header 1122 ! enactfiles ENACT input observation file names 1123 ! coriofiles Coriolis input observation file name 1124 ! ! profbfiles: Profile feedback input observation file name 1125 profbfiles = 'profiles_01.nc' 1126 ! ln_profb_enatim Enact feedback input time setting switch 1127 ! slafilesact Active SLA input observation file name 1128 ! slafilespas Passive SLA input observation file name 1129 ! ! slafbfiles: Feedback SLA input observation file name 1130 slafbfiles = 'sla_01.nc' 1131 ! sstfiles GHRSST input observation file name 1132 ! ! sstfbfiles: Feedback SST input observation file name 1133 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 1134 ! seaicefiles Sea Ice input observation file name 1135 ! velavcurfiles Vel. cur. daily av. input file name 1136 ! velhvcurfiles Vel. cur. high freq. input file name 1137 ! velavadcpfiles Vel. ADCP daily av. input file name 1138 ! velhvadcpfiles Vel. ADCP high freq. input file name 1139 ! velfbfiles Vel. feedback input observation file name 1140 ! dobsini Initial date in window YYYYMMDD.HHMMSS 1141 ! dobsend Final date in window YYYYMMDD.HHMMSS 1142 ! n1dint Type of vertical interpolation method 1143 ! n2dint Type of horizontal interpolation method 1144 ! ln_nea Rejection of observations near land switch 1145 nmsshc = 0 ! MSSH correction scheme 1146 ! mdtcorr MDT correction 1147 ! mdtcutoff MDT cutoff for computed correction 1148 ln_altbias = .false. ! Logical switch for alt bias 1149 ln_ignmis = .true. ! Logical switch for ignoring missing files 1150 ! endailyavtypes ENACT daily average types 1151 ln_grid_global = .true. 1152 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' 1153 1222 / 1154 1223 !----------------------------------------------------------------------- … … 1171 1240 / 1172 1241 !----------------------------------------------------------------------- 1173 &namsbc_wave ! External fields from wave model 1174 !----------------------------------------------------------------------- 1175 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 1176 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 1177 sn_cdg = 'cdg_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' , '' 1178 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1179 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1180 sn_wn = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' , '' , '' , '' 1181 ! 1182 cn_dir_cdg = './' ! root directory for the location of drag coefficient files 1183 / 1184 !----------------------------------------------------------------------- 1185 &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) 1186 !----------------------------------------------------------------------- 1187 ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 1188 ln_neptsimp = .false. ! yes/no use simplified neptune 1189 1190 ln_smooth_neptvel = .false. ! yes/no smooth zunep, zvnep 1191 rn_tslse = 1.2e4 ! value of lengthscale L at the equator 1192 rn_tslsp = 3.0e3 ! value of lengthscale L at the pole 1193 ! Specify whether to ramp down the Neptune velocity in shallow 1194 ! water, and if so the depth range controlling such ramping down 1195 ln_neptramp = .true. ! ramp down Neptune velocity in shallow water 1196 rn_htrmin = 100.0 ! min. depth of transition range 1197 rn_htrmax = 200.0 ! max. depth of transition range 1198 / 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.