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 6497 for trunk/NEMOGCM/CONFIG – NEMO

Changeset 6497 for trunk/NEMOGCM/CONFIG


Ignore:
Timestamp:
2016-04-27T09:33:46+02:00 (8 years ago)
Author:
gm
Message:

#1720 - trunk: add Casimir tidal parameterization

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMOGCM/CONFIG/SHARED/namelist_ref

    r6489 r6497  
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!! NEMO/OPA  :  1 - run manager      (namrun) 
    5 !! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd) 
     5!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
    66!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    77!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    … …  
    5959!!====================================================================== 
    6060!!   namcfg       parameters of the configuration 
    61 !!   namzgr       vertical coordinate 
     61!!   namzgr       vertical coordinate                                   (default: NO selection) 
    6262!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    6363!!   namdom       space and time domain (bathymetry, mesh, timestep) 
     64!!   namwad       Wetting and drying                                    (default F) 
     65!!   namtsd       data: temperature & salinity 
    6466!!   namcrs       coarsened grid (for outputs and/or TOP)               ("key_crs") 
    6567!!   namc1d       1D configuration options                              ("key_c1d") 
     68!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
    6669!!   namc1d_uvd   1D data (currents)                                    ("key_c1d") 
    67 !!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
    68 !!   namtsd       data: temperature & salinity 
    6970!!====================================================================== 
    7071! 
    … …  
    101102/ 
    102103!----------------------------------------------------------------------- 
    103 &namzgr_sco    !   s-coordinate or hybrid z-s-coordinate 
     104&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default F) 
    104105!----------------------------------------------------------------------- 
    105106   ln_s_sh94   = .false.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
    … …  
    165166/ 
    166167!----------------------------------------------------------------------- 
     168&namwad        !   Wetting and drying                                   (default F) 
     169!----------------------------------------------------------------------- 
     170   ln_wd       = .false.   !  T/F activation of wetting and drying 
     171   rn_wdmin1   =  0.1      !  Minimum wet depth on dried cells 
     172   rn_wdmin2   =  0.01     !  Tolerance of min wet depth on dried cells 
     173   rn_wdld     =  20.0     !  Land elevation below which wetting/drying is allowed 
     174   nn_wdit     =  10       !  Max iterations for W/D limiter 
     175/ 
     176!----------------------------------------------------------------------- 
     177&namtsd        !   data : Temperature  & Salinity 
     178!----------------------------------------------------------------------- 
     179!              !  file name                 ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     180!              !                            !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     181   sn_tem = 'data_1m_potential_temperature_nomask',     -1      ,'votemper',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     182   sn_sal = 'data_1m_salinity_nomask'             ,     -1      ,'vosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     183   ! 
     184   cn_dir      = './'      !  root directory for the location of the runoff files 
     185   ln_tsd_init = .true.    !  Initialisation of ocean T & S with T & S input data (T) or not (F) 
     186   ln_tsd_tradmp = .true.  !  damping of ocean T & S toward T & S input data (T) or not (F) 
     187/ 
     188!----------------------------------------------------------------------- 
    167189&namcrs        !   coarsened grid (for outputs and/or TOP)              ("key_crs") 
    168190!----------------------------------------------------------------------- 
    … …  
    202224   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F) 
    203225   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F) 
    204 / 
    205 !----------------------------------------------------------------------- 
    206 &namtsd    !   data : Temperature  & Salinity 
    207 !----------------------------------------------------------------------- 
    208 !          !  file name                            ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    209 !          !                                       !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    210    sn_tem  = 'data_1m_potential_temperature_nomask',         -1        ,'votemper' ,    .true.    , .true. , 'yearly'   , ''       ,   ''    ,    '' 
    211    sn_sal  = 'data_1m_salinity_nomask'             ,         -1        ,'vosaline' ,    .true.    , .true. , 'yearly'   , ''       ,   ''    ,    '' 
    212    ! 
    213    cn_dir        = './'     !  root directory for the location of the runoff files 
    214    ln_tsd_init   = .true.   !  Initialisation of ocean T & S with T &S input data (T) or not (F) 
    215    ln_tsd_tradmp = .true.   !  damping of ocean T & S toward T &S input data (T) or not (F) 
    216226/ 
    217227 
    … …  
    278288   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    279289   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf) 
    280    ln_wave = .false.       !  coupling with surface wave                (T => fill namsbc_wave) 
    281    nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
     290   ln_wave     = .false.   !  coupling with surface wave                (T => fill namsbc_wave) 
     291   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
    282292                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
    283293/ 
    … …  
    382392   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    383393! 
    384    nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
    385    ln_usecplmask = .false.   !  use a coupling mask file to merge data received from several models 
    386                              !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
     394   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
     395   ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models 
     396   !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    387397/ 
    388398!----------------------------------------------------------------------- 
    389399&namsbc_sas    !   analytical surface boundary condition 
    390400!----------------------------------------------------------------------- 
    391 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    392 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    393    sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    394    sn_vsp      = 'sas_grid_V' ,    120           , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    395    sn_tem      = 'sas_grid_T' ,    120           , 'sosstsst' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    396    sn_sal      = 'sas_grid_T' ,    120           , 'sosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    397    sn_ssh      = 'sas_grid_T' ,    120           , 'sossheig' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    398    sn_e3t      = 'sas_grid_T' ,    120           , 'e3t_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
    399    sn_frq      = 'sas_grid_T' ,    120           , 'frq_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     401!              !  file name  ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     402!              !             !  (if <0  months)  !   name    !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     403   sn_usp      = 'sas_grid_U',     120           , 'vozocrtx',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     404   sn_vsp      = 'sas_grid_V',     120           , 'vomecrty',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     405   sn_tem      = 'sas_grid_T',     120           , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     406   sn_sal      = 'sas_grid_T',     120           , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     407   sn_ssh      = 'sas_grid_T',     120           , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     408   sn_e3t      = 'sas_grid_T',     120           , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     409   sn_frq      = 'sas_grid_T',     120           , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    400410 
    401411   ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field 
    402    ln_read_frq = .false.    !  specify whether we must read frq or not 
     412   ln_read_frq = .false.   !  specify whether we must read frq or not 
    403413   cn_dir      = './'      !  root directory for the location of the bulk files are 
    404414/ 
    … …  
    431441   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    432442 
    433    cn_dir       = './'      !  root directory for the location of the runoff files 
    434    ln_rnf_mouth = .true.    !  specific treatment at rivers mouths 
    435    rn_hrnf      =  15.e0    !  depth over which enhanced vertical mixing is used 
    436    rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] 
    437    rn_rfact     =   1.e0    !  multiplicative factor for runoff 
    438    ln_rnf_depth = .false.   !  read in depth information for runoff 
    439    ln_rnf_tem   = .false.   !  read in temperature information for runoff 
    440    ln_rnf_sal   = .false.   !  read in salinity information for runoff 
    441    ln_rnf_depth_ini = .false.  ! compute depth at initialisation from runoff file 
    442    rn_rnf_max   = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 
    443    rn_dep_max   = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
    444    nn_rnf_depth_file = 0    !  create (=1) a runoff depth file or not (=0) 
     443   cn_dir      = './'      !  root directory for the location of the runoff files 
     444   ln_rnf_mouth= .true.    !  specific treatment at rivers mouths 
     445      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T) 
     446      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 
     447   rn_rfact    =   1.e0    !  multiplicative factor for runoff 
     448   ln_rnf_depth= .false.   !  read in depth information for runoff 
     449   ln_rnf_tem  = .false.   !  read in temperature information for runoff 
     450   ln_rnf_sal  = .false.   !  read in salinity information for runoff 
     451   ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file 
     452      rn_rnf_max  = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 
     453      rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
     454      nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0) 
    445455/ 
    446456!----------------------------------------------------------------------- 
    447457&namsbc_isf    !  Top boundary layer (ISF)                              (nn_isf >0) 
    448458!----------------------------------------------------------------------- 
    449 !              ! file name ! frequency (hours) ! variable ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    450 !              !           !  (if <0  months)  !   name   !  (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     459!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     460!              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    451461! nn_isf == 4 
    452    sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
     462   sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    453463! nn_isf == 3 
    454    sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
     464   sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    455465! nn_isf == 2 and 3 
    456    sn_depmax_isf='rnfisf'  ,         -12       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    457    sn_depmin_isf='rnfisf'  ,         -12       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
     466   sn_depmax_isf='rnfisf'  ,         -12       ,'sozisfmax',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
     467   sn_depmin_isf='rnfisf'  ,         -12       ,'sozisfmin',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    458468! nn_isf == 2 
    459    sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
     469   sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    460470! 
    461471! for all case 
    … …  
    465475                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    466476! only for nn_isf = 1 or 2 
    467    rn_gammat0  = 1.e-4    ! gammat coefficient used in blk formula 
    468    rn_gammas0  = 1.e-4    ! gammas coefficient used in blk formula 
     477   rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
     478   rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    469479! only for nn_isf = 1 or 4 
    470480   rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    471                           ! 0 => thickness of the tbl = thickness of the first wet cell 
     481   !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    472482! only for nn_isf = 1 
    473    nn_isfblk   = 1        ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    474                           ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    475    nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s) 
    476                           ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    477                           ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
     483   nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
     484   !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
     485   nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
     486   !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     487   !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    478488/ 
    479489!----------------------------------------------------------------------- 
    480490&namsbc_iscpl  !   land ice / ocean coupling option                      
    481491!----------------------------------------------------------------------- 
    482    nn_drown  = 10       ! number of iteration of the extrapolation loop (fill the new wet cells) 
    483    ln_hsb    = .false.  ! activate conservation module (conservation exact after a time of rn_fiscpl) 
    484    nn_fiscpl = 43800    ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
     492   nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     493   ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
     494   nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
    485495/ 
    486496!----------------------------------------------------------------------- 
    487497&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk 
    488498!----------------------------------------------------------------------- 
    489 !              !  file name ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    490 !              !            !  (if <0  months)  !   name    !  (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    491    sn_apr      = 'patm'     ,         -1        ,'somslpre',    .true.     , .true. , 'yearly'  ,  ''      ,   ''     , '' 
    492  
    493    cn_dir      = './'       !  root directory for the location of the bulk files 
    494    rn_pref     = 101000.    !  reference atmospheric pressure   [N/m2]/ 
    495    ln_ref_apr  = .false.    !  ref. pressure: global mean Patm (T) or a constant (F) 
    496    ln_apr_obc  = .false.    !  inverse barometer added to OBC ssh data 
     499!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     500!              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     501   sn_apr      = 'patm'    ,         -1        ,'somslpre',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,      '' 
     502 
     503   cn_dir      = './'      !  root directory for the location of the bulk files 
     504   rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/ 
     505   ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F) 
     506   ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data 
    497507/ 
    498508!----------------------------------------------------------------------- 
    499509&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr=T) 
    500510!----------------------------------------------------------------------- 
    501 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    502 !              !             !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    503    sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , ''       , '' 
    504    sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , ''       , '' 
     511!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     512!              !           !  (if <0  months)  !   name   !   (logical) !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     513   sn_sst      = 'sst_data',        24         ,  'sst'   ,    .false.  , .false., 'yearly'  ,    ''    ,    ''    ,    '' 
     514   sn_sss      = 'sss_data',        -1         ,  'sss'   ,    .true.   , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    505515 
    506516   cn_dir      = './'      !  root directory for the location of the runoff files 
    … …  
    510520   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K] 
    511521   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
    512    ln_sssr_bnd =   .true.  !  flag to bound erp term (associated with nn_sssr=2) 
     522   ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    513523   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
    514524/ 
    … …  
    516526&namsbc_alb    !   albedo parameters 
    517527!----------------------------------------------------------------------- 
    518    nn_ice_alb  =    0   !  parameterization of ice/snow albedo 
    519                         !     0: Shine & Henderson-Sellers (JGR 1985) 
    520                         !     1: "home made" based on Brandt et al. (J. Climate 2005) 
    521                         !                         and Grenfell & Perovich (JGR 2004) 
    522    rn_albice   =  0.53  !  albedo of bare puddled ice (values from 0.49 to 0.58) 
    523                         !     0.53 (default) => if nn_ice_alb=0 
    524                         !     0.50 (default) => if nn_ice_alb=1 
     528   nn_ice_alb  =    0      !  parameterization of ice/snow albedo 
     529                           !     0: Shine & Henderson-Sellers (JGR 1985) 
     530                           !     1: "home made" based on Brandt et al. (J. Climate 2005) 
     531                           !                         and Grenfell & Perovich (JGR 2004) 
     532   rn_albice   =  0.53     !  albedo of bare puddled ice (values from 0.49 to 0.58) 
     533                           !     0.53 (default) => if nn_ice_alb=0 
     534                           !     0.50 (default) => if nn_ice_alb=1 
    525535/ 
    526536!----------------------------------------------------------------------- 
    527537&namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
    528538!----------------------------------------------------------------------- 
    529 !              !  file name  ! frequency (hours) ! variable    ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    530 !              !             !  (if <0  months)  !   name      !  (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    531    sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff',     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
    532    sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'    ,     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
    533    sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'    ,     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
    534    sn_wn       =  'sdw_wave' ,        1          , 'wave_num'  ,     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
     539!              ! file name ! frequency (hours) ! variable    ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     540!              !           !  (if <0  months)  !   name      !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     541   sn_cdg      = 'cdg_wave',        1          , 'drag_coeff',   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
     542   sn_usd      = 'sdw_wave',        1          , 'u_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
     543   sn_vsd      = 'sdw_wave',        1          , 'v_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
     544   sn_wn       = 'sdw_wave',        1          , 'wave_num'  ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    535545! 
    536546   cn_dir_cdg  = './'      !  root directory for the location of drag coefficient files 
    537    ln_cdgw = .false.       !  Neutral drag coefficient read from wave model 
    538    ln_sdw  = .false.       !  Computation of 3D stokes drift                
     547   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model 
     548   ln_sdw      = .false.   !  Computation of 3D stokes drift                
    539549/ 
    540550!----------------------------------------------------------------------- 
    … …  
    566576      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
    567577 
    568 !            ! file name ! frequency (hours) !   variable   ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    569 !            !           !  (if <0  months)  !     name     !   (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    570       sn_icb =  'calving',       -1          , 'calvingmask',  .true.      , .true.  , 'yearly'  , ''       , ''       , '' 
     578!            ! file name ! frequency (hours) !   variable   ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     579!            !           !  (if <0  months)  !     name     !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     580      sn_icb =  'calving',       -1          , 'calvingmask',   .true.    , .true.  , 'yearly'  ,    ''    ,    ''    ,    '' 
    571581 
    572582      cn_dir = './' 
    … …  
    577587!!====================================================================== 
    578588!!   namlbc        lateral momentum boundary condition 
    579 !!   namobc        open boundaries parameters                           ("key_obc") 
    580589!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
     590!!   nam_tide      Tidal forcing  
    581591!!   nambdy        Unstructured open boundaries                         ("key_bdy") 
    582 !!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides") 
     592!!   nambdy_dta    Unstructured open boundaries - external data         ("key_bdy") 
     593!!   nambdy_tide   tidal forcing at open boundaries                     ("key_bdy_tides") 
    583594!!====================================================================== 
    584595! 
    … …  
    602613&nam_tide      !   tide parameters                                      ("key_tide") 
    603614!----------------------------------------------------------------------- 
    604    ln_tide_pot   = .true.   !  use tidal potential forcing 
    605    ln_tide_ramp  = .false.  ! 
    606    rdttideramp   =    0.    ! 
    607    clname(1)     = 'DUMMY'  !  name of constituent - all tidal components must be set in namelist_cfg 
     615   ln_tide_pot = .true.    !  use tidal potential forcing 
     616   ln_tide_ramp= .false.   ! 
     617   rdttideramp =    0.     ! 
     618   clname(1)   = 'DUMMY'   !  name of constituent - all tidal components must be set in namelist_cfg 
    608619/ 
    609620!----------------------------------------------------------------------- 
    … …  
    644655&nambdy_dta    !  open boundaries - external data                       ("key_bdy") 
    645656!----------------------------------------------------------------------- 
    646 !              !  file name      ! frequency (hours) ! variable  ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    647 !              !                 !  (if <0  months)  !   name    !  (logical)   !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    648    bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    649    bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    650    bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    651    bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    652    bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    653    bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    654    bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     657!              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     658!              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     659   bn_ssh      = 'amm12_bdyT_u2d',         24        , 'sossheig',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     660   bn_u2d      = 'amm12_bdyU_u2d',         24        , 'vobtcrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     661   bn_v2d      = 'amm12_bdyV_u2d',         24        , 'vobtcrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     662   bn_u3d      = 'amm12_bdyU_u3d',         24        , 'vozocrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     663   bn_v3d      = 'amm12_bdyV_u3d',         24        , 'vomecrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     664   bn_tem      = 'amm12_bdyT_tra',         24        , 'votemper',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     665   bn_sal      = 'amm12_bdyT_tra',         24        , 'vosaline',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
    655666! for lim2 
    656 !   bn_frld  =   'amm12_bdyT_ice' ,         24        , 'ileadfra',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    657 !   bn_hicif =   'amm12_bdyT_ice' ,         24        , 'iicethic',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    658 !   bn_hsnif =   'amm12_bdyT_ice' ,         24        , 'isnowthi',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     667!   bn_frld    = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     668!   bn_hicif   = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     669!   bn_hsnif   = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
    659670! for lim3 
    660 !   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    661 !   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    662 !   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    663  
    664    cn_dir      =    'bdydta/'  !  root directory for the location of the bulk files 
    665    ln_full_vel = .false.        !   
    666 / 
    667 !----------------------------------------------------------------------- 
    668 &nambdy_tide     ! tidal forcing at open boundaries 
     671!   bn_a_i     = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     672!   bn_ht_i    = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     673!   bn_ht_s    = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,    '' 
     674 
     675   cn_dir      = 'bdydta/' !  root directory for the location of the bulk files 
     676   ln_full_vel = .false.   !   
     677/ 
     678!----------------------------------------------------------------------- 
     679&nambdy_tide   ! tidal forcing at open boundaries 
    669680!----------------------------------------------------------------------- 
    670681   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
    … …  
    672683   ln_bdytide_conj  = .false.                   !  
    673684/ 
     685 
    674686!!====================================================================== 
    675687!!                 ***  Bottom boundary condition  *** 
    … …  
    687699   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case) 
    688700   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
    689    rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T) 
     701   rn_bfri2_max=    1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T) 
    690702   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    691703   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T 
    … …  
    694706   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case) 
    695707   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
    696    rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T) 
     708   rn_tfri2_max=    1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T) 
    697709   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2) 
    698710   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T 
    699711   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file ) 
    700    rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T) 
     712   rn_tfrien   =   50.     !  local multiplying factor of tfr (ln_tfr2d=T) 
    701713 
    702714   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true) 
    … …  
    720732&nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    721733!----------------------------------------------------------------------- 
    722    nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0) 
    723    nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0) 
    724    rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s] 
    725    rn_gambbl   =  10.    !  advective bbl coefficient                 [s] 
     734   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
     735   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     736   rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
     737   rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
    726738/ 
    727739 
    … …  
    744756   ln_seos     = .false.         !  = Use simplified equation of state (S-EOS) 
    745757                                 ! 
    746    !                     ! S-EOS coefficients : 
    747                                  !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     758   !                     ! S-EOS coefficients (ln_seos=T): 
     759   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    748760   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
    749761   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
    … …  
    757769&namtra_adv    !   advection scheme for tracer                          (default: NO advection) 
    758770!----------------------------------------------------------------------- 
    759    ln_traadv_cen =  .false. !  2nd order centered scheme 
    760       nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
    761       nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
    762    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  
    765       nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
    766       !                             !        (number of sub-timestep = nn_fct_zts) 
    767    ln_traadv_mus =  .false. !  MUSCL scheme 
    768       ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
    769    ln_traadv_ubs =  .false. !  UBS scheme 
    770       nn_ubs_v   =  2               !  =2  , vertical 2nd order FCT / COMPACT 4th order 
    771    ln_traadv_qck =  .false. !  QUICKEST scheme 
     771   ln_traadv_cen = .false. !  2nd order centered scheme 
     772      nn_cen_h   =  4            !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     773      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
     774   ln_traadv_fct = .false. !  FCT scheme 
     775      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order  
     776      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order  
     777      nn_fct_zts =  0            !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     778      !                          !        (number of sub-timestep = nn_fct_zts) 
     779   ln_traadv_mus = .false. !  MUSCL scheme 
     780      ln_mus_ups = .false.       !  use upstream scheme near river mouths 
     781   ln_traadv_ubs = .false. !  UBS scheme 
     782      nn_ubs_v   =  2            !  =2  , vertical 2nd order FCT / COMPACT 4th order 
     783   ln_traadv_qck = .false. !  QUICKEST scheme 
    772784/ 
    773785!----------------------------------------------------------------------- 
    774786&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) 
    775787!----------------------------------------------------------------------- 
    776    ln_mle    = .false.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
    777    rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
    778    nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
    779    rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
    780    rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
    781    rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
    782    nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
    783    nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE 
    784    rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK 
     788   ln_mle      = .false.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     789   rn_ce       = 0.06      ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     790   nn_mle      = 1         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     791   rn_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
     792   rn_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
     793   rn_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     794   nn_mld_uv   = 0         ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
     795   nn_conv     = 0         ! =1 no MLE in case of convection ; =0 always MLE 
     796   rn_rho_c_mle= 0.01      ! delta rho criterion used to calculate MLD for FK 
    785797/ 
    786798!----------------------------------------------------------------------- 
    … …  
    791803   ln_traldf_lap   =  .false.  !    laplacian operator 
    792804   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     805   ! 
    793806   !                       !  Direction of action: 
    794807   ln_traldf_lev   =  .false.  !  iso-level 
    … …  
    866879   ln_vvl_layer  = .false.          !  full layer vertical coordinate 
    867880   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar 
    868    ln_vvl_zstar_at_eqtor = .false. !  ztilde near the equator 
     881   ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator 
    869882   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient 
    870883   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days] 
    … …  
    874887/ 
    875888!----------------------------------------------------------------------- 
    876 &namdyn_vor    !   option of physics/algorithm                          (default: NO) 
     889&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO) 
    877890!----------------------------------------------------------------------- 
    878891   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
    … …  
    953966   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0) 
    954967   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F) 
    955    nn_evdm     =    0      ! evd apply on tracer (=0) or on tracer and momentum (=1) 
    956    rn_avevd    =  100.     !  evd mixing coefficient [m2/s] 
     968      nn_evdm     =    0        ! evd apply on tracer (=0) or on tracer and momentum (=1) 
     969      rn_avevd    =  100.       !  evd mixing coefficient [m2/s] 
    957970   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F) 
    958    nn_npc      =    1            ! frequency of application of npc 
    959    nn_npcp     =  365            ! npc control print frequency 
     971      nn_npc      =    1        ! frequency of application of npc 
     972      nn_npcp     =  365        ! npc control print frequency 
    960973   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping 
    961    nn_zdfexp   =    3            ! number of sub-timestep for ln_zdfexp=T 
     974      nn_zdfexp   =    3        ! number of sub-timestep for ln_zdfexp=T 
    962975/ 
    963976!----------------------------------------------------------------------- 
    964977&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
    965978!----------------------------------------------------------------------- 
    966    rn_avmri    = 100.e-4   !  maximum value of the vertical viscosity 
    967    rn_alp      =   5.      !  coefficient of the parameterization 
    968    nn_ric      =   2       !  coefficient of the parameterization 
    969    rn_ekmfc    =   0.7     !  Factor in the Ekman depth Equation 
    970    rn_mldmin   =   1.0     !  minimum allowable mixed-layer depth estimate (m) 
    971    rn_mldmax   =1000.0     !  maximum allowable mixed-layer depth estimate (m) 
    972    rn_wtmix    =  10.0     !  vertical eddy viscosity coeff [m2/s] in the mixed-layer 
    973    rn_wvmix    =  10.0     !  vertical eddy diffusion coeff [m2/s] in the mixed-layer 
    974    ln_mldw     = .true.    !  Flag to use or not the mixed layer depth param. 
     979   rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity 
     980   rn_alp      =    5.     !  coefficient of the parameterization 
     981   nn_ric      =    2      !  coefficient of the parameterization 
     982   rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation 
     983   rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m) 
     984   rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m) 
     985   rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer 
     986   rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer 
     987   ln_mldw     =  .true.   !  Flag to use or not the mixed layer depth param. 
    975988/ 
    976989!----------------------------------------------------------------------- 
    … …  
    9921005   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
    9931006   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
    994    nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
     1007   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to near intertial waves 
    9951008                           !        = 0 no penetration 
    9961009                           !        = 1 add a tke source below the ML 
    9971010                           !        = 2 add a tke source just at the base of the ML 
    998                            !        = 3 as = 1 applied on HF part of the stress    ("key_oasis3") 
     1011                           !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T) 
    9991012   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 
    10001013   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML 
    … …  
    10031016/ 
    10041017!----------------------------------------------------------------------- 
    1005 &namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     1018&namzdf_gls    !   GLS vertical diffusion                               ("key_zdfgls") 
    10061019!----------------------------------------------------------------------- 
    10071020   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2] 
    … …  
    10361049   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency 
    10371050/ 
     1051!----------------------------------------------------------------------- 
     1052&namzdf_tmx_new !   internal wave-driven mixing parameterization        ("key_zdftmx_new" & "key_zdfddm") 
     1053!----------------------------------------------------------------------- 
     1054   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) 
     1055   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
     1056   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
     1057/ 
     1058 
    10381059 
    10391060!!====================================================================== 
    … …  
    10731094/ 
    10741095!----------------------------------------------------------------------- 
    1075 &namsto       ! Stochastic parametrization of EOS                       (default: NO) 
    1076 !----------------------------------------------------------------------- 
    1077    ln_sto_eos   = .false.  ! stochastic equation of state 
    1078    nn_sto_eos   = 1        ! number of independent random walks 
    1079    rn_eos_stdxy = 1.4      ! random walk horz. standard deviation (in grid points) 
    1080    rn_eos_stdz  = 0.7      ! random walk vert. standard deviation (in grid points) 
    1081    rn_eos_tcor  = 1440.    ! random walk time correlation (in timesteps) 
    1082    nn_eos_ord   = 1        ! order of autoregressive processes 
    1083    nn_eos_flt   = 0        ! passes of Laplacian filter 
    1084    rn_eos_lim   = 2.0      ! limitation factor (default = 3.0) 
    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 
     1096&namsto        ! Stochastic parametrization of EOS                      (default: NO) 
     1097!----------------------------------------------------------------------- 
     1098   ln_sto_eos  = .false.   ! stochastic equation of state 
     1099   nn_sto_eos  = 1         ! number of independent random walks 
     1100   rn_eos_stdxy= 1.4       ! random walk horz. standard deviation (in grid points) 
     1101   rn_eos_stdz = 0.7       ! random walk vert. standard deviation (in grid points) 
     1102   rn_eos_tcor = 1440.     ! random walk time correlation (in timesteps) 
     1103   nn_eos_ord  = 1         ! order of autoregressive processes 
     1104   nn_eos_flt  = 0         ! passes of Laplacian filter 
     1105   rn_eos_lim  = 2.0       ! limitation factor (default = 3.0) 
     1106   ln_rststo   = .false.   ! start from mean parameter (F) or from restart file (T) 
     1107   ln_rstseed  = .true.    ! read seed of RNG from restart file 
    10871108   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input) 
    10881109   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output) 
    … …  
    10921113!!                  ***  Diagnostics namelists  *** 
    10931114!!====================================================================== 
    1094 !!   namtrd       dynamics and/or tracer trends 
    1095 !!   namptr       Poleward Transport Diagnostics 
    1096 !!   namhsb       Heat and salt budgets 
     1115!!   namtrd       dynamics and/or tracer trends                         (default F) 
     1116!!   namptr       Poleward Transport Diagnostics                        (default F) 
     1117!!   namhsb       Heat and salt budgets                                 (default F) 
     1118!!   namdiu       Cool skin and warm layer models                       (default F) 
    10971119!!   namflo       float parameters                                      ("key_float") 
    1098 !!   nam_diaharm  Harmonic analysis of tidal constituents               ('key_diaharm') 
    1099 !!   namdct       transports through some sections 
     1120!!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
     1121!!   namdct       transports through some sections                      ("key_diadct") 
     1122!!   nam_diatmb   Top Middle Bottom Output                              (default F) 
     1123!!   nam_dia25h   25h Mean Output                                       (default F) 
    11001124!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    11011125!!====================================================================== 
    11021126! 
    11031127!----------------------------------------------------------------------- 
    1104 &namtrd        !   diagnostics on dynamics and/or tracer trends         (default F) 
    1105 !              !   and/or mixed-layer trends and/or barotropic vorticity 
     1128&namtrd        !   trend diagnostics                                    (default F) 
    11061129!----------------------------------------------------------------------- 
    11071130   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
    … …  
    11231146!!gm 
    11241147!----------------------------------------------------------------------- 
    1125 &namptr       !   Poleward Transport Diagnostic                         (default F) 
    1126 !----------------------------------------------------------------------- 
    1127    ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
    1128    ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    1129 / 
    1130 !----------------------------------------------------------------------- 
    1131 &namhsb       !  Heat and salt budgets                                  (default F) 
    1132 !----------------------------------------------------------------------- 
    1133    ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F) 
    1134 / 
    1135 !----------------------------------------------------------------------- 
    1136 &namflo       !   float parameters                                      ("key_float") 
    1137 !----------------------------------------------------------------------- 
    1138    jpnfl         = 1          !  total number of floats during the run 
    1139    jpnnewflo     = 0          !  number of floats for the restart 
    1140    ln_rstflo     = .false.    !  float restart (T) or not (F) 
    1141    nn_writefl    =      75    !  frequency of writing in float output file 
    1142    nn_stockfl    =    5475    !  frequency of creation of the float restart file 
    1143    ln_argo       = .false.    !  Argo type floats (stay at the surface each 10 days) 
    1144    ln_flork4     = .false.    !  trajectories computed with a 4th order Runge-Kutta (T) 
    1145                               !  or computed with Blanke' scheme (F) 
    1146    ln_ariane     = .true.     !  Input with Ariane tool convention(T) 
    1147    ln_flo_ascii  = .true.     !  Output with Ariane tool netcdf convention(F) or ascii file (T) 
    1148 / 
    1149 !----------------------------------------------------------------------- 
    1150 &nam_diaharm   !   Harmonic analysis of tidal constituents              ('key_diaharm') 
     1148&namptr        !   Poleward Transport Diagnostic                         (default F) 
     1149!----------------------------------------------------------------------- 
     1150   ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F) 
     1151   ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not 
     1152/ 
     1153!----------------------------------------------------------------------- 
     1154&namhsb        !  Heat and salt budgets                                  (default F) 
     1155!----------------------------------------------------------------------- 
     1156   ln_diahsb   = .false.   !  check the heat and salt budgets (T) or not (F) 
     1157/ 
     1158!----------------------------------------------------------------------- 
     1159&namdiu        !   Cool skin and warm layer models                       (default F) 
     1160!----------------------------------------------------------------------- 
     1161   ln_diurnal      = .false.   !  
     1162   ln_diurnal_only = .false.   ! 
     1163/ 
     1164!----------------------------------------------------------------------- 
     1165&namflo        !   float parameters                                      ("key_float") 
     1166!----------------------------------------------------------------------- 
     1167   jpnfl       = 1         !  total number of floats during the run 
     1168   jpnnewflo   = 0         !  number of floats for the restart 
     1169   ln_rstflo   = .false.   !  float restart (T) or not (F) 
     1170   nn_writefl  =      75   !  frequency of writing in float output file 
     1171   nn_stockfl  =    5475   !  frequency of creation of the float restart file 
     1172   ln_argo     = .false.   !  Argo type floats (stay at the surface each 10 days) 
     1173   ln_flork4   = .false.   !  trajectories computed with a 4th order Runge-Kutta (T) 
     1174   !                       !  or computed with Blanke' scheme (F) 
     1175   ln_ariane   = .true.    !  Input with Ariane tool convention(T) 
     1176   ln_flo_ascii= .true.    !  Output with Ariane tool netcdf convention(F) or ascii file (T) 
     1177/ 
     1178!----------------------------------------------------------------------- 
     1179&nam_diaharm   !   Harmonic analysis of tidal constituents               ("key_diaharm") 
    11511180!----------------------------------------------------------------------- 
    11521181    nit000_han = 1         ! First time step used for harmonic analysis 
    … …  
    11571186/ 
    11581187!----------------------------------------------------------------------- 
    1159 &namdct        ! transports through some sections 
    1160 !----------------------------------------------------------------------- 
    1161     nn_dct      = 15       !  time step frequency for transports computing 
    1162     nn_dctwri   = 15       !  time step frequency for transports writing 
    1163     nn_secdebug = 112      !      0 : no section to debug 
    1164                            !     -1 : debug all section 
    1165                            !  0 < n : debug section number n 
    1166 / 
    1167 !----------------------------------------------------------------------- 
    1168 &namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
     1188&namdct        ! transports through some sections                        ("key_diadct") 
     1189!----------------------------------------------------------------------- 
     1190    nn_dct     = 15        !  time step frequency for transports computing 
     1191    nn_dctwri  = 15        !  time step frequency for transports writing 
     1192    nn_secdebug= 112       !      0 : no section to debug 
     1193    !                      !     -1 : debug all section 
     1194    !                      !  0 < n : debug section number n 
     1195/ 
     1196!----------------------------------------------------------------------- 
     1197&nam_diatmb    !  Top Middle Bottom Output                               (default F) 
     1198!----------------------------------------------------------------------- 
     1199   ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
     1200/ 
     1201!----------------------------------------------------------------------- 
     1202&nam_dia25h    !  25h Mean Output                                        (default F) 
     1203!----------------------------------------------------------------------- 
     1204   ln_dia25h   = .false.   ! Choose 25h mean output or not 
     1205/ 
     1206!----------------------------------------------------------------------- 
     1207&namnc4        !   netcdf4 chunking and compression settings             ("key_netcdf4") 
    11691208!----------------------------------------------------------------------- 
    11701209   nn_nchunks_i=   4       !  number of chunks in i-dimension 
    11711210   nn_nchunks_j=   4       !  number of chunks in j-dimension 
    11721211   nn_nchunks_k=   31      !  number of chunks in k-dimension 
    1173                            !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
    1174                            !  is optimal for postprocessing which works exclusively with horizontal slabs 
     1212   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
     1213   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs 
    11751214   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
    1176                            !  (F) ignore chunking information and produce netcdf3-compatible files 
     1215   !                       !  (F) ignore chunking information and produce netcdf3-compatible files 
    11771216/ 
    11781217 
    … …  
    11851224! 
    11861225!----------------------------------------------------------------------- 
    1187 &namobs       !  observation usage switch 
    1188 !----------------------------------------------------------------------- 
    1189    ln_diaobs  = .false.             ! Logical switch for the observation operator 
    1190    ln_t3d     = .false.             ! Logical switch for T profile observations 
    1191    ln_s3d     = .false.             ! Logical switch for S profile observations 
    1192    ln_sla     = .false.             ! Logical switch for SLA observations 
    1193    ln_sst     = .false.             ! Logical switch for SST observations 
    1194    ln_sic     = .false.             ! Logical switch for Sea Ice observations 
    1195    ln_vel3d   = .false.             ! Logical switch for velocity observations 
    1196    ln_altbias = .false.             ! Logical switch for altimeter bias correction 
    1197    ln_nea     = .false.             ! Logical switch for rejection of observations near land 
    1198    ln_grid_global = .true.          ! Logical switch for global distribution of observations 
    1199    ln_grid_search_lookup = .false.  ! Logical switch for obs grid search w/lookup table 
    1200    ln_ignmis  = .true.              ! Logical switch for ignoring missing files 
    1201    ln_s_at_t  = .false.             ! Logical switch for computing model S at T obs if not there 
    1202    ln_sstnight = .false.            ! Logical switch for calculating night-time average for SST obs 
     1226&namobs        !  observation usage switch 
     1227!----------------------------------------------------------------------- 
     1228   ln_diaobs   = .false.             ! Logical switch for the observation operator 
     1229   ln_t3d      = .false.             ! Logical switch for T profile observations 
     1230   ln_s3d      = .false.             ! Logical switch for S profile observations 
     1231   ln_sla      = .false.             ! Logical switch for SLA observations 
     1232   ln_sst      = .false.             ! Logical switch for SST observations 
     1233   ln_sic      = .false.             ! Logical switch for Sea Ice observations 
     1234   ln_vel3d    = .false.             ! Logical switch for velocity observations 
     1235   ln_altbias  = .false.             ! Logical switch for altimeter bias correction 
     1236   ln_nea      = .false.             ! Logical switch for rejection of observations near land 
     1237   ln_grid_global = .true.           ! Logical switch for global distribution of observations 
     1238   ln_grid_search_lookup = .false.   ! Logical switch for obs grid search w/lookup table 
     1239   ln_ignmis   = .true.              ! Logical switch for ignoring missing files 
     1240   ln_s_at_t   = .false.             ! Logical switch for computing model S at T obs if not there 
     1241   ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs 
    12031242! All of the *files* variables below are arrays. Use namelist_cfg to add more files 
    1204    cn_profbfiles = 'profiles_01.nc'    ! Profile feedback input observation file names 
    1205    cn_slafbfiles = 'sla_01.nc'         ! SLA feedback input observation file names 
    1206    cn_sstfbfiles = 'sst_01.nc'         ! SST feedback input observation file names 
    1207    cn_sicfbfiles = 'sic_01.nc'         ! SIC feedback input observation file names 
    1208    cn_velfbfiles = 'vel_01.nc'         ! Velocity feedback input observation file names 
    1209    cn_altbiasfile = 'altbias.nc'       ! Altimeter bias input file name 
    1210    cn_gridsearchfile = 'gridsearch.nc' ! Grid search file name 
    1211    rn_gridsearchres = 0.5              ! Grid search resolution 
    1212    rn_dobsini = 00010101.000000        ! Initial date in window YYYYMMDD.HHMMSS 
    1213    rn_dobsend = 00010102.000000        ! Final date in window YYYYMMDD.HHMMSS 
    1214    nn_1dint = 0                        ! Type of vertical interpolation method 
    1215    nn_2dint = 0                        ! Type of horizontal interpolation method 
    1216    nn_msshc = 0                        ! MSSH correction scheme 
    1217    rn_mdtcorr = 1.61                   ! MDT  correction 
    1218    rn_mdtcutoff = 65.0                 ! MDT cutoff for computed correction 
    1219    nn_profdavtypes = -1                ! Profile daily average types - array 
    1220    ln_sstbias = .false.  
    1221    cn_sstbias_files = 'sstbias.nc' 
    1222 / 
    1223 !----------------------------------------------------------------------- 
    1224 &nam_asminc   !   assimilation increments                               ('key_asminc') 
    1225 !----------------------------------------------------------------------- 
    1226     ln_bkgwri = .false.    !  Logical switch for writing out background state 
    1227     ln_trainc = .false.    !  Logical switch for applying tracer increments 
    1228     ln_dyninc = .false.    !  Logical switch for applying velocity increments 
    1229     ln_sshinc = .false.    !  Logical switch for applying SSH increments 
    1230     ln_asmdin = .false.    !  Logical switch for Direct Initialization (DI) 
    1231     ln_asmiau = .false.    !  Logical switch for Incremental Analysis Updating (IAU) 
    1232     nitbkg    = 0          !  Timestep of background in [0,nitend-nit000-1] 
    1233     nitdin    = 0          !  Timestep of background for DI in [0,nitend-nit000-1] 
    1234     nitiaustr = 1          !  Timestep of start of IAU interval in [0,nitend-nit000-1] 
    1235     nitiaufin = 15         !  Timestep of end of IAU interval in [0,nitend-nit000-1] 
    1236     niaufn    = 0          !  Type of IAU weighting function 
    1237     ln_salfix = .false.    !  Logical switch for ensuring that the sa > salfixmin 
    1238     salfixmin = -9999      !  Minimum salinity after applying the increments 
    1239     nn_divdmp = 0          !  Number of iterations of divergence damping operator 
    1240 / 
    1241 !----------------------------------------------------------------------- 
    1242 &namdiu !   Cool skin and warm layer models 
    1243 !----------------------------------------------------------------------- 
    1244    ln_diurnal      = .false.   !  
    1245    ln_diurnal_only = .false.   ! 
    1246 / 
    1247 !----------------------------------------------------------------------- 
    1248 &nam_diatmb  !  Top Middle Bottom Output 
    1249 !----------------------------------------------------------------------- 
    1250    ln_diatmb  = .false.    !  Choose Top Middle and Bottom output or not 
    1251 / 
    1252 !----------------------------------------------------------------------- 
    1253 &namwad  !   Wetting and drying 
    1254 !----------------------------------------------------------------------- 
    1255    ln_wd             = .false.  ! T/F activation of wetting and drying 
    1256    rn_wdmin1         =  0.1     ! Minimum wet depth on dried cells 
    1257    rn_wdmin2         =  0.01    ! Tolerance of min wet depth on dried cells 
    1258    rn_wdld           =  20.0    ! Land elevation below which wetting/drying is allowed 
    1259    nn_wdit           =  10      ! Max iterations for W/D limiter 
    1260 / 
    1261 !----------------------------------------------------------------------- 
    1262 &nam_dia25h  !  25h Mean Output 
    1263 !----------------------------------------------------------------------- 
    1264    ln_dia25h  = .false.    ! Choose 25h mean output or not 
    1265 / 
     1243   cn_profbfiles = 'profiles_01.nc'  ! Profile feedback input observation file names 
     1244   cn_slafbfiles = 'sla_01.nc'       ! SLA feedback input observation file names 
     1245   cn_sstfbfiles = 'sst_01.nc'       ! SST feedback input observation file names 
     1246   cn_sicfbfiles = 'sic_01.nc'       ! SIC feedback input observation file names 
     1247   cn_velfbfiles = 'vel_01.nc'       ! Velocity feedback input observation file names 
     1248   cn_altbiasfile = 'altbias.nc'     ! Altimeter bias input file name 
     1249   cn_gridsearchfile='gridsearch.nc' ! Grid search file name 
     1250   rn_gridsearchres = 0.5            ! Grid search resolution 
     1251   rn_dobsini  = 00010101.000000     ! Initial date in window YYYYMMDD.HHMMSS 
     1252   rn_dobsend  = 00010102.000000     ! Final date in window YYYYMMDD.HHMMSS 
     1253   nn_1dint    = 0                   ! Type of vertical interpolation method 
     1254   nn_2dint    = 0                   ! Type of horizontal interpolation method 
     1255   nn_msshc    = 0                   ! MSSH correction scheme 
     1256   rn_mdtcorr  = 1.61                ! MDT  correction 
     1257   rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction 
     1258   nn_profdavtypes = -1              ! Profile daily average types - array 
     1259   ln_sstbias  = .false.             ! 
     1260   cn_sstbias_files = 'sstbias.nc'   ! 
     1261/ 
     1262!----------------------------------------------------------------------- 
     1263&nam_asminc    !   assimilation increments                              ('key_asminc') 
     1264!----------------------------------------------------------------------- 
     1265    ln_bkgwri  = .false.   !  Logical switch for writing out background state 
     1266    ln_trainc  = .false.   !  Logical switch for applying tracer increments 
     1267    ln_dyninc  = .false.   !  Logical switch for applying velocity increments 
     1268    ln_sshinc  = .false.   !  Logical switch for applying SSH increments 
     1269    ln_asmdin  = .false.   !  Logical switch for Direct Initialization (DI) 
     1270    ln_asmiau  = .false.   !  Logical switch for Incremental Analysis Updating (IAU) 
     1271    nitbkg     = 0         !  Timestep of background in [0,nitend-nit000-1] 
     1272    nitdin     = 0         !  Timestep of background for DI in [0,nitend-nit000-1] 
     1273    nitiaustr  = 1         !  Timestep of start of IAU interval in [0,nitend-nit000-1] 
     1274    nitiaufin  = 15        !  Timestep of end of IAU interval in [0,nitend-nit000-1] 
     1275    niaufn     = 0         !  Type of IAU weighting function 
     1276    ln_salfix  = .false.   !  Logical switch for ensuring that the sa > salfixmin 
     1277    salfixmin  = -9999     !  Minimum salinity after applying the increments 
     1278    nn_divdmp  = 0         !  Number of iterations of divergence damping operator 
     1279/ 
Note: See TracChangeset for help on using the changeset viewer.