- Timestamp:
- 2019-12-10T15:44:23+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/cfgs/SHARED/namelist_ref
r12113 r12154 5 5 !! namelists 2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 6 6 !! namsbc_sas, namtra_qsr, namsbc_rnf, 7 !! namsbc_isf, namsbc_iscpl, namsbc_apr, 7 !! namsbc_isf, namsbc_iscpl, namsbc_apr, 8 8 !! namsbc_ssr, namsbc_wave, namberg) 9 9 !! 3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) … … 65 65 ln_clobber = .true. ! clobber (overwrite) an existing file 66 66 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 67 ln_xios_read = . FALSE. ! use XIOS to read restart file (only for a single file restart)67 ln_xios_read = .false. ! use XIOS to read restart file (only for a single file restart) 68 68 nn_wxios = 0 ! use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 69 69 / … … 88 88 cn_domcfg = "domain_cfg" ! domain configuration filename 89 89 ! 90 ln_closea = .false. ! T => keep closed seas (defined by closea_mask field) in the 90 ln_closea = .false. ! T => keep closed seas (defined by closea_mask field) in the 91 91 ! ! domain and apply special treatment of freshwater fluxes. 92 ! ! F => suppress closed seas (defined by closea_mask field) 92 ! ! F => suppress closed seas (defined by closea_mask field) 93 93 ! ! from the bathymetry at runtime. 94 94 ! ! If closea_mask field doesn't exist in the domain_cfg file … … 106 106 ln_tsd_init = .false. ! ocean initialisation 107 107 ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp) 108 108 109 109 cn_dir = './' ! root directory for the T-S data location 110 110 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! … … 195 195 nn_fsbc = 2 ! frequency of SBC module call 196 196 ! ! (control sea-ice & iceberg model call) 197 ! Type of air-sea fluxes 197 ! Type of air-sea fluxes 198 198 ln_usr = .false. ! user defined formulation (T => check usrdef_sbc) 199 199 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 200 200 ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk ) 201 ln_abl = .false. ! ABL formulation (T => fill namsbc_abl ) 201 202 ! ! Type of coupling (Ocean/Ice/Atmosphere) : 202 203 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) … … 205 206 ! ! =0 no opa-sas OASIS coupling: default single executable config. 206 207 ! ! =1 opa-sas OASIS coupling: multi executable config., OPA component 207 ! ! =2 opa-sas OASIS coupling: multi executable config., SAS component 208 ! ! =2 opa-sas OASIS coupling: multi executable config., SAS component 208 209 ! Sea-ice : 209 nn_ice = 0 ! =0 no ice boundary condition 210 nn_ice = 0 ! =0 no ice boundary condition 210 211 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 211 212 ! ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") … … 213 214 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 214 215 ! ! =F levitating ice (no pressure, mass and salt exchanges) 215 ! Misc. options of sbc : 216 ! Misc. options of sbc : 216 217 ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr) 217 218 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave … … 225 226 ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave) 226 227 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 227 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 228 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 228 229 nn_sdrift = 0 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 229 230 ! ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] … … 250 251 / 251 252 !----------------------------------------------------------------------- 252 &namsbc_blk ! namsbc_blk generic Bulk formula 253 &namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T) 253 254 !----------------------------------------------------------------------- 254 255 ! ! bulk algorithm : 255 ln_NCAR = .false.! "NCAR" algorithm (Large and Yeager 2008)256 ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) 256 257 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 257 ln_COARE_3p 5 = .false. ! "COARE 3.5" algorithm (Edson et al. 2013)258 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31)258 ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013) 259 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 45r1) 259 260 ! 260 rn_zqt = 10. ! Air temperature & humidity reference height (m) 261 rn_zu = 10. ! Wind vector reference height (m) 262 ln_Cd_L12 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2012) 263 ln_Cd_L15 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2015) 264 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 265 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 266 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) 267 rn_vfac = 0. ! multiplicative factor for ocean & ice velocity used to 268 ! ! calculate the wind stress (0.=absolute or 1.=relative winds) 269 261 rn_zqt = 10. ! Air temperature & humidity reference height (m) 262 rn_zu = 10. ! Wind vector reference height (m) 263 ln_Cd_L12 = .false. ! air-ice drags = F(ice conc.) (Lupkes et al. 2012) 264 ln_Cd_L15 = .false. ! air-ice drags = F(ice conc.) (Lupkes et al. 2015) 265 ! ! - module of the mean stress" data 266 rn_pfac = 1. ! multipl. factor for precipitation (total & snow) 267 rn_efac = 1. ! multipl. factor for evaporation (0. or 1.) 268 rn_vfac = 0. ! multipl. factor for ocean & ice velocity 269 ! ! used to calculate the wind stress 270 ! ! (0. => absolute or 1. => relative winds) 271 ln_skin_cs = .false. ! use the cool-skin parameterization 272 ln_skin_wl = .false. ! use the warm-layer parameterization 273 ! ! ==> only available in ECMWF and COARE algorithms 274 ln_humi_sph = .true. ! humidity "sn_humi" is specific humidity [kg/kg] 275 ln_humi_dpt = .false. ! humidity "sn_humi" is dew-point temperature [K] 276 ln_humi_rlh = .false. ! humidity "sn_humi" is relative humidity [%] 277 ! 270 278 cn_dir = './' ! root directory for the bulk data location 271 279 !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! … … 278 286 sn_tair = 't_10.15JUNE2009_fill' , 6. , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 279 287 sn_humi = 'q_10.15JUNE2009_fill' , 6. , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 288 sn_hpgi = 'NONE' , 24. , 'uhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG' , '' 289 sn_hpgj = 'NONE' , 24. , 'vhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG' , '' 280 290 sn_prec = 'ncar_precip.15JUNE2009_fill', -1. , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 281 291 sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 282 292 sn_slp = 'slp.15JUNE2009_fill' , 6. , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 283 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 293 / 294 !----------------------------------------------------------------------- 295 &namsbc_abl ! Atmospheric Boundary Layer formulation (ln_abl = T) 296 !----------------------------------------------------------------------- 297 cn_dir = './' ! root directory for the location of the ABL grid file 298 cn_dom = 'dom_cfg_abl.nc' 299 300 cn_ablrst_in = "restart_abl" ! suffix of abl restart name (input) 301 cn_ablrst_out = "restart_abl" ! suffix of abl restart name (output) 302 cn_ablrst_indir = "." ! directory to read input abl restarts 303 cn_ablrst_outdir = "." ! directory to write output abl restarts 304 305 ln_hpgls_frc = .false. 306 ln_geos_winds = .false. 307 nn_dyn_restore = 2 ! restoring option for dynamical ABL variables: = 0 no restoring 308 ! = 1 equatorial restoring 309 ! = 2 global restoring 310 rn_ldyn_min = 4.5 ! magnitude of the nudging on ABL dynamics at the bottom of the ABL [hour] 311 rn_ldyn_max = 1.5 ! magnitude of the nudging on ABL dynamics at the top of the ABL [hour] 312 rn_ltra_min = 4.5 ! magnitude of the nudging on ABL tracers at the bottom of the ABL [hour] 313 rn_ltra_max = 1.5 ! magnitude of the nudging on ABL tracers at the top of the ABL [hour] 314 nn_amxl = 0 ! mixing length: = 0 Deardorff 80 length-scale 315 ! = 1 length-scale based on the distance to the PBL height 316 ! = 2 Bougeault & Lacarrere 89 length-scale 317 rn_Cm = 0.0667 ! 0.126 in MesoNH 318 rn_Ct = 0.1667 ! 0.143 in MesoNH 319 rn_Ce = 0.4 ! 0.4 in MesoNH 320 rn_Ceps = 0.7 ! 0.85 in MesoNH 321 rn_Rod = 0.15 ! c0 in RMCA17 mixing length formulation (not yet implemented) 322 rn_Ric = 0.139 ! Critical Richardson number (to compute PBL height and diffusivities) 284 323 / 285 324 !----------------------------------------------------------------------- … … 375 414 nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) 376 415 rn_si1 = 23.0 ! 2BD : longest depth of extinction 377 416 378 417 cn_dir = './' ! root directory for the chlorophyl data location 379 418 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! … … 443 482 / 444 483 !----------------------------------------------------------------------- 445 &namsbc_isf ! Top boundary layer (ISF) (ln_isfcav =T : read (ln_read_cfg=T) 484 &namsbc_isf ! Top boundary layer (ISF) (ln_isfcav =T : read (ln_read_cfg=T) 446 485 !----------------------------------------------------------------------- or set or usr_def_zgr ) 447 ! ! type of top boundary layer 486 ! ! type of top boundary layer 448 487 nn_isf = 1 ! ice shelf melting/freezing 449 ! 1 = presence of ISF ; 2 = bg03 parametrisation 488 ! 1 = presence of ISF ; 2 = bg03 parametrisation 450 489 ! 3 = rnf file for ISF ; 4 = ISF specified freshwater flux 451 490 ! options 1 and 4 need ln_isfcav = .true. (domzgr) … … 470 509 !* nn_isf = 3 case 471 510 sn_rnfisf = 'rnfisf' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' 472 !* nn_isf = 2 and 3 cases 511 !* nn_isf = 2 and 3 cases 473 512 sn_depmax_isf ='rnfisf' , -12. ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 474 513 sn_depmin_isf ='rnfisf' , -12. ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' … … 477 516 / 478 517 !----------------------------------------------------------------------- 479 &namsbc_iscpl ! land ice / ocean coupling option (ln_isfcav =T : read (ln_read_cfg=T) 518 &namsbc_iscpl ! land ice / ocean coupling option (ln_isfcav =T : read (ln_read_cfg=T) 480 519 !----------------------------------------------------------------------- or set or usr_def_zgr ) 481 520 nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells) … … 572 611 !----------------------------------------------------------------------- 573 612 ln_tide = .false. ! Activate tides 574 ln_tide_pot = . true.! use tidal potential forcing613 ln_tide_pot = .false. ! use tidal potential forcing 575 614 ln_scal_load = .false. ! Use scalar approximation for 576 615 rn_scal_load = 0.094 ! load potential 577 616 ln_read_load = .false. ! Or read load potential from file 578 617 cn_tide_load = 'tide_LOAD_grid_T.nc' ! filename for load potential 579 ! 618 ! 580 619 ln_tide_ramp = .false. ! Use linear ramp for tides at startup 581 620 rdttideramp = 0. ! ramp duration in days … … 656 695 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 657 696 ln_bdytide_2ddta = .false. ! 658 ln_bdytide_conj = .false. ! 697 ln_bdytide_conj = .false. ! 659 698 / 660 699 … … 683 722 !----------------------------------------------------------------------- 684 723 rn_Cd0 = 1.e-3 ! drag coefficient [-] 685 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 724 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 686 725 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 687 726 rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) … … 694 733 !----------------------------------------------------------------------- 695 734 rn_Cd0 = 1.e-3 ! drag coefficient [-] 696 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 735 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 697 736 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 698 737 rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) … … 761 800 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 762 801 ln_traadv_fct = .false. ! FCT scheme 763 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 764 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 802 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 803 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 765 804 ln_traadv_mus = .false. ! MUSCL scheme 766 805 ln_mus_ups = .false. ! use upstream scheme near river mouths … … 783 822 ln_traldf_triad = .false. ! iso-neutral (triad operator) 784 823 ! 785 ! ! iso-neutral options: 824 ! ! iso-neutral options: 786 825 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 787 826 rn_slpmax = 0.01 ! slope limit (both operators) … … 793 832 nn_aht_ijk_t = 0 ! space/time variation of eddy coefficient: 794 833 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 795 ! ! = 0 constant 796 ! ! = 10 F(k) =ldf_c1d 797 ! ! = 20 F(i,j) =ldf_c2d 834 ! ! = 0 constant 835 ! ! = 10 F(k) =ldf_c1d 836 ! ! = 20 F(i,j) =ldf_c2d 798 837 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 799 838 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 800 839 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 801 ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) 840 ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) 802 841 ! ! or = 1/12 Ud*Ld^3 (blp case) 803 842 rn_Ud = 0.01 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) … … 825 864 nn_aei_ijk_t = 0 ! space/time variation of eddy coefficient: 826 865 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 827 ! ! = 0 constant 828 ! ! = 10 F(k) =ldf_c1d 829 ! ! = 20 F(i,j) =ldf_c2d 866 ! ! = 0 constant 867 ! ! = 10 F(k) =ldf_c1d 868 ! ! = 20 F(i,j) =ldf_c2d 830 869 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 831 870 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 832 ! ! time invariant coefficients: aei0 = 1/2 Ue*Le 871 ! ! time invariant coefficients: aei0 = 1/2 Ue*Le 833 872 rn_Ue = 0.02 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 834 873 rn_Le = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10) … … 870 909 rn_lf_cutoff = 5.0 ! cutoff frequency for low-pass filter [days] 871 910 rn_zdef_max = 0.9 ! maximum fractional e3t deformation 872 ln_vvl_dbg = . true.! debug prints (T/F)911 ln_vvl_dbg = .false. ! debug prints (T/F) 873 912 / 874 913 !----------------------------------------------------------------------- … … 890 929 ln_dynvor_eeT = .false. ! energy conserving scheme (een using e3t) 891 930 ln_dynvor_een = .false. ! energy & enstrophy scheme 892 nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4 931 nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4 893 932 ! ! =1 e3f = mi(mj(e3t))/mi(mj( tmask)) 894 933 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) ==>>> PLEASE DO NOT ACTIVATE … … 935 974 ! ! =-30 read in eddy_viscosity_3D.nc file 936 975 ! ! =-20 read in eddy_viscosity_2D.nc file 937 ! ! = 0 constant 976 ! ! = 0 constant 938 977 ! ! = 10 F(k)=c1d 939 978 ! ! = 20 F(i,j)=F(grid spacing)=c2d … … 941 980 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 942 981 ! ! = 32 F(i,j,k)=F(local gridscale and deformation rate) 943 ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) 982 ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) 944 983 ! ! or = 1/12 Uv*Lv^3 (blp case) 945 984 rn_Uv = 0.1 ! lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) … … 1065 1104 ! = 0 constant 10 m length scale 1066 1105 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 1067 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1106 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1068 1107 / 1069 1108 !----------------------------------------------------------------------- … … 1323 1362 ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 1324 1363 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 1325 sn_cfctl%l_runstat = . FALSE.! switches and which areas produce reports with the proc integer settings.1364 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1326 1365 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1327 1366 sn_cfctl%l_oceout = .FALSE. ! that all areas report.
Note: See TracChangeset
for help on using the changeset viewer.