Changeset 6041 for branches/2015/dev_r5776_UKMO2_OBS_efficiency_improvs/NEMOGCM/CONFIG/SHARED/namelist_ref
- Timestamp:
- 2015-12-14T10:06:06+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5776_UKMO2_OBS_efficiency_improvs/NEMOGCM/CONFIG/SHARED/namelist_ref
r6018 r6041 1 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 !! namelist_ref 1 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 2 4 !! NEMO/OPA : 1 - run manager (namrun) … … 5 7 !! namsbc_cpl, namtra_qsr, namsbc_rnf, 6 8 !! namsbc_apr, namsbc_ssr, namsbc_alb) 7 !! 4 - lateral boundary (namlbc, nam cla, namobc, namagrif, nambdy, nambdy_tide)9 !! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 8 10 !! 5 - bottom boundary (nambfr, nambbc, nambbl) 9 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ dmp)11 !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 10 12 !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 11 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ kpp, namzdf_ddm, namzdf_tmx)13 !! 8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 12 14 !! 9 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto) 13 !! 10 - miscellaneous (nam sol, nammpp, namctl)15 !! 10 - miscellaneous (nammpp, namctl) 14 16 !! 11 - Obs & Assim (namobs, nam_asminc) 15 17 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 161 163 / 162 164 !----------------------------------------------------------------------- 163 &namsplit ! time splitting parameters ("key_dynspg_ts")164 !-----------------------------------------------------------------------165 ln_bt_fw = .TRUE. ! Forward integration of barotropic equations166 ln_bt_av = .TRUE. ! Time filtering of barotropic variables167 ln_bt_nn_auto = .TRUE. ! Set nn_baro automatically to be just below168 ! a user defined maximum courant number (rn_bt_cmax)169 nn_baro = 30 ! Number of iterations of barotropic mode170 ! during rn_rdt seconds. Only used if ln_bt_nn_auto=F171 rn_bt_cmax = 0.8 ! Maximum courant number allowed if ln_bt_nn_auto=T172 nn_bt_flt = 1 ! Time filter choice173 ! = 0 None174 ! = 1 Boxcar over nn_baro barotropic steps175 ! = 2 Boxcar over 2*nn_baro " "176 /177 !-----------------------------------------------------------------------178 165 &namcrs ! Grid coarsening for dynamics output and/or 179 166 ! passive tracer coarsened online simulations … … 200 187 !----------------------------------------------------------------------- 201 188 &namtsd ! data : Temperature & Salinity 202 !-----------------------------------------------------------------------203 189 !----------------------------------------------------------------------- 204 190 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! … … 545 531 !!====================================================================== 546 532 !! namlbc lateral momentum boundary condition 547 !! namcla cross land advection548 533 !! namobc open boundaries parameters ("key_obc") 549 534 !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") … … 558 543 ! free slip ! partial slip ! no slip ! strong slip 559 544 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical eqs. 560 /561 !-----------------------------------------------------------------------562 &namcla ! cross land advection563 !-----------------------------------------------------------------------564 nn_cla = 0 ! advection between 2 ocean pts separates by land565 /566 !-----------------------------------------------------------------------567 &namobc ! open boundaries parameters ("key_obc")568 !-----------------------------------------------------------------------569 ln_obc_clim = .false. ! climatological obc data files (T) or not (F)570 ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F)571 ln_obc_fla = .false. ! Flather open boundary condition572 nn_obcdta = 1 ! = 0 the obc data are equal to the initial state573 ! = 1 the obc data are read in 'obc.dta' files574 cn_obcdta = 'annual' ! set to annual if obc datafile hold 1 year of data575 ! set to monthly if obc datafile hold 1 month of data576 rn_dpein = 1. ! damping time scale for inflow at east open boundary577 rn_dpwin = 1. ! - - - west - -578 rn_dpnin = 1. ! - - - north - -579 rn_dpsin = 1. ! - - - south - -580 rn_dpeob = 3000. ! time relaxation (days) for the east open boundary581 rn_dpwob = 15. ! - - - west - -582 rn_dpnob = 3000. ! - - - north - -583 rn_dpsob = 15. ! - - - south - -584 rn_volemp = 1. ! = 0 the total volume change with the surface flux (E-P-R)585 ! = 1 the total volume remains constant586 545 / 587 546 !----------------------------------------------------------------------- … … 711 670 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) 712 671 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] 713 714 672 / 715 673 !----------------------------------------------------------------------- … … 725 683 !! Tracer (T & S ) namelists 726 684 !!====================================================================== 727 !! nameos equation of state728 !! namtra_adv advection scheme685 !! nameos equation of state 686 !! namtra_adv advection scheme 729 687 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) 730 !! namtra_ldf lateral diffusion scheme 731 !! namtra_dmp T & S newtonian damping 688 !! namtra_ldf lateral diffusion scheme 689 !! namtra_ldfeiv eddy induced velocity param. 690 !! namtra_dmp T & S newtonian damping 732 691 !!====================================================================== 733 692 ! … … 740 699 ! = 1, S-EOS (simplified eos) 741 700 ln_useCT = .true. ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 742 !!701 ! 743 702 ! ! S-EOS coefficients : 744 !! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS703 ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 745 704 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) 746 705 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) … … 754 713 &namtra_adv ! advection scheme for tracer 755 714 !----------------------------------------------------------------------- 756 ln_traadv_cen2 = .false. ! 2nd order centered scheme 757 ln_traadv_tvd = .true. ! TVD scheme 758 ln_traadv_muscl = .false. ! MUSCL scheme 759 ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries 760 ln_traadv_ubs = .false. ! UBS scheme 761 ln_traadv_qck = .false. ! QUICKEST scheme 762 ln_traadv_msc_ups= .false. ! use upstream scheme within muscl 763 ln_traadv_tvd_zts= .false. ! TVD scheme with sub-timestepping of vertical tracer advection 715 ln_traadv_cen = .false. ! 2nd order centered scheme 716 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN 717 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 718 ln_traadv_fct = .false. ! FCT scheme 719 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 720 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 721 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping 722 ! ! (number of sub-timestep = nn_fct_zts) 723 ln_traadv_mus = .false. ! MUSCL scheme 724 ln_mus_ups = .false. ! use upstream scheme near river mouths 725 ln_traadv_ubs = .false. ! UBS scheme 726 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT 727 ln_traadv_qck = .false. ! QUICKEST scheme 764 728 / 765 729 !----------------------------------------------------------------------- 766 730 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 767 731 !----------------------------------------------------------------------- 768 ln_mle = . true. ! (T) use the Mixed Layer Eddy (MLE) parameterisation732 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 769 733 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) 770 734 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation … … 780 744 !---------------------------------------------------------------------------------- 781 745 ! ! Operator type: 782 ln_traldf_lap = .true. ! laplacian operator 783 ln_traldf_bilap = .false. ! bilaplacian operator 746 ! ! no diffusion: set ln_traldf_lap=..._blp=F 747 ln_traldf_lap = .false. ! laplacian operator 748 ln_traldf_blp = .false. ! bilaplacian operator 784 749 ! ! Direction of action: 785 ln_traldf_level = .false. ! iso-level 786 ln_traldf_hor = .false. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) 787 ln_traldf_iso = .true. ! iso-neutral (needs "key_ldfslp") 788 ! ! Griffies parameters (all need "key_ldfslp") 789 ln_traldf_grif = .false. ! use griffies triads 790 ln_traldf_gdia = .false. ! output griffies eddy velocities 791 ln_triad_iso = .false. ! pure lateral mixing in ML 792 ln_botmix_grif = .false. ! lateral mixing on bottom 793 ! ! Coefficients 794 ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 795 ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 796 ! (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 797 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 798 rn_aht_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] 799 rn_ahtb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] 800 ! (normally=0; not used with Griffies) 801 rn_slpmax = 0.01 ! slope limit 802 rn_chsmag = 1. ! multiplicative factor in Smagorinsky diffusivity 803 rn_smsh = 1. ! Smagorinsky diffusivity: = 0 - use only sheer 804 rn_aht_m = 2000. ! upper limit or stability criteria for lateral eddy diffusivity (m2/s) 750 ln_traldf_lev = .false. ! iso-level 751 ln_traldf_hor = .false. ! horizontal (geopotential) 752 ln_traldf_iso = .false. ! iso-neutral (standard operator) 753 ln_traldf_triad = .false. ! iso-neutral (triad operator) 754 ! 755 ! ! iso-neutral options: 756 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 757 rn_slpmax = 0.01 ! slope limit (both operators) 758 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) 759 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) 760 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) 761 ! 762 ! ! Coefficients: 763 nn_aht_ijk_t = 0 ! space/time variation of eddy coef 764 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 765 ! ! = 0 constant 766 ! ! = 10 F(k) =ldf_c1d 767 ! ! = 20 F(i,j) =ldf_c2d 768 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 769 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 770 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 771 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s] 772 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] 773 / 774 !---------------------------------------------------------------------------------- 775 &namtra_ldfeiv ! eddy induced velocity param. 776 !---------------------------------------------------------------------------------- 777 ln_ldfeiv =.false. ! use eddy induced velocity parameterization 778 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities 779 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s] 780 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient 781 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 782 ! ! = 0 constant 783 ! ! = 10 F(k) =ldf_c1d 784 ! ! = 20 F(i,j) =ldf_c2d 785 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 786 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d 805 787 / 806 788 !----------------------------------------------------------------------- … … 820 802 !! namdyn_vor advection scheme 821 803 !! namdyn_hpg hydrostatic pressure gradient 822 !! namdyn_spg surface pressure gradient (CPP key only)804 !! namdyn_spg surface pressure gradient 823 805 !! namdyn_ldf lateral diffusion scheme 824 806 !!====================================================================== … … 853 835 ln_dynvor_ens = .false. ! energy conserving scheme 854 836 ln_dynvor_mix = .false. ! mixed scheme 855 ln_dynvor_een = .true. ! energy & enstrophy scheme 856 ln_dynvor_een_old = .false. ! energy & enstrophy scheme - original formulation 837 ln_dynvor_een = .false. ! energy & enstrophy scheme 838 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 839 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT USE 857 840 / 858 841 !----------------------------------------------------------------------- … … 865 848 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) 866 849 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) 867 ln_dynhpg_imp = .false. ! time stepping: semi-implicit time scheme (T) 868 ! centered time scheme (F) 869 / 870 !----------------------------------------------------------------------- 871 !namdyn_spg ! surface pressure gradient (CPP key only) 872 !----------------------------------------------------------------------- 873 ! ! explicit free surface ("key_dynspg_exp") 874 ! ! filtered free surface ("key_dynspg_flt") 875 ! ! split-explicit free surface ("key_dynspg_ts") 876 850 / 851 !----------------------------------------------------------------------- 852 &namdyn_spg ! surface pressure gradient 853 !----------------------------------------------------------------------- 854 ln_dynspg_exp = .false. ! explicit free surface 855 ln_dynspg_ts = .false. ! split-explicit free surface 856 ln_bt_fw = .true. ! Forward integration of barotropic Eqs. 857 ln_bt_av = .true. ! Time filtering of barotropic variables 858 nn_bt_flt = 1 ! Time filter choice = 0 None 859 ! ! = 1 Boxcar over nn_baro sub-steps 860 ! ! = 2 Boxcar over 2*nn_baro " " 861 ln_bt_auto = .true. ! Number of sub-step defined from: 862 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed 863 nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds 864 / 877 865 !----------------------------------------------------------------------- 878 866 &namdyn_ldf ! lateral diffusion on momentum 879 867 !----------------------------------------------------------------------- 880 868 ! ! Type of the operator : 881 ln_dynldf_lap = .true. ! laplacian operator 882 ln_dynldf_bilap = .false. ! bilaplacian operator 869 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 870 ln_dynldf_lap = .false. ! laplacian operator 871 ln_dynldf_blp = .false. ! bilaplacian operator 883 872 ! ! Direction of action : 884 ln_dynldf_lev el = .false.! iso-level885 ln_dynldf_hor = .true. ! horizontal (geopotential) (require "key_ldfslp" in s-coord.)886 ln_dynldf_iso = .false. ! iso-neutral (require "key_ldfslp")873 ln_dynldf_lev = .false. ! iso-level 874 ln_dynldf_hor = .false. ! horizontal (geopotential) 875 ln_dynldf_iso = .false. ! iso-neutral 887 876 ! ! Coefficient 888 rn_ahm_0_lap = 40000. ! horizontal laplacian eddy viscosity [m2/s] 889 rn_ahmb_0 = 0. ! background eddy viscosity for ldf_iso [m2/s] 890 rn_ahm_0_blp = 0. ! horizontal bilaplacian eddy viscosity [m4/s] 891 rn_cmsmag_1 = 3. ! constant in laplacian Smagorinsky viscosity 892 rn_cmsmag_2 = 3 ! constant in bilaplacian Smagorinsky viscosity 893 rn_cmsh = 1. ! 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 894 rn_ahm_m_blp = -1.e12 ! upper limit for bilap abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 895 rn_ahm_m_lap = 40000. ! upper limit for lap ahm < min(dx^2/16rdt, rn_ahm_m_lap) 877 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef 878 ! ! =-30 read in eddy_viscosity_3D.nc file 879 ! ! =-20 read in eddy_viscosity_2D.nc file 880 ! ! = 0 constant 881 ! ! = 10 F(k)=c1d 882 ! ! = 20 F(i,j)=F(grid spacing)=c2d 883 ! ! = 30 F(i,j,k)=c2d*c1d 884 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 885 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] 886 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] 887 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] 888 ! 889 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 896 890 / 897 891 … … 902 896 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric") 903 897 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke") 904 !! namzdf_kpp KPP dependent vertical mixing ("key_zdfkpp")905 898 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm") 906 899 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx") … … 964 957 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 965 958 / 966 !------------------------------------------------------------------------967 &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally:968 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")969 ln_kpprimix = .true. ! shear instability mixing970 rn_difmiw = 1.0e-04 ! constant internal wave viscosity [m2/s]971 rn_difsiw = 0.1e-04 ! constant internal wave diffusivity [m2/s]972 rn_riinfty = 0.8 ! local Richardson Number limit for shear instability973 rn_difri = 0.0050 ! maximum shear mixing at Rig = 0 [m2/s]974 rn_bvsqcon = -0.01e-07 ! Brunt-Vaisala squared for maximum convection [1/s2]975 rn_difcon = 1. ! maximum mixing in interior convection [m2/s]976 nn_avb = 0 ! horizontal averaged (=1) or not (=0) on avt and amv977 nn_ave = 1 ! constant (=0) or profile (=1) background on avt978 /979 959 !----------------------------------------------------------------------- 980 960 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") … … 1015 995 !! *** Miscellaneous namelists *** 1016 996 !!====================================================================== 1017 !! namsol elliptic solver / island / free surface1018 997 !! nammpp Massively Parallel Processing ("key_mpp_mpi) 1019 998 !! namctl Control prints & Benchmark … … 1024 1003 !!====================================================================== 1025 1004 ! 1026 !-----------------------------------------------------------------------1027 &namsol ! elliptic solver / island / free surface1028 !-----------------------------------------------------------------------1029 nn_solv = 1 ! elliptic solver: =1 preconditioned conjugate gradient (pcg)1030 ! =2 successive-over-relaxation (sor)1031 nn_sol_arp = 0 ! absolute/relative (0/1) precision convergence test1032 rn_eps = 1.e-6 ! absolute precision of the solver1033 nn_nmin = 300 ! minimum of iterations for the SOR solver1034 nn_nmax = 800 ! maximum of iterations for the SOR solver1035 nn_nmod = 10 ! frequency of test for the SOR solver1036 rn_resmax = 1.e-10 ! absolute precision for the SOR solver1037 rn_sor = 1.92 ! optimal coefficient for SOR solver (to be adjusted with the domain)1038 /1039 1005 !----------------------------------------------------------------------- 1040 1006 &nammpp ! Massively Parallel Processing ("key_mpp_mpi)
Note: See TracChangeset
for help on using the changeset viewer.