New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 5980 for branches/2014/dev_r4650_UKMO7_STARTHOUR/NEMOGCM/CONFIG/SHARED/namelist_ref – NEMO

Ignore:
Timestamp:
2015-12-02T16:20:47+01:00 (8 years ago)
Author:
timgraham
Message:

Upgraded to v3.6 revision of trunk (r5518)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4650_UKMO7_STARTHOUR/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5872 r5980  
    1010!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    1111!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
    12 !!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namptr, namhsb) 
     12!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    1313!!             10 - miscellaneous    (namsol, nammpp, namctl) 
    1414!!             11 - Obs & Assim      (namobs, nam_asminc) 
     
    3838                           !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
    3939   cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input) 
     40   cn_ocerst_indir = "."       !  directory from which to read input ocean restarts 
    4041   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
     42   cn_ocerst_outdir = "."      !  directory in which to write output ocean restarts 
    4143   nn_istate   =       0   !  output the initial state (1) or not (0) 
     44   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
    4245   nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1) 
     46   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 
    4347   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    4448   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4549   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
     50   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard 
    4651   ln_clobber  = .false.   !  clobber (overwrite) an existing file 
    4752   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
     
    5156!!                      ***  Domain namelists  *** 
    5257!!====================================================================== 
    53 !!   namcfg       parameters of the configuration       
     58!!   namcfg       parameters of the configuration 
    5459!!   namzgr       vertical coordinate 
    5560!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
     
    5964! 
    6065!----------------------------------------------------------------------- 
    61 &namcfg     !   parameters of the configuration       
     66&namcfg     !   parameters of the configuration 
    6267!----------------------------------------------------------------------- 
    6368   cp_cfg      =  "default"            !  name of the configuration 
     
    7378   jperio      =       0               !  lateral cond. type (between 0 and 6) 
    7479                                       !  = 0 closed                 ;   = 1 cyclic East-West 
    75                                        !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot  
     80                                       !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
    7681                                       !  = 4 cyclic East-West AND North fold T-point pivot 
    7782                                       !  = 5 North fold F-point pivot 
    7883                                       !  = 6 cyclic East-West AND North fold F-point pivot 
     84   ln_use_jattr = .false.              !  use (T) the file attribute: open_ocean_jstart, if present 
     85                                       !  in netcdf input files, as the start j-row for reading 
    7986/ 
    8087!----------------------------------------------------------------------- 
     
    8491   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    8592   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
    86    ln_isfcav   = .false.   !  ice shelf cavity 
     93   ln_isfcav   = .false.   !  ice shelf cavity               (T/F) 
    8794/ 
    8895!----------------------------------------------------------------------- 
     
    100107                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
    101108   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
    102    rn_bb       =    0.8    !  stretching with SH94 s-sigma    
     109   rn_bb       =    0.8    !  stretching with SH94 s-sigma 
    103110                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
    104111   rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
     
    109116   rn_zb_b     =   -0.2    !  offset for calculating Zb 
    110117                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
    111    rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)  
     118   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1) 
    112119/ 
    113120!----------------------------------------------------------------------- 
     
    117124   rn_bathy    =    0.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
    118125   nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    119    nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
     126   nn_msh      =    1      !  create (=1) a mesh file or not (=0) 
    120127   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
    121128   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
     
    163170   nn_baro       =    30               !  Number of iterations of barotropic mode 
    164171                                       !  during rn_rdt seconds. Only used if ln_bt_nn_auto=F 
    165    rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T  
     172   rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T 
    166173   nn_bt_flt     =    1                !  Time filter choice 
    167174                                       !  = 0 None 
    168175                                       !  = 1 Boxcar over   nn_baro barotropic steps 
    169                                        !  = 2 Boxcar over 2*nn_baro     "        "   
     176                                       !  = 2 Boxcar over 2*nn_baro     "        " 
    170177/ 
    171178!----------------------------------------------------------------------- 
     
    214221!!   namsbc_core     CORE bulk formulae formulation 
    215222!!   namsbc_mfs      MFS  bulk formulae formulation 
    216 !!   namsbc_cpl      CouPLed            formulation                     ("key_coupled") 
     223!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3") 
    217224!!   namsbc_sas      StAndalone Surface module 
    218225!!   namtra_qsr      penetrative solar radiation 
     
    234241   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    235242   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
     243   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
     244   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 ) 
     245   nn_components = 0       !  configuration of the opa-sas OASIS coupling 
     246                           !  =0 no opa-sas OASIS coupling: default single executable configuration 
     247                           !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component 
     248                           !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component  
    236249   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    237250   nn_ice      = 2         !  =0 no ice boundary condition   , 
     
    244257   ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf) 
    245258   nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf) 
    246                            !  0 =no isf                  1 = presence of ISF  
    247                            !  2 = bg03 parametrisation   3 = rnf file for isf    
     259                           !  0 =no isf                  1 = presence of ISF 
     260                           !  2 = bg03 parametrisation   3 = rnf file for isf 
    248261                           !  4 = ISF fwf specified 
    249262                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
     
    276289&namsbc_flx    !   surface boundary condition : flux formulation 
    277290!----------------------------------------------------------------------- 
    278 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !  
     291!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    279292!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    280293   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
     
    319332   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    320333   rn_zqt      = 10.        !  Air temperature and humidity reference height (m) 
    321    rn_zu       = 10.        !  Wind vector reference height (m)                  
     334   rn_zu       = 10.        !  Wind vector reference height (m) 
    322335   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    323336   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    324    rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity  
     337   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity 
    325338                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    326339/ 
     
    341354/ 
    342355!----------------------------------------------------------------------- 
    343 &namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_coupled") 
     356&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
    344357!----------------------------------------------------------------------- 
    345358!                    !     description       !  multiple  !    vector   !      vector          ! vector ! 
     
    372385!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    373386!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    374    sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''       
     387   sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    375388   sn_vsp      = 'sas_grid_V' ,    120           , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    376389   sn_tem      = 'sas_grid_T' ,    120           , 'sosstsst' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    377390   sn_sal      = 'sas_grid_T' ,    120           , 'sosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    378391   sn_ssh      = 'sas_grid_T' ,    120           , 'sossheig' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    379  
    380    ln_3d_uv    = .true.    !  specify whether we are supplying a 3D u,v field 
     392   sn_e3t      = 'sas_grid_T' ,    120           , 'e3t_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     393   sn_frq      = 'sas_grid_T' ,    120           , 'frq_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     394 
     395   ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field 
     396   ln_read_frq = .false.    !  specify whether we must read frq or not 
    381397   cn_dir      = './'      !  root directory for the location of the bulk files are 
    382398/ 
     
    411427 
    412428   cn_dir       = './'      !  root directory for the location of the runoff files 
    413    ln_rnf_emp   = .false.   !  runoffs included into precipitation field (T) or into a file (F) 
    414429   ln_rnf_mouth = .true.    !  specific treatment at rivers mouths 
    415430   rn_hrnf      =  15.e0    !  depth over which enhanced vertical mixing is used 
     
    419434   ln_rnf_tem   = .false.   !  read in temperature information for runoff 
    420435   ln_rnf_sal   = .false.   !  read in salinity information for runoff 
    421 / 
    422 !----------------------------------------------------------------------- 
    423 &namsbc_isf    !  Top boundary layer (ISF)  
     436   ln_rnf_depth_ini = .false.  ! compute depth at initialisation from runoff file 
     437   rn_rnf_max   = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 
     438   rn_dep_max   = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
     439   nn_rnf_depth_file = 0    !  create (=1) a runoff depth file or not (=0) 
     440/ 
     441!----------------------------------------------------------------------- 
     442&namsbc_isf    !  Top boundary layer (ISF) 
    424443!----------------------------------------------------------------------- 
    425444!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
     
    458477 
    459478   cn_dir      = './'       !  root directory for the location of the bulk files 
    460    rn_pref     = 101000._wp !  reference atmospheric pressure   [N/m2]/ 
     479   rn_pref     = 101000.    !  reference atmospheric pressure   [N/m2]/ 
    461480   ln_ref_apr  = .false.    !  ref. pressure: global mean Patm (T) or a constant (F) 
    462481   ln_apr_obc  = .false.    !  inverse barometer added to OBC ssh data 
     
    498517                                                      ! Initial mass required for an iceberg of each class 
    499518      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
    500                                                       ! Proportion of calving mass to apportion to each class   
     519                                                      ! Proportion of calving mass to apportion to each class 
    501520      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    502521                                                      ! Ratio between effective and real iceberg mass (non-dim) 
    503                                                       ! i.e. number of icebergs represented at a point          
     522                                                      ! i.e. number of icebergs represented at a point 
    504523      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    505524                                                      ! thickness of newly calved bergs (m) 
     
    510529      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
    511530      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
    512       ln_passive_mode          = .false.              ! iceberg - ocean decoupling    
     531      ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
    513532      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
    514533                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
    515534      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
    516       rn_speed_limit           = 0.                   ! CFL speed limit for a berg    
     535      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
    517536 
    518537!              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    519538!              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    520539      sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , '' 
    521     
    522       cn_dir = './'  
     540 
     541      cn_dir = './' 
    523542/ 
    524543 
     
    581600   ln_tide_ramp  = .false.  ! 
    582601   rdttideramp   =    0.    ! 
    583    clname(1)     =   'M2'   !  name of constituent 
    584    clname(2)     =   'S2' 
    585    clname(3)     =   'N2' 
    586    clname(4)     =   'K1' 
    587    clname(5)     =   'O1' 
    588    clname(6)     =   'Q1' 
    589    clname(7)     =   'M4' 
    590    clname(8)     =   'K2' 
    591    clname(9)     =   'P1' 
    592    clname(10)    =   'Mf' 
    593    clname(11)    =   'Mm' 
     602   clname(1)     = 'DUMMY'  !  name of constituent - all tidal components must be set in namelist_cfg 
    594603/ 
    595604!----------------------------------------------------------------------- 
     
    606615                                          !  = 2, use tidal harmonic forcing data from files 
    607616                                          !  = 3, use external data AND tidal harmonic forcing 
    608     cn_dyn3d      =  'none'               !   
     617    cn_dyn3d      =  'none'               ! 
    609618    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state 
    610619                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    611     cn_tra        =  'none'               !  
     620    cn_tra        =  'none'               ! 
    612621    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state 
    613622                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    614     cn_ice_lim      =  'none'             !   
     623    cn_ice_lim      =  'none'             ! 
    615624    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state 
    616625                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     
    621630    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers 
    622631    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities 
    623     rn_time_dmp   =  1.                   ! Damping time scale in days  
     632    rn_time_dmp   =  1.                   ! Damping time scale in days 
    624633    rn_time_dmp_out =  1.                 ! Outflow damping time scale 
    625634    nn_rimwidth   = 10                    !  width of the relaxation zone 
     
    674683   rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T) 
    675684   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    676    rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T  
     685   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T 
    677686   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    678687   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
     
    691700&nambbc        !   bottom temperature boundary condition 
    692701!----------------------------------------------------------------------- 
     702!              !                              !  (if <0  months)  !   
     703!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     704!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     705   sn_qgh      ='geothermal_heating.nc',  -12.  , 'heatflow'      ,   .false.      , .true.  , 'yearly'  , ''       , ''       , '' 
     706   ! 
     707   cn_dir      = './'      !  root directory for the location of the runoff files 
    693708   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    694709   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux 
     
    696711                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    697712   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
     713 
    698714/ 
    699715!----------------------------------------------------------------------- 
     
    720736!----------------------------------------------------------------------- 
    721737   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency 
    722                                  !  =-1, TEOS-10  
    723                                  !  = 0, EOS-80  
     738                                 !  =-1, TEOS-10 
     739                                 !  = 0, EOS-80 
    724740                                 !  = 1, S-EOS   (simplified eos) 
    725741   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     
    792808!----------------------------------------------------------------------- 
    793809   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F) 
    794    nn_hdmp     =   -1      !  horizontal shape =-1, damping in Med and Red Seas only 
    795                            !                   =XX, damping poleward of XX degrees (XX>0) 
    796                            !                      + F(distance-to-coast) + Red and Med Seas 
    797810   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column 
    798811                           !                   =1 no damping in the mixing layer (kz  criteria) 
    799812                           !                   =2 no damping in the mixed  layer (rho crieria) 
    800    rn_surf     =   50.     !  surface time scale of damping   [days] 
    801    rn_bot      =  360.     !  bottom  time scale of damping   [days] 
    802    rn_dep      =  800.     !  depth of transition between rn_surf and rn_bot [meters] 
    803    nn_file     =    0      !  create a damping.coeff NetCDF file (=1) or not (=0) 
     813   cn_resto    = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this) 
    804814/ 
    805815 
     
    818828!----------------------------------------------------------------------- 
    819829   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
     830   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction 
    820831   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    821832   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
     
    825836&nam_vvl    !   vertical coordinate options 
    826837!----------------------------------------------------------------------- 
    827    ln_vvl_zstar  = .true.           !  zstar vertical coordinate                    
     838   ln_vvl_zstar  = .true.           !  zstar vertical coordinate 
    828839   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations 
    829840   ln_vvl_layer  = .false.          !  full layer vertical coordinate 
     
    851862   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
    852863   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
     864   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf 
    853865   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    854866   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
     
    946958                           !        = 1 add a tke source below the ML 
    947959                           !        = 2 add a tke source just at the base of the ML 
    948                            !        = 3 as = 1 applied on HF part of the stress    ("key_coupled") 
     960                           !        = 3 as = 1 applied on HF part of the stress    ("key_oasis3") 
    949961   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 
    950962   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML 
     
    968980&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    969981!----------------------------------------------------------------------- 
    970    rn_emin       = 1.e-6   !  minimum value of e   [m2/s2] 
     982   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2] 
    971983   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3] 
    972984   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988) 
    973    rn_clim_galp  = 0.53    !  galperin limit 
    974    ln_crban      = .true.  !  Use Craig & Banner (1994) surface wave mixing parametrisation 
     985   rn_clim_galp  = 0.267   !  galperin limit 
    975986   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case 
    976987   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux 
    977988   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length 
    978    nn_tkebc_surf =   1     !  surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 
    979    nn_tkebc_bot  =   1     !  bottom tke condition (0/1=Dir/Neum) 
    980    nn_psibc_surf =   1     !  surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg) 
    981    nn_psibc_bot  =   1     !  bottom psi condition (0/1=Dir/Neum) 
    982    nn_stab_func  =   2     !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 
    983    nn_clos       =   1     !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 
     989   rn_hsro       =  0.02   !  Minimum surface roughness 
     990   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2) 
     991   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2) 
     992   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum) 
     993   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum) 
     994   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 
     995   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 
    984996/ 
    985997!----------------------------------------------------------------------- 
     
    10091021!!   namc1d_uvd        data: U & V currents                             ("key_c1d") 
    10101022!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d") 
     1023!!   namsto            Stochastic parametrization of EOS 
    10111024!!====================================================================== 
    10121025! 
     
    10671080   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    10681081/ 
     1082!----------------------------------------------------------------------- 
     1083&namsto       ! Stochastic parametrization of EOS 
     1084!----------------------------------------------------------------------- 
     1085   ln_rststo = .false.           ! start from mean parameter (F) or from restart file (T) 
     1086   ln_rstseed = .true.           ! read seed of RNG from restart file 
     1087   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input) 
     1088   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output) 
     1089 
     1090   ln_sto_eos = .false.          ! stochastic equation of state 
     1091   nn_sto_eos = 1                ! number of independent random walks 
     1092   rn_eos_stdxy = 1.4            ! random walk horz. standard deviation (in grid points) 
     1093   rn_eos_stdz  = 0.7            ! random walk vert. standard deviation (in grid points) 
     1094   rn_eos_tcor  = 1440.0         ! random walk time correlation (in timesteps) 
     1095   nn_eos_ord  = 1               ! order of autoregressive processes 
     1096   nn_eos_flt  = 0               ! passes of Laplacian filter 
     1097   rn_eos_lim  = 2.0             ! limitation factor (default = 3.0) 
     1098/ 
    10691099 
    10701100!!====================================================================== 
     
    10731103!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    10741104!!   namtrd       dynamics and/or tracer trends 
     1105!!   namptr       Poleward Transport Diagnostics 
    10751106!!   namflo       float parameters                                      ("key_float") 
    1076 !!   namptr       Poleward Transport Diagnostics 
    10771107!!   namhsb       Heat and salt budgets 
    10781108!!====================================================================== 
     
    11281158!----------------------------------------------------------------------- 
    11291159   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
    1130    ln_diaznl  = .true.     !  Add zonal means and meridional stream functions 
    1131    ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    1132                            !  (orca configuration only, need input basins mask file named "subbasins.nc" 
    1133    ln_ptrcomp = .true.     !  Add decomposition : overturning 
    1134    nn_fptr    =  1         !  Frequency of ptr computation [time step] 
    1135    nn_fwri    =  15        !  Frequency of ptr outputs [time step] 
     1160   ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    11361161/ 
    11371162!----------------------------------------------------------------------- 
     
    11831208   ln_sst     = .false.     ! Logical switch for SST observations 
    11841209   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations 
    1185    ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations       
     1210   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations 
    11861211 
    11871212   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
     
    12101235   sstfbfiles = 'sst_01.nc' 
    12111236                           !     seaicefiles             Sea Ice input observation file names 
    1212    seaicefiles = 'seaice_01.nc'   
     1237   seaicefiles = 'seaice_01.nc' 
    12131238                           !     velavcurfiles           Vel. cur. daily av. input file name 
    12141239                           !     velhvcurfiles           Vel. cur. high freq. input file name 
Note: See TracChangeset for help on using the changeset viewer.