Changeset 907 for trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist
- Timestamp:
- 2008-04-25T09:35:26+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist
r889 r907 1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 ! OPA namelist : model option and parameter input 3 ! ------------- 4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 5 ! 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : 1 - run manager (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 3 !! namelists 2 - Domain (nam_zgr, nam_zgr_sco, namdom) 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 5 !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 6 !! 4 - lateral boundary (namlbc, namcla, namobc) 7 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 8 !! 6 - Tracer (nameos, nam_traadv, nam_traldf, namtdp) 9 !! 7 - dynamics (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) 10 !! 8 - Verical physics (namzdf, namnpc, namric, namtke, namkpp, namddm) 11 !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) 12 !! 9 - miscellaneous (namsol) 13 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 14 ! CAUTION: some scripts does not support CAPITALs for logical definition 15 ! ******* use .true. or .false. and NOT .TRUE. or .FALSE. 16 17 !!====================================================================== 18 !! *** Run management namelists *** 19 !!====================================================================== 20 !! namrun parameters of the run 21 !! nam_mpp Massively Parallel Processing 22 !! nam_mpp_dyndist 23 !! nam_ctl Control prints & Benchmark 24 !!====================================================================== 25 6 26 !----------------------------------------------------------------------- 7 27 ! namrun parameters of the run … … 10 30 ! cexper experience name for vairmer format 11 31 ! ln_rstart boolean term for restart (true or false) 12 ! nrstdt control of the restart timestep: 13 ! = 0 restart, do not control nit000 in the restart file. 14 ! = 1 restart, control nit000 in the restart file. Do not 15 ! use the date in the restart file (use ndate0 in namelist) 16 ! = 2 restart, control nit000 in the restart file, use the date 17 ! in the restart file. ndate0 in the namelist is ignored. 32 ! nrstdt restart control = 0 restart, do not control nit000 in the restart file. 33 ! ! = 1 restart, control nit000 in the restart file. Do not 34 ! ! use the date in the restart file (use ndate0 in namelist) 35 ! ! = 2 restart, control nit000 in the restart file, use the date 36 ! ! in the restart file. ndate0 in the namelist is ignored. 18 37 ! nit000 number of the first time step 19 38 ! nitend number of the last time step … … 24 43 ! nwrite frequency of OUTPUT file 25 44 ! ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 26 !27 ! CAUTION: for usual run scripts, logical value of28 ! ******* ln_rstart must be .true. or .false.29 ! and NOT .TRUE. or .FALSE.30 45 &namrun 31 46 no = 0 … … 38 53 nleapy = 0 39 54 ninist = 0 40 nstock = 10950 41 nwrite = 10950 42 ln_dimgnnn = .FALSE. 55 nstock = 10950 56 nwrite = 10950 57 ln_dimgnnn = .false. 58 / 59 !----------------------------------------------------------------------- 60 ! nam_mpp Massively Parallel Processing 61 !----------------------------------------------------------------------- 62 ! c_mpi_send mpi send/recieve type 63 ! = 'S' : standard blocking send 64 ! = 'B' : buffer blocking send 65 ! = 'I' : immediate non-blocking send 66 &nam_mpp 67 c_mpi_send = 'S' 68 / 69 !----------------------------------------------------------------------- 70 ! nam_mpp_dyndist Massively Parallel Distribution ("key_agrif" && "key_mpp_dyndist") 71 !----------------------------------------------------------------------- 72 ! jpni number of processors following i 73 ! jpnj number of processors following j 74 ! jpnij number of local domains 75 &nam_mpp_dyndist 76 jpni = 1 77 jpnj = 1 78 jpnij = 1 43 79 / 44 80 !----------------------------------------------------------------------- … … 69 105 nbit_cmp = 0 70 106 / 71 !----------------------------------------------------------------------- 72 ! nam_mpp Massively Parallel Processing 73 !----------------------------------------------------------------------- 74 ! c_mpi_send mpi send/recieve type 75 ! = 'S' : standard blocking send 76 ! = 'B' : buffer blocking send 77 ! = 'I' : immediate non-blocking send 78 &nam_mpp 79 c_mpi_send = 'S' 80 / 81 !----------------------------------------------------------------------- 82 ! nam_mpp_dyndist Massively Parallel Distribution 83 ! #ifdef ( key_agrif && key_mpp_dyndist ) 84 !----------------------------------------------------------------------- 85 ! jpni number of processors following i 86 ! jpnj number of processors following j 87 ! jpnij number of local domains 88 &nam_mpp_dyndist 89 jpni = 1 90 jpnj = 1 91 jpnij = 1 92 / 107 108 !!====================================================================== 109 !! *** Domain namelists *** 110 !!====================================================================== 111 !! nam_zgr vertical coordinate 112 !! nam_zgr_sco s-coordinate or hybrid z-s-coordinate 113 !! namdom space and time domain (bathymetry, mesh, timestep) 114 !!====================================================================== 115 93 116 !----------------------------------------------------------------------- 94 117 ! nam_zgr vertical coordinate 95 118 !----------------------------------------------------------------------- 96 ! ln_zco z-coordinate - full steps (T/F) 119 ! ln_zco z-coordinate - full steps (T/F) ("key_zco" may also be defined) 97 120 ! ln_zps z-coordinate - partial steps (T/F) 98 121 ! ln_sco s- or hybrid z-s-coordinate (T/F) … … 118 141 / 119 142 !----------------------------------------------------------------------- 120 ! nam_traadv advection scheme for tracer (option not control by CPP keys)121 !-----------------------------------------------------------------------122 ! ln_traadv_cen2 2nd order centered scheme (default T)123 ! ln_traadv_tvd TVD scheme (default F)124 ! ln_traadv_muscl MUSCL scheme (default F)125 ! ln_traadv_muscl2 MUSCL2 scheme (default F)126 ! ln_traadv_ubs UBS scheme (default F)127 ! ln_traadv_qck QUICKEST scheme (default F)128 &nam_traadv129 ln_traadv_cen2 = .false.130 ln_traadv_tvd = .true.131 ln_traadv_muscl = .false.132 ln_traadv_muscl2 = .false.133 ln_traadv_ubs = .false.134 ln_traadv_qck = .false.135 /136 !-----------------------------------------------------------------------137 ! nam_traldf lateral diffusion scheme for tracer (option not control by CPP keys)138 !-----------------------------------------------------------------------139 ! Type of the operator :140 ! ln_traldf_lap laplacian operator (default T)141 ! ln_traldf_bilap bilaplacian operator (default F)142 ! Direction of action :143 ! ln_traldf_level iso-level (default F)144 ! ln_traldf_hor horizontal (geopotential) (default F)^**145 ! ln_traldf_iso iso-neutral (default T)^*146 ! Coefficient147 ! aht0 horizontal eddy diffusivity for tracers (m2/s)148 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s)149 ! aeiv0 eddy induced velocity coefficient (m2/s)150 ! ^* require key_ldfslp to compute the direction of the lateral diffusion151 ! ^** require key_ldfslp in s-coordinate152 &nam_traldf153 ln_traldf_lap = .true.154 ln_traldf_bilap = .false.155 ln_traldf_level = .false.156 ln_traldf_hor = .false.157 ln_traldf_iso = .true.158 aht0 = 1000.159 ahtb0 = 0.160 aeiv0 = 0.161 /162 !-----------------------------------------------------------------------163 ! nam_dynldf lateral diffusion on momentum164 !-----------------------------------------------------------------------165 ! Type of the operator :166 ! ln_dynldf_lap laplacian operator (default T)167 ! ln_dynldf_bilap bilaplacian operator (default F)168 ! Direction of action :169 ! ln_dynldf_level iso-level (default F)170 ! ln_dynldf_hor horizontal (geopotential) (default F)^**171 ! ln_dynldf_iso iso-neutral (default T)^*172 ! Coefficient173 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s)174 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s)175 &nam_dynldf176 ln_dynldf_lap = .false.177 ln_dynldf_bilap = .true.178 ln_dynldf_level = .false.179 ln_dynldf_hor = .true.180 ln_dynldf_iso = .false.181 ahm0 = -8.5e+11182 ahmb0 = 0.183 /184 !-----------------------------------------------------------------------185 ! namflg algorithm flags (algorithm not control by CPP keys)186 !-----------------------------------------------------------------------187 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T)188 ! centered time scheme (F)189 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0)190 &namflg191 ln_dynhpg_imp = .false.192 nn_dynhpg_rst = 0193 /194 !-----------------------------------------------------------------------195 ! nam_dynhpg Hydrostatic pressure gradient option196 !-----------------------------------------------------------------------197 ! type of pressure gradient scheme (choose one only!)198 ! ln_hpg_zco z-coordinate - full steps (default T)199 ! ln_hpg_zps z-coordinate - partial steps (interpolation)200 ! ln_hpg_sco s-coordinate (standard jacobian formulation)201 ! ln_hpg_hel s-coordinate (helsinki modification)202 ! ln_hpg_wdj s-coordinate (weighted density jacobian)203 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial)204 ! ln_hpg_rot s-coordinate (ROTated axes scheme)205 ! parameters206 ! gamm weighting coefficient (wdj scheme)207 &nam_dynhpg208 ln_hpg_zco = .false.209 ln_hpg_zps = .true.210 ln_hpg_sco = .false.211 ln_hpg_hel = .false.212 ln_hpg_wdj = .false.213 ln_hpg_djc = .false.214 ln_hpg_rot = .false.215 gamm = 0.e0216 /217 !-----------------------------------------------------------------------218 ! nam_dynadv option of physics/algorithm (not control by CPP keys)219 !-----------------------------------------------------------------------220 ! ln_dynadv_vec vector form flag221 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T)222 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F)223 &nam_dynadv224 ln_dynadv_vec = .TRUE.225 ln_dynadv_cen2 = .FALSE.226 ln_dynadv_ubs = .FALSE.227 /228 !-----------------------------------------------------------------------229 ! nam_dynvor option of physics/algorithm (not control by CPP keys)230 !-----------------------------------------------------------------------231 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T)232 ! ln_dynvor_ene " " : energy conserving scheme (default F)233 ! ln_dynvor_mix " " : mixed scheme (default F)234 ! ln_dynvor_een " " : energy & enstrophy scheme (default F)235 &nam_dynvor236 ln_dynvor_ene = .FALSE.237 ln_dynvor_ens = .FALSE.238 ln_dynvor_mix = .FALSE.239 ln_dynvor_een = .TRUE.240 /241 !-----------------------------------------------------------------------242 ! namalb albedo parameters243 !-----------------------------------------------------------------------244 ! cgren correction of the snow or ice albedo to take into account245 ! albice albedo of melting ice in the arctic and antarctic246 ! alphd coefficients for linear interpolation used to compute albedo247 ! between two extremes values (Pyane, 1972)248 ! alphc " "249 ! alphdi " "250 &namalb251 cgren = 0.06252 albice = 0.5253 alphd = 0.80254 alphc = 0.65255 alphdi = 0.72256 /257 !-----------------------------------------------------------------------258 143 ! namdom space and time domain (bathymetry, mesh, timestep) 259 144 !----------------------------------------------------------------------- 260 ! ntopo = 0/1 ,compute/read the bathymetry file 261 ! (mbathy, nb of T-ocean levels) 145 ! ntopo = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 262 146 ! e3zps_min the thickness of the partial step is set larger than the 263 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t 264 ! (N.B. 0<e3zps_rat<1) 265 ! (coordinates, scale factors) 147 ! e3zps_rat the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 266 148 ! nmsh =1 create a mesh file (coordinates, scale factors, masks) 267 ! nacc the acceleration of convergence method 268 ! = 0, no acceleration, rdt = rdttra 269 ! = 1, acceleration used, rdt < rdttra(k) 149 ! nacc acceleration of convergence method = 0, not used, rdt = rdttra 150 ! = 1, used, rdt < rdttra(k) 270 151 ! atfp asselin time filter parameter 271 152 ! rdt time step for the dynamics (and tracer if nacc=0) … … 274 155 ! rdth depth variation of tracer time step 275 156 ! rdtbt barotropic time step (for the time splitting algorithm) 276 ! nclosea 0 no closed sea277 ! 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)157 ! nclosea =0 remove the closed sea from the global domain (orca configuration) 158 ! =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 278 159 &namdom 279 160 ntopo = 1 … … 290 171 nclosea = 0 291 172 / 292 !----------------------------------------------------------------------- 293 ! namfwb freshwater budget correction 294 !----------------------------------------------------------------------- 295 ! ln_fwb logical flag for freshwater budget correction (0 annual mean) 296 &namfwb 297 ln_fwb = .false. 298 / 299 !----------------------------------------------------------------------- 300 ! namptr Poleward Transport Diagnostic 301 !----------------------------------------------------------------------- 302 ! ln_diaptr logical flag for Poleward transport computation 303 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 304 ! need input basins mask file named "subbasins.nc" 305 ! nf_ptr Frequency of computation 306 &namptr 307 ln_diaptr = .false. 308 ln_subbas = .false. 309 nf_ptr = 15 310 / 311 !----------------------------------------------------------------------- 312 ! namcro cross land advection 313 !----------------------------------------------------------------------- 314 ! n_cla advection between 2 ocean pts separates by land 315 &namcla 316 n_cla = 0 317 / 318 !----------------------------------------------------------------------- 319 ! namzdf vertical physics 320 !----------------------------------------------------------------------- 321 ! ln_zdfevd enhanced vertical diffusion (default T) 322 ! ln_zdfnpc Non-Penetrative Convection (default T) 323 ! avm0 vertical eddy viscosity for the dynamic (m2/s) 324 ! avt0 vertical eddy diffusivity for tracers (m2/s) 325 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 326 ! nevdm = 0 apply enhanced mixing on tracer only 327 ! = 1 apply enhanced mixing on both tracer and momentum 328 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 329 ! (=F) euler backward (F) 330 ! n_zdfexp number of sub-timestep for time splitting scheme 331 &namzdf 332 ln_zdfevd = .true. 333 ln_zdfnpc = .false. 334 avm0 = 1.2e-4 335 avt0 = 1.2e-5 336 avevd = 100. 337 n_evdm = 1 338 ln_zdfexp = .false. 339 n_zdfexp = 3 340 / 341 !----------------------------------------------------------------------- 342 ! namnpc vnon penetrative convection 343 !----------------------------------------------------------------------- 344 ! nnpc1 non penetrative convective scheme frequency 345 ! nnpc2 non penetrative convective scheme print frequency 346 &namnpc 347 nnpc1 = 1 348 nnpc2 = 365 349 / 350 !----------------------------------------------------------------------- 351 ! nambbl bottom boundary layer scheme 352 !----------------------------------------------------------------------- 353 ! atrbbl lateral tracer coeff. for bottom boundary layer scheme(m2/s) 354 &nambbl 355 atrbbl = 10000. 356 / 357 !----------------------------------------------------------------------- 358 ! namric richardson number dependent vertical diffusion 359 ! ( #ifdef "key_zdfrichardson" ) 360 !----------------------------------------------------------------------- 361 ! avmri maximum value of the vertical viscosity 362 ! alp coefficient of the parameterization 363 ! nric coefficient of the parameterization 364 &namric 365 avmri = 100.e-4 366 alp = 5. 367 nric = 2 368 / 369 !----------------------------------------------------------------------- 370 ! namtke turbulent eddy kinetic dependent vertical diffusion 371 ! ( #ifdef "key_zdftke" ) 372 !----------------------------------------------------------------------- 373 ! ln_rstke flag to restart with tke from a run without tke (default F) 374 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 375 ! ediss coef. of the Kolmogoroff dissipation 376 ! ebb coef. of the surface input of tke 377 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 378 ! emin minimum value of tke (m^2/s^2) 379 ! emin0 surface minimum value of tke (m^2/s^2) 380 ! nitke number of restart iterative loops 381 ! ri_c critic richardson number 382 ! nmxl flag on mixing length used 383 ! = 0 bounded by the distance to surface and bottom 384 ! = 1 bounded by the local vertical scale factor 385 ! = 2 first vertical derivative of mixing length bounded by 1 386 ! npdl flag on prandtl number 387 ! = 0 no vertical prandtl number (avt=avm) 388 ! = 1 prandtl number function of richarson number (avt=pdl*avm) 389 ! = 2 same as = 1 but a shapiro filter is applied on pdl 390 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 391 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 392 &namtke 393 ln_rstke = .false. 394 ediff = 0.1 395 ediss = 0.7 396 ebb = 60. 397 efave = 1. 398 emin = 1.e-6 399 emin0 = 1.e-4 400 nitke = 50 401 nmxl = 2 402 npdl = 1 403 navb = 0 404 / 405 !----------------------------------------------------------------------- 406 ! namkpp K-Profile Parameterization dependent vertical diffusion 407 ! ( #ifdef "key_zdfkpp" ) 408 !----------------------------------------------------------------------- 409 ! ln_kpprimix shear instability mixing (default T) 410 ! difmiw constant internal wave viscosity (m2/s) 411 ! difsiw constant internal wave diffusivity (m2/s) 412 ! Riinfty local Richardson Number limit for shear instability 413 ! difri maximum shear mixing at Rig = 0 (m2/s) 414 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 415 ! difcon maximum mixing in interior convection (m2/s) 416 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 417 ! navb = 0/1 flag for constant or profile background avt 418 &namkpp 419 ln_kpprimix = .true. 420 difmiw = 1.e-04 421 difsiw = 0.1e-04 422 Riinfty = 0.8 423 difri = 0.0050 424 bvsqcon = -0.01e-07 425 difcon = 1. 426 navb = 0 427 nave = 1 428 / 429 !----------------------------------------------------------------------- 430 ! namddm double diffusive mixing parameterization 431 !----------------------------------------------------------------------- 432 ! avts maximum avs for dd mixing 433 ! hsbfr heat/salt buoyancy flux ratio 434 &namddm 435 avts = 1.e-4 436 hsbfr = 1.6 437 / 438 !----------------------------------------------------------------------- 439 ! namlbc lateral momentum boundary condition 440 !----------------------------------------------------------------------- 441 ! shlat lateral boundary condition on velocity 442 ! shlat = 0 , free slip 443 ! 0 < shlat < 2 , partial slip 444 ! shlat = 2 , no slip 445 ! 2 < shlat , strong slip 446 &namlbc 447 shlat = 2. 448 / 449 !----------------------------------------------------------------------- 450 ! nambfr bottom friction 451 !----------------------------------------------------------------------- 452 ! nbotfr type of bottom friction 453 ! nbotfr = 0 , no slip 454 ! nbotfr = 1 , linear friction 455 ! nbotfr = 2 , nonlinear friction 456 ! nbotfr = 3 , free slip 457 ! bfri1 bottom drag coefficient (linear case) 458 ! bfri2 bottom drag coefficient (non linear case) 459 ! bfeb2 bottom turbulent kinetic energy (m^2/s^2) 460 &nambfr 461 nbotfr = 1 462 bfri1 = 4.e-4 463 bfri2 = 1.e-3 464 bfeb2 = 2.5e-3 465 / 466 !----------------------------------------------------------------------- 467 ! nambbc bottom temperature boundary condition 468 !----------------------------------------------------------------------- 469 ! ngeo_flux = 0 no geothermal heat flux 470 ! = 1 constant geothermal heat flux 471 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc) 472 ! ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 473 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 474 &nambbc 475 ngeo_flux = 2 476 ngeo_flux_const = 86.4e-3 173 174 !!====================================================================== 175 !! *** Surface Boundary Condition namelists *** 176 !!====================================================================== 177 !! namsbc surface boundary condition 178 !! namsbc_ana analytical formulation 179 !! namsbc_flx flux formulation 180 !! namsbc_clio CLIO bulk formulea formulation 181 !! namsbc_core CORE bulk formulea formulation 182 !! namsbc_cpl CouPLed formulation ("key_coupled") 183 !! namqsr penetrative solar radiation 184 !! namsbc_rnf river runoffs 185 !! namsbc_ssr sea surface restoring term (for T and/or S) 186 !! namalb albedo parameters 187 !!====================================================================== 188 189 !----------------------------------------------------------------------- 190 ! namsbc Surface Boundary Condition (surface module) 191 !----------------------------------------------------------------------- 192 ! nn_fsbc frequency of surface boundary condition computation 193 ! (= the frequency of sea-ice model call) 194 ! ln_ana analytical formulation (fill namsbc_ana ) 195 ! ln_flx flux formulation (fill namsbc_flx ) 196 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) 197 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) 198 ! ln_cpl Coupled formulation (fill namsbc_cpl ) 199 ! nn_ice =0 no ice boundary condition , 200 ! =1 observed ice-cover , 201 ! =2 ice-model used ("key_lim3" or "key_lim2) 202 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 203 ! ln_ssr Sea Surface Restoring on T and/or S (fill namsbc_ssr) 204 ! nn_fwb FreshWater Budget: =0 no control , 205 ! =1 annual global mean of e-p-r set to zero , 206 ! =2 global mean of e-p-r set to zero at each nn_fsbc time step 207 &namsbc 208 nn_fsbc = 5 209 ln_ana = .false. 210 ln_flx = .false. 211 ln_blk_clio = .true. 212 ln_blk_core = .false. 213 ln_cpl = .false. 214 nn_ice = 2 215 ln_dm2dc = .false. 216 ln_rnf = .false. 217 ln_ssr = .false. 218 nn_fwb = 0 219 / 220 !----------------------------------------------------------------------- 221 ! namsbc_ana analytical surface boundary condition 222 !----------------------------------------------------------------------- 223 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 224 ! rn_utau0 uniform value for the i-stress 225 ! rn_vtau0 uniform value for the j-stress 226 ! rn_q0 uniform value for the total heat flux 227 ! rn_qsr0 uniform value for the solar radiation 228 ! rn_emp0 uniform value for the freswater budget (E-P) 229 &namsbc_ana 230 nn_tau000 = 0 231 rn_utau0 = 0.5 232 rn_vtau0 = 0.e0 233 rn_q0 = 0.e0 234 rn_qsr0 = 0.e0 235 rn_emp0 = 0.e0 236 / 237 !----------------------------------------------------------------------- 238 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 239 !----------------------------------------------------------------------- 240 ! cn_dir directory in which the model is executed 241 ! 242 ! THE ORDER OF THE FILES MATTER: 243 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 244 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 245 &namsbc_flx 246 cn_dir = './' 247 sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 248 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 249 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 250 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 251 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 252 / 253 !----------------------------------------------------------------------- 254 ! namsbc_clio CLIO bulk formulea 255 !----------------------------------------------------------------------- 256 ! cn_dir directory in which the model is executed 257 ! 258 ! THE ORDER OF THE FILES MATTER: 259 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 260 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 261 &namsbc_clio 262 cn_dir = './' 263 sn_utau = 'taux_1m' , -12. , 'sozotaux', .false. , 1 , 0 264 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .false. , 1 , 0 265 sn_wndm = 'flx' , -12. , 'socliowi', .false. , 1 , 0 266 sn_tair = 'flx' , -12. , 'socliot1', .false. , 1 , 0 267 sn_humi = 'flx' , -12. , 'socliohu', .false. , 1 , 0 268 sn_ccov = 'flx' , -12. , 'socliocl', .false. , 1 , 0 269 sn_prec = 'flx' , -12. , 'socliopl', .false. , 1 , 0 270 / 271 !----------------------------------------------------------------------- 272 ! namsbc_core CORE bulk formulea 273 !----------------------------------------------------------------------- 274 ! cn_dir directory in which the model is executed 275 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 276 ! alpha_precip multiplication factor for precipitation (total & snow) 277 ! 278 ! THE ORDER OF THE FILES MATTER: 279 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 280 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 281 &namsbc_core 282 cn_dir = './' 283 ln_2m = .false. 284 alpha_precip = 1. 285 sn_wndi = 'u10' , 24. , 'U_10_MOD', .false. , 1 , 0 286 sn_wndj = 'v10' , 24. , 'V_10_MOD', .false. , 1 , 0 287 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .false. , 1 , 0 288 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .false. , 1 , 0 289 sn_tair = 't10' , 24. , 'T_10_MOD', .false. , 1 , 0 290 sn_humi = 'q10' , 24. , 'Q_10_MOD', .false. , 1 , 0 291 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 292 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 293 / 294 !----------------------------------------------------------------------- 295 ! namsbc_cpl coupled ocean/atmosphere model ("key_coupled") 296 !----------------------------------------------------------------------- 297 &namsbc_cpl 477 298 / 478 299 !----------------------------------------------------------------------- … … 490 311 / 491 312 !----------------------------------------------------------------------- 492 ! namtdp tracer newtonian damping ('key_tradmp') 493 !----------------------------------------------------------------------- 494 ! ndmp type of damping in temperature and salinity 495 ! (='latitude', damping poleward of 'ndmp' degrees and function 496 ! of the distance-to-coast. Red and Med Seas as ndmp=-1) 497 ! (=-1 damping only in Med and Red Seas) 498 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 499 ! nmldmp type of damping in the mixed layer 500 ! (=0 damping throughout the water column) 501 ! (=1 no damping in the mixed layer defined by avt >5cm2/s ) 502 ! (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 503 ! sdmp surface time scale for internal damping (days) 504 ! bdmp bottom time scale for internal damping (days) 505 ! hdmp depth of transition between sdmp and bdmp (meters) 506 &namtdp 507 ndmp = 90 508 ndmpf = 1 509 nmldmp = 1 510 sdmp = 50. 511 bdmp = 360. 512 hdmp = 800. 513 / 514 !----------------------------------------------------------------------- 515 ! nameos ocean physical parameters 516 !----------------------------------------------------------------------- 517 ! neos type of equation of state and Brunt-Vaisala frequency 518 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) 519 ! and of McDougall (1987) ) 520 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 521 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 522 ! with rau0=1020 set in parcst routine 523 ! ralpha thermal expension coefficient (linear equation of state) 524 ! rbeta saline expension coefficient (linear equation of state) 525 &nameos 526 neos = 0 527 ralpha = 2.e-4 528 rbeta = 0.001 529 / 530 !----------------------------------------------------------------------- 531 ! namsol elliptic solver / island / free surface 532 !----------------------------------------------------------------------- 533 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 534 ! (=2 successive-over-relaxation: sor) 535 ! (=3 FETI: fet, all require "key_feti" defined) 536 ! nsol_arp absolute/relative (0/1) precision convergence test 537 ! nmin minimum of iterations for the SOR solver 538 ! nmax maximum of iterations for the SOR solver 539 ! nmod frequency of test for the SOR solver 540 ! eps absolute precision of the solver 541 ! resmax absolute precision for the SOR solver 542 ! sor optimal coefficient for SOR solver 543 ! epsisl absolute precision on stream function solver 544 ! nmisl maximum pcg iterations for island 545 ! rnu strength of the additional force used in free surface b.c. 546 &namsol 547 nsolv = 1 548 nsol_arp = 0 549 nmin = 300 550 nmax = 800 551 nmod = 10 552 eps = 1.E-6 553 resmax = 1.E-10 554 sor = 1.92 555 epsisl = 1.e-10 556 nmisl = 4000 557 rnu = 1. 558 / 559 !======================================================================= 560 ! Diagnostics namelists 561 ! namtrd dynamics and/or tracer trends 562 ! namgap level mean model-data gap 563 ! namznl zonal mean heat & freshwater fluxes computation 564 ! namspr surface pressure in rigid-lid 565 !======================================================================= 566 !----------------------------------------------------------------------- 567 ! namtrd diagnostics on dynamics and/or tracer trends 568 ! ('key_trdyn' and/or 'key_trdtra') 569 ! or mixed-layer trends ('key_trdmld') 570 !----------------------------------------------------------------------- 571 ! ntrd time step frequency dynamics and tracers trends 572 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 573 ! ln_trdmld_restart restart for ML diagnostics 574 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 575 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 576 &namtrd 577 ntrd = 365 578 nctls = 0 579 ln_trdmld_restart = .false. 580 ucf = 1. 581 ln_trdmld_instant = .false. 582 / 583 !----------------------------------------------------------------------- 584 ! namgap level mean model-data gap ('key_diagap') 585 !----------------------------------------------------------------------- 586 ! ngap time-step frequency of model-data gap computation 587 ! nprg time-step frequency of gap print in model output 588 &namgap 589 ngap = 15 590 nprg = 10 591 / 592 !----------------------------------------------------------------------- 593 ! namznl zonal mean heat & freshwater fluxes computation 594 ! (#ifdef "key_diaznl") 595 !----------------------------------------------------------------------- 596 ! nfznl time-step frequency of zonal mean fluxes computation 597 &namznl 598 nfznl = 15 599 / 600 !----------------------------------------------------------------------- 601 ! namspr surface pressure diagnostic 602 !----------------------------------------------------------------------- 603 ! nmaxp maximum of iterations for the solver 604 ! epsp absolute precision of the solver 605 ! niterp number of iteration done by the solver 606 &namspr 607 nmaxp = 1000 608 epsp = 1.e-3 609 niterp = 400 313 ! namsbc_rnf runoffs namelist surface boundary condition 314 !----------------------------------------------------------------------- 315 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 316 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 317 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 318 !cn_dir Root directory for location of ssr files 319 !sn_rnf informations about the runoff file to be read 320 !sn_cnf informations about the runoff mouth file to be read 321 ! THE ORDER OF THE FILES MATTER: 322 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 323 ! ! ! (if <0 months) ! name ! (logical) ! 324 &namsbc_rnf 325 nn_runoff = 2 326 rn_hrnf = 0.e0 327 rn_avt_rnf = 1.e-3 328 cn_dir = './' 329 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .true. , 1 , 0 330 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .false. , 1 , 0 331 / 332 !----------------------------------------------------------------------- 333 ! namsbc_ssr surface boundary condition : sea surface restoring 334 !----------------------------------------------------------------------- 335 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 336 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 337 ! (=1, salt flux, concentration/dillution emps only) 338 ! (=2, volume flux, both emp and emps are updated) 339 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 340 ! deds magnitude of the damping on salinity [mm/day] 341 ! THE ORDER OF THE FILES MATTER: 342 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 343 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 344 &namsbc_ssr 345 nn_sstr = 0 346 nn_sssr = 0 347 dqdt = -40. 348 deds = 27.7 349 sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 350 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 351 / 352 !----------------------------------------------------------------------- 353 ! namalb albedo parameters 354 !----------------------------------------------------------------------- 355 ! cgren correction of the snow or ice albedo to take into account 356 ! albice albedo of melting ice in the arctic and antarctic 357 ! alphd coefficients for linear interpolation used to compute albedo 358 ! between two extremes values (Pyane, 1972) 359 ! alphc " " 360 ! alphdi " " 361 &namalb 362 cgren = 0.06 363 albice = 0.5 364 alphd = 0.80 365 alphc = 0.65 366 alphdi = 0.72 367 / 368 369 !!====================================================================== 370 !! *** Lateral boundary condition *** 371 !!====================================================================== 372 !! namlbc lateral momentum boundary condition 373 !! namcla cross land advection 374 !! namobc open boundaries parameters ("key_obc") 375 !!====================================================================== 376 377 !----------------------------------------------------------------------- 378 ! namlbc lateral momentum boundary condition 379 !----------------------------------------------------------------------- 380 ! shlat shlat = 0 , free slip ; 0 < shlat < 2 , partial slip 381 ! shlat = 2 , no slip ; 2 < shlat , strong slip 382 &namlbc 383 shlat = 2. 384 / 385 !----------------------------------------------------------------------- 386 ! namcla cross land advection 387 !----------------------------------------------------------------------- 388 ! n_cla advection between 2 ocean pts separates by land 389 &namcla 390 n_cla = 0 610 391 / 611 392 !----------------------------------------------------------------------- … … 639 420 ln_vol_cst = .false. 640 421 / 641 !----------------------------------------------------------------------- 642 ! namflo float parameters (#ifdef key_float) 422 423 !!====================================================================== 424 !! *** Bottom boundary condition *** 425 !!====================================================================== 426 !! nambfr bottom friction 427 !! nambbc bottom temperature boundary condition ("key_trabbc") 428 !! nambbl bottom boundary layer scheme ("key_trabbl_dif","key_trabbl_adv") 429 !!====================================================================== 430 !----------------------------------------------------------------------- 431 ! nambfr bottom friction 432 !----------------------------------------------------------------------- 433 ! nbotfr type of bottom friction : = 0 : no slip ; = 1 : linear friction 434 ! = 3 : free slip ; = 2 : non linear friction 435 ! bfri1 bottom drag coefficient (linear case) 436 ! bfri2 bottom drag coefficient (non linear case) 437 ! bfeb2 bottom turbulent kinetic energy (m2/s2) 438 &nambfr 439 nbotfr = 1 440 bfri1 = 4.e-4 441 bfri2 = 1.e-3 442 bfeb2 = 2.5e-3 443 / 444 !----------------------------------------------------------------------- 445 ! nambbc bottom temperature boundary condition 446 !----------------------------------------------------------------------- 447 ! ngeo_flux = 0 no geothermal heat flux 448 ! = 1 constant geothermal heat flux 449 ! = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 450 ! ngeo_flux_const Constant value of geothermal heat flux (W/m2) 451 &nambbc 452 ngeo_flux = 2 453 ngeo_flux_const = 86.4e-3 454 / 455 !----------------------------------------------------------------------- 456 ! nambbl bottom boundary layer scheme 457 !----------------------------------------------------------------------- 458 ! ! diffusive bbl ("key_trabbl") 459 ! ! advective bbl ("key_trabbl_adv") 460 ! atrbbl mixing coefficient of the bottom boundary layer scheme (m2/s) 461 &nambbl 462 atrbbl = 10000. 463 / 464 465 !!====================================================================== 466 !! Tracer (T & S ) namelists 467 !!====================================================================== 468 !! nameos equation of state 469 !! nam_traadv advection scheme 470 !! nam_traldf lateral diffusion scheme 471 !! namtdp tracer newtonian damping ("key_tradmp") 472 !!====================================================================== 473 474 !----------------------------------------------------------------------- 475 ! nameos ocean physical parameters 476 !----------------------------------------------------------------------- 477 ! neos type of equation of state and Brunt-Vaisala frequency 478 ! ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 479 ! ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 480 ! ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 481 ! ralpha thermal expension coefficient (neos= 1 or 2) 482 ! rbeta saline expension coefficient (neos= 2) 483 &nameos 484 neos = 0 485 ralpha = 2.e-4 486 rbeta = 0.001 487 / 488 !----------------------------------------------------------------------- 489 ! nam_traadv advection scheme for tracer 490 !----------------------------------------------------------------------- 491 ! ln_traadv_cen2 2nd order centered scheme (default T) 492 ! ln_traadv_tvd TVD scheme (default F) 493 ! ln_traadv_muscl MUSCL scheme (default F) 494 ! ln_traadv_muscl2 MUSCL scheme + cen2 at boundaries (default F) 495 ! ln_traadv_ubs UBS scheme (default F) 496 ! ln_traadv_qck QUICKEST scheme (default F) 497 &nam_traadv 498 ln_traadv_cen2 = .false. 499 ln_traadv_tvd = .true. 500 ln_traadv_muscl = .false. 501 ln_traadv_muscl2 = .false. 502 ln_traadv_ubs = .false. 503 ln_traadv_qck = .false. 504 / 505 !----------------------------------------------------------------------- 506 ! nam_traldf lateral diffusion scheme for tracer 507 !----------------------------------------------------------------------- 508 ! Type of the operator : 509 ! ln_traldf_lap laplacian operator (default T) 510 ! ln_traldf_bilap bilaplacian operator (default F) 511 ! Direction of action : 512 ! ln_traldf_level iso-level (default F) 513 ! ln_traldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" when ln_sco=T) 514 ! ln_traldf_iso iso-neutral (default T) (require "key_ldfslp") 515 ! Coefficient 516 ! aht0 horizontal eddy diffusivity for tracers (m2/s) 517 ! ahtb0 background eddy diffusivity for isopycnal diffusion (m2/s) 518 ! aeiv0 eddy induced velocity coefficient (m2/s) (require "key_traldf_eiv") 519 &nam_traldf 520 ln_traldf_lap = .true. 521 ln_traldf_bilap = .false. 522 ln_traldf_level = .false. 523 ln_traldf_hor = .false. 524 ln_traldf_iso = .true. 525 aht0 = 1000. 526 ahtb0 = 0. 527 aeiv0 = 0. 528 / 529 !----------------------------------------------------------------------- 530 ! namtdp tracer newtonian damping ('key_tradmp') 531 !----------------------------------------------------------------------- 532 ! ndmp type of damping in temperature and salinity 533 ! ! ='latitude', damping poleward of 'ndmp' degrees and function 534 ! ! of the distance-to-coast. Red and Med Seas as ndmp=-1 535 ! ! =-1 damping only in Med and Red Seas 536 ! ndmpf =1 create a damping.coeff NetCDF file (the 3D damping array) 537 ! nmldmp type of damping in the mixed layer 538 ! ! =0 damping throughout the water column 539 ! ! =1 no damping in the mixed layer defined by avt >5cm2/s 540 ! ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 541 ! sdmp surface time scale for internal damping (days) 542 ! bdmp bottom time scale for internal damping (days) 543 ! hdmp depth of transition between sdmp and bdmp (meters) 544 &namtdp 545 ndmp = 90 546 ndmpf = 1 547 nmldmp = 1 548 sdmp = 50. 549 bdmp = 360. 550 hdmp = 800. 551 / 552 553 !!====================================================================== 554 !! *** Dynamics namelists *** 555 !!====================================================================== 556 !! nam_dynadv formulation of the momentum advection 557 !! nam_dynvor advection scheme 558 !! nam_dynhpg hydrostatic pressure gradient 559 !! namflg hydrostatic pressure gradient time stepping 560 !! nam_dynspg surface pressure gradient (CPP key only) 561 !! nam_dynldf lateral diffusion scheme 562 !!====================================================================== 563 564 !----------------------------------------------------------------------- 565 ! nam_dynadv formulation of the momentum advection 566 !----------------------------------------------------------------------- 567 ! ln_dynadv_vec vector form (T) or flux form (F) (default T) 568 ! ln_dynadv_cen2 flux form - 2nd order centered scheme (default T) 569 ! ln_dynadv_ubs flux form - 3rd order UBS scheme (default F) 570 &nam_dynadv 571 ln_dynadv_vec = .true. 572 ln_dynadv_cen2 = .false. 573 ln_dynadv_ubs = .false. 574 / 575 !----------------------------------------------------------------------- 576 ! nam_dynvor option of physics/algorithm (not control by CPP keys) 577 !----------------------------------------------------------------------- 578 ! ln_dynvor_ens vorticity trends: enstrophy conserving scheme (default T) 579 ! ln_dynvor_ene " " : energy conserving scheme (default F) 580 ! ln_dynvor_mix " " : mixed scheme (default F) 581 ! ln_dynvor_een " " : energy & enstrophy scheme (default F) 582 &nam_dynvor 583 ln_dynvor_ene = .false. 584 ln_dynvor_ens = .false. 585 ln_dynvor_mix = .false. 586 ln_dynvor_een = .true. 587 / 588 !----------------------------------------------------------------------- 589 ! nam_dynhpg Hydrostatic pressure gradient option 590 !----------------------------------------------------------------------- 591 ! type of pressure gradient scheme (choose one only!) 592 ! ln_hpg_zco z-coordinate - full steps (default T) 593 ! ln_hpg_zps z-coordinate - partial steps (interpolation) 594 ! ln_hpg_sco s-coordinate (standard jacobian formulation) 595 ! ln_hpg_hel s-coordinate (helsinki modification) 596 ! ln_hpg_wdj s-coordinate (weighted density jacobian) 597 ! ln_hpg_djc s-coordinate (Density Jacobian with Cubic polynomial) 598 ! ln_hpg_rot s-coordinate (ROTated axes scheme) 599 ! parameters 600 ! gamm weighting coefficient (wdj scheme) 601 &nam_dynhpg 602 ln_hpg_zco = .false. 603 ln_hpg_zps = .true. 604 ln_hpg_sco = .false. 605 ln_hpg_hel = .false. 606 ln_hpg_wdj = .false. 607 ln_hpg_djc = .false. 608 ln_hpg_rot = .false. 609 gamm = 0.e0 610 / 611 !----------------------------------------------------------------------- 612 ! namflg algorithm flags (algorithm not control by CPP keys) 613 !----------------------------------------------------------------------- 614 ! ln_dynhpg_imp hydrostatic pressure gradient: semi-implicit time scheme (T) 615 ! centered time scheme (F) 616 ! nn_dynhpg_rst add dynhpg implicit variables in restart ot not (1/0) 617 &namflg 618 ln_dynhpg_imp = .false. 619 nn_dynhpg_rst = 0 620 / 621 !----------------------------------------------------------------------- 622 ! nam_dynspg surface pressure gradient (CPP key only) 623 !----------------------------------------------------------------------- 624 ! ! explicit free surface ("key_dynspg_exp") 625 ! ! filtered free surface ("key_dynspg_flt") 626 ! ! split-explicit free surface ("key_dynspg_ts") 627 ! ! rigid-lid ("key_dynspg_rl") 628 629 !----------------------------------------------------------------------- 630 ! nam_dynldf lateral diffusion on momentum 631 !----------------------------------------------------------------------- 632 ! Type of the operator : 633 ! ln_dynldf_lap laplacian operator (default T) 634 ! ln_dynldf_bilap bilaplacian operator (default F) 635 ! Direction of action : 636 ! ln_dynldf_level iso-level (default F) 637 ! ln_dynldf_hor horizontal (geopotential) (default F) (require "key_ldfslp" if ln_sco=T) 638 ! ln_dynldf_iso iso-neutral (default T) (require "key_ldfslp") 639 ! Coefficient 640 ! ahm0 horizontal eddy viscosity for the dynamics (m2/s) 641 ! ahmb0 background eddy viscosity for isopycnal diffusion (m2/s) 642 &nam_dynldf 643 ln_dynldf_lap = .false. 644 ln_dynldf_bilap = .true. 645 ln_dynldf_level = .false. 646 ln_dynldf_hor = .true. 647 ln_dynldf_iso = .false. 648 ahm0 = -8.5e+11 649 ahmb0 = 0. 650 / 651 652 !!====================================================================== 653 !! Tracers & Dynamics vertical physics namelists 654 !!====================================================================== 655 !! namzdf vertical physics 656 !! namnpc non penetrative convection 657 !! namric richardson number dependent vertical mixing ("key_zdfric" ) 658 !! namtke TKE dependent vertical mixing ("key_zdftke" ) 659 !! namkpp KPP dependent vertical mixing ("key_zdfkpp" ) 660 !! namddm double diffusive mixing parameterization ("key_zdfddm" ) 661 !!====================================================================== 662 663 !----------------------------------------------------------------------- 664 ! namzdf vertical physics 665 !----------------------------------------------------------------------- 666 ! ln_zdfevd enhanced vertical diffusion (default T) 667 ! ln_zdfnpc Non-Penetrative Convection (default T) 668 ! avm0 vertical eddy viscosity for the dynamic (m2/s) (also background Kz if not "key_zdfcst") 669 ! avt0 vertical eddy diffusivity for tracers (m2/s) (also background Kz if not "key_zdfcst") 670 ! avevd vertical coefficient for enhanced diffusion scheme (m2/s) 671 ! n_evdm = 0 apply enhanced mixing on tracer only 672 ! = 1 apply enhanced mixing on both tracer and momentum 673 ! ln_zdfexp vertical physics: (=T) time splitting (T) (Default=F) 674 ! (=F) euler backward (F) 675 ! n_zdfexp number of sub-timestep for time splitting scheme 676 &namzdf 677 ln_zdfevd = .true. 678 ln_zdfnpc = .false. 679 avm0 = 1.2e-4 680 avt0 = 1.2e-5 681 avevd = 100. 682 n_evdm = 1 683 ln_zdfexp = .false. 684 n_zdfexp = 3 685 / 686 !----------------------------------------------------------------------- 687 ! namnpc non penetrative convection 688 !----------------------------------------------------------------------- 689 ! nnpc1 non penetrative convective scheme frequency 690 ! nnpc2 non penetrative convective scheme print frequency 691 &namnpc 692 nnpc1 = 1 693 nnpc2 = 365 694 / 695 !----------------------------------------------------------------------- 696 ! namric richardson number dependent vertical diffusion ("key_zdfric" ) 697 !----------------------------------------------------------------------- 698 ! avmri maximum value of the vertical viscosity 699 ! alp coefficient of the parameterization 700 ! nric coefficient of the parameterization 701 &namric 702 avmri = 100.e-4 703 alp = 5. 704 nric = 2 705 / 706 !----------------------------------------------------------------------- 707 ! namtke turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 708 !----------------------------------------------------------------------- 709 ! ln_rstke flag to restart with tke from a run without tke (default F) 710 ! ediff coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 711 ! ediss coef. of the Kolmogoroff dissipation 712 ! ebb coef. of the surface input of tke 713 ! efave coef. to applied to the tke diffusion ( avtke=efave*avm ) 714 ! emin minimum value of tke (m^2/s^2) 715 ! emin0 surface minimum value of tke (m^2/s^2) 716 ! nitke number of restart iterative loops 717 ! ri_c critic richardson number = 2/9 = 0.22222222 (hard coded) 718 ! nmxl length used = 0 bounded by the distance to surface and bottom 719 ! ! = 1 bounded by the local vertical scale factor 720 ! ! = 2 first vertical derivative of mixing length bounded by 1 721 ! npdl Prandtl number = 0 no vertical prandtl number (avt=avm) 722 ! ! = 1 prandtl number function of richarson number (avt=pdl*avm) 723 ! ! = 2 same as = 1 but a shapiro filter is applied on pdl 724 ! nave = horizontal averaged (=1) or not (=0) of avt (default =1) 725 ! navb = 0 cst background avt0, avm0 / =1 profile used on avtb 726 &namtke 727 ln_rstke = .false. 728 ediff = 0.1 729 ediss = 0.7 730 ebb = 60. 731 efave = 1. 732 emin = 1.e-6 733 emin0 = 1.e-4 734 nitke = 50 735 nmxl = 2 736 npdl = 1 737 navb = 0 738 / 739 !----------------------------------------------------------------------- 740 ! namkpp K-Profile Parameterization dependent vertical mixing ("key_zdfkpp" ) 741 !----------------------------------------------------------------------- 742 ! ln_kpprimix shear instability mixing (default T) 743 ! difmiw constant internal wave viscosity (m2/s) 744 ! difsiw constant internal wave diffusivity (m2/s) 745 ! Riinfty local Richardson Number limit for shear instability 746 ! difri maximum shear mixing at Rig = 0 (m2/s) 747 ! bvsqcon Brunt-Vaisala squared (1/s**2) for maximum convection 748 ! difcon maximum mixing in interior convection (m2/s) 749 ! nave = 0/1 flag for horizontal average on avt, avmu, avmv 750 ! navb = 0/1 flag for constant or profile background avt 751 &namkpp 752 ln_kpprimix = .true. 753 difmiw = 1.e-04 754 difsiw = 0.1e-04 755 Riinfty = 0.8 756 difri = 0.0050 757 bvsqcon = -0.01e-07 758 difcon = 1. 759 navb = 0 760 nave = 1 761 / 762 !----------------------------------------------------------------------- 763 ! namddm double diffusive mixing parameterization ("key_zdfddm") 764 !----------------------------------------------------------------------- 765 ! avts maximum avs for dd mixing 766 ! hsbfr heat/salt buoyancy flux ratio 767 &namddm 768 avts = 1.e-4 769 hsbfr = 1.6 770 / 771 772 !!====================================================================== 773 !! *** Miscelaneous namelists *** 774 !!====================================================================== 775 !! namsol elliptic solver / island / free surface 776 !!====================================================================== 777 778 !----------------------------------------------------------------------- 779 ! namsol elliptic solver / island / free surface 780 !----------------------------------------------------------------------- 781 ! nsolv elliptic solver (=1 preconditioned conjugate gradient: pcg) 782 ! (=2 successive-over-relaxation: sor) 783 ! (=3 FETI currently it does not work! ("key_feti") 784 ! nsol_arp absolute/relative (0/1) precision convergence test 785 ! nmin minimum of iterations for the SOR solver 786 ! nmax maximum of iterations for the SOR solver 787 ! nmod frequency of test for the SOR solver 788 ! eps absolute precision of the solver 789 ! resmax absolute precision for the SOR solver 790 ! sor optimal coefficient for SOR solver 791 ! epsisl absolute precision on stream function solver 792 ! nmisl maximum pcg iterations for island 793 ! rnu strength of the additional force used in free surface b.c. 794 &namsol 795 nsolv = 1 796 nsol_arp = 0 797 nmin = 300 798 nmax = 800 799 nmod = 10 800 eps = 1.E-6 801 resmax = 1.E-10 802 sor = 1.92 803 epsisl = 1.e-10 804 nmisl = 4000 805 rnu = 1. 806 / 807 !!====================================================================== 808 !! *** Diagnostics namelists *** 809 !!====================================================================== 810 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 811 !! namgap level mean model-data gap ("key_diagap") 812 !! namspr surface pressure diagnosed in rigid-lid ("key_diaspr") 813 !! namflo float parameters ("key_float") 814 !! namptr Poleward Transport Diagnostics 815 !!====================================================================== 816 817 !----------------------------------------------------------------------- 818 ! namtrd diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 819 ! or mixed-layer trends ('key_trdmld') 820 ! or barotropic vorticity ("key_trdvor") 821 !----------------------------------------------------------------------- 822 ! ntrd time step frequency dynamics and tracers trends 823 ! nctls control surface type in mixed-layer trends (0,1 or n<jpk) 824 ! ln_trdmld_restart restart for ML diagnostics 825 ! ucf unit conversion factor (=1 -> /seconds | =86400. -> /day) 826 ! ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 827 &namtrd 828 ntrd = 365 829 nctls = 0 830 ln_trdmld_restart = .false. 831 ucf = 1. 832 ln_trdmld_instant = .false. 833 / 834 !----------------------------------------------------------------------- 835 ! namgap level mean model-data gap ('key_diagap') 836 !----------------------------------------------------------------------- 837 ! ngap time-step frequency of model-data gap computation 838 ! nprg time-step frequency of gap print in model output 839 &namgap 840 ngap = 15 841 nprg = 10 842 / 843 !----------------------------------------------------------------------- 844 ! namspr surface pressure diagnostic 845 !----------------------------------------------------------------------- 846 ! nmaxp maximum of iterations for the solver 847 ! epsp absolute precision of the solver 848 ! niterp number of iteration done by the solver 849 &namspr 850 nmaxp = 1000 851 epsp = 1.e-3 852 niterp = 400 853 / 854 !----------------------------------------------------------------------- 855 ! namflo float parameters ("key_float") 643 856 !----------------------------------------------------------------------- 644 857 ! ln_rstflo boolean term for float restart (true or false) … … 655 868 ln_flork4 = .false. 656 869 / 657 !!====================================================================== 658 !! Surface Boundary Condition namelists 659 !! 660 !! namsbc surface boundary condition 661 !! namsbc_ana analytical formulation of sbc 662 !! namsbc_flx flux formulation of sbc 663 !! namsbc_core CORE bulk formulea formulation of sbc 664 !! namsbc_clio CLIO bulk formulea formulation of sbc 665 !! namsbc_cpl CouPLed formulation of sbc 666 !! namsbc_ssr sea surface restoring term (for T and/or S) 667 !!====================================================================== 668 !----------------------------------------------------------------------- 669 ! namsbc Surface Boundary Condition (surface module) 670 !----------------------------------------------------------------------- 671 ! nn_fsbc frequency of surface boundary condition computation 672 ! (also correspond to the frequency of sea-ice model call) 673 ! 674 ! ln_ana analytical formulation (fill namsbc_ana ) 675 ! ln_flx flux formulation (fill namsbc_flx ) (overwritten by key_sbc_flux ) 676 ! ln_blk_clio CLIO bulk formulation (fill namsbc_core) (overwritten by key_sbc_bulk_clio) 677 ! ln_blk_core CORE bulk formulation (fill namsbc_clio) (overwritten by key_sbc_bulk_core) 678 ! ln_cpl Coupled formulation (fill namsbc_cpl ) (overwritten by key_sbc_coupled ) 679 ! 680 ! nn_ice =0 no ice boundary condition , 681 ! =1 observed ice-cover , 682 ! =2 ice-model ice-cover (overwritten by key_ice_lim) 683 ! 684 ! ln_dm2dc Daily mean to Diurnal Cycle short wave (qsr) 685 ! 686 ! ln_ssr Sea Surface restoring on SST and/or SSS (fill namsbc_ssr) 687 ! (overwritten by key_dta_sst and/or key_dtasss) 688 ! 689 ! nn_fwb FreshWater Budget: =0 no control , 690 ! =1 annual global mean of emp set to zero , 691 ! =2 global mean of emp set to zero at each nn_fsbc time step 692 &namsbc 693 nn_fsbc = 5 694 ln_ana = .false. 695 ln_flx = .false. 696 ln_blk_clio = .false. 697 ln_blk_core = .true. 698 ln_cpl = .false. 699 nn_ice = 2 700 ln_dm2dc = .false. 701 ln_rnf = .false. 702 ln_ssr = .false. 703 nn_fwb = 0 704 / 705 !----------------------------------------------------------------------- 706 ! namsbc_rnf runoffs namelist surface boundary condition 707 !----------------------------------------------------------------------- 708 !nn_runoff =0 no, 1 runoff, 2 runoff+river mouth ups adv 709 !rn_hrnf runoffs, depth over which enhanced vertical mixing is used 710 !rn_avt_rnf runoffs, value of the additional vertical mixing coef. [m2/s] 711 !cn_dir Root directory for location of ssr files 712 !sn_rnf informations about the runoff file to be read 713 !sn_cnf informations about the runoff mouth file to be read 714 ! THE ORDER OF THE FILES MATTER: 715 ! ! file name ! frequency (hours) ! variable ! time interpolation ! 716 ! ! ! (if <0 months) ! name ! (logical) ! 717 &namsbc_rnf 718 nn_runoff = 2 719 rn_hrnf = 0.e0 720 rn_avt_rnf = 1.e-3 721 cn_dir = './' 722 sn_rnf = 'runoff_1m_nomask.nc' , -12. , 'sorunoff', .TRUE. , 1 , 0 723 sn_cnf = 'runoff_1m_nomask.nc' , 0. , 'socoefr' , .FALSE. , 1 , 0 724 / 725 !----------------------------------------------------------------------- 726 ! namsbc_ana analytical surface boundary condition 727 !----------------------------------------------------------------------- 728 ! nn_tau000 gently increase the stress over the first ntau_rst time-steps 729 ! rn_utau0 uniform value used as default surface heat flux 730 ! rn_vtau0 uniform value used as default solar radiation flux 731 ! rn_q0 uniform value used as default surface heat flux 732 ! rn_qsr0 uniform value used as default solar radiation flux 733 ! rn_emp0 uniform value used as default surface freswater budget (E-P) 734 &namsbc_ana 735 nn_tau000 = 0 736 rn_utau0 = 0.5 737 rn_vtau0 = 0.e0 738 rn_q0 = 0.e0 739 rn_qsr0 = 0.e0 740 rn_emp0 = 0.e0 741 / 742 !----------------------------------------------------------------------- 743 ! namsbc_flx surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 744 !----------------------------------------------------------------------- 745 ! cn_dir directory in which the model is executed 746 ! 747 ! THE ORDER OF THE FILES MATTER: 748 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 749 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 750 &namsbc_flx 751 cn_dir = './' 752 sn_utau = 'utau.nc' , 24. , 'utau' , .FALSE. , 0 , 0 753 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .FALSE. , 0 , 0 754 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .FALSE. , 0 , 0 755 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .FALSE. , 0 , 0 756 sn_emp = 'emp.nc' , 24. , 'emp' , .FALSE. , 0 , 0 757 758 / 759 !----------------------------------------------------------------------- 760 ! namsbc_clio CLIO bulk formulea 761 !----------------------------------------------------------------------- 762 ! cn_dir directory in which the model is executed 763 ! 764 ! THE ORDER OF THE FILES MATTER: 765 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 766 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 767 &namsbc_clio 768 cn_dir = './' 769 sn_utau = 'taux_1m' , -12. , 'sozotaux', .FALSE. , 1 , 0 770 sn_vtau = 'tauy_1m' , -12. , 'sometauy', .FALSE. , 1 , 0 771 sn_wndm = 'flx' , -12. , 'socliowi', .FALSE. , 1 , 0 772 sn_tair = 'flx' , -12. , 'socliot1', .FALSE. , 1 , 0 773 sn_humi = 'flx' , -12. , 'socliohu', .FALSE. , 1 , 0 774 sn_ccov = 'flx' , -12. , 'socliocl', .FALSE. , 1 , 0 775 sn_prec = 'flx' , -12. , 'socliopl', .FALSE. , 1 , 0 776 / 777 !----------------------------------------------------------------------- 778 ! namsbc_core CORE bulk formulea 779 !----------------------------------------------------------------------- 780 ! cn_dir directory in which the model is executed 781 ! ln_2m logical flag to use air temp. and hum referenced at 2m instead 10m 782 ! alpha_precip multiplication factor for precipitation (total & snow) 783 ! 784 ! THE ORDER OF THE FILES MATTER: 785 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 786 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 787 &namsbc_core 788 cn_dir = './' 789 ln_2m = .FALSE. 790 alpha_precip = 1. 791 sn_wndi = 'u10' , 24. , 'U_10_MOD', .FALSE. , 1 , 0 792 sn_wndj = 'v10' , 24. , 'V_10_MOD', .FALSE. , 1 , 0 793 sn_qsr = 'rad' , 24. , 'SWDN_MOD', .FALSE. , 1 , 0 794 sn_qlw = 'rad' , 24. , 'LWDN_MOD', .FALSE. , 1 , 0 795 sn_tair = 't10' , 24. , 'T_10_MOD', .FALSE. , 1 , 0 796 sn_humi = 'q10' , 24. , 'Q_10_MOD', .FALSE. , 1 , 0 797 sn_prec = 'precip' , -12. , 'RAIN' , .FALSE. , 1 , 0 798 sn_snow = 'precip' , -12. , 'SNOW' , .FALSE. , 1 , 0 799 / 800 !----------------------------------------------------------------------- 801 ! namsbc_cpl coupled ocean/atmosphere model (#ifdef "key_coupled") 802 !----------------------------------------------------------------------- 803 &namsbc_cpl 804 / 805 !----------------------------------------------------------------------- 806 ! namsbc_ssr surface boundary condition : sea surface restoring 807 !----------------------------------------------------------------------- 808 ! nn_sstr =0/1 add a retroaction term in the surface heat flux 809 ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux 810 ! (=1, salt flux, concentration/dillution emps only) 811 ! (=2, volume flux, both emp and emps are updated) 812 ! dqdt magnitude of the retroaction on temperature [W/m2/K] 813 ! deds magnitude of the damping on salinity [mm/day] 814 ! THE ORDER OF THE FILES MATTER: 815 ! ! file name ! frequency (hours) ! variable ! time intepolation ! clim ! starting ! 816 ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! 817 &namsbc_ssr 818 nn_sstr = 0 819 nn_sssr = 0 820 dqdt = -40. 821 deds = 27.7 822 sn_sst = 'sst_data.nc' , 24. , 'sst' , .FALSE. , 0 , 0 823 sn_sss = 'sss_data.nc' , -12. , 'sss' , .TRUE. , 0 , 0 824 / 825 !----------------------------------------------------------------------- 826 ! namagrif agrif parameters (#ifdef key_agrif) 827 !----------------------------------------------------------------------- 828 ! nbclineupdate baroclinic update frequency 829 ! ln_spc_dyn use 0 as special value for dynamics 830 ! visc_tra viscosity coeeficient for tracers sponge layer 831 ! visc_dyn viscosity coeeficient for dynamics sponge layer 832 &namagrif 833 nbclineupdate = 3 834 ln_spc_dyn = .TRUE. 835 visc_tra = 2880. 836 visc_dyn = 2880. 837 / 870 !----------------------------------------------------------------------- 871 ! namptr Poleward Transport Diagnostic 872 !----------------------------------------------------------------------- 873 ! ln_diaptr logical flag for Poleward transport computation 874 ! ln_subbas logical flag for Atlantic/Pacific/Indian basins computation 875 ! need input basins mask file named "subbasins.nc" 876 ! nf_ptr Frequency of computation 877 &namptr 878 ln_diaptr = .false. 879 ln_subbas = .false. 880 nf_ptr = 15 881 /
Note: See TracChangeset
for help on using the changeset viewer.