- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/CONFIG/SHARED/namelist_ref
- Property svn:keywords set to Id
r7350 r7351 3 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 4 4 !! NEMO/OPA : 1 - run manager (namrun) 5 !! namelists 2 - Domain (namcfg, namzgr, namzgr_sco, namdom, namtsd )5 !! namelists 2 - Domain (namcfg, namzgr, namzgr_sco, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 6 6 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 7 7 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 8 !! namsbc_apr, namsbc_ssr, namsbc_alb )8 !! namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave) 9 9 !! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 10 10 !! 5 - bottom boundary (nambfr, nambbc, nambbl) … … 31 31 nn_itend = 5475 ! last time step (std 5475) 32 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 33 34 nn_leapy = 0 ! Leap year calendar (1) or not (0) 34 35 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 35 nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T 36 nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T 37 ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 38 ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 39 ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart 40 cn_ocerst_in = "restart" ! suffix of ocean restart name (input) 41 cn_ocerst_indir = "." ! directory from which to read input ocean restarts 42 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 43 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 44 46 nn_istate = 0 ! output the initial state (1) or not (0) 45 47 ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) … … 47 49 nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 48 50 nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000) 49 ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T)50 51 ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) 51 52 ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard 52 ln_clobber = .true. ! clobber (overwrite) an existing file53 ln_clobber = .true. ! clobber (overwrite) an existing file 53 54 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 54 55 / … … 58 59 !!====================================================================== 59 60 !! namcfg parameters of the configuration 60 !! namzgr vertical coordinate 61 !! namzgr vertical coordinate (default: NO selection) 61 62 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 62 63 !! namdom space and time domain (bathymetry, mesh, timestep) 64 !! namwad Wetting and drying (default F) 63 65 !! namtsd data: temperature & salinity 64 !!====================================================================== 65 ! 66 !----------------------------------------------------------------------- 67 &namcfg ! parameters of the configuration 68 !----------------------------------------------------------------------- 69 cp_cfg = "default" ! name of the configuration 70 cp_cfz = "no zoom" ! name of the zoom of configuration 71 jp_cfg = 0 ! resolution of the configuration 72 jpidta = 10 ! 1st lateral dimension ( >= jpi ) 73 jpjdta = 12 ! 2nd " " ( >= jpj ) 74 jpkdta = 31 ! number of levels ( >= jpk ) 75 jpiglo = 10 ! 1st dimension of global domain --> i =jpidta 76 jpjglo = 12 ! 2nd - - --> j =jpjdta 77 jpizoom = 1 ! left bottom (i,j) indices of the zoom 78 jpjzoom = 1 ! in data domain indices 79 jperio = 0 ! lateral cond. type (between 0 and 6) 80 ! = 0 closed ; = 1 cyclic East-West 81 ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 82 ! = 4 cyclic East-West AND North fold T-point pivot 83 ! = 5 North fold F-point pivot 84 ! = 6 cyclic East-West AND North fold F-point pivot 85 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 86 ! in netcdf input files, as the start j-row for reading 87 / 88 !----------------------------------------------------------------------- 89 &namzgr ! vertical coordinate 90 !----------------------------------------------------------------------- 91 ln_zco = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) 92 ln_zps = .true. ! z-coordinate - partial steps (T/F) 93 ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F) 94 ln_isfcav = .false. ! ice shelf cavity (T/F) 95 / 96 !----------------------------------------------------------------------- 97 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate 98 !----------------------------------------------------------------------- 99 ln_s_sh94 = .true. ! Song & Haidvogel 1994 hybrid S-sigma (T)| 66 !! namcrs coarsened grid (for outputs and/or TOP) ("key_crs") 67 !! namc1d 1D configuration options ("key_c1d") 68 !! namc1d_dyndmp 1D newtonian damping applied on currents ("key_c1d") 69 !! namc1d_uvd 1D data (currents) ("key_c1d") 70 !!====================================================================== 71 ! 72 !----------------------------------------------------------------------- 73 &namcfg ! parameters of the configuration 74 !----------------------------------------------------------------------- 75 cp_cfg = "default" ! name of the configuration 76 cp_cfz = "no zoom" ! name of the zoom of configuration 77 jp_cfg = 0 ! resolution of the configuration 78 jpidta = 10 ! 1st lateral dimension ( >= jpi ) 79 jpjdta = 12 ! 2nd " " ( >= jpj ) 80 jpkdta = 31 ! number of levels ( >= jpk ) 81 jpiglo = 10 ! 1st dimension of global domain --> i =jpidta 82 jpjglo = 12 ! 2nd - - --> j =jpjdta 83 jpizoom = 1 ! left bottom (i,j) indices of the zoom 84 jpjzoom = 1 ! in data domain indices 85 jperio = 0 ! lateral cond. type (between 0 and 6) 86 ! = 0 closed ; = 1 cyclic East-West 87 ! = 2 equatorial symmetric ; = 3 North fold T-point pivot 88 ! = 4 cyclic East-West AND North fold T-point pivot 89 ! = 5 North fold F-point pivot 90 ! = 6 cyclic East-West AND North fold F-point pivot 91 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 92 ! in netcdf input files, as the start j-row for reading 93 / 94 !----------------------------------------------------------------------- 95 &namzgr ! vertical coordinate (default: NO selection) 96 !----------------------------------------------------------------------- 97 ln_zco = .false. ! z-coordinate - full steps 98 ln_zps = .false. ! z-coordinate - partial steps 99 ln_sco = .false. ! s- or hybrid z-s-coordinate 100 ln_isfcav = .false. ! ice shelf cavity 101 ln_linssh = .false. ! linear free surface 102 / 103 !----------------------------------------------------------------------- 104 &namzgr_sco ! s-coordinate or hybrid z-s-coordinate (default F) 105 !----------------------------------------------------------------------- 106 ln_s_sh94 = .false. ! Song & Haidvogel 1994 hybrid S-sigma (T)| 100 107 ln_s_sf12 = .false. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 101 108 ln_sigcrit = .false. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch … … 127 134 nn_msh = 1 ! create (=1) a mesh file or not (=0) 128 135 rn_hmin = -3. ! min depth of the ocean (>0) or min number of ocean level (<0) 136 rn_isfhmin = 1.00 ! treshold (m) to discriminate grounding ice to floating ice 129 137 rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of 130 138 rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 … … 132 140 rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0) 133 141 rn_atfp = 0.1 ! asselin time filter parameter 134 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k)135 ! =0, not used, rdt = rdttra136 rn_rdtmin = 28800. ! minimum time step on tracers (used if nn_acc=1)137 rn_rdtmax = 28800. ! maximum time step on tracers (used if nn_acc=1)138 rn_rdth = 800. ! depth variation of tracer time step (used if nn_acc=1)139 142 ln_crs = .false. ! Logical switch for coarsening module 140 143 jphgr_msh = 0 ! type of horizontal mesh … … 163 166 / 164 167 !----------------------------------------------------------------------- 165 &namcrs ! Grid coarsening for dynamics output and/or 166 ! passive tracer coarsened online simulations 168 &namwad ! Wetting and drying (default F) 169 !----------------------------------------------------------------------- 170 ln_wd = .false. ! T/F activation of wetting and drying 171 rn_wdmin1 = 0.1 ! Minimum wet depth on dried cells 172 rn_wdmin2 = 0.01 ! Tolerance of min wet depth on dried cells 173 rn_wdld = 20.0 ! Land elevation below which wetting/drying is allowed 174 nn_wdit = 10 ! Max iterations for W/D limiter 175 / 176 !----------------------------------------------------------------------- 177 &namtsd ! data : Temperature & Salinity 178 !----------------------------------------------------------------------- 179 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 180 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 181 sn_tem = 'data_1m_potential_temperature_nomask', -1 ,'votemper', .true. , .true. , 'yearly' , '' , '' , '' 182 sn_sal = 'data_1m_salinity_nomask' , -1 ,'vosaline', .true. , .true. , 'yearly' , '' , '' , '' 183 ! 184 cn_dir = './' ! root directory for the location of the runoff files 185 ln_tsd_init = .true. ! Initialisation of ocean T & S with T & S input data (T) or not (F) 186 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T & S input data (T) or not (F) 187 / 188 !----------------------------------------------------------------------- 189 &namcrs ! coarsened grid (for outputs and/or TOP) ("key_crs") 167 190 !----------------------------------------------------------------------- 168 191 nn_factx = 3 ! Reduction factor of x-direction … … 186 209 / 187 210 !----------------------------------------------------------------------- 188 &namtsd ! data : Temperature & Salinity 189 !----------------------------------------------------------------------- 190 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 191 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 192 sn_tem = 'data_1m_potential_temperature_nomask', -1 ,'votemper' , .true. , .true. , 'yearly' , '' , '' , '' 193 sn_sal = 'data_1m_salinity_nomask' , -1 ,'vosaline' , .true. , .true. , 'yearly' , '' , '' , '' 194 ! 195 cn_dir = './' ! root directory for the location of the runoff files 196 ln_tsd_init = .true. ! Initialisation of ocean T & S with T &S input data (T) or not (F) 197 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T &S input data (T) or not (F) 198 / 211 &namc1d_dyndmp ! U & V newtonian damping ("key_c1d") 212 !----------------------------------------------------------------------- 213 ln_dyndmp = .false. ! add a damping term (T) or not (F) 214 / 215 !----------------------------------------------------------------------- 216 &namc1d_uvd ! data: U & V currents ("key_c1d") 217 !----------------------------------------------------------------------- 218 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 219 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 220 sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' 221 sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' 222 ! 223 cn_dir = './' ! root directory for the location of the files 224 ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F) 225 ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F) 226 / 227 199 228 !!====================================================================== 200 229 !! *** Surface Boundary Condition namelists *** 201 230 !!====================================================================== 202 231 !! namsbc surface boundary condition 203 !! namsbc_ana analytical formulation 204 !! namsbc_flx flux formulation 205 !! namsbc_clio CLIO bulk formulae formulation 206 !! namsbc_core CORE bulk formulae formulation 207 !! namsbc_mfs MFS bulk formulae formulation 208 !! namsbc_cpl CouPLed formulation ("key_oasis3" )232 !! namsbc_ana analytical formulation (ln_ana =T) 233 !! namsbc_flx flux formulation (ln_flx =T) 234 !! namsbc_clio CLIO bulk formulae formulation (ln_blk_clio=T) 235 !! namsbc_core CORE bulk formulae formulation (ln_blk_core=T) 236 !! namsbc_mfs MFS bulk formulae formulation (ln_blk_mfs =T) 237 !! namsbc_cpl CouPLed formulation ("key_oasis3" ) 209 238 !! namsbc_sas StAndalone Surface module 210 !! namtra_qsr penetrative solar radiation 211 !! namsbc_rnf river runoffs 212 !! namsbc_isf ice shelf melting/freezing 213 !! namsbc_apr Atmospheric Pressure 214 !! namsbc_ssr sea surface restoring term (for T and/or S) 239 !! namtra_qsr penetrative solar radiation (ln_traqsr =T) 240 !! namsbc_rnf river runoffs (ln_rnf =T) 241 !! namsbc_isf ice shelf melting/freezing (nn_isf >0) 242 !! namsbc_iscpl coupling option between land ice model and ocean 243 !! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T) 244 !! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T) 215 245 !! namsbc_alb albedo parameters 246 !! namsbc_wave external fields from wave model (ln_wave =T) 247 !! namberg iceberg floats (ln_icebergs=T) 216 248 !!====================================================================== 217 249 ! … … 220 252 !----------------------------------------------------------------------- 221 253 nn_fsbc = 5 ! frequency of surface boundary condition computation 222 ! (also = the frequency of sea-ice model call) 254 ! (also = the frequency of sea-ice & iceberg model call) 255 ! Type of air-sea fluxes 223 256 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana ) 224 257 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) … … 226 259 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core) 227 260 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 261 ! Type of coupling (Ocean/Ice/Atmosphere) : 228 262 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) 229 263 ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 ) … … 232 266 ! =1 opa-sas OASIS coupling: multi executable configuration, OPA component 233 267 ! =2 opa-sas OASIS coupling: multi executable configuration, SAS component 234 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr )235 nn_ice = 2 ! =0 no ice boundary condition ,236 ! =1 use observed ice-cover ,237 ! =2 ice-model used ("key_lim3" or "key_lim2")238 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect)239 ! =1 levitating ice with mass and salt exchange but no presure effect240 ! =2 embedded sea-ice (full salt and mass exchanges and pressure)241 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave242 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf)243 nn_isf = 0 ! ice shelf melting/freezing (/=0 => fill namsbc_isf)244 ! 0 =no isf 1 = presence of ISF245 ! 2 = bg03 parametrisation 3 = rnf file for isf246 ! 4 = ISF fwf specified247 ! option 1 and 4 need ln_isfcav = .true. (domzgr)248 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)249 nn_fwb = 2 ! FreshWater Budget: =0 unchecked250 ! =1 global mean of e-p-r set to zero at each time step251 ! =2 annual global mean of e-p-r set to zero252 ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave)253 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)254 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)255 ln_tauoc= .false. ! Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)256 ln_stcor= .false. ! Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)257 nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,258 ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)259 268 nn_limflx = -1 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 260 269 ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled … … 262 271 ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 263 272 ! = 2 Redistribute a single flux over categories (coupled mode only) 273 ! Sea-ice : 274 nn_ice = 2 ! =0 no ice boundary condition , 275 ! =1 use observed ice-cover , 276 ! =2 ice-model used ("key_lim3", "key_lim2", "key_cice") 277 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 278 ! =1 levitating ice with mass and salt exchange but no presure effect 279 ! =2 embedded sea-ice (full salt and mass exchanges and pressure) 280 ! Misc. options of sbc : 281 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr ) 282 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 283 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 284 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 285 nn_fwb = 2 ! FreshWater Budget: =0 unchecked 286 ! =1 global mean of e-p-r set to zero at each time step 287 ! =2 annual global mean of e-p-r set to zero 288 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 289 ln_isf = .false. ! ice shelf (T => fill namsbc_isf) 290 ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave) 291 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 292 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 293 ln_tauoc = .false. ! Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave) 294 ln_stcor = .false. ! Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) 295 nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 296 ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 264 297 / 265 298 !----------------------------------------------------------------------- … … 304 337 &namsbc_core ! namsbc_core CORE bulk formulae 305 338 !----------------------------------------------------------------------- 306 ! ! file name 307 ! ! 308 sn_wndi = 'u_10.15JUNE2009_fill' 309 sn_wndj = 'v_10.15JUNE2009_fill' 310 sn_qsr = 'ncar_rad.15JUNE2009_fill' 311 sn_qlw = 'ncar_rad.15JUNE2009_fill' 312 sn_tair = 't_10.15JUNE2009_fill' 313 sn_humi = 'q_10.15JUNE2009_fill' 314 sn_prec = 'ncar_precip.15JUNE2009_fill' 315 sn_snow = 'ncar_precip.15JUNE2009_fill' 316 sn_tdif = 'taudif_core' 339 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 340 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 341 sn_wndi = 'u_10.15JUNE2009_fill' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Uwnd' , '' 342 sn_wndj = 'v_10.15JUNE2009_fill' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Vwnd' , '' 343 sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 344 sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 345 sn_tair = 't_10.15JUNE2009_fill' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 346 sn_humi = 'q_10.15JUNE2009_fill' , 6 , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 347 sn_prec = 'ncar_precip.15JUNE2009_fill', -1 , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 348 sn_snow = 'ncar_precip.15JUNE2009_fill', -1 , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 349 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 317 350 318 351 cn_dir = './' ! root directory for the location of the bulk files 319 352 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 320 rn_zqt = 10. 321 rn_zu = 10. 353 rn_zqt = 10. ! Air temperature and humidity reference height (m) 354 rn_zu = 10. ! Wind vector reference height (m) 322 355 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 323 356 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) … … 328 361 &namsbc_mfs ! namsbc_mfs MFS bulk formulae 329 362 !----------------------------------------------------------------------- 330 ! ! file name ! frequency (hours) ! variable 331 ! ! ! (if <0 months) ! name 332 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''333 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''334 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false. , 'daily' ,'bilinear.nc', '' ,''335 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''336 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false. , 'daily' ,'bicubic.nc' , '' ,''337 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false. , 'daily' ,'bilinear.nc', '' ,''338 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' ,''363 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 364 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 365 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 366 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 367 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false., 'daily' ,'bilinear.nc', '' , '' 368 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 369 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false., 'daily' ,'bicubic.nc' , '' , '' 370 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false., 'daily' ,'bilinear.nc', '' , '' 371 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' , '' 339 372 340 373 cn_dir = './ECMWF/' ! root directory for the location of the bulk files … … 343 376 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 344 377 !----------------------------------------------------------------------- 345 ! ! description 346 ! ! 378 ! ! description ! multiple ! vector ! vector ! vector ! 379 ! ! ! categories ! reference ! orientation ! grids ! 347 380 ! send 348 sn_snd_temp = 349 sn_snd_alb = 350 sn_snd_thick = 'none' , 'no', '' , '' , ''351 sn_snd_crt = 352 sn_snd_co2 = 381 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' 382 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' 383 sn_snd_thick = 'none' , 'no' , '' , '' , '' 384 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 385 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 353 386 sn_snd_crtw = 'none' , 'no' , '' , '' , 'U,V' 354 387 sn_snd_ifrac = 'none' , 'no' , '' , '' , '' 355 388 sn_snd_wlev = 'coupled' , 'no' , '' , '' , '' 356 389 ! receive 357 sn_rcv_w10m = 358 sn_rcv_taumod = 359 sn_rcv_tau = 360 sn_rcv_dqnsdt = 361 sn_rcv_qsr = 362 sn_rcv_qns = 363 sn_rcv_emp = 364 sn_rcv_rnf = 365 sn_rcv_cal = 366 sn_rcv_co2 = 390 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 391 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 392 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 393 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 394 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 395 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 396 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 397 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 398 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 399 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 367 400 sn_rcv_hsig = 'none' , 'no' , '' , '' , '' 368 401 sn_rcv_iceflx = 'none' , 'no' , '' , '' , '' … … 376 409 sn_rcv_wdrag = 'none' , 'no' , '' , '' , '' 377 410 ! 378 nn_cplmodel = 1 379 ln_usecplmask = .false. 380 411 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentialy sending/receiving data 412 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 413 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 381 414 / 382 415 !----------------------------------------------------------------------- 383 416 &namsbc_sas ! analytical surface boundary condition 384 417 !----------------------------------------------------------------------- 385 ! ! file name ! frequency (hours) ! variable ! time interp. 386 ! ! ! (if <0 months) ! name ! 387 sn_usp = 'sas_grid_U' , 120 , 'vozocrtx' , .true. , .true. , 'yearly' , '' , '' ,''388 sn_vsp = 'sas_grid_V' , 120 , 'vomecrty' , .true. , .true. , 'yearly' , '' , '' ,''389 sn_tem = 'sas_grid_T' , 120 , 'sosstsst' , .true. , .true. , 'yearly' , '' , '' ,''390 sn_sal = 'sas_grid_T' , 120 , 'sosaline' , .true. , .true. , 'yearly' , '' , '' ,''391 sn_ssh = 'sas_grid_T' , 120 , 'sossheig' , .true. , .true. , 'yearly' , '' , '' ,''392 sn_e3t = 'sas_grid_T' , 120 , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' ,''393 sn_frq = 'sas_grid_T' , 120 , 'frq_m' , .true. , .true. , 'yearly' , '' , '' ,''418 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 419 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 420 sn_usp = 'sas_grid_U', 120 , 'vozocrtx', .true. , .true. , 'yearly' , '' , '' , '' 421 sn_vsp = 'sas_grid_V', 120 , 'vomecrty', .true. , .true. , 'yearly' , '' , '' , '' 422 sn_tem = 'sas_grid_T', 120 , 'sosstsst', .true. , .true. , 'yearly' , '' , '' , '' 423 sn_sal = 'sas_grid_T', 120 , 'sosaline', .true. , .true. , 'yearly' , '' , '' , '' 424 sn_ssh = 'sas_grid_T', 120 , 'sossheig', .true. , .true. , 'yearly' , '' , '' , '' 425 sn_e3t = 'sas_grid_T', 120 , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , '' 426 sn_frq = 'sas_grid_T', 120 , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , '' 394 427 395 428 ln_3d_uve = .true. ! specify whether we are supplying a 3D u,v and e3 field 396 ln_read_frq = .false. 429 ln_read_frq = .false. ! specify whether we must read frq or not 397 430 cn_dir = './' ! root directory for the location of the bulk files are 398 431 / 399 432 !----------------------------------------------------------------------- 400 &namtra_qsr ! penetrative solar radiation 433 &namtra_qsr ! penetrative solar radiation (ln_traqsr=T) 401 434 !----------------------------------------------------------------------- 402 435 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 405 438 406 439 cn_dir = './' ! root directory for the location of the runoff files 407 ln_traqsr = .true. ! Light penetration (T) or not (F)408 440 ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration 409 441 ln_qsr_2bd = .false. ! 2 bands light penetration … … 416 448 / 417 449 !----------------------------------------------------------------------- 418 &namsbc_rnf ! runoffs namelist surface boundary condition 450 &namsbc_rnf ! runoffs namelist surface boundary condition (ln_rnf=T) 419 451 !----------------------------------------------------------------------- 420 452 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 426 458 sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , '' 427 459 428 cn_dir 429 ln_rnf_mouth 430 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used431 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s]432 rn_rfact 433 ln_rnf_depth 434 ln_rnf_tem 435 ln_rnf_sal 436 ln_rnf_depth_ini = .false. 437 rn_rnf_max= 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )438 rn_dep_max= 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true )439 nn_rnf_depth_file = 0! create (=1) a runoff depth file or not (=0)440 / 441 !----------------------------------------------------------------------- 442 &namsbc_isf ! Top boundary layer (ISF) 443 !----------------------------------------------------------------------- 444 ! ! file name ! frequency (hours) ! variable ! time interp ol. ! clim ! 'yearly'/ ! weights ! rotation!445 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing!460 cn_dir = './' ! root directory for the location of the runoff files 461 ln_rnf_mouth= .true. ! specific treatment at rivers mouths 462 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T) 463 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 464 rn_rfact = 1.e0 ! multiplicative factor for runoff 465 ln_rnf_depth= .false. ! read in depth information for runoff 466 ln_rnf_tem = .false. ! read in temperature information for runoff 467 ln_rnf_sal = .false. ! read in salinity information for runoff 468 ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file 469 rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 470 rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 471 nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0) 472 / 473 !----------------------------------------------------------------------- 474 &namsbc_isf ! Top boundary layer (ISF) (nn_isf >0) 475 !----------------------------------------------------------------------- 476 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 477 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 446 478 ! nn_isf == 4 447 sn_qisf = 'rnfisf' , -12 ,'sohflisf', .false. , .true. , 'yearly' , '' , '' 448 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' 479 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' , '' 449 480 ! nn_isf == 3 450 sn_rnfisf = 'runoffs' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' ,''481 sn_rnfisf = 'rnfisf' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' , '' 451 482 ! nn_isf == 2 and 3 452 sn_depmax_isf = 'runoffs' , -12 ,'sozisfmax' , .false. , .true. , 'yearly' , '' ,''453 sn_depmin_isf = 'runoffs' , -12 ,'sozisfmin' , .false. , .true. , 'yearly' , '' ,''483 sn_depmax_isf='rnfisf' , -12 ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 484 sn_depmin_isf='rnfisf' , -12 ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' 454 485 ! nn_isf == 2 455 sn_Leff_isf = 'rnfisf' , 0 ,'Leff' , .false. , .true. , 'yearly' , '' , '' 486 sn_Leff_isf = 'rnfisf' , -12 ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 487 ! 456 488 ! for all case 457 ln_divisf = .true. ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) 489 nn_isf = 1 ! ice shelf melting/freezing 490 ! 1 = presence of ISF 2 = bg03 parametrisation 491 ! 3 = rnf file for isf 4 = ISF fwf specified 492 ! option 1 and 4 need ln_isfcav = .true. (domzgr) 458 493 ! only for nn_isf = 1 or 2 459 rn_gammat0 = 1.0e-4 ! gammat coefficient used in blk formula 460 rn_gammas0 = 1.0e-4 ! gammas coefficient used in blk formula 494 rn_gammat0 = 1.e-4 ! gammat coefficient used in blk formula 495 rn_gammas0 = 1.e-4 ! gammas coefficient used in blk formula 496 ! only for nn_isf = 1 or 4 497 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 498 ! ! 0 => thickness of the tbl = thickness of the first wet cell 461 499 ! only for nn_isf = 1 462 nn_isfblk = 1 ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 463 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 464 ! 0 => thickness of the tbl = thickness of the first wet cell 465 ln_conserve = .true. ! conservative case (take into account meltwater advection) 466 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) 467 ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 468 ! if you want to keep the cd as in global config, adjust rn_gammat0 to compensate 469 ! 2 = velocity and stability dependent Gamma Holland et al. 1999 500 nn_isfblk = 1 ! 1 ISOMIP like: 2 equations formulation (Hunter et al., 2006) 501 ! ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 502 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) 503 ! ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 504 ! ! 2 = velocity and stability dependent Gamma (Holland et al. 1999) 505 / 506 !----------------------------------------------------------------------- 507 &namsbc_iscpl ! land ice / ocean coupling option 508 !----------------------------------------------------------------------- 509 nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells) 510 ln_hsb = .false. ! activate conservation module (conservation exact after a time of rn_fiscpl) 511 nn_fiscpl = 43800 ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 470 512 / 471 513 !----------------------------------------------------------------------- 472 514 &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk 473 515 !----------------------------------------------------------------------- 474 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !475 ! ! ! (if <0 months) ! name !(logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !476 sn_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' ,''477 478 cn_dir = './' 479 rn_pref = 101000. 480 ln_ref_apr = .false. 481 ln_apr_obc = .false. 482 / 483 !----------------------------------------------------------------------- 484 &namsbc_ssr ! surface boundary condition : sea surface restoring 485 !----------------------------------------------------------------------- 486 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !487 ! ! ! (if <0 months) ! name ! (logical)! (T/F) ! 'monthly' ! filename ! pairing ! filename !488 sn_sst = 'sst_data' , 24 , 'sst' , .false. , .false., 'yearly' , '' , '' ,''489 sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' ,''516 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 517 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 518 sn_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' , '' 519 520 cn_dir = './' ! root directory for the location of the bulk files 521 rn_pref = 101000. ! reference atmospheric pressure [N/m2]/ 522 ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F) 523 ln_apr_obc = .false. ! inverse barometer added to OBC ssh data 524 / 525 !----------------------------------------------------------------------- 526 &namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr=T) 527 !----------------------------------------------------------------------- 528 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 529 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 530 sn_sst = 'sst_data', 24 , 'sst' , .false. , .false., 'yearly' , '' , '' , '' 531 sn_sss = 'sss_data', -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' , '' 490 532 491 533 cn_dir = './' ! root directory for the location of the runoff files … … 495 537 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] 496 538 rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] 497 ln_sssr_bnd = .true.! flag to bound erp term (associated with nn_sssr=2)539 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 498 540 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 499 541 / … … 501 543 &namsbc_alb ! albedo parameters 502 544 !----------------------------------------------------------------------- 503 rn_cloud = 0.06 ! cloud correction to snow and ice albedo 504 rn_albice = 0.53 ! albedo of melting ice in the arctic and antarctic 505 rn_alphd = 0.80 ! coefficients for linear interpolation used to 506 rn_alphc = 0.65 ! compute albedo between two extremes values 507 rn_alphdi = 0.72 ! (Pyane, 1972) 508 / 509 !----------------------------------------------------------------------- 510 &namberg ! iceberg parameters 511 !----------------------------------------------------------------------- 512 ln_icebergs = .false. 545 nn_ice_alb = 0 ! parameterization of ice/snow albedo 546 ! 0: Shine & Henderson-Sellers (JGR 1985) 547 ! 1: "home made" based on Brandt et al. (J. Climate 2005) 548 ! and Grenfell & Perovich (JGR 2004) 549 rn_albice = 0.53 ! albedo of bare puddled ice (values from 0.49 to 0.58) 550 ! 0.53 (default) => if nn_ice_alb=0 551 ! 0.50 (default) => if nn_ice_alb=1 552 / 553 !----------------------------------------------------------------------- 554 &namsbc_wave ! External fields from wave model 555 !----------------------------------------------------------------------- 556 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 557 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 558 sn_cdg = 'sdw_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' , '' 559 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 560 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 561 sn_swh = 'sdw_wave' , 1 , 'hs' , .true. , .false. , 'daily' , '' , '' , '' 562 sn_wmp = 'sdw_wave' , 1 , 'wmp' , .true. , .false. , 'daily' , '' , '' , '' 563 sn_wnum = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' , '' , '' , '' 564 sn_tauoc = 'sdw_wave' , 1 , 'wave_stress', .true. , .false. , 'daily' , '' , '' , '' 565 ! 566 cn_dir = './' ! root directory for the location of drag coefficient files 567 / 568 !----------------------------------------------------------------------- 569 &namberg ! iceberg parameters (default: No iceberg) 570 !----------------------------------------------------------------------- 571 ln_icebergs = .false. ! iceberg floats or not 513 572 ln_bergdia = .true. ! Calculate budgets 514 573 nn_verbose_level = 1 ! Turn on more verbose output if level > 0 … … 535 594 rn_speed_limit = 0. ! CFL speed limit for a berg 536 595 537 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !538 ! ! ! (if <0 months) ! name ! (logical)! (T/F ) ! 'monthly' ! filename ! pairing ! filename !539 sn_icb = 'calving' , -1 , 'calvingmask', .true. , .true. , 'yearly' , '' , '' ,''596 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 597 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 598 sn_icb = 'calving', -1 , 'calvingmask', .true. , .true. , 'yearly' , '' , '' , '' 540 599 541 600 cn_dir = './' … … 546 605 !!====================================================================== 547 606 !! namlbc lateral momentum boundary condition 548 !! namobc open boundaries parameters ("key_obc")549 607 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") 608 !! nam_tide Tidal forcing 550 609 !! nambdy Unstructured open boundaries ("key_bdy") 551 !! namtide Tidal forcing at open boundaries ("key_bdy_tides") 610 !! nambdy_dta Unstructured open boundaries - external data ("key_bdy") 611 !! nambdy_tide tidal forcing at open boundaries ("key_bdy_tides") 552 612 !!====================================================================== 553 613 ! … … 555 615 &namlbc ! lateral momentum boundary condition 556 616 !----------------------------------------------------------------------- 617 ! ! free slip ! partial slip ! no slip ! strong slip 557 618 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. 619 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs. 560 620 / 561 621 !----------------------------------------------------------------------- … … 569 629 / 570 630 !----------------------------------------------------------------------- 571 &nam_tide ! tide parameters (#ifdef key_tide)572 !----------------------------------------------------------------------- 573 ln_tide_pot = .true.! use tidal potential forcing574 ln_tide_ramp = .false.!575 rdttideramp = 0.!576 clname(1) = 'DUMMY'! name of constituent - all tidal components must be set in namelist_cfg631 &nam_tide ! tide parameters ("key_tide") 632 !----------------------------------------------------------------------- 633 ln_tide_pot = .true. ! use tidal potential forcing 634 ln_tide_ramp= .false. ! 635 rdttideramp = 0. ! 636 clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg 577 637 / 578 638 !----------------------------------------------------------------------- … … 611 671 / 612 672 !----------------------------------------------------------------------- 613 &nambdy_dta ! open boundaries - external data("key_bdy")614 !----------------------------------------------------------------------- 615 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !616 ! ! ! (if <0 months) ! name ! (logical)! (T/F ) ! 'monthly' ! filename ! pairing ! filename !617 bn_ssh = 'amm12_bdyT_u2d' , 24 , 'sossheig' , .true. , .false. , 'daily' , '' , '' ,''618 bn_u2d = 'amm12_bdyU_u2d' , 24 , 'vobtcrtx' , .true. , .false. , 'daily' , '' , '' ,''619 bn_v2d = 'amm12_bdyV_u2d' , 24 , 'vobtcrty' , .true. , .false. , 'daily' , '' , '' ,''620 bn_u3d = 'amm12_bdyU_u3d' , 24 , 'vozocrtx' , .true. , .false. , 'daily' , '' , '' ,''621 bn_v3d = 'amm12_bdyV_u3d' , 24 , 'vomecrty' , .true. , .false. , 'daily' , '' , '' ,''622 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true. , .false. , 'daily' , '' , '' ,''623 bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline' , .true. , .false. , 'daily' , '' , '' ,''673 &nambdy_dta ! open boundaries - external data ("key_bdy") 674 !----------------------------------------------------------------------- 675 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 676 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 677 bn_ssh = 'amm12_bdyT_u2d', 24 , 'sossheig', .true. , .false. , 'daily' , '' , '' , '' 678 bn_u2d = 'amm12_bdyU_u2d', 24 , 'vobtcrtx', .true. , .false. , 'daily' , '' , '' , '' 679 bn_v2d = 'amm12_bdyV_u2d', 24 , 'vobtcrty', .true. , .false. , 'daily' , '' , '' , '' 680 bn_u3d = 'amm12_bdyU_u3d', 24 , 'vozocrtx', .true. , .false. , 'daily' , '' , '' , '' 681 bn_v3d = 'amm12_bdyV_u3d', 24 , 'vomecrty', .true. , .false. , 'daily' , '' , '' , '' 682 bn_tem = 'amm12_bdyT_tra', 24 , 'votemper', .true. , .false. , 'daily' , '' , '' , '' 683 bn_sal = 'amm12_bdyT_tra', 24 , 'vosaline', .true. , .false. , 'daily' , '' , '' , '' 624 684 ! for lim2 625 ! bn_frld = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' ,''626 ! bn_hicif = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' ,''627 ! bn_hsnif = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' ,''685 ! bn_frld = 'amm12_bdyT_ice', 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , '' 686 ! bn_hicif = 'amm12_bdyT_ice', 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , '' 687 ! bn_hsnif = 'amm12_bdyT_ice', 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , '' 628 688 ! for lim3 629 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 630 ! bn_ht_i = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 631 ! bn_ht_s = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 632 cn_dir = 'bdydta/' 633 ln_full_vel = .false. 634 / 635 !----------------------------------------------------------------------- 636 &nambdy_tide ! tidal forcing at open boundaries 637 !----------------------------------------------------------------------- 638 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 639 ln_bdytide_2ddta = .false. 640 ln_bdytide_conj = .false. 641 / 689 ! bn_a_i = 'amm12_bdyT_ice', 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , '' 690 ! bn_ht_i = 'amm12_bdyT_ice', 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , '' 691 ! bn_ht_s = 'amm12_bdyT_ice', 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , '' 692 693 cn_dir = 'bdydta/' ! root directory for the location of the bulk files 694 ln_full_vel = .false. ! 695 / 696 !----------------------------------------------------------------------- 697 &nambdy_tide ! tidal forcing at open boundaries 698 !----------------------------------------------------------------------- 699 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 700 ln_bdytide_2ddta = .false. ! 701 ln_bdytide_conj = .false. ! 702 / 703 642 704 !!====================================================================== 643 705 !! *** Bottom boundary condition *** … … 649 711 ! 650 712 !----------------------------------------------------------------------- 651 &nambfr ! bottom friction 713 &nambfr ! bottom friction (default: linear) 652 714 !----------------------------------------------------------------------- 653 715 nn_bfr = 1 ! type of bottom friction : = 0 : free slip, = 1 : linear friction … … 655 717 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 656 718 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 657 rn_bfri2_max =1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T)719 rn_bfri2_max= 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T) 658 720 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 659 721 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T … … 662 724 rn_tfri1 = 4.e-4 ! top drag coefficient (linear case) 663 725 rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 664 rn_tfri2_max =1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T)726 rn_tfri2_max= 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T) 665 727 rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2) 666 728 rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T 667 729 ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file ) 668 rn_tfrien = 50.! local multiplying factor of tfr (ln_tfr2d=T)730 rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T) 669 731 670 732 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) … … 672 734 / 673 735 !----------------------------------------------------------------------- 674 &nambbc ! bottom temperature boundary condition 675 !----------------------------------------------------------------------- 676 ! ! ! (if <0 months) ! 677 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 678 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 679 sn_qgh ='geothermal_heating.nc', -12. , 'heatflow' , .false. , .true. , 'yearly' , '' , '' , '' 736 &nambbc ! bottom temperature boundary condition (default: NO) 737 !----------------------------------------------------------------------- 738 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 739 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! 740 sn_qgh ='geothermal_heating.nc', -12. , 'heatflow', .false. , .true. , 'yearly' , '' , '' , '' 680 741 ! 681 cn_dir = './' ! root directory for the location of the runoff files 682 ln_trabbc = .true. ! Apply a geothermal heating at the ocean bottom 742 ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom 683 743 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux 684 744 ! = 1 constant flux 685 745 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 686 746 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 687 / 688 !----------------------------------------------------------------------- 689 &nambbl ! bottom boundary layer scheme 690 !----------------------------------------------------------------------- 691 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) 692 nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0) 693 rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s] 694 rn_gambbl = 10. ! advective bbl coefficient [s] 747 cn_dir = './' ! root directory for the location of the runoff files 748 / 749 !----------------------------------------------------------------------- 750 &nambbl ! bottom boundary layer scheme ("key_trabbl") 751 !----------------------------------------------------------------------- 752 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) 753 nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0) 754 rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s] 755 rn_gambbl = 10. ! advective bbl coefficient [s] 695 756 / 696 757 … … 709 770 &nameos ! ocean physical parameters 710 771 !----------------------------------------------------------------------- 711 nn_eos = -1 ! type of equation of state and Brunt-Vaisala frequency 712 ! =-1, TEOS-10 713 ! = 0, EOS-80 714 ! = 1, S-EOS (simplified eos) 715 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 772 ln_teos10 = .false. ! = Use TEOS-10 equation of state 773 ln_eos80 = .false. ! = Use EOS80 equation of state 774 ln_seos = .false. ! = Use simplified equation of state (S-EOS) 716 775 ! 717 ! ! S-EOS coefficients :718 776 ! ! S-EOS coefficients (ln_seos=T): 777 ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 719 778 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 720 779 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 726 785 / 727 786 !----------------------------------------------------------------------- 728 &namtra_adv ! advection scheme for tracer 729 !----------------------------------------------------------------------- 730 ln_traadv_cen = .false.! 2nd order centered scheme731 nn_cen_h = 4 732 nn_cen_v = 4 733 ln_traadv_fct = .false.! FCT scheme734 nn_fct_h = 2 735 nn_fct_v = 2 736 nn_fct_zts = 0 737 ! 738 ln_traadv_mus = .false.! MUSCL scheme739 ln_mus_ups = .false.! use upstream scheme near river mouths740 ln_traadv_ubs = .false.! UBS scheme741 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT742 ln_traadv_qck = .false.! QUICKEST scheme743 / 744 !----------------------------------------------------------------------- 745 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 746 !----------------------------------------------------------------------- 747 ln_mle = .false.! (T) use the Mixed Layer Eddy (MLE) parameterisation748 rn_ce = 0.06! magnitude of the MLE (typical value: 0.06 to 0.08)749 nn_mle = 1! MLE type: =0 standard Fox-Kemper ; =1 new formulation750 rn_lf = 5.e+3! typical scale of mixed layer front (meters) (case rn_mle=0)751 rn_time = 172800.! time scale for mixing momentum across the mixed layer (seconds) (case rn_mle=0)752 rn_lat = 20.! reference latitude (degrees) of MLE coef. (case rn_mle=1)753 nn_mld_uv = 0! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)754 nn_conv = 0! =1 no MLE in case of convection ; =0 always MLE755 rn_rho_c_mle = 0.01! delta rho criterion used to calculate MLD for FK756 / 757 !----------------------------------------------------------------------- -----------758 &namtra_ldf ! lateral diffusion scheme for tracers 759 !----------------------------------------------------------------------- -----------787 &namtra_adv ! advection scheme for tracer (default: NO advection) 788 !----------------------------------------------------------------------- 789 ln_traadv_cen = .false. ! 2nd order centered scheme 790 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 791 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 792 ln_traadv_fct = .false. ! FCT scheme 793 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 794 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 795 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 796 ! ! (number of sub-timestep = nn_fct_zts) 797 ln_traadv_mus = .false. ! MUSCL scheme 798 ln_mus_ups = .false. ! use upstream scheme near river mouths 799 ln_traadv_ubs = .false. ! UBS scheme 800 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT / COMPACT 4th order 801 ln_traadv_qck = .false. ! QUICKEST scheme 802 / 803 !----------------------------------------------------------------------- 804 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) 805 !----------------------------------------------------------------------- 806 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 807 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 808 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 809 rn_lf = 5.e+3 ! typical scale of mixed layer front (meters) (case rn_mle=0) 810 rn_time = 172800. ! time scale for mixing momentum across the mixed layer (seconds) (case rn_mle=0) 811 rn_lat = 20. ! reference latitude (degrees) of MLE coef. (case rn_mle=1) 812 nn_mld_uv = 0 ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 813 nn_conv = 0 ! =1 no MLE in case of convection ; =0 always MLE 814 rn_rho_c_mle= 0.01 ! delta rho criterion used to calculate MLD for FK 815 / 816 !----------------------------------------------------------------------- 817 &namtra_ldf ! lateral diffusion scheme for tracers (default: NO diffusion) 818 !----------------------------------------------------------------------- 760 819 ! ! Operator type: 761 820 ! ! no diffusion: set ln_traldf_lap=..._blp=F 762 821 ln_traldf_lap = .false. ! laplacian operator 763 822 ln_traldf_blp = .false. ! bilaplacian operator 823 ! 764 824 ! ! Direction of action: 765 825 ln_traldf_lev = .false. ! iso-level … … 787 847 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 788 848 / 789 !----------------------------------------------------------------------- -----------790 &namtra_ldfeiv ! eddy induced velocity param. 791 !----------------------------------------------------------------------- -----------792 ln_ldfeiv =.false. 793 ln_ldfeiv_dia =.false. 794 rn_aeiv_0 = 2000. 795 nn_aei_ijk_t = 21 849 !----------------------------------------------------------------------- 850 &namtra_ldfeiv ! eddy induced velocity param. (default: NO) 851 !----------------------------------------------------------------------- 852 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 853 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 854 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 855 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 796 856 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 797 857 ! ! = 0 constant … … 802 862 / 803 863 !----------------------------------------------------------------------- 804 &namtra_dmp ! tracer: T & S newtonian damping 864 &namtra_dmp ! tracer: T & S newtonian damping (default: NO) 805 865 !----------------------------------------------------------------------- 806 866 ln_tradmp = .true. ! add a damping termn (T) or not (F) … … 808 868 ! =1 no damping in the mixing layer (kz criteria) 809 869 ! =2 no damping in the mixed layer (rho crieria) 810 cn_resto = 'resto.nc' ! Name of file containing restoration coefficientfield (use dmp_tools to create this)870 cn_resto ='resto.nc' ! Name of file containing restoration coeff. field (use dmp_tools to create this) 811 871 / 812 872 … … 822 882 ! 823 883 !----------------------------------------------------------------------- 824 &namdyn_adv ! formulation of the momentum advection 884 &namdyn_adv ! formulation of the momentum advection (default: vector form) 825 885 !----------------------------------------------------------------------- 826 886 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) … … 831 891 / 832 892 !----------------------------------------------------------------------- 833 &nam_vvl ! vertical coordinate options 893 &nam_vvl ! vertical coordinate options (default: zstar) 834 894 !----------------------------------------------------------------------- 835 895 ln_vvl_zstar = .true. ! zstar vertical coordinate … … 837 897 ln_vvl_layer = .false. ! full layer vertical coordinate 838 898 ln_vvl_ztilde_as_zstar = .false. ! ztilde vertical coordinate emulating zstar 839 ln_vvl_zstar_at_eqtor = .false.! ztilde near the equator899 ln_vvl_zstar_at_eqtor = .false. ! ztilde near the equator 840 900 rn_ahe3 = 0.0e0 ! thickness diffusion coefficient 841 901 rn_rst_e3t = 30.e0 ! ztilde to zstar restoration timescale [days] … … 845 905 / 846 906 !----------------------------------------------------------------------- 847 &namdyn_vor ! option of physics/algorithm (not control by CPP keys)907 &namdyn_vor ! Vorticity / Coriolis scheme (default: NO) 848 908 !----------------------------------------------------------------------- 849 909 ln_dynvor_ene = .false. ! enstrophy conserving scheme … … 851 911 ln_dynvor_mix = .false. ! mixed scheme 852 912 ln_dynvor_een = .false. ! energy & enstrophy scheme 853 nn_een_e3f = 1 !e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)854 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT USE855 / 856 !----------------------------------------------------------------------- 857 &namdyn_hpg ! Hydrostatic pressure gradient option 913 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 914 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT ACTIVATE 915 / 916 !----------------------------------------------------------------------- 917 &namdyn_hpg ! Hydrostatic pressure gradient option (default: zps) 858 918 !----------------------------------------------------------------------- 859 919 ln_hpg_zco = .false. ! z-coordinate - full steps … … 865 925 / 866 926 !----------------------------------------------------------------------- 867 &namdyn_spg ! surface pressure gradient 927 &namdyn_spg ! surface pressure gradient (default: NO) 868 928 !----------------------------------------------------------------------- 869 929 ln_dynspg_exp = .false. ! explicit free surface 870 930 ln_dynspg_ts = .false. ! split-explicit free surface 871 ln_bt_fw = .true. ! Forward integration of barotropic Eqs.872 ln_bt_av = .true. ! Time filtering of barotropic variables873 nn_bt_flt = 1! Time filter choice = 0 None874 ! ! = 1 Boxcar over nn_baro sub-steps875 ! ! = 2 Boxcar over 2*nn_baro " "876 ln_bt_auto = .true. ! Number of sub-step defined from:877 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed878 nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds879 / 880 !----------------------------------------------------------------------- 881 &namdyn_ldf ! lateral diffusion on momentum 931 ln_bt_fw = .true. ! Forward integration of barotropic Eqs. 932 ln_bt_av = .true. ! Time filtering of barotropic variables 933 nn_bt_flt = 1 ! Time filter choice = 0 None 934 ! ! = 1 Boxcar over nn_baro sub-steps 935 ! ! = 2 Boxcar over 2*nn_baro " " 936 ln_bt_auto = .true. ! Number of sub-step defined from: 937 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed 938 nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds 939 / 940 !----------------------------------------------------------------------- 941 &namdyn_ldf ! lateral diffusion on momentum (default: NO) 882 942 !----------------------------------------------------------------------- 883 943 ! ! Type of the operator : … … 911 971 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 912 972 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 973 !! namzdf_gls GLS vertical mixing ("key_zdfgls") 913 974 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 914 975 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 923 984 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 924 985 ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F) 925 nn_evdm = 0 !evd apply on tracer (=0) or on tracer and momentum (=1)926 rn_avevd = 100.! evd mixing coefficient [m2/s]986 nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1) 987 rn_avevd = 100. ! evd mixing coefficient [m2/s] 927 988 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F) 928 nn_npc = 1 !frequency of application of npc929 nn_npcp = 365 !npc control print frequency989 nn_npc = 1 ! frequency of application of npc 990 nn_npcp = 365 ! npc control print frequency 930 991 ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping 931 nn_zdfexp = 3 !number of sub-timestep for ln_zdfexp=T992 nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 932 993 ln_zdfqiao = .false. ! Enhanced wave vertical mixing Qiao (2010) (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) 933 994 / … … 935 996 &namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" ) 936 997 !----------------------------------------------------------------------- 937 rn_avmri = 100.e-4! maximum value of the vertical viscosity938 rn_alp = 5.! coefficient of the parameterization939 nn_ric = 2! coefficient of the parameterization940 rn_ekmfc = 0.7! Factor in the Ekman depth Equation941 rn_mldmin = 1.0! minimum allowable mixed-layer depth estimate (m)942 rn_mldmax = 1000.0! maximum allowable mixed-layer depth estimate (m)943 rn_wtmix = 10.0! vertical eddy viscosity coeff [m2/s] in the mixed-layer944 rn_wvmix = 10.0! vertical eddy diffusion coeff [m2/s] in the mixed-layer945 ln_mldw = .true. ! Flag to use or not the mized layer depth param.998 rn_avmri = 100.e-4 ! maximum value of the vertical viscosity 999 rn_alp = 5. ! coefficient of the parameterization 1000 nn_ric = 2 ! coefficient of the parameterization 1001 rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation 1002 rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m) 1003 rn_mldmax = 1000.0 ! maximum allowable mixed-layer depth estimate (m) 1004 rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer 1005 rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer 1006 ln_mldw = .true. ! Flag to use or not the mixed layer depth param. 946 1007 / 947 1008 !----------------------------------------------------------------------- … … 963 1024 ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002) 964 1025 rn_lc = 0.15 ! coef. associated to Langmuir cells 965 nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to internal &intertial waves1026 nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to near intertial waves 966 1027 ! = 0 no penetration 967 1028 ! = 1 add a tke source below the ML 968 1029 ! = 2 add a tke source just at the base of the ML 969 ! = 3 as = 1 applied on HF part of the stress ("key_oasis3")1030 ! = 3 as = 1 applied on HF part of the stress (ln_cpl=T) 970 1031 rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 971 1032 nn_htau = 1 ! type of exponential decrease of tke penetration below the ML … … 974 1035 / 975 1036 !----------------------------------------------------------------------- 976 &namzdf_gls ! GLS vertical diffusion("key_zdfgls")1037 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") 977 1038 !----------------------------------------------------------------------- 978 1039 rn_emin = 1.e-7 ! minimum value of e [m2/s2] … … 1007 1068 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency 1008 1069 / 1070 !----------------------------------------------------------------------- 1071 &namzdf_tmx_new ! internal wave-driven mixing parameterization ("key_zdftmx_new" & "key_zdfddm") 1072 !----------------------------------------------------------------------- 1073 nn_zpyc = 1 ! pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) 1074 ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency 1075 ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F) 1076 / 1077 1009 1078 1010 1079 !!====================================================================== … … 1013 1082 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 1014 1083 !! namctl Control prints & Benchmark 1015 !! namc1d 1D configuration options ("key_c1d")1016 !! namc1d_uvd data: U & V currents ("key_c1d")1017 !! namc1d_dyndmp U & V newtonian damping ("key_c1d")1018 1084 !! namsto Stochastic parametrization of EOS 1019 1085 !!====================================================================== … … 1044 1110 ! (no physical validity of the results) 1045 1111 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 1046 / 1047 !----------------------------------------------------------------------- 1048 &namc1d_uvd ! data: U & V currents ("key_c1d") 1049 !----------------------------------------------------------------------- 1050 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 1051 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 1052 sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' 1053 sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' 1054 ! 1055 cn_dir = './' ! root directory for the location of the files 1056 ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F) 1057 ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F) 1058 / 1059 !----------------------------------------------------------------------- 1060 &namc1d_dyndmp ! U & V newtonian damping ("key_c1d") 1061 !----------------------------------------------------------------------- 1062 ln_dyndmp = .false. ! add a damping term (T) or not (F) 1063 / 1064 !----------------------------------------------------------------------- 1065 &namsto ! Stochastic parametrization of EOS 1066 !----------------------------------------------------------------------- 1067 ln_rststo = .false. ! start from mean parameter (F) or from restart file (T) 1068 ln_rstseed = .true. ! read seed of RNG from restart file 1112 nn_diacfl = 0 ! Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 1113 / 1114 !----------------------------------------------------------------------- 1115 &namsto ! Stochastic parametrization of EOS (default: NO) 1116 !----------------------------------------------------------------------- 1117 ln_sto_eos = .false. ! stochastic equation of state 1118 nn_sto_eos = 1 ! number of independent random walks 1119 rn_eos_stdxy= 1.4 ! random walk horz. standard deviation (in grid points) 1120 rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points) 1121 rn_eos_tcor = 1440. ! random walk time correlation (in timesteps) 1122 nn_eos_ord = 1 ! order of autoregressive processes 1123 nn_eos_flt = 0 ! passes of Laplacian filter 1124 rn_eos_lim = 2.0 ! limitation factor (default = 3.0) 1125 ln_rststo = .false. ! start from mean parameter (F) or from restart file (T) 1126 ln_rstseed = .true. ! read seed of RNG from restart file 1069 1127 cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input) 1070 1128 cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output) 1071 1072 ln_sto_eos = .false. ! stochastic equation of state1073 nn_sto_eos = 1 ! number of independent random walks1074 rn_eos_stdxy = 1.4 ! random walk horz. standard deviation (in grid points)1075 rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points)1076 rn_eos_tcor = 1440.0 ! random walk time correlation (in timesteps)1077 nn_eos_ord = 1 ! order of autoregressive processes1078 nn_eos_flt = 0 ! passes of Laplacian filter1079 rn_eos_lim = 2.0 ! limitation factor (default = 3.0)1080 1129 / 1081 1130 … … 1083 1132 !! *** Diagnostics namelists *** 1084 1133 !!====================================================================== 1134 !! namtrd dynamics and/or tracer trends (default F) 1135 !! namptr Poleward Transport Diagnostics (default F) 1136 !! namhsb Heat and salt budgets (default F) 1137 !! namdiu Cool skin and warm layer models (default F) 1138 !! namflo float parameters ("key_float") 1139 !! nam_diaharm Harmonic analysis of tidal constituents ("key_diaharm") 1140 !! namdct transports through some sections ("key_diadct") 1141 !! nam_diatmb Top Middle Bottom Output (default F) 1142 !! nam_dia25h 25h Mean Output (default F) 1085 1143 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 1086 !! namtrd dynamics and/or tracer trends 1087 !! namptr Poleward Transport Diagnostics 1088 !! namflo float parameters ("key_float") 1089 !! namhsb Heat and salt budgets 1090 !!====================================================================== 1091 ! 1092 !----------------------------------------------------------------------- 1093 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 1094 !----------------------------------------------------------------------- 1095 nn_nchunks_i= 4 ! number of chunks in i-dimension 1096 nn_nchunks_j= 4 ! number of chunks in j-dimension 1097 nn_nchunks_k= 31 ! number of chunks in k-dimension 1098 ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 1099 ! is optimal for postprocessing which works exclusively with horizontal slabs 1100 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 1101 ! (F) ignore chunking information and produce netcdf3-compatible files 1102 / 1103 !----------------------------------------------------------------------- 1104 &namtrd ! diagnostics on dynamics and/or tracer trends 1105 ! ! and/or mixed-layer trends and/or barotropic vorticity 1144 !!====================================================================== 1145 ! 1146 !----------------------------------------------------------------------- 1147 &namtrd ! trend diagnostics (default F) 1106 1148 !----------------------------------------------------------------------- 1107 1149 ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE 1108 1150 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)1151 ln_dyn_mxl = .false. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 1152 ln_vor_trd = .false. ! (T) 2D barotropic vorticity trends (not coded yet) 1111 1153 ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends 1112 1154 ln_PE_trd = .false. ! (T) 3D Potential Energy trends 1113 ln_tra_trd = . FALSE. ! (T) 3D tracer trend output1155 ln_tra_trd = .false. ! (T) 3D tracer trend output 1114 1156 ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 1115 1157 nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step) … … 1123 1165 !!gm 1124 1166 !----------------------------------------------------------------------- 1125 &namflo ! float parameters ("key_float") 1126 !----------------------------------------------------------------------- 1127 jpnfl = 1 ! total number of floats during the run 1128 jpnnewflo = 0 ! number of floats for the restart 1129 ln_rstflo = .false. ! float restart (T) or not (F) 1130 nn_writefl = 75 ! frequency of writing in float output file 1131 nn_stockfl = 5475 ! frequency of creation of the float restart file 1132 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 1133 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 1134 ! or computed with Blanke' scheme (F) 1135 ln_ariane = .true. ! Input with Ariane tool convention(T) 1136 ln_flo_ascii = .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 1137 / 1138 !----------------------------------------------------------------------- 1139 &namptr ! Poleward Transport Diagnostic 1140 !----------------------------------------------------------------------- 1141 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1142 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 1143 / 1144 !----------------------------------------------------------------------- 1145 &namhsb ! Heat and salt budgets 1146 !----------------------------------------------------------------------- 1147 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 1148 / 1149 !----------------------------------------------------------------------- 1150 &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') 1167 &namptr ! Poleward Transport Diagnostic (default F) 1168 !----------------------------------------------------------------------- 1169 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1170 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 1171 / 1172 !----------------------------------------------------------------------- 1173 &namhsb ! Heat and salt budgets (default F) 1174 !----------------------------------------------------------------------- 1175 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 1176 / 1177 !----------------------------------------------------------------------- 1178 &namdiu ! Cool skin and warm layer models (default F) 1179 !----------------------------------------------------------------------- 1180 ln_diurnal = .false. ! 1181 ln_diurnal_only = .false. ! 1182 / 1183 !----------------------------------------------------------------------- 1184 &namflo ! float parameters ("key_float") 1185 !----------------------------------------------------------------------- 1186 jpnfl = 1 ! total number of floats during the run 1187 jpnnewflo = 0 ! number of floats for the restart 1188 ln_rstflo = .false. ! float restart (T) or not (F) 1189 nn_writefl = 75 ! frequency of writing in float output file 1190 nn_stockfl = 5475 ! frequency of creation of the float restart file 1191 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 1192 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 1193 ! ! or computed with Blanke' scheme (F) 1194 ln_ariane = .true. ! Input with Ariane tool convention(T) 1195 ln_flo_ascii= .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 1196 / 1197 !----------------------------------------------------------------------- 1198 &nam_diaharm ! Harmonic analysis of tidal constituents ("key_diaharm") 1151 1199 !----------------------------------------------------------------------- 1152 1200 nit000_han = 1 ! First time step used for harmonic analysis … … 1157 1205 / 1158 1206 !----------------------------------------------------------------------- 1159 &namdct ! transports through sections 1160 !----------------------------------------------------------------------- 1161 nn_dct = 15 ! time step frequency for transports computing 1162 nn_dctwri = 15 ! time step frequency for transports writing 1163 nn_secdebug = 112 ! 0 : no section to debug 1164 ! -1 : debug all section 1165 ! 0 < n : debug section number n 1166 / 1167 1168 !!====================================================================== 1169 !! *** Observation & Assimilation namelists *** 1170 !!====================================================================== 1171 !! namobs observation and model comparison ('key_diaobs') 1207 &namdct ! transports through some sections ("key_diadct") 1208 !----------------------------------------------------------------------- 1209 nn_dct = 15 ! time step frequency for transports computing 1210 nn_dctwri = 15 ! time step frequency for transports writing 1211 nn_secdebug= 112 ! 0 : no section to debug 1212 ! ! -1 : debug all section 1213 ! ! 0 < n : debug section number n 1214 / 1215 !----------------------------------------------------------------------- 1216 &nam_diatmb ! Top Middle Bottom Output (default F) 1217 !----------------------------------------------------------------------- 1218 ln_diatmb = .false. ! Choose Top Middle and Bottom output or not 1219 / 1220 !----------------------------------------------------------------------- 1221 &nam_dia25h ! 25h Mean Output (default F) 1222 !----------------------------------------------------------------------- 1223 ln_dia25h = .false. ! Choose 25h mean output or not 1224 / 1225 !----------------------------------------------------------------------- 1226 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 1227 !----------------------------------------------------------------------- 1228 nn_nchunks_i= 4 ! number of chunks in i-dimension 1229 nn_nchunks_j= 4 ! number of chunks in j-dimension 1230 nn_nchunks_k= 31 ! number of chunks in k-dimension 1231 ! ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 1232 ! ! is optimal for postprocessing which works exclusively with horizontal slabs 1233 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 1234 ! ! (F) ignore chunking information and produce netcdf3-compatible files 1235 / 1236 1237 !!====================================================================== 1238 !! *** Observation & Assimilation *** 1239 !!====================================================================== 1240 !! namobs observation and model comparison 1172 1241 !! nam_asminc assimilation increments ('key_asminc') 1173 1242 !!====================================================================== 1174 1243 ! 1175 1244 !----------------------------------------------------------------------- 1176 &namobs ! observation usage switch ('key_diaobs') 1177 !----------------------------------------------------------------------- 1178 ln_t3d = .false. ! Logical switch for T profile observations 1179 ln_s3d = .false. ! Logical switch for S profile observations 1180 ln_ena = .false. ! Logical switch for ENACT insitu data set 1181 ln_cor = .false. ! Logical switch for Coriolis insitu data set 1182 ln_profb = .false. ! Logical switch for feedback insitu data set 1183 ln_sla = .false. ! Logical switch for SLA observations 1184 ln_sladt = .false. ! Logical switch for AVISO SLA data 1185 ln_slafb = .false. ! Logical switch for feedback SLA data 1186 ln_ssh = .false. ! Logical switch for SSH observations 1187 ln_sst = .false. ! Logical switch for SST observations 1188 ln_reysst = .false. ! Logical switch for Reynolds observations 1189 ln_ghrsst = .false. ! Logical switch for GHRSST observations 1190 ln_sstfb = .false. ! Logical switch for feedback SST data 1191 ln_sss = .false. ! Logical switch for SSS observations 1192 ln_seaice = .false. ! Logical switch for Sea Ice observations 1193 ln_vel3d = .false. ! Logical switch for velocity observations 1194 ln_velavcur= .false ! Logical switch for velocity daily av. cur. 1195 ln_velhrcur= .false ! Logical switch for velocity high freq. cur. 1196 ln_velavadcp = .false. ! Logical switch for velocity daily av. ADCP 1197 ln_velhradcp = .false. ! Logical switch for velocity high freq. ADCP 1198 ln_velfb = .false. ! Logical switch for feedback velocity data 1199 ln_grid_global = .false. ! Global distribtion of observations 1200 ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table 1201 grid_search_file = 'grid_search' ! Grid search lookup file header 1245 &namobs ! observation usage switch 1246 !----------------------------------------------------------------------- 1247 ln_diaobs = .false. ! Logical switch for the observation operator 1248 ln_t3d = .false. ! Logical switch for T profile observations 1249 ln_s3d = .false. ! Logical switch for S profile observations 1250 ln_sla = .false. ! Logical switch for SLA observations 1251 ln_sst = .false. ! Logical switch for SST observations 1252 ln_sic = .false. ! Logical switch for Sea Ice observations 1253 ln_vel3d = .false. ! Logical switch for velocity observations 1254 ln_altbias = .false. ! Logical switch for altimeter bias correction 1255 ln_nea = .false. ! Logical switch for rejection of observations near land 1256 ln_grid_global = .true. ! Logical switch for global distribution of observations 1257 ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table 1258 ln_ignmis = .true. ! Logical switch for ignoring missing files 1259 ln_s_at_t = .false. ! Logical switch for computing model S at T obs if not there 1260 ln_sstnight = .false. ! Logical switch for calculating night-time average for SST obs 1202 1261 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files 1203 enactfiles = 'enact.nc' ! ENACT input observation file names (specify full array in namelist_cfg) 1204 coriofiles = 'corio.nc' ! Coriolis input observation file name 1205 profbfiles = 'profiles_01.nc' ! Profile feedback input observation file name 1206 ln_profb_enatim = .false ! Enact feedback input time setting switch 1207 slafilesact = 'sla_act.nc' ! Active SLA input observation file names 1208 slafilespas = 'sla_pass.nc' ! Passive SLA input observation file names 1209 slafbfiles = 'sla_01.nc' ! slafbfiles: Feedback SLA input observation file names 1210 sstfiles = 'ghrsst.nc' ! GHRSST input observation file names 1211 sstfbfiles = 'sst_01.nc' ! Feedback SST input observation file names 1212 seaicefiles = 'seaice_01.nc' ! Sea Ice input observation file names 1213 velavcurfiles = 'velavcurfile.nc' ! Vel. cur. daily av. input file name 1214 velhrcurfiles = 'velhrcurfile.nc' ! Vel. cur. high freq. input file name 1215 velavadcpfiles = 'velavadcpfile.nc' ! Vel. ADCP daily av. input file name 1216 velhradcpfiles = 'velhradcpfile.nc' ! Vel. ADCP high freq. input file name 1217 velfbfiles = 'velfbfile.nc' ! Vel. feedback input observation file name 1218 dobsini = 20000101.000000 ! Initial date in window YYYYMMDD.HHMMSS 1219 dobsend = 20010101.000000 ! Final date in window YYYYMMDD.HHMMSS 1220 n1dint = 0 ! Type of vertical interpolation method 1221 n2dint = 0 ! Type of horizontal interpolation method 1222 ln_nea = .false. ! Rejection of observations near land switch 1223 nmsshc = 0 ! MSSH correction scheme 1224 mdtcorr = 1.61 ! MDT correction 1225 mdtcutoff = 65.0 ! MDT cutoff for computed correction 1226 ln_altbias = .false. ! Logical switch for alt bias 1227 ln_ignmis = .true. ! Logical switch for ignoring missing files 1228 endailyavtypes = 820 ! ENACT daily average types - array (use namelist_cfg to set more values) 1229 ln_grid_global = .true. 1230 ln_grid_search_lookup = .false. 1231 / 1232 !----------------------------------------------------------------------- 1233 &nam_asminc ! assimilation increments ('key_asminc') 1234 !----------------------------------------------------------------------- 1235 ln_bkgwri = .false. ! Logical switch for writing out background state 1236 ln_trainc = .false. ! Logical switch for applying tracer increments 1237 ln_dyninc = .false. ! Logical switch for applying velocity increments 1238 ln_sshinc = .false. ! Logical switch for applying SSH increments 1239 ln_asmdin = .false. ! Logical switch for Direct Initialization (DI) 1240 ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) 1241 nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1] 1242 nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1] 1243 nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1] 1244 nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1] 1245 niaufn = 0 ! Type of IAU weighting function 1246 ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin 1247 salfixmin = -9999 ! Minimum salinity after applying the increments 1248 nn_divdmp = 0 ! Number of iterations of divergence damping operator 1249 / 1250 !----------------------------------------------------------------------- 1251 &namsbc_wave ! External fields from wave model 1252 !----------------------------------------------------------------------- 1253 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 1254 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 1255 sn_cdg = 'sdw_wave' , 1 , 'drag_coeff' , .true. , .false. , 'daily' , '' , '' , '' 1256 sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1257 sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false. , 'daily' , '' , '' , '' 1258 sn_swh = 'sdw_wave' , 1 , 'hs' , .true. , .false. , 'daily' , '' , '' , '' 1259 sn_wmp = 'sdw_wave' , 1 , 'wmp' , .true. , .false. , 'daily' , '' , '' , '' 1260 sn_wnum = 'sdw_wave' , 1 , 'wave_num' , .true. , .false. , 'daily' , '' , '' , '' 1261 sn_tauoc = 'sdw_wave' , 1 , 'wave_stress', .true. , .false. , 'daily' , '' , '' , '' 1262 ! 1263 cn_dir = './' ! root directory for the location of drag coefficient files 1264 / 1265 !----------------------------------------------------------------------- 1266 &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) 1267 !----------------------------------------------------------------------- 1268 ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 1269 ln_neptsimp = .false. ! yes/no use simplified neptune 1270 1271 ln_smooth_neptvel = .false. ! yes/no smooth zunep, zvnep 1272 rn_tslse = 1.2e4 ! value of lengthscale L at the equator 1273 rn_tslsp = 3.0e3 ! value of lengthscale L at the pole 1274 ! Specify whether to ramp down the Neptune velocity in shallow 1275 ! water, and if so the depth range controlling such ramping down 1276 ln_neptramp = .true. ! ramp down Neptune velocity in shallow water 1277 rn_htrmin = 100.0 ! min. depth of transition range 1278 rn_htrmax = 200.0 ! max. depth of transition range 1279 / 1262 cn_profbfiles = 'profiles_01.nc' ! Profile feedback input observation file names 1263 cn_slafbfiles = 'sla_01.nc' ! SLA feedback input observation file names 1264 cn_sstfbfiles = 'sst_01.nc' ! SST feedback input observation file names 1265 cn_sicfbfiles = 'sic_01.nc' ! SIC feedback input observation file names 1266 cn_velfbfiles = 'vel_01.nc' ! Velocity feedback input observation file names 1267 cn_altbiasfile = 'altbias.nc' ! Altimeter bias input file name 1268 cn_gridsearchfile='gridsearch.nc' ! Grid search file name 1269 rn_gridsearchres = 0.5 ! Grid search resolution 1270 rn_dobsini = 00010101.000000 ! Initial date in window YYYYMMDD.HHMMSS 1271 rn_dobsend = 00010102.000000 ! Final date in window YYYYMMDD.HHMMSS 1272 nn_1dint = 0 ! Type of vertical interpolation method 1273 nn_2dint = 0 ! Type of horizontal interpolation method 1274 nn_msshc = 0 ! MSSH correction scheme 1275 rn_mdtcorr = 1.61 ! MDT correction 1276 rn_mdtcutoff = 65.0 ! MDT cutoff for computed correction 1277 nn_profdavtypes = -1 ! Profile daily average types - array 1278 ln_sstbias = .false. ! 1279 cn_sstbias_files = 'sstbias.nc' ! 1280 / 1281 !----------------------------------------------------------------------- 1282 &nam_asminc ! assimilation increments ('key_asminc') 1283 !----------------------------------------------------------------------- 1284 ln_bkgwri = .false. ! Logical switch for writing out background state 1285 ln_trainc = .false. ! Logical switch for applying tracer increments 1286 ln_dyninc = .false. ! Logical switch for applying velocity increments 1287 ln_sshinc = .false. ! Logical switch for applying SSH increments 1288 ln_asmdin = .false. ! Logical switch for Direct Initialization (DI) 1289 ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) 1290 nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1] 1291 nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1] 1292 nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1] 1293 nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1] 1294 niaufn = 0 ! Type of IAU weighting function 1295 ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin 1296 salfixmin = -9999 ! Minimum salinity after applying the increments 1297 nn_divdmp = 0 ! Number of iterations of divergence damping operator 1298 /
Note: See TracChangeset
for help on using the changeset viewer.