Changeset 5260 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/SHARED/namelist_ref
- Timestamp:
- 2015-05-12T12:37:15+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/CONFIG/SHARED/namelist_ref
r4756 r5260 1 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! NEMO/OPA : 1 - run manager (namrun , namcfg)3 !! namelists 2 - Domain (nam zgr, namzgr_sco, namdom, namtsd)2 !! NEMO/OPA : 1 - run manager (namrun) 3 !! namelists 2 - Domain (namcfg, namzgr, namzgr_sco, namdom, namtsd) 4 4 !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 5 5 !! namsbc_cpl, namtra_qsr, namsbc_rnf, … … 10 10 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 11 11 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 12 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, nam ptr, namhsb)12 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb) 13 13 !! 10 - miscellaneous (namsol, nammpp, namctl) 14 14 !! 11 - Obs & Assim (namobs, nam_asminc) … … 18 18 !! *** Run management namelists *** 19 19 !!====================================================================== 20 !! namrun 20 !! namrun parameters of the run 21 21 !!====================================================================== 22 22 ! … … 31 31 nn_leapy = 0 ! Leap year calendar (1) or not (0) 32 32 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 33 nn_euler = 1 ! = 0 : start with forward time step if ln_rstart= .true.34 nn_rstctl = 0 ! restart control = > activated only if ln_rstart =T33 nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T 34 nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T 35 35 ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 36 36 ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart … … 47 47 / 48 48 ! 49 !----------------------------------------------------------------------- 50 &namcfg ! default parameters of the configuration 49 !!====================================================================== 50 !! *** Domain namelists *** 51 !!====================================================================== 52 !! namcfg parameters of the configuration 53 !! namzgr vertical coordinate 54 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 55 !! namdom space and time domain (bathymetry, mesh, timestep) 56 !! namtsd data: temperature & salinity 57 !!====================================================================== 58 ! 59 !----------------------------------------------------------------------- 60 &namcfg ! parameters of the configuration 51 61 !----------------------------------------------------------------------- 52 62 cp_cfg = "default" ! name of the configuration 53 cp_cfz = ''! name of the zoom of configuration63 cp_cfz = "no zoom" ! name of the zoom of configuration 54 64 jp_cfg = 0 ! resolution of the configuration 55 65 jpidta = 10 ! 1st lateral dimension ( >= jpi ) … … 57 67 jpkdta = 31 ! number of levels ( >= jpk ) 58 68 jpiglo = 10 ! 1st dimension of global domain --> i =jpidta 59 jpjglo = 12 ! 2nd - - --> j 69 jpjglo = 12 ! 2nd - - --> j =jpjdta 60 70 jpizoom = 1 ! left bottom (i,j) indices of the zoom 61 71 jpjzoom = 1 ! in data domain indices … … 66 76 ! = 5 North fold F-point pivot 67 77 ! = 6 cyclic East-West AND North fold F-point pivot 68 / 69 !!====================================================================== 70 !! *** Domain namelists *** 71 !!====================================================================== 72 !! namzgr vertical coordinate 73 !! namzgr_sco s-coordinate or hybrid z-s-coordinate 74 !! namdom space and time domain (bathymetry, mesh, timestep) 75 !! namtsd data: temperature & salinity 76 !!====================================================================== 77 ! 78 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 79 ! in netcdf input files, as the start j-row for reading 80 / 78 81 !----------------------------------------------------------------------- 79 82 &namzgr ! vertical coordinate … … 82 85 ln_zps = .true. ! z-coordinate - partial steps (T/F) 83 86 ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F) 87 ln_isfcav = .false. ! ice shelf cavity (T/F) 84 88 / 85 89 !----------------------------------------------------------------------- … … 183 187 / 184 188 !----------------------------------------------------------------------- 189 &namc1d ! 1D configuration options ("key_c1d") 190 !----------------------------------------------------------------------- 191 rn_lat1d = 50 ! Column latitude (default at PAPA station) 192 rn_lon1d = -145 ! Column longitude (default at PAPA station) 193 ln_c1d_locpt= .true. ! Localization of 1D config in a grid (T) or independant point (F) 194 / 195 !----------------------------------------------------------------------- 185 196 &namtsd ! data : Temperature & Salinity 186 197 !----------------------------------------------------------------------- … … 208 219 !! namtra_qsr penetrative solar radiation 209 220 !! namsbc_rnf river runoffs 221 !! namsbc_isf ice shelf melting/freezing 210 222 !! namsbc_apr Atmospheric Pressure 211 223 !! namsbc_ssr sea surface restoring term (for T and/or S) … … 223 235 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core) 224 236 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs ) 225 ln_cpl = .false. ! Coupled formulation (T => fill namsbc_cpl )226 237 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 227 238 nn_ice = 2 ! =0 no ice boundary condition , 228 239 ! =1 use observed ice-cover , 229 ! =2 ice-model used ("key_lim3" or "key_lim2 )240 ! =2 ice-model used ("key_lim3" or "key_lim2") 230 241 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect) 231 242 ! =1 levitating ice with mass and salt exchange but no presure effect 232 243 ! =2 embedded sea-ice (full salt and mass exchanges and pressure) 233 244 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 234 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 245 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 246 nn_isf = 0 ! ice shelf melting/freezing (/=0 => fill namsbc_isf) 247 ! 0 =no isf 1 = presence of ISF 248 ! 2 = bg03 parametrisation 3 = rnf file for isf 249 ! 4 = ISF fwf specified 250 ! option 1 and 4 need ln_isfcav = .true. (domzgr) 235 251 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 236 252 nn_fwb = 2 ! FreshWater Budget: =0 unchecked … … 240 256 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave) 241 257 ln_sdw = .false. ! Computation of 3D stokes drift (T => fill namsbc_wave) 242 cn_iceflx = 'linear' ! redistribution of solar input into ice categories during coupling ice/atm. 258 nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 259 ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 260 nn_limflx = -1 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 261 ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 262 ! = 0 Average per-category fluxes (forced and coupled mode) 263 ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 264 ! = 2 Redistribute a single flux over categories (coupled mode only) 243 265 / 244 266 !----------------------------------------------------------------------- … … 296 318 297 319 cn_dir = './' ! root directory for the location of the bulk files 298 ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F)299 320 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 300 ln_bulk2z = .false. ! Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 301 rn_zqt = 3. ! Air temperature and humidity reference height (m) (ln_bulk2z) 302 rn_zu = 4. ! Wind vector reference height (m) (ln_bulk2z) 321 rn_zqt = 10. ! Air temperature and humidity reference height (m) 322 rn_zu = 10. ! Wind vector reference height (m) 303 323 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 304 324 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) … … 327 347 ! ! ! categories ! reference ! orientation ! grids ! 328 348 ! send 329 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , ''330 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , ''331 sn_snd_thick = 'none' , 'no' , '' , '' , ''332 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T'333 sn_snd_co2 = 'coupled' , 'no' , '' , '' , ''349 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' 350 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' 351 sn_snd_thick = 'none' , 'no' , '' , '' , '' 352 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 353 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 334 354 ! receive 335 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 336 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 337 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 338 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 339 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 340 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 341 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 342 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 343 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 344 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 355 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' 356 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' 357 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V' 358 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' 359 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' 360 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' 361 sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' 362 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' 363 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 364 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 365 ! 366 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentialy sending/receiving data 367 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 368 ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 345 369 / 346 370 !----------------------------------------------------------------------- … … 396 420 ln_rnf_tem = .false. ! read in temperature information for runoff 397 421 ln_rnf_sal = .false. ! read in salinity information for runoff 422 / 423 !----------------------------------------------------------------------- 424 &namsbc_isf ! Top boundary layer (ISF) 425 !----------------------------------------------------------------------- 426 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 427 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 428 ! nn_isf == 4 429 sn_qisf = 'rnfisf' , -12 ,'sohflisf', .false. , .true. , 'yearly' , '' , '' 430 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' 431 ! nn_isf == 3 432 sn_rnfisf = 'runoffs' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' 433 ! nn_isf == 2 and 3 434 sn_depmax_isf = 'runoffs' , -12 ,'sozisfmax' , .false. , .true. , 'yearly' , '' , '' 435 sn_depmin_isf = 'runoffs' , -12 ,'sozisfmin' , .false. , .true. , 'yearly' , '' , '' 436 ! nn_isf == 2 437 sn_Leff_isf = 'rnfisf' , 0 ,'Leff' , .false. , .true. , 'yearly' , '' , '' 438 ! for all case 439 ln_divisf = .true. ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) 440 ! only for nn_isf = 1 or 2 441 rn_gammat0 = 1.0e-4 ! gammat coefficient used in blk formula 442 rn_gammas0 = 1.0e-4 ! gammas coefficient used in blk formula 443 ! only for nn_isf = 1 444 nn_isfblk = 1 ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 445 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 446 ! 0 => thickness of the tbl = thickness of the first wet cell 447 ln_conserve = .true. ! conservative case (take into account meltwater advection) 448 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) 449 ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 450 ! if you want to keep the cd as in global config, adjust rn_gammat0 to compensate 451 ! 2 = velocity and stability dependent Gamma Holland et al. 1999 398 452 / 399 453 !----------------------------------------------------------------------- … … 528 582 ln_tide_ramp = .false. ! 529 583 rdttideramp = 0. ! 530 clname(1) = 'M2' ! name of constituent 531 clname(2) = 'S2' 532 clname(3) = 'N2' 533 clname(4) = 'K1' 534 clname(5) = 'O1' 535 clname(6) = 'Q1' 536 clname(7) = 'M4' 537 clname(8) = 'K2' 538 clname(9) = 'P1' 539 clname(10) = 'Mf' 540 clname(11) = 'Mm' 584 clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg 541 585 / 542 586 !----------------------------------------------------------------------- … … 559 603 nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state 560 604 ! = 1, bdy data are read in 'bdydata .nc' files 605 cn_ice_lim = 'none' ! 606 nn_ice_lim_dta = 0 ! = 0, bdy data are equal to the initial state 607 ! = 1, bdy data are read in 'bdydata .nc' files 608 rn_ice_tem = 270. ! lim3 only: arbitrary temperature of incoming sea ice 609 rn_ice_sal = 10. ! lim3 only: -- salinity -- 610 rn_ice_age = 30. ! lim3 only: -- age -- 611 561 612 ln_tra_dmp =.false. ! open boudaries conditions for tracers 562 613 ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities … … 579 630 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true. , .false. , 'daily' , '' , '' , '' 580 631 bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline' , .true. , .false. , 'daily' , '' , '' , '' 632 ! for lim2 633 ! bn_frld = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 634 ! bn_hicif = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 635 ! bn_hsnif = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 636 ! for lim3 637 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 638 ! bn_ht_i = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 639 ! bn_ht_s = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 581 640 cn_dir = 'bdydta/' 582 641 ln_full_vel = .false. … … 607 666 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 608 667 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T 609 ln_loglayer = .false. ! logarithmic formulation (non linear case)610 668 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 611 669 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) 670 rn_tfri1 = 4.e-4 ! top drag coefficient (linear case) 671 rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 672 rn_tfri2_max = 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T) 673 rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2) 674 rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T 675 ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file ) 676 rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T) 677 612 678 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) 679 ln_loglayer = .false. ! logarithmic formulation (non linear case) 613 680 / 614 681 !----------------------------------------------------------------------- … … 635 702 !! nameos equation of state 636 703 !! namtra_adv advection scheme 704 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 637 705 !! namtra_ldf lateral diffusion scheme 638 706 !! namtra_dmp T & S newtonian damping … … 642 710 &nameos ! ocean physical parameters 643 711 !----------------------------------------------------------------------- 644 nn_eos = 0 ! type of equation of state and Brunt-Vaisala frequency 645 ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 646 ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) 647 ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 648 rn_alpha = 2.0e-4 ! thermal expension coefficient (nn_eos= 1 or 2) 649 rn_beta = 7.7e-4 ! saline expension coefficient (nn_eos= 2) 712 nn_eos = -1 ! type of equation of state and Brunt-Vaisala frequency 713 ! =-1, TEOS-10 714 ! = 0, EOS-80 715 ! = 1, S-EOS (simplified eos) 716 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 717 ! ! 718 ! ! S-EOS coefficients : 719 ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 720 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 721 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) 722 rn_lambda1 = 5.9520e-2 ! cabbeling coeff in T^2 (=0 for linear eos) 723 rn_lambda2 = 7.4914e-4 ! cabbeling coeff in S^2 (=0 for linear eos) 724 rn_mu1 = 1.4970e-4 ! thermobaric coeff. in T (=0 for linear eos) 725 rn_mu2 = 1.1090e-5 ! thermobaric coeff. in S (=0 for linear eos) 726 rn_nu = 2.4341e-3 ! cabbeling coeff in T*S (=0 for linear eos) 650 727 / 651 728 !----------------------------------------------------------------------- 652 729 &namtra_adv ! advection scheme for tracer 653 730 !----------------------------------------------------------------------- 654 ln_traadv_cen2 = .false. ! 2nd order centered scheme 655 ln_traadv_tvd = .true. ! TVD scheme 656 ln_traadv_muscl = .false. ! MUSCL scheme 657 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 658 ln_traadv_ubs = .false. ! UBS scheme 659 ln_traadv_qck = .false. ! QUICKEST scheme 660 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 731 ln_traadv_cen2 = .false. ! 2nd order centered scheme 732 ln_traadv_tvd = .true. ! TVD scheme 733 ln_traadv_muscl = .false. ! MUSCL scheme 734 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 735 ln_traadv_ubs = .false. ! UBS scheme 736 ln_traadv_qck = .false. ! QUICKEST scheme 737 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 738 ln_traadv_tvd_zts= .false. ! TVD scheme with sub-timestepping of vertical tracer advection 661 739 / 662 740 !----------------------------------------------------------------------- … … 705 783 !----------------------------------------------------------------------- 706 784 ln_tradmp = .true. ! add a damping termn (T) or not (F) 707 nn_hdmp = -1 ! horizontal shape =-1, damping in Med and Red Seas only708 ! =XX, damping poleward of XX degrees (XX>0)709 ! + F(distance-to-coast) + Red and Med Seas710 785 nn_zdmp = 0 ! vertical shape =0 damping throughout the water column 711 786 ! =1 no damping in the mixing layer (kz criteria) 712 787 ! =2 no damping in the mixed layer (rho crieria) 713 rn_surf = 50. ! surface time scale of damping [days] 714 rn_bot = 360. ! bottom time scale of damping [days] 715 rn_dep = 800. ! depth of transition between rn_surf and rn_bot [meters] 716 nn_file = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) 788 cn_resto = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this) 717 789 / 718 790 … … 733 805 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme 734 806 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 807 ln_dynzad_zts = .false. ! Use (T) sub timestepping for vertical momentum advection 735 808 / 736 809 !----------------------------------------------------------------------- … … 755 828 ln_dynvor_mix = .false. ! mixed scheme 756 829 ln_dynvor_een = .true. ! energy & enstrophy scheme 830 ln_dynvor_een_old = .false. ! energy & enstrophy scheme - original formulation 757 831 / 758 832 !----------------------------------------------------------------------- … … 762 836 ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation) 763 837 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) 838 ln_hpg_isf = .false. ! s-coordinate (sco ) adapted to isf 764 839 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 765 840 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) … … 879 954 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") 880 955 !----------------------------------------------------------------------- 881 rn_emin = 1.e- 6! minimum value of e [m2/s2]956 rn_emin = 1.e-7 ! minimum value of e [m2/s2] 882 957 rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] 883 958 ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) 884 rn_clim_galp = 0.53 ! galperin limit 885 ln_crban = .true. ! Use Craig & Banner (1994) surface wave mixing parametrisation 959 rn_clim_galp = 0.267 ! galperin limit 886 960 ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case 887 961 rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux 888 962 rn_charn = 70000. ! Charnock constant for wb induced roughness length 889 nn_tkebc_surf = 1 ! surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 890 nn_tkebc_bot = 1 ! bottom tke condition (0/1=Dir/Neum) 891 nn_psibc_surf = 1 ! surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 892 nn_psibc_bot = 1 ! bottom psi condition (0/1=Dir/Neum) 893 nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 894 nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 963 rn_hsro = 0.02 ! Minimum surface roughness 964 rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met=2) 965 nn_z0_met = 2 ! Method for surface roughness computation (0/1/2) 966 nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum) 967 nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum) 968 nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 969 nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 895 970 / 896 971 !----------------------------------------------------------------------- … … 914 989 !! *** Miscellaneous namelists *** 915 990 !!====================================================================== 991 !! namsol elliptic solver / island / free surface 916 992 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 917 993 !! namctl Control prints & Benchmark 918 !! namsol elliptic solver / island / free surface 994 !! namc1d 1D configuration options ("key_c1d") 995 !! namc1d_uvd data: U & V currents ("key_c1d") 996 !! namc1d_dyndmp U & V newtonian damping ("key_c1d") 919 997 !!====================================================================== 920 998 ! … … 959 1037 / 960 1038 !----------------------------------------------------------------------- 961 &namc1d ! 1D configuration options ("key_c1d")962 !-----------------------------------------------------------------------963 rn_lat = 50 ! Column latitude964 rn_lon = -145 ! Column longitude965 /966 !-----------------------------------------------------------------------967 1039 &namc1d_uvd ! data: U & V currents ("key_c1d") 968 1040 !----------------------------------------------------------------------- … … 981 1053 ln_dyndmp = .false. ! add a damping term (T) or not (F) 982 1054 / 1055 983 1056 !!====================================================================== 984 1057 !! *** Diagnostics namelists *** 985 1058 !!====================================================================== 986 1059 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 987 !! namtrd dynamics and/or tracer trends ("key_trddyn","key_trdtra","key_trdmld") 1060 !! namtrd dynamics and/or tracer trends 1061 !! namptr Poleward Transport Diagnostics 988 1062 !! namflo float parameters ("key_float") 989 !! namptr Poleward Transport Diagnostics990 1063 !! namhsb Heat and salt budgets 991 1064 !!====================================================================== … … 1003 1076 / 1004 1077 !----------------------------------------------------------------------- 1005 &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") 1006 ! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor") 1007 !----------------------------------------------------------------------- 1008 nn_trd = 365 ! time step frequency dynamics and tracers trends 1009 nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 1010 rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 1011 cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 1012 cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 1013 ln_trdmld_restart = .false. ! restart for ML diagnostics 1014 ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 1015 / 1078 &namtrd ! diagnostics on dynamics and/or tracer trends 1079 ! ! and/or mixed-layer trends and/or barotropic vorticity 1080 !----------------------------------------------------------------------- 1081 ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE 1082 ln_dyn_trd = .false. ! (T) 3D momentum trend output 1083 ln_dyn_mxl = .FALSE. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 1084 ln_vor_trd = .FALSE. ! (T) 2D barotropic vorticity trends (not coded yet) 1085 ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends 1086 ln_PE_trd = .false. ! (T) 3D Potential Energy trends 1087 ln_tra_trd = .FALSE. ! (T) 3D tracer trend output 1088 ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 1089 nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step) 1090 / 1091 !!gm nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 1092 !!gm rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 1093 !!gm cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 1094 !!gm cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 1095 !!gm ln_trdmld_restart = .false. ! restart for ML diagnostics 1096 !!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 1097 !!gm 1016 1098 !----------------------------------------------------------------------- 1017 1099 &namflo ! float parameters ("key_float") … … 1032 1114 !----------------------------------------------------------------------- 1033 1115 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1034 ln_diaznl = .true. ! Add zonal means and meridional stream functions 1035 ln_subbas = .true. ! Atlantic/Pacific/Indian basins computation (T) or not 1036 ! (orca configuration only, need input basins mask file named "subbasins.nc" 1037 ln_ptrcomp = .true. ! Add decomposition : overturning 1038 nn_fptr = 1 ! Frequency of ptr computation [time step] 1039 nn_fwri = 15 ! Frequency of ptr outputs [time step] 1116 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 1040 1117 / 1041 1118 !----------------------------------------------------------------------- … … 1085 1162 ! ln_ssh Logical switch for SSH observations 1086 1163 1087 ln_sst = . true. ! Logical switch for SST observations1088 ln_reysst = . true. ! ln_reysst Logical switch for Reynolds observations1164 ln_sst = .false. ! Logical switch for SST observations 1165 ln_reysst = .false. ! ln_reysst Logical switch for Reynolds observations 1089 1166 ln_ghrsst = .false. ! ln_ghrsst Logical switch for GHRSST observations 1090 1167 1091 1168 ln_sstfb = .false. ! Logical switch for feedback SST data 1092 1169 ! ln_sss Logical switch for SSS observations 1093 ! ln_seaiceLogical switch for Sea Ice observations1170 ln_seaice = .false. ! Logical switch for Sea Ice observations 1094 1171 ! ln_vel3d Logical switch for velocity observations 1095 1172 ! ln_velavcur Logical switch for velocity daily av. cur. … … 1112 1189 ! sstfiles GHRSST input observation file name 1113 1190 ! ! sstfbfiles: Feedback SST input observation file name 1114 sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 1115 ! seaicefiles Sea Ice input observation file name 1191 sstfbfiles = 'sst_01.nc' 1192 ! seaicefiles Sea Ice input observation file names 1193 seaicefiles = 'seaice_01.nc' 1116 1194 ! velavcurfiles Vel. cur. daily av. input file name 1117 1195 ! velhvcurfiles Vel. cur. high freq. input file name
Note: See TracChangeset
for help on using the changeset viewer.