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 9490 for branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED/namelist_ref – NEMO

Ignore:
Timestamp:
2018-04-23T10:44:07+02:00 (6 years ago)
Author:
gm
Message:

#2075 - dev_merge_2017: scale-aware setting of lateral viscous and diffusive coefficient

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED/namelist_ref

    r9367 r9490  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !!                            namelist_ref                            !! 
     2!! NEMO/OCE :   Reference namelist_ref                                !! 
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!! NEMO/OPA  :  1 - run manager      (namrun) 
     
    99!!                                    namsbc_ssr, namsbc_wave, namberg) 
    1010!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    11 !!              5 - bottom boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl) 
    12 !!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
     11!!              5 - top/bot boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl) 
     12!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_eiv, namtra_dmp) 
    1313!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    14 !!              8 - Verical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm) 
     14!!              8 - Vertical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm) 
    1515!!              9 - miscellaneous    (nammpp, namctl) 
    1616!!             10 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
     
    1919 
    2020!!====================================================================== 
    21 !!                   ***  Run management namelists  ***               !! 
    22 !!====================================================================== 
     21!!              ***  Domain & Run management namelists  ***           !! 
     22!!                                                                    !! 
    2323!!   namrun       parameters of the run 
     24!!   namdom       space and time domain 
     25!!   namcfg       parameters of the configuration                       (default: user defined GYRE) 
     26!!   namwad       Wetting and drying                                    (default NO) 
     27!!   namtsd       data: temperature & salinity                          (default NO) 
     28!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T) 
     29!!   namc1d       1D configuration options                              ("key_c1d") 
     30!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
     31!!   namc1d_uvd   1D data (currents)                                    ("key_c1d") 
    2432!!====================================================================== 
    2533! 
     
    2735&namrun        !   parameters of the run 
    2836!----------------------------------------------------------------------- 
    29    nn_no       =       0   !  job number (no more used...) 
     37   nn_no       =       0   !  Assimilation cycle index 
    3038   cn_exp      =  "ORCA2"  !  experience name 
    3139   nn_it000    =       1   !  first time step 
     
    3543   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    3644   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    37       nn_euler    =    1            !  = 0 : start with forward time step if ln_rstart=T 
    38       nn_rstctl   =    0            !  restart control ==> activated only if ln_rstart=T 
    39       !                             !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    40       !                             !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
    41       !                             !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
     45      nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
     46      nn_rstctl   =    0      !  restart control ==> activated only if ln_rstart=T 
     47      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
     48      !                          !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     49      !                          !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
    4250      cn_ocerst_in    = "restart"   !  suffix of ocean restart name (input) 
    4351      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts 
     
    5765   nn_wxios = 0      !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 
    5866/ 
    59 ! 
    60 !!====================================================================== 
    61 !!                      ***  Domain namelists  *** 
    62 !!====================================================================== 
    63 !!   namcfg       parameters of the configuration 
    64 !!   namdom       space and time domain (bathymetry, mesh, timestep) 
    65 !!   namwad       Wetting and drying                                    (default F) 
    66 !!   namtsd       data: temperature & salinity 
    67 !!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T) 
    68 !!   namc1d       1D configuration options                              ("key_c1d") 
    69 !!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
    70 !!   namc1d_uvd   1D data (currents)                                    ("key_c1d") 
    71 !!====================================================================== 
    72 ! 
    73 !----------------------------------------------------------------------- 
    74 &namcfg        !   parameters of the configuration                      (default: user defined GYRE) 
     67!----------------------------------------------------------------------- 
     68&namdom        !   time and space domain                                 
     69!----------------------------------------------------------------------- 
     70   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
     71   rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
     72   ! 
     73   rn_rdt      = 5760.     !  time step for the dynamics and tracer 
     74   rn_atfp     =    0.1    !  asselin time filter parameter 
     75   ! 
     76   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
     77   ! 
     78   ln_meshmask = .false.   !  =T create a mesh file 
     79/ 
     80!----------------------------------------------------------------------- 
     81&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg) 
    7582!----------------------------------------------------------------------- 
    7683   ln_read_cfg = .false.   !  (=T) read the domain configuration file 
    77       !                    !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules 
     84      !                    !  (=F) user defined configuration           (F => create/check namusr_def) 
    7885      cn_domcfg = "domain_cfg"  ! domain configuration filename 
    7986      ! 
     
    9198/ 
    9299!----------------------------------------------------------------------- 
    93 &namdom        !   time and space domain 
    94 !----------------------------------------------------------------------- 
    95    ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    96    rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    97    ! 
    98    rn_rdt      = 5760.     !  time step for the dynamics and tracer 
    99    rn_atfp     =    0.1    !  asselin time filter parameter 
    100    ! 
    101    ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    102    ! 
    103    ln_meshmask = .false.   !  =T create a mesh file 
    104 / 
    105 !----------------------------------------------------------------------- 
    106 &namtsd        !   data : Temperature  & Salinity 
    107 !----------------------------------------------------------------------- 
    108 !              !  file name                 ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    109 !              !                            !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    110    sn_tem = 'data_1m_potential_temperature_nomask',     -1      ,'votemper',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    111    sn_sal = 'data_1m_salinity_nomask'             ,     -1      ,'vosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    112    ! 
    113    cn_dir      = './'      !  root directory for the location of the runoff files 
    114    ln_tsd_init = .true.    !  Initialisation of ocean T & S with T & S input data (T) or not (F) 
    115    ln_tsd_tradmp = .true.  !  damping of ocean T & S toward T & S input data (T) or not (F) 
    116 / 
    117 !----------------------------------------------------------------------- 
    118 &namwad        !   Wetting and drying  default is no WAD 
    119 !----------------------------------------------------------------------- 
    120    ln_wd_il          = .false   ! T/F activation of iterative limiter for  wetting and drying scheme 
    121    ln_wd_dl          = .false.   ! T/F activation of directional llimiter for wetting drying scheme 
    122    ln_wd_dl_bc       = .false.   ! T/F Directional limiteer Baroclinic option 
    123    ln_wd_dl_rmp      = .false.   ! T/F Turn on directional limiter ramp 
    124    rn_wdmin0         =  0.30    ! dpoth at which wetting/drying starts 
    125    rn_wdmin1         =  0.2     ! Minimum wet depth on dried cells 
    126    rn_wdmin2         =  0.0001  ! Tolerance of min wet depth on dried cells 
    127    rn_wdld           =  2.5     ! Land elevation below which wetting/drying is allowed 
    128    nn_wdit           =   20     ! Max iterations for W/D limiter 
     100&namtsd        !    Temperature & Salinity Data  (init/dmp)             (default: NO) 
     101!----------------------------------------------------------------------- 
     102   !                       ! =T  read T-S fields for: 
     103   ln_tsd_init = .false.         !  ocean initialisation 
     104   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp) 
     105    
     106   cn_dir      = './'      !  root directory for the T-S data location 
     107   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     108   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     109   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     110   sn_tem = 'data_1m_potential_temperature_nomask',  -1      , 'votemper',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     111   sn_sal = 'data_1m_salinity_nomask'             ,  -1      , 'vosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     112/ 
     113!----------------------------------------------------------------------- 
     114&namwad        !   Wetting and Drying (WaD)                             (default: NO) 
     115!----------------------------------------------------------------------- 
     116   ln_wd_il    = .false    !  T/F activation of iterative   limiter 
     117   ln_wd_dl    = .false.   !  T/F activation of directional limiter 
     118   ln_wd_dl_bc = .false.   !  T/F Directional limiteer Baroclinic option 
     119   ln_wd_dl_rmp= .false.   !  T/F Turn on directional limiter ramp 
     120   rn_wdmin0   =  0.30     !  depth at which WaD starts 
     121   rn_wdmin1   =  0.2      !  Minimum wet depth on dried cells 
     122   rn_wdmin2   =  0.0001   !  Tolerance of min wet depth on dried cells 
     123   rn_wdld     =  2.5      !  Land elevation below which WaD is allowed 
     124   nn_wdit     =   20      !  Max iterations for WaD limiter 
    129125/ 
    130126!----------------------------------------------------------------------- 
     
    134130   nn_facty    = 3         !  Reduction factor of y-direction 
    135131   nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT 
    136                            !  0, coarse grid is binned with preferential treatment of the north fold 
    137                            !  1, coarse grid is binned with centering at the equator 
    138                            !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty. 
     132      !                    !  0, coarse grid is binned with preferential treatment of the north fold 
     133      !                    !  1, coarse grid is binned with centering at the equator 
     134      !                    !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty. 
    139135   ln_msh_crs  = .false.   ! =T create a mesh & mask file 
    140136   nn_crs_kz   = 0         ! 0, MEAN of volume boxes 
    141                            ! 1, MAX of boxes 
    142                            ! 2, MIN of boxes 
     137      !                    ! 1, MAX of boxes 
     138      !                    ! 2, MIN of boxes 
    143139   ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F ) 
    144140/ 
    145141!----------------------------------------------------------------------- 
    146 &namc1d        !   1D configuration options                             ("key_c1d") 
     142&namc1d        !   1D configuration options                             ("key_c1d" default: PAPA station) 
    147143!----------------------------------------------------------------------- 
    148144   rn_lat1d    =      50   !  Column latitude (default at PAPA station) 
     
    151147/ 
    152148!----------------------------------------------------------------------- 
    153 &namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d") 
     149&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: NO) 
    154150!----------------------------------------------------------------------- 
    155151   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    156152/ 
    157153!----------------------------------------------------------------------- 
    158 &namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
    159 !----------------------------------------------------------------------- 
    160 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    161 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    162    sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , '' 
    163    sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , '' 
    164 ! 
    165    cn_dir        = './'    !  root directory for the location of the files 
    166    ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F) 
    167    ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F) 
    168 / 
    169  
    170 !!====================================================================== 
    171 !!            ***  Surface Boundary Condition namelists  *** 
    172 !!====================================================================== 
    173 !!   namsbc          surface boundary condition 
     154&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: NO) 
     155!-----------------------------------------------------------------------    
     156   !                       !  =T read U-V fields for: 
     157   ln_uvd_init   = .false.       !  ocean initialisation 
     158   ln_uvd_dyndmp = .false.       !  U-V restoring 
     159 
     160   cn_dir      = './'      !  root directory for the U-V data location 
     161   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     162   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     163   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     164   sn_ucur     = 'ucurrent'              ,         -1        ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , '' 
     165   sn_vcur     = 'vcurrent'              ,         -1        ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , '' 
     166/ 
     167 
     168!!====================================================================== 
     169!!            ***  Surface Boundary Condition namelists  ***          !! 
     170!!                                                                    !! 
     171!!   namsbc          surface boundary condition manager                 (default: NO selection) 
    174172!!   namsbc_flx      flux               formulation                     (ln_flx     =T) 
    175173!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T) 
    176174!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" ) 
    177 !!   namsbc_sas      Stand-Alone Surface module 
     175!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only) 
     176!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   ) 
    178177!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T) 
     178!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
    179179!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    180 !!   namsbc_isf      ice shelf melting/freezing                         (nn_isf     >0) 
    181 !!   namsbc_iscpl    coupling option between land ice model and ocean 
    182180!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    183 !!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
     181!!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
     182!!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    184183!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    185184!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     
    187186! 
    188187!----------------------------------------------------------------------- 
    189 &namsbc        !   Surface Boundary Condition (surface module) 
    190 !----------------------------------------------------------------------- 
    191    nn_fsbc     = 5         !  frequency of surface boundary condition computation 
    192                            !     (also = the frequency of sea-ice & iceberg model call) 
     188&namsbc        !   Surface Boundary Condition manager                   (default: NO selection) 
     189!----------------------------------------------------------------------- 
     190   nn_fsbc     = 5         !  frequency of SBC module call 
     191      !                    !  (control sea-ice & iceberg model call) 
    193192                     ! Type of air-sea fluxes  
    194193   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc) 
    195194   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    196    ln_blk      = .false.    !  Bulk formulation                         (T => fill namsbc_blk ) 
    197                      ! Type of coupling (Ocean/Ice/Atmosphere) : 
     195   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk ) 
     196      !              ! Type of coupling (Ocean/Ice/Atmosphere) : 
    198197   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
    199198   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 ) 
    200199   nn_components = 0       !  configuration of the opa-sas OASIS coupling 
    201                            !  =0 no opa-sas OASIS coupling: default single executable config. 
    202                            !  =1 opa-sas OASIS coupling: multi executable config., OPA component 
    203                            !  =2 opa-sas OASIS coupling: multi executable config., SAS component  
     200      !                    !  =0 no opa-sas OASIS coupling: default single executable config. 
     201      !                    !  =1 opa-sas OASIS coupling: multi executable config., OPA component 
     202      !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component  
    204203                     ! Sea-ice : 
    205    nn_ice      = 2         !  =0 no ice boundary condition     
    206                            !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
    207                            !  =2 or 3 automatically for LIM3 or CICE    ("key_lim3" or "key_cice") 
    208                            !          except in AGRIF zoom where it has to be specified 
     204   nn_ice      = 0         !  =0 no ice boundary condition     
     205      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
     206      !                    !  =2 or 3 automatically for LIM3 or CICE    ("key_lim3" or "key_cice") 
     207      !                    !          except in AGRIF zoom where it has to be specified 
    209208   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
    210                            !  =F levitating ice (no pressure, mass and salt exchanges) 
     209      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
    211210                     ! Misc. options of sbc :  
    212    ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
     211   ln_traqsr   = .false.   !  Light penetration in the ocean            (T => fill namtra_qsr) 
    213212   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    214    ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
    215    ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    216    nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
    217                            !     =1 global mean of e-p-r set to zero at each time step 
    218                            !     =2 annual global mean of e-p-r set to zero 
     213   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
     214   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
     215      !                    !     =1 global mean of e-p-r set to zero at each time step 
     216      !                    !     =2 annual global mean of e-p-r set to zero 
     217   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    219218   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    220    ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf) 
     219   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf & namsbc_iscpl) 
    221220   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
    222221   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
    223222   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)  
    224223   nn_sdrift   =  0        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 
    225                            !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] 
    226                            !   = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))] 
    227                            !   = 2 Phillips as (1) but using the wave frequency from a wave model 
     224      !                    !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] 
     225      !                    !   = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))] 
     226      !                    !   = 2 Phillips as (1) but using the wave frequency from a wave model 
    228227   ln_tauwoc   = .false.   !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave) 
    229228   ln_tauw     = .false.   !  Activate ocean stress components from wave model 
     
    233232/ 
    234233!----------------------------------------------------------------------- 
    235 &namsbc_flx    !   surface boundary condition : flux formulation 
    236 !----------------------------------------------------------------------- 
    237 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    238 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    239    sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    240    sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    241    sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    242    sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    243    sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    244    ! 
    245    cn_dir      = './'      !  root directory for the location of the flux files 
    246 / 
    247 !----------------------------------------------------------------------- 
    248 &namsbc_blk   !   namsbc_blk  generic Bulk formula                      (ln_blk =T) 
    249 !----------------------------------------------------------------------- 
    250 !              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                              ! rotation ! land/sea mask ! 
    251 !              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                             ! pairing  ! filename      ! 
    252    sn_wndi     = 'u_10.15JUNE2009_fill'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , '' 
    253    sn_wndj     = 'v_10.15JUNE2009_fill'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , '' 
    254    sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    255    sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    256    sn_tair     = 't_10.15JUNE2009_fill'       ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    257    sn_humi     = 'q_10.15JUNE2009_fill'       ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    258    sn_prec     = 'ncar_precip.15JUNE2009_fill',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    259    sn_snow     = 'ncar_precip.15JUNE2009_fill',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    260    sn_slp      = 'slp.15JUNE2009_fill'        ,         6         , 'SLP'     ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    261    sn_tdif     = 'taudif_core'                ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     234&namsbc_flx    !   surface boundary condition : flux formulation        (ln_flx =T) 
     235!----------------------------------------------------------------------- 
     236   cn_dir      = './'      !  root directory for the fluxes data location 
     237   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     238   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     239   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     240   sn_utau     = 'utau'                  ,        24         , 'utau'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     241   sn_vtau     = 'vtau'                  ,        24         , 'vtau'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     242   sn_qtot     = 'qtot'                  ,        24         , 'qtot'    , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     243   sn_qsr      = 'qsr'                   ,        24         , 'qsr'     , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     244   sn_emp      = 'emp'                   ,        24         , 'emp'     , .false.     , .false., 'yearly'  , ''               , ''       , '' 
     245/ 
     246!----------------------------------------------------------------------- 
     247&namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
     248!----------------------------------------------------------------------- 
    262249   !                    !  bulk algorithm : 
    263250   ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
     
    265252   ln_COARE_3p5= .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    266253   ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
    267    ! 
    268    cn_dir      = './'      !  root directory for the location of the bulk files 
    269    ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    270    rn_zqt      = 10.       !  Air temperature and humidity reference height (m) 
    271    rn_zu       = 10.       !  Wind vector reference height (m) 
    272    rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    273    rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    274    rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    275    !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    276    ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    277    ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
     254      ! 
     255      rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
     256      rn_zu       = 10.       !  Wind vector reference height (m) 
     257      ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
     258      ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
     259      ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
     260      rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
     261      rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     262      rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
     263      !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
     264 
     265   cn_dir      = './'      !  root directory for the bulk data location 
     266   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
     267   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask ! 
     268   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   ! 
     269   sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6         , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Uwnd'   , '' 
     270   sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6         , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Vwnd'   , '' 
     271   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24         , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     272   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24         , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     273   sn_tair     = 't_10.15JUNE2009_fill'       ,    6         , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     274   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6         , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     275   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1         , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     276   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1         , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     277   sn_slp      = 'slp.15JUNE2009_fill'        ,    6         , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     278   sn_tdif     = 'taudif_core'                ,   24         , 'taudif'  ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    278279/ 
    279280!----------------------------------------------------------------------- 
    280281&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
    281282!----------------------------------------------------------------------- 
    282 !                    !     description      !  multiple  !    vector   !      vector          ! vector ! 
    283 !                    !                      ! categories !  reference  !    orientation       ! grids  ! 
    284 ! send 
     283   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data 
     284   ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models 
     285   !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
     286   nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
     287 
     288   !_____________!__________________________!____________!_____________!______________________!________! 
     289   !             !        description       !  multiple  !    vector   !       vector         ! vector ! 
     290   !             !                          ! categories !  reference  !     orientation      ! grids  ! 
     291!***   send    *** 
    285292   sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    286293   sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     
    296303   sn_snd_sstfrz =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    297304   sn_snd_ttilyr =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    298 ! receive 
     305!***  receive  *** 
    299306   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    300307   sn_rcv_taumod =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     
    323330   sn_rcv_tauw   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    324331   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    325 ! 
    326    nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
    327    ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models 
    328    !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    329    nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1) 
    330 / 
    331 !----------------------------------------------------------------------- 
    332 &namsbc_sas    !   Stand-Alone Surface boundary condition 
    333 !----------------------------------------------------------------------- 
    334 !              !  file name  ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    335 !              !             !  (if <0  months)  !   name    !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    336    sn_usp      = 'sas_grid_U',     120           , 'uos'     ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    337    sn_vsp      = 'sas_grid_V',     120           , 'vos'     ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    338    sn_tem      = 'sas_grid_T',     120           , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    339    sn_sal      = 'sas_grid_T',     120           , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    340    sn_ssh      = 'sas_grid_T',     120           , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    341    sn_e3t      = 'sas_grid_T',     120           , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    342    sn_frq      = 'sas_grid_T',     120           , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    343  
    344    l_sasread   = .true.    !  =T Read the above fields in a file, =F initialize to 0. in sbcssm.F90 
    345    ln_3d_uve   = .false.   !  specify whether we are supplying a 3D u,v and e3 field 
    346    ln_read_frq = .false.   !  specify whether we must read frq or not 
    347    cn_dir      = './'      !  root directory for the location of the bulk files are 
     332/ 
     333!----------------------------------------------------------------------- 
     334&namsbc_sas    !   Stand-Alone Surface module: ocean data               (SAS_SRC  only) 
     335!----------------------------------------------------------------------- 
     336   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm) 
     337      ln_3d_uve   = .false.   !  specify whether we are supplying a 3D u,v and e3 field 
     338      ln_read_frq = .false.   !  specify whether we must read frq or not 
     339 
     340   cn_dir      = './'      !  root directory for the ocean data location 
     341   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     342   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     343   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     344   sn_usp      = 'sas_grid_U'            ,       120         , 'uos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     345   sn_vsp      = 'sas_grid_V'            ,       120         , 'vos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     346   sn_tem      = 'sas_grid_T'            ,       120         , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     347   sn_sal      = 'sas_grid_T'            ,       120         , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     348   sn_ssh      = 'sas_grid_T'            ,       120         , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     349   sn_e3t      = 'sas_grid_T'            ,       120         , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
     350   sn_frq      = 'sas_grid_T'            ,       120         , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    '' 
    348351/ 
    349352!----------------------------------------------------------------------- 
    350353&namsbc_iif    !   Ice-IF : use observed ice cover                      (nn_ice = 1) 
    351354!----------------------------------------------------------------------- 
    352 !              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    353 !              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    354    sn_ice      ='ice_cover_clim.nc',     -12.        ,'ice_cover',  .true.     , .true.  , 'yearly'  ,   ''     ,   ''     ,   '' 
    355    ! 
    356    cn_dir      = './'      !  root directory for the location of the runoff files  
     355   cn_dir      = './'      !  root directory for the ice cover data location 
     356   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     357   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     358   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     359   sn_ice      ='ice_cover_clim.nc'      ,        -12.       ,'ice_cover',   .true.    , .true.  , 'yearly'  ,   ''            ,   ''     ,   '' 
    357360/ 
    358361!----------------------------------------------------------------------- 
    359362&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T) 
    360363!----------------------------------------------------------------------- 
    361 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    362 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    363    sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    364  
    365    cn_dir      = './'      !  root directory for the location of the runoff files 
    366    ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration 
    367    ln_qsr_2bd  = .false.   !  2 bands              light penetration 
    368    ln_qsr_bio  = .false.   !  bio-model light penetration 
    369    nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0) 
    370    rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1) 
    371    rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction 
    372    rn_si1      =   23.0    !  2 bands: longest depth of extinction 
    373 / 
    374 !----------------------------------------------------------------------- 
    375 &namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf =T) 
    376 !----------------------------------------------------------------------- 
    377 !              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    378 !              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    379    sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    380    sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    381    sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    382    sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    383    sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    384  
    385    cn_dir      = './'      !  root directory for the location of the runoff files 
    386    ln_rnf_mouth= .true.    !  specific treatment at rivers mouths 
     364   !                       !  type of penetration                        (default: NO selection) 
     365   ln_qsr_rgb  = .false.      !  RGB light penetration (Red-Green-Blue) 
     366   ln_qsr_2bd  = .false.      !  2BD light penetration (two bands) 
     367   ln_qsr_bio  = .false.      !  bio-model light penetration 
     368   !                       !  RGB & 2BD choices: 
     369   rn_abs      =   0.58       !  RGB & 2BD: fraction absorbed in the very near surface 
     370   rn_si0      =   0.35       !  RGB & 2BD: shortess depth of extinction 
     371   nn_chldta   =      0       !  RGB : Chl data (=1) or cst value (=0) 
     372   rn_si1      =   23.0       !  2BD : longest depth of extinction 
     373    
     374   cn_dir      = './'      !  root directory for the chlorophyl data location 
     375   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     376   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     377   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     378   sn_chl      ='chlorophyll'            ,        -1         , 'CHLA'    ,   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     379/ 
     380!----------------------------------------------------------------------- 
     381&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
     382!----------------------------------------------------------------------- 
     383   nn_sstr     =     0     !  add a retroaction term to the surface heat flux (=1) or not (=0) 
     384      rn_dqdt     = -40.      !  magnitude of the retroaction on temperature   [W/m2/K] 
     385   nn_sssr     =     0     !  add a damping term to the surface freshwater flux (=2) 
     386      !                    !  or to SSS only (=1) or no damping term (=0) 
     387      rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
     388      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
     389      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
     390 
     391   cn_dir      = './'      !  root directory for the SST/SSS data location 
     392   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     393   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     394   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     395   sn_sst      = 'sst_data'              ,        24         ,  'sst'    ,    .false.  , .false., 'yearly'  ,    ''            ,    ''    ,     '' 
     396   sn_sss      = 'sss_data'              ,        -1         ,  'sss'    ,    .true.   , .true. , 'yearly'  ,    ''            ,    ''    ,     '' 
     397/ 
     398!----------------------------------------------------------------------- 
     399&namsbc_rnf    !   runoffs                                              (ln_rnf =T) 
     400!----------------------------------------------------------------------- 
     401   ln_rnf_mouth= .false.   !  specific treatment at rivers mouths 
    387402      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T) 
    388403      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 
     
    395410      rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
    396411      nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0) 
    397 / 
    398 !----------------------------------------------------------------------- 
    399 &namsbc_isf    !  Top boundary layer (ISF)                              (nn_isf >0) 
    400 !----------------------------------------------------------------------- 
    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 ! nn_isf == 4 
    404    sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    405 ! nn_isf == 3 
    406    sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    407 ! nn_isf == 2 and 3 
    408    sn_depmax_isf='rnfisf'  ,         -12       ,'sozisfmax',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    409    sn_depmin_isf='rnfisf'  ,         -12       ,'sozisfmin',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    410 ! nn_isf == 2 
    411    sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    412 ! 
    413 ! for all case 
     412 
     413   cn_dir      = './'      !  root directory for the runoff data location 
     414   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     415   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     416   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     417   sn_rnf      = 'runoff_core_monthly'   ,        -1         , 'sorunoff',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     418   sn_cnf      = 'runoff_core_monthly'   ,         0         , 'socoefr0',   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
     419   sn_s_rnf    = 'runoffs'               ,        24         , 'rosaline',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     420   sn_t_rnf    = 'runoffs'               ,        24         , 'rotemper',   .true.    , .true. , 'yearly'  , ''               , ''       , '' 
     421   sn_dep_rnf  = 'runoffs'               ,         0         , 'rodepth' ,   .false.   , .true. , 'yearly'  , ''               , ''       , '' 
     422/ 
     423!----------------------------------------------------------------------- 
     424&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T) 
     425!----------------------------------------------------------------------- 
     426   rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/ 
     427   ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F) 
     428   ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data 
     429 
     430   cn_dir = './'        !  root directory for the Patm data location 
     431   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     432   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     433   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     434   sn_apr      = 'patm'                  ,         -1        ,'somslpre' ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,      '' 
     435/ 
     436!----------------------------------------------------------------------- 
     437&namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)  
     438!-----------------------------------------------------------------------             or set or usr_def_zgr ) 
     439   !                 ! type of top boundary layer  
    414440   nn_isf      = 1         !  ice shelf melting/freezing 
    415                            !  1 = presence of ISF    2 = bg03 parametrisation  
    416                            !  3 = rnf file for isf   4 = ISF fwf specified 
     441                           !  1 = presence of ISF   ; 2 = bg03 parametrisation  
     442                           !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux 
    417443                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    418 ! only for nn_isf = 1 or 2 
    419    rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    420    rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    421 ! only for nn_isf = 1 or 4 
    422    rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    423    !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    424 ! only for nn_isf = 1 
    425    nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    426    !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    427    nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    428    !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    429    !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    430 / 
    431 !----------------------------------------------------------------------- 
    432 &namsbc_iscpl  !   land ice / ocean coupling option                      
    433 !----------------------------------------------------------------------- 
     444      !              !  nn_isf = 1 or 2 cases: 
     445      rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
     446      rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
     447      !              !  nn_isf = 1 or 4 cases: 
     448      rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     449      !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     450      !              ! nn_isf = 1 case 
     451      nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
     452      !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
     453      nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
     454      !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     455      !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
     456 
     457   !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     458   !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     459   !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     460!* nn_isf = 4 case 
     461   sn_fwfisf   = 'rnfisf'    ,         -12       ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     462!* nn_isf = 3 case 
     463   sn_rnfisf   = 'rnfisf'    ,         -12       ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     464!* nn_isf = 2 and 3 cases  
     465   sn_depmax_isf='rnfisf'    ,         -12       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     466   sn_depmin_isf='rnfisf'    ,         -12       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     467!* nn_isf = 2 case 
     468   sn_Leff_isf = 'rnfisf'    ,         -12       ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     469/ 
     470!----------------------------------------------------------------------- 
     471&namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T)  
     472!-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    434473   nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    435474   ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
     
    437476/ 
    438477!----------------------------------------------------------------------- 
    439 &namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T) 
    440 !----------------------------------------------------------------------- 
    441 !              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    442 !              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    443    sn_apr      = 'patm'    ,         -1        ,'somslpre',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,      '' 
    444  
    445    cn_dir      = './'      !  root directory for the location of the bulk files 
    446    rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/ 
    447    ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F) 
    448    ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data 
    449 / 
    450 !----------------------------------------------------------------------- 
    451 &namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T) 
    452 !----------------------------------------------------------------------- 
    453 !              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    454 !              !           !  (if <0  months)  !   name   !   (logical) !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    455    sn_sst      = 'sst_data',        24         ,  'sst'   ,    .false.  , .false., 'yearly'  ,    ''    ,    ''    ,     '' 
    456    sn_sss      = 'sss_data',        -1         ,  'sss'   ,    .true.   , .true. , 'yearly'  ,    ''    ,    ''    ,     '' 
    457    cn_dir      = './'      !  root directory for the location of the runoff files 
    458    ! 
    459    nn_sstr     =     0     !  add a retroaction term in the surface heat flux (=1) or not (=0) 
    460       rn_dqdt     = -40.      !  magnitude of the retroaction on temperature   [W/m2/K] 
    461    nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2) 
    462       !                    !  or to SSS only (=1) or no damping term (=0) 
    463       rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
    464       ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    465       rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
    466 / 
    467 !----------------------------------------------------------------------- 
    468478&namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
    469479!----------------------------------------------------------------------- 
    470 !              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    471 !              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    472    sn_cdg      =  'sdw_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    473    sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    474    sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    475    sn_hsw      =  'sdw_wave' ,        1          , 'hs'         ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    476    sn_wmp      =  'sdw_wave' ,        1          , 'wmp'        ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    477    sn_wfr      =  'sdw_wave' ,        1          , 'wfr'        ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    478    sn_wnum     =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    479    sn_tauwoc   =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    480    sn_tauwx    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    481    sn_tauwy    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     480   cn_dir      = './'      !  root directory for the waves data location 
     481   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     482   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     483   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     484   sn_cdg      =  'sdw_wave'             ,        1          , 'drag_coeff' ,  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     485   sn_usd      =  'sdw_wave'             ,        1          , 'u_sd2d'     ,  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     486   sn_vsd      =  'sdw_wave'             ,        1          , 'v_sd2d'     ,  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     487   sn_hsw      =  'sdw_wave'             ,        1          , 'hs'         ,  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     488   sn_wmp      =  'sdw_wave'             ,        1          , 'wmp'        ,  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     489   sn_wfr      =  'sdw_wave'             ,        1          , 'wfr'        ,  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     490   sn_wnum     =  'sdw_wave'             ,        1          , 'wave_num'   ,  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     491   sn_tauwoc   =  'sdw_wave'             ,        1          , 'wave_stress',  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     492   sn_tauwx    =  'sdw_wave'             ,        1          , 'wave_stress',  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     493   sn_tauwy    =  'sdw_wave'             ,        1          , 'wave_stress',  .true.  , .false., 'daily'   ,  ''             , ''       , '' 
     494/ 
     495!----------------------------------------------------------------------- 
     496&namberg       !   iceberg parameters                                   (default: NO) 
     497!----------------------------------------------------------------------- 
     498   ln_icebergs = .false.      ! activate iceberg floats (force =F with "key_agrif") 
     499   ! 
     500   !                          ! diagnostics: 
     501   ln_bergdia        = .true.       ! Calculate budgets 
     502   nn_verbose_level  = 1            ! Turn on more verbose output if level > 0 
     503   nn_verbose_write  = 15           ! Timesteps between verbose messages 
     504   nn_sample_rate    = 1            ! Timesteps between sampling for trajectory storage 
     505   ! 
     506   !                          ! iceberg setting: 
     507   !                                ! Initial mass required for an iceberg of each class 
     508   rn_initial_mass   = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
     509   !                                ! Proportion of calving mass to apportion to each class 
     510   rn_distribution   = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
     511   !                                ! Ratio between effective and real iceberg mass (non-dim) 
     512   !                                ! i.e. number of icebergs represented at a point 
     513   rn_mass_scaling   = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
     514                                    ! thickness of newly calved bergs (m) 
     515   rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
     516   ! 
     517   rn_rho_bergs            = 850.   ! Density of icebergs 
     518   rn_LoW_ratio            = 1.5    ! Initial ratio L/W for newly calved icebergs 
     519   ln_operator_splitting   = .true. ! Use first order operator splitting for thermodynamics 
     520   rn_bits_erosion_fraction = 0.    ! Fraction of erosion melt flux to divert to bergy bits 
     521   rn_sicn_shift           = 0.     ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
     522   ln_passive_mode         = .false.! iceberg - ocean decoupling 
     523   nn_test_icebergs        =  10    ! Create test icebergs of this class (-1 = no) 
     524   !                                ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
     525   rn_test_box             = 108.0,  116.0, -66.0, -58.0 
     526   rn_speed_limit          = 0.     ! CFL speed limit for a berg 
     527 
     528   cn_dir      = './'      !  root directory for the calving data location 
     529   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     530   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     531   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     532   sn_icb     =  'calving'              ,         -1         ,'calvingmask',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     533/ 
     534 
     535!!====================================================================== 
     536!!               ***  Lateral boundary condition  ***                 !! 
     537!!                                                                    !! 
     538!!   namlbc        lateral momentum boundary condition                  (default: no slip) 
     539!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
     540!!   nam_tide      Tidal forcing                                        (default: NO) 
     541!!   nambdy        Unstructured open boundaries                         (default: NO) 
     542!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy) 
     543!!   nambdy_tide   tidal forcing at open boundaries                     (default: NO) 
     544!!====================================================================== 
    482545! 
    483    cn_dir  = './'  !  root directory for the location of drag coefficient files 
    484 / 
    485 !----------------------------------------------------------------------- 
    486 &namberg       !   iceberg parameters                                   (default: No iceberg) 
    487 !----------------------------------------------------------------------- 
    488    ln_icebergs              = .false.              ! iceberg floats or not 
    489    ln_bergdia               = .true.               ! Calculate budgets 
    490    nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
    491    nn_verbose_write         = 15                   ! Timesteps between verbose messages 
    492    nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
    493                                                    ! Initial mass required for an iceberg of each class 
    494    rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
    495                                                    ! Proportion of calving mass to apportion to each class 
    496    rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    497                                                    ! Ratio between effective and real iceberg mass (non-dim) 
    498                                                    ! i.e. number of icebergs represented at a point 
    499    rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    500                                                    ! thickness of newly calved bergs (m) 
    501    rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
    502    rn_rho_bergs             = 850.                 ! Density of icebergs 
    503    rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
    504    ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
    505    rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
    506    rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
    507    ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
    508    nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
    509                                                    ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
    510    rn_test_box              = 108.0,  116.0, -66.0, -58.0 
    511    rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
    512  
    513 !         ! file name ! frequency (hours) !   variable   ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    514 !         !           !  (if <0  months)  !     name     !   (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    515    sn_icb =  'calving',       -1          , 'calvingmask',  .true.      , .true.  , 'yearly'  , ''       , ''       , '' 
    516  
    517    cn_dir = './' 
    518 / 
    519  
    520 !!====================================================================== 
    521 !!               ***  Lateral boundary condition  ***                 !! 
    522 !!====================================================================== 
    523 !!   namlbc        lateral momentum boundary condition 
    524 !!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
    525 !!   nam_tide      Tidal forcing  
    526 !!   nambdy        Unstructured open boundaries                          
    527 !!   nambdy_dta    Unstructured open boundaries - external data          
    528 !!   nambdy_tide   tidal forcing at open boundaries                      
    529 !!====================================================================== 
    530 ! 
    531 !----------------------------------------------------------------------- 
    532 &namlbc        !   lateral momentum boundary condition 
     546!----------------------------------------------------------------------- 
     547&namlbc        !   lateral momentum boundary condition                  (default: no slip) 
    533548!----------------------------------------------------------------------- 
    534549   !                       !  free slip  !   partial slip  !   no slip   ! strong slip 
     
    542557   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    543558   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
    544    ln_chk_bathy  = .FALSE. ! 
    545 / 
    546 !----------------------------------------------------------------------- 
    547 &nam_tide      !   tide parameters 
    548 !----------------------------------------------------------------------- 
    549    ln_tide       = .false.                ! Activate tides 
     559   ln_chk_bathy  = .false. !  =T  check the parent bathymetry 
     560/ 
     561!----------------------------------------------------------------------- 
     562&nam_tide      !   tide parameters                                      (default: NO) 
     563!----------------------------------------------------------------------- 
     564   ln_tide     = .false.      ! Activate tides 
    550565      ln_tide_pot   = .true.                !  use tidal potential forcing 
    551566         ln_scal_load  = .false.               ! Use scalar approximation for 
     
    559574/ 
    560575!----------------------------------------------------------------------- 
    561 &nambdy        !  unstructured open boundaries                           
    562 !----------------------------------------------------------------------- 
    563     ln_bdy         = .false.              !  Use unstructured open boundaries 
    564     nb_bdy         = 0                    !  number of open boundary sets 
    565     ln_coords_file = .true.               !  =T : read bdy coordinates from file 
    566     cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files 
    567     ln_mask_file   = .false.              !  =T : read mask from file 
    568     cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.) 
    569     cn_dyn2d       = 'none'               ! 
    570     nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state 
    571                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    572                                           !  = 2, use tidal harmonic forcing data from files 
    573                                           !  = 3, use external data AND tidal harmonic forcing 
    574     cn_dyn3d      =  'none'               ! 
    575     nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state 
    576                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    577     cn_tra        =  'none'               ! 
    578     nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state 
    579                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    580     cn_ice_lim      =  'none'             ! 
    581     nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state 
    582                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    583     rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice 
    584     rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           -- 
    585     rn_ice_age      = 30.                 !  lim3 only:      --   age                -- 
    586  
    587     ln_tra_dmp    =.false.                !  open boudaries conditions for tracers 
    588     ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities 
    589     rn_time_dmp   =  1.                   ! Damping time scale in days 
    590     rn_time_dmp_out =  1.                 ! Outflow damping time scale 
    591     nn_rimwidth   = 10                    !  width of the relaxation zone 
    592     ln_vol        = .false.               !  total volume correction (see nn_volctl parameter) 
    593     nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero 
    594     nb_jpk_bdy    = -1                    ! number of levels in the bdy data (set < 0 if consistent with planned run) 
    595 / 
    596 !----------------------------------------------------------------------- 
    597 &nambdy_dta    !  open boundaries - external data                        
    598 !----------------------------------------------------------------------- 
    599 !              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    600 !              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    601    bn_ssh      = 'amm12_bdyT_u2d',         24        , 'sossheig',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    602    bn_u2d      = 'amm12_bdyU_u2d',         24        , 'vobtcrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    603    bn_v2d      = 'amm12_bdyV_u2d',         24        , 'vobtcrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    604    bn_u3d      = 'amm12_bdyU_u3d',         24        , 'vozocrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    605    bn_v3d      = 'amm12_bdyV_u3d',         24        , 'vomecrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    606    bn_tem      = 'amm12_bdyT_tra',         24        , 'votemper',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    607    bn_sal      = 'amm12_bdyT_tra',         24        , 'vosaline',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    608 ! for lim3 
    609 !   bn_a_i     = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    610 !   bn_h_i     = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    611 !   bn_h_s     = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    612  
    613    cn_dir      = 'bdydta/' !  root directory for the location of the bulk files 
    614    ln_full_vel = .false.   !   
    615 / 
    616 !----------------------------------------------------------------------- 
    617 &nambdy_tide   !  tidal forcing at open boundaries 
     576&nambdy        !  unstructured open boundaries                          (default: NO) 
     577!----------------------------------------------------------------------- 
     578   ln_bdy         = .false.   !  Use unstructured open boundaries 
     579   nb_bdy         = 0         !  number of open boundary sets 
     580   ln_coords_file = .true.    !  =T : read bdy coordinates from file 
     581      cn_coords_file = 'coordinates.bdy.nc'  !  bdy coordinates files 
     582   ln_mask_file   = .false.   !  =T : read mask from file 
     583      cn_mask_file= ''        !  name of mask file (if ln_mask_file=.TRUE.) 
     584   cn_dyn2d    = 'none'       ! 
     585   nn_dyn2d_dta   =  0        !  = 0, bdy data are equal to the initial state 
     586      !                       !  = 1, bdy data are read in 'bdydata   .nc' files 
     587      !                       !  = 2, use tidal harmonic forcing data from files 
     588      !                       !  = 3, use external data AND tidal harmonic forcing 
     589   cn_dyn3d      =  'none'    ! 
     590   nn_dyn3d_dta  =  0         !  = 0, bdy data are equal to the initial state 
     591   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     592   cn_tra        =  'none'    ! 
     593   nn_tra_dta    =  0         !  = 0, bdy data are equal to the initial state 
     594   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     595   cn_ice_lim    =  'none'    ! 
     596   nn_ice_lim_dta=  0         !  = 0, bdy data are equal to the initial state 
     597   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     598   rn_ice_tem    = 270.       !  lim3 only: arbitrary temperature of incoming sea ice 
     599   rn_ice_sal    = 10.        !  lim3 only:      --   salinity           -- 
     600   rn_ice_age    = 30.        !  lim3 only:      --   age                -- 
     601   ! 
     602   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers 
     603   ln_dyn3d_dmp  =.false.     !  open boundary condition for baroclinic velocities 
     604   rn_time_dmp   =  1.        !  Damping time scale in days 
     605   rn_time_dmp_out= 1.        !  Outflow damping time scale 
     606   nn_rimwidth   = 10         !  width of the relaxation zone 
     607   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter) 
     608   nn_volctl     =  1         !  = 0, the total water flux across open boundaries is zero 
     609   nb_jpk_bdy    = -1         ! number of levels in the bdy data (set < 0 if consistent with planned run) 
     610/ 
     611!----------------------------------------------------------------------- 
     612&nambdy_dta    !  open boundaries - external data                       (see nam_bdy) 
     613!----------------------------------------------------------------------- 
     614   ln_full_vel = .false.      !  ??? 
     615 
     616   cn_dir      = 'bdydta/'    !  root directory for the BDY data location 
     617   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     618   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     619   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     620   bn_ssh      = 'amm12_bdyT_u2d'        ,         24        , 'sossheig',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     621   bn_u2d      = 'amm12_bdyU_u2d'        ,         24        , 'vobtcrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     622   bn_v2d      = 'amm12_bdyV_u2d'        ,         24        , 'vobtcrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     623   bn_u3d      = 'amm12_bdyU_u3d'        ,         24        , 'vozocrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     624   bn_v3d      = 'amm12_bdyV_u3d'        ,         24        , 'vomecrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     625   bn_tem      = 'amm12_bdyT_tra'        ,         24        , 'votemper',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     626   bn_sal      = 'amm12_bdyT_tra'        ,         24        , 'vosaline',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     627!* for lim3 
     628!   bn_a_i     = 'amm12_bdyT_ice'        ,         24        , 'ileadfra',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     629!   bn_h_i     = 'amm12_bdyT_ice'        ,         24        , 'iicethic',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     630!   bn_h_s     = 'amm12_bdyT_ice'        ,         24        , 'isnowthi',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     '' 
     631/ 
     632!----------------------------------------------------------------------- 
     633&nambdy_tide   !  tidal forcing at open boundaries                      (default: NO) 
    618634!----------------------------------------------------------------------- 
    619635   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
     
    624640!!====================================================================== 
    625641!!                ***  Top/Bottom boundary condition  ***             !! 
    626 !!====================================================================== 
     642!!                                                                    !! 
    627643!!   namdrg        top/bottom drag coefficient                          (default: NO selection) 
    628 !!   namdrg_top    top    friction                                      (ln_isfcav=T) 
    629 !!   namdrg_bot    bottom friction                                       
     644!!   namdrg_top    top    friction                                      (ln_NONE=F & ln_isfcav=T) 
     645!!   namdrg_bot    bottom friction                                      (ln_NONE=F) 
    630646!!   nambbc        bottom temperature boundary condition                (default: NO) 
    631647!!   nambbl        bottom boundary layer scheme                         (default: NO) 
     
    643659/ 
    644660!----------------------------------------------------------------------- 
    645 &namdrg_top    !   TOP friction                                         (ln_isfcav=T) 
     661&namdrg_top    !   TOP friction                                         (ln_NONE=F & ln_isfcav=T) 
    646662!----------------------------------------------------------------------- 
    647663   rn_Cd0     =  1.e-3     !  drag coefficient [-] 
     
    654670/ 
    655671!----------------------------------------------------------------------- 
    656 &namdrg_bot    !   BOTTOM friction                                   
     672&namdrg_bot    !   BOTTOM friction                                      (ln_NONE=F) 
    657673!----------------------------------------------------------------------- 
    658674   rn_Cd0     =  1.e-3    !  drag coefficient [-] 
     
    667683&nambbc        !   bottom temperature boundary condition                (default: NO) 
    668684!----------------------------------------------------------------------- 
    669 !              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    670 !              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    671    sn_qgh      ='geothermal_heating.nc',  -12.       , 'heatflow',   .false.   , .true.  , 'yearly'  ,   ''     ,   ''     ,   '' 
    672    ! 
    673685   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    674       nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux 
    675       !                       !     = 1 constant flux 
    676       !                       !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    677       rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    678       cn_dir      = './'      !  root directory for the location of the runoff files  
     686      nn_geoflx   =    2      !  geothermal heat flux: = 1 constant flux 
     687      !                       !                        = 2 read variable flux [mW/m2] 
     688      rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux       [mW/m2] 
     689 
     690   cn_dir      = './'      !  root directory for the geothermal data location 
     691   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     692   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     693   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     694   sn_qgh      ='geothermal_heating.nc'  ,       -12.        , 'heatflow',   .false.   , .true. , 'yearly'  ,   ''             ,   ''     ,   '' 
    679695/ 
    680696!----------------------------------------------------------------------- 
     
    689705 
    690706!!====================================================================== 
    691 !!                        Tracer (T & S) namelists 
    692 !!====================================================================== 
    693 !!   nameos           equation of state 
    694 !!   namtra_adv       advection scheme 
    695 !!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    696 !!   namtra_ldf       lateral diffusion scheme 
    697 !!   namtra_ldfeiv    eddy induced velocity param. 
    698 !!   namtra_dmp       T & S newtonian damping 
     707!!                        Tracer (T & S) namelists                    !! 
     708!!                                                                    !! 
     709!!   nameos           equation of state                                 (default: NO selection) 
     710!!   namtra_adv       advection scheme                                  (default: NO selection) 
     711!!   namtra_ldf       lateral diffusion scheme                          (default: NO selection) 
     712!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.)       (default: NO) 
     713!!   namtra_eiv       eddy induced velocity param.                      (default: NO) 
     714!!   namtra_dmp       T & S newtonian damping                           (default: NO) 
    699715!!====================================================================== 
    700716! 
    701717!----------------------------------------------------------------------- 
    702 &nameos        !   ocean Equation Of Seawater                           (default: NO) 
     718&nameos        !   ocean Equation Of Seawater                           (default: NO selection) 
    703719!----------------------------------------------------------------------- 
    704720   ln_teos10   = .false.         !  = Use TEOS-10 
     
    733749/ 
    734750!----------------------------------------------------------------------- 
    735 &namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) 
     751&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection) 
     752!----------------------------------------------------------------------- 
     753   !                       !  Operator type: 
     754   ln_traldf_NONE  = .false.   !  No explicit diffusion 
     755   ln_traldf_lap   = .false.   !    laplacian operator 
     756   ln_traldf_blp   = .false.   !  bilaplacian operator 
     757   ! 
     758   !                       !  Direction of action: 
     759   ln_traldf_lev   = .false.   !  iso-level 
     760   ln_traldf_hor   = .false.   !  horizontal  (geopotential) 
     761   ln_traldf_iso   = .false.   !  iso-neutral (standard operator) 
     762   ln_traldf_triad = .false.   !  iso-neutral (triad    operator) 
     763   ! 
     764   !                       !  iso-neutral options:         
     765   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators) 
     766   rn_slpmax       =  0.01     !  slope limit                           (both operators) 
     767   ln_triad_iso    = .false.   !  pure horizontal mixing in ML              (triad only) 
     768   rn_sw_triad     = 1         !  =1 switching triad ; =0 all 4 triads used (triad only) 
     769   ln_botmix_triad = .false.   !  lateral mixing on bottom                  (triad only) 
     770   ! 
     771   !                       !  Coefficients: 
     772   nn_aht_ijk_t    = 0         !  space/time variation of eddy coefficient: 
     773      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     774      !                             !   =  0           constant  
     775      !                             !   = 10 F(k)      =ldf_c1d  
     776      !                             !   = 20 F(i,j)    =ldf_c2d  
     777      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     778      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     779      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     780      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)  
     781      !                             !                           or   = 1/12 Ud*Ld^3 (blp case) 
     782      rn_Ud        = 0.01           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
     783      rn_Ld        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10) 
     784/ 
     785!----------------------------------------------------------------------- 
     786&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper)       (default: NO) 
    736787!----------------------------------------------------------------------- 
    737788   ln_mle      = .false.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     
    746797/ 
    747798!----------------------------------------------------------------------- 
    748 &namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection) 
    749 !----------------------------------------------------------------------- 
    750    !                       !  Operator type: 
    751    ln_traldf_NONE  =  .false.  !  No explicit diffusion 
    752    ln_traldf_lap   =  .false.  !    laplacian operator 
    753    ln_traldf_blp   =  .false.  !  bilaplacian operator 
    754    ! 
    755    !                       !  Direction of action: 
    756    ln_traldf_lev   =  .false.  !  iso-level 
    757    ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    758    ln_traldf_iso   =  .false.  !  iso-neutral (standard operator) 
    759    ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
    760    ! 
    761    !                       !  iso-neutral options:         
    762    ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
    763    rn_slpmax       =   0.01    !  slope limit                      (both operators) 
    764    ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
    765    rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
    766    ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
    767    ! 
    768    !                       !  Coefficients: 
    769    nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
    770    !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    771    !                                !   =  0           constant  
    772    !                                !   = 10 F(k)      =ldf_c1d  
    773    !                                !   = 20 F(i,j)    =ldf_c2d  
    774    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    775    !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    776    !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    777    rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
    778    rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    779 / 
    780 !----------------------------------------------------------------------- 
    781 &namtra_ldfeiv !   eddy induced velocity param.                         (default: NO) 
    782 !----------------------------------------------------------------------- 
    783    ln_ldfeiv     = .false. ! use eddy induced velocity parameterization 
    784       rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
    785       nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     799&namtra_eiv    !   eddy induced velocity param.                         (default: NO) 
     800!----------------------------------------------------------------------- 
     801   ln_ldfeiv   = .false.   ! use eddy induced velocity parameterization 
     802      ! 
     803      !                        !  Coefficients: 
     804      nn_aei_ijk_t    = 0           !  space/time variation of eddy coefficient: 
    786805      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    787806      !                             !   =  0           constant  
     
    789808      !                             !   = 20 F(i,j)    =ldf_c2d  
    790809      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    791       !                             !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    792       ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities 
     810      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     811      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le  
     812      rn_Ue        = 0.02           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
     813      rn_Le        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10) 
     814      ! 
     815      ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
    793816/ 
    794817!----------------------------------------------------------------------- 
    795818&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO) 
    796819!----------------------------------------------------------------------- 
    797    ln_tradmp   =  .true.   !  add a damping term 
     820   ln_tradmp   =  .false.  !  add a damping term (using resto.nc coef.) 
    798821      nn_zdmp  =    0         !  vertical shape =0    damping throughout the water column 
    799822      !                       !                 =1 no damping in the mixing layer (kz  criteria) 
     
    804827!!====================================================================== 
    805828!!                      ***  Dynamics namelists  ***                  !! 
    806 !!====================================================================== 
    807 !!   nam_vvl       vertical coordinate options 
    808 !!   namdyn_adv    formulation of the momentum advection 
    809 !!   namdyn_vor    advection scheme 
    810 !!   namdyn_hpg    hydrostatic pressure gradient 
    811 !!   namdyn_spg    surface pressure gradient 
    812 !!   namdyn_ldf    lateral diffusion scheme 
     829!!                                                                    !! 
     830!!   nam_vvl       vertical coordinate options                          (default: z-star) 
     831!!   namdyn_adv    formulation of the momentum advection                (default: NO selection) 
     832!!   namdyn_vor    advection scheme                                     (default: NO selection) 
     833!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection) 
     834!!   namdyn_spg    surface pressure gradient                            (default: NO selection) 
     835!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection) 
     836!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only) 
    813837!!====================================================================== 
    814838! 
     
    837861/ 
    838862!----------------------------------------------------------------------- 
    839 &namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO) 
    840 !----------------------------------------------------------------------- 
    841    ln_dynvor_ene = .false. !  enstrophy conserving scheme 
    842    ln_dynvor_ens = .false. !  energy conserving scheme 
     863&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection) 
     864!----------------------------------------------------------------------- 
     865   ln_dynvor_ene = .false. !  energy    conserving scheme 
     866   ln_dynvor_ens = .false. !  enstrophy conserving scheme 
    843867   ln_dynvor_mix = .false. !  mixed scheme 
    844868   ln_dynvor_een = .false. !  energy & enstrophy scheme 
    845       nn_een_e3f = 1          ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     869      nn_een_e3f = 1          ! =0   e3f = mean masked e3t divided by 4 
     870      !                       ! =1   e3f = mean masked e3t divided by the sum of mask 
    846871   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes)  ! PLEASE DO NOT ACTIVATE 
    847872/ 
     
    857882/ 
    858883!----------------------------------------------------------------------- 
    859 &namdyn_spg    !   surface pressure gradient                            (default: NO) 
     884&namdyn_spg    !   surface pressure gradient                            (default: NO selection) 
    860885!----------------------------------------------------------------------- 
    861886   ln_dynspg_exp  = .false.   ! explicit free surface 
     
    875900!----------------------------------------------------------------------- 
    876901   !                       !  Type of the operator : 
    877    ln_dynldf_NONE=  .false.    !  No operator (i.e. no explicit diffusion) 
    878    ln_dynldf_lap =  .false.    !    laplacian operator 
    879    ln_dynldf_blp =  .false.    !  bilaplacian operator 
     902   ln_dynldf_NONE= .false.     !  No operator (i.e. no explicit diffusion) 
     903   ln_dynldf_lap = .false.     !    laplacian operator 
     904   ln_dynldf_blp = .false.     !  bilaplacian operator 
    880905   !                       !  Direction of action  : 
    881    ln_dynldf_lev =  .false.    !  iso-level 
    882    ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
    883    ln_dynldf_iso =  .false.    !  iso-neutral 
     906   ln_dynldf_lev = .false.     !  iso-level 
     907   ln_dynldf_hor = .false.     !  horizontal (geopotential) 
     908   ln_dynldf_iso = .false.     !  iso-neutral (lap only) 
    884909   !                       !  Coefficient 
    885    nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
    886    !                                !  =-30  read in eddy_viscosity_3D.nc file 
    887    !                                !  =-20  read in eddy_viscosity_2D.nc file 
    888    !                                !  =  0  constant  
    889    !                                !  = 10  F(k)=c1d 
    890    !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
    891    !                                !  = 30  F(i,j,k)=c2d*c1d 
    892    !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
    893    !                                !  = 32  F(i,j,k)=F(local gridscale and deformation rate) 
    894    ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    895    rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
    896    rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
    897    rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    898    !                       !  Smagorinsky settings (nn_ahm_ijk_t  = 32) : 
    899    rn_csmc       = 3.5         !  Smagorinsky constant of proportionality 
    900    rn_minfac     = 1.0         !  multiplier of theorectical lower limit 
    901    rn_maxfac     = 1.0         !  multiplier of theorectical upper limit 
    902 / 
    903 ! 
     910   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coefficient : 
     911      !                             !  =-30  read in eddy_viscosity_3D.nc file 
     912      !                             !  =-20  read in eddy_viscosity_2D.nc file 
     913      !                             !  =  0  constant  
     914      !                             !  = 10  F(k)=c1d 
     915      !                             !  = 20  F(i,j)=F(grid spacing)=c2d 
     916      !                             !  = 30  F(i,j,k)=c2d*c1d 
     917      !                             !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     918      !                             !  = 32  F(i,j,k)=F(local gridscale and deformation rate) 
     919      !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case)  
     920      !                             !                            or  = 1/12 Uv*Lv^3 (blp case) 
     921      rn_Uv      = 0.1              !  lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) 
     922      rn_Lv      = 10.e+3           !  lateral viscous length   [m]   (nn_ahm_ijk_t= 0, 10) 
     923      !                       !  Smagorinsky settings  (nn_ahm_ijk_t= 32) : 
     924      rn_csmc       = 3.5         !  Smagorinsky constant of proportionality 
     925      rn_minfac     = 1.0         !  multiplier of theorectical lower limit 
     926      rn_maxfac     = 1.0         !  multiplier of theorectical upper limit 
     927      !                       !  iso-neutral laplacian operator (ln_dynldf_iso=T) : 
     928      rn_ahm_b      = 0.0         !  background eddy viscosity  [m2/s] 
     929/ 
     930!----------------------------------------------------------------------- 
     931&namdta_dyn    !   offline ocean input files                            (OFF_SRC only) 
     932!----------------------------------------------------------------------- 
     933   ln_dynrnf       =  .false.    !  runoffs option enabled (T) or not (F) 
     934   ln_dynrnf_depth =  .false.    !  runoffs is spread in vertical (T) or not (F) 
     935!   fwbcorr        = 3.786e-06   !  annual global mean of empmr for ssh correction 
     936 
     937   cn_dir      = './'      !  root directory for the ocean data location 
     938   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     939   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
     940   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
     941   sn_tem      = 'dyna_grid_T'           ,       120         , 'votemper'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     942   sn_sal      = 'dyna_grid_T'           ,       120         , 'vosaline'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     943   sn_mld      = 'dyna_grid_T'           ,       120         , 'somixhgt'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     944   sn_emp      = 'dyna_grid_T'           ,       120         , 'sowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     945   sn_fmf      = 'dyna_grid_T'           ,       120         , 'iowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     946   sn_ice      = 'dyna_grid_T'           ,       120         , 'soicecov'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     947   sn_qsr      = 'dyna_grid_T'           ,       120         , 'soshfldo'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     948   sn_wnd      = 'dyna_grid_T'           ,       120         , 'sowindsp'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     949   sn_uwd      = 'dyna_grid_U'           ,       120         , 'uocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     950   sn_vwd      = 'dyna_grid_V'           ,       120         , 'vocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     951   sn_wwd      = 'dyna_grid_W'           ,       120         , 'wocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     952   sn_avt      = 'dyna_grid_W'           ,       120         , 'voddmavs'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     953   sn_ubl      = 'dyna_grid_U'           ,       120         , 'sobblcox'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     954   sn_vbl      = 'dyna_grid_V'           ,       120         , 'sobblcoy'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , '' 
     955/ 
     956 
    904957!!====================================================================== 
    905958!!                     vertical physics namelists                     !! 
    906 !!====================================================================== 
    907 !!    namzdf        vertical physics 
     959!!                                                                    !! 
     960!!    namzdf        vertical physics manager                            (default: NO selection) 
    908961!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T) 
    909962!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T) 
    910963!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T) 
     964!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T) 
    911965!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T) 
    912966!!====================================================================== 
    913967! 
    914968!----------------------------------------------------------------------- 
    915 &namzdf        !   vertical physics                                     (default: NO selection) 
     969&namzdf        !   vertical physics manager                             (default: NO selection) 
    916970!----------------------------------------------------------------------- 
    917971   !                       ! type of vertical closure (required) 
     
    10311085/ 
    10321086!!====================================================================== 
    1033 !!                  ***  Miscellaneous namelists  *** 
    1034 !!====================================================================== 
     1087!!                  ***  Diagnostics namelists  ***                   !! 
     1088!!                                                                    !! 
     1089!!   namtrd       dynamics and/or tracer trends                         (default NO) 
     1090!!   namptr       Poleward Transport Diagnostics                        (default NO) 
     1091!!   namhsb       Heat and salt budgets                                 (default NO) 
     1092!!   namdiu       Cool skin and warm layer models                       (default NO) 
     1093!!   namdiu       Cool skin and warm layer models                       (default NO) 
     1094!!   namflo       float parameters                                      ("key_float") 
     1095!!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
     1096!!   namdct       transports through some sections                      ("key_diadct") 
     1097!!   nam_diatmb   Top Middle Bottom Output                              (default NO) 
     1098!!   nam_dia25h   25h Mean Output                                       (default NO) 
     1099!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
     1100!!====================================================================== 
     1101! 
     1102!----------------------------------------------------------------------- 
     1103&namtrd        !   trend diagnostics                                    (default NO) 
     1104!----------------------------------------------------------------------- 
     1105   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
     1106   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
     1107   ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
     1108   ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet) 
     1109   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
     1110   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
     1111   ln_tra_trd  = .false.   ! (T) 3D tracer trend output 
     1112   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
     1113   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
     1114/ 
     1115!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     1116!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     1117!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
     1118!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
     1119!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics 
     1120!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
     1121!!gm 
     1122!----------------------------------------------------------------------- 
     1123&namptr        !   Poleward Transport Diagnostic                        (default NO) 
     1124!----------------------------------------------------------------------- 
     1125   ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F) 
     1126   ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not 
     1127/ 
     1128!----------------------------------------------------------------------- 
     1129&namhsb        !  Heat and salt budgets                                 (default NO) 
     1130!----------------------------------------------------------------------- 
     1131   ln_diahsb   = .false.   !  check the heat and salt budgets (T) or not (F) 
     1132/ 
     1133!----------------------------------------------------------------------- 
     1134&namdiu        !   Cool skin and warm layer models                      (default NO) 
     1135!----------------------------------------------------------------------- 
     1136   ln_diurnal      = .false.   ! 
     1137   ln_diurnal_only = .false.   ! 
     1138/ 
     1139!----------------------------------------------------------------------- 
     1140&namflo        !   float parameters                                     ("key_float") 
     1141!----------------------------------------------------------------------- 
     1142   jpnfl       = 1         !  total number of floats during the run 
     1143   jpnnewflo   = 0         !  number of floats for the restart 
     1144   ln_rstflo   = .false.   !  float restart (T) or not (F) 
     1145   nn_writefl  =      75   !  frequency of writing in float output file 
     1146   nn_stockfl  =    5475   !  frequency of creation of the float restart file 
     1147   ln_argo     = .false.   !  Argo type floats (stay at the surface each 10 days) 
     1148   ln_flork4   = .false.   !  trajectories computed with a 4th order Runge-Kutta (T) 
     1149   !                       !  or computed with Blanke' scheme (F) 
     1150   ln_ariane   = .true.    !  Input with Ariane tool convention(T) 
     1151   ln_flo_ascii= .true.    !  Output with Ariane tool netcdf convention(F) or ascii file (T) 
     1152/ 
     1153!----------------------------------------------------------------------- 
     1154&nam_diaharm   !   Harmonic analysis of tidal constituents              ("key_diaharm") 
     1155!----------------------------------------------------------------------- 
     1156    nit000_han = 1         ! First time step used for harmonic analysis 
     1157    nitend_han = 75        ! Last time step used for harmonic analysis 
     1158    nstep_han  = 15        ! Time step frequency for harmonic analysis 
     1159    tname(1)   = 'M2'      ! Name of tidal constituents 
     1160    tname(2)   = 'K1' 
     1161/ 
     1162!----------------------------------------------------------------------- 
     1163&namdct        ! transports through some sections                       ("key_diadct") 
     1164!----------------------------------------------------------------------- 
     1165    nn_dct     = 15        !  time step frequency for transports computing 
     1166    nn_dctwri  = 15        !  time step frequency for transports writing 
     1167    nn_secdebug= 112       !      0 : no section to debug 
     1168    !                      !     -1 : debug all section 
     1169    !                      !  0 < n : debug section number n 
     1170/ 
     1171!----------------------------------------------------------------------- 
     1172&nam_diatmb    !  Top Middle Bottom Output                              (default NO) 
     1173!----------------------------------------------------------------------- 
     1174   ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
     1175/ 
     1176!----------------------------------------------------------------------- 
     1177&nam_dia25h    !  25h Mean Output                                       (default NO) 
     1178!----------------------------------------------------------------------- 
     1179   ln_dia25h   = .false.   ! Choose 25h mean output or not 
     1180/ 
     1181!----------------------------------------------------------------------- 
     1182&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
     1183!----------------------------------------------------------------------- 
     1184   nn_nchunks_i=   4       !  number of chunks in i-dimension 
     1185   nn_nchunks_j=   4       !  number of chunks in j-dimension 
     1186   nn_nchunks_k=   31      !  number of chunks in k-dimension 
     1187   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
     1188   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs 
     1189   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
     1190   !                       !  (F) ignore chunking information and produce netcdf3-compatible files 
     1191/ 
     1192 
     1193!!====================================================================== 
     1194!!                  ***  Miscellaneous namelists  ***                 !! 
     1195!!                                                                    !! 
    10351196!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi") 
    1036 !!   namctl            Control prints  
    1037 !!   namsto            Stochastic parametrization of EOS 
     1197!!   namctl            Control prints                                   (default NO) 
     1198!!   namsto            Stochastic parametrization of EOS                (default NO) 
    10381199!!====================================================================== 
    10391200! 
     
    10501211/ 
    10511212!----------------------------------------------------------------------- 
    1052 &namctl        !   Control prints  
     1213&namctl        !   Control prints                                       (default: NO) 
    10531214!----------------------------------------------------------------------- 
    10541215   ln_ctl      = .false.   !  trends control print (expensive!) 
     
    10811242 
    10821243!!====================================================================== 
    1083 !!                  ***  Diagnostics namelists  *** 
    1084 !!====================================================================== 
    1085 !!   namtrd       dynamics and/or tracer trends                         (default F) 
    1086 !!   namptr       Poleward Transport Diagnostics                        (default F) 
    1087 !!   namhsb       Heat and salt budgets                                 (default F) 
    1088 !!   namdiu       Cool skin and warm layer models                       (default F) 
    1089 !!   namdiu       Cool skin and warm layer models                       (default F) 
    1090 !!   namflo       float parameters                                      ("key_float") 
    1091 !!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
    1092 !!   namdct       transports through some sections                      ("key_diadct") 
    1093 !!   nam_diatmb   Top Middle Bottom Output                              (default F) 
    1094 !!   nam_dia25h   25h Mean Output                                       (default F) 
    1095 !!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
     1244!!               ***  Observation & Assimilation  ***                 !! 
     1245!!                                                                    !! 
     1246!!   namobs       observation and model comparison                      (default: NO) 
     1247!!   nam_asminc   assimilation increments                               ('key_asminc') 
    10961248!!====================================================================== 
    10971249! 
    10981250!----------------------------------------------------------------------- 
    1099 &namtrd        !   trend diagnostics                                    (default F) 
    1100 !----------------------------------------------------------------------- 
    1101    ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
    1102    ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
    1103    ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
    1104    ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet) 
    1105    ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
    1106    ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
    1107    ln_tra_trd  = .false.   ! (T) 3D tracer trend output 
    1108    ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
    1109    nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
    1110 / 
    1111 !!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
    1112 !!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
    1113 !!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
    1114 !!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
    1115 !!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics 
    1116 !!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    1117 !!gm 
    1118 !----------------------------------------------------------------------- 
    1119 &namptr        !   Poleward Transport Diagnostic                        (default F) 
    1120 !----------------------------------------------------------------------- 
    1121    ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F) 
    1122    ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not 
    1123 / 
    1124 !----------------------------------------------------------------------- 
    1125 &namhsb        !  Heat and salt budgets                                 (default F) 
    1126 !----------------------------------------------------------------------- 
    1127    ln_diahsb   = .false.   !  check the heat and salt budgets (T) or not (F) 
    1128 / 
    1129 !----------------------------------------------------------------------- 
    1130 &namdiu        !   Cool skin and warm layer models                      (default F) 
    1131 !----------------------------------------------------------------------- 
    1132    ln_diurnal      = .false.   ! 
    1133    ln_diurnal_only = .false.   ! 
    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") 
    1151 !----------------------------------------------------------------------- 
    1152     nit000_han = 1         ! First time step used for harmonic analysis 
    1153     nitend_han = 75        ! Last time step used for harmonic analysis 
    1154     nstep_han  = 15        ! Time step frequency for harmonic analysis 
    1155     tname(1)   = 'M2'      ! Name of tidal constituents 
    1156     tname(2)   = 'K1' 
    1157 / 
    1158 !----------------------------------------------------------------------- 
    1159 &namdct        ! transports through some sections                       ("key_diadct") 
    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 &nam_diatmb    !  Top Middle Bottom Output                              (default F) 
    1169 !----------------------------------------------------------------------- 
    1170    ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
    1171 / 
    1172 !----------------------------------------------------------------------- 
    1173 &nam_dia25h    !  25h Mean Output                                       (default F) 
    1174 !----------------------------------------------------------------------- 
    1175    ln_dia25h   = .false.   ! Choose 25h mean output or not 
    1176 / 
    1177 !----------------------------------------------------------------------- 
    1178 &namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
    1179 !----------------------------------------------------------------------- 
    1180    nn_nchunks_i=   4       !  number of chunks in i-dimension 
    1181    nn_nchunks_j=   4       !  number of chunks in j-dimension 
    1182    nn_nchunks_k=   31      !  number of chunks in k-dimension 
    1183    !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
    1184    !                       !  is optimal for postprocessing which works exclusively with horizontal slabs 
    1185    ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
    1186    !                       !  (F) ignore chunking information and produce netcdf3-compatible files 
    1187 / 
    1188  
    1189 !!====================================================================== 
    1190 !!               ***  Observation & Assimilation  *** 
    1191 !!====================================================================== 
    1192 !!   namobs       observation and model comparison 
    1193 !!   nam_asminc   assimilation increments                               ('key_asminc') 
    1194 !!====================================================================== 
    1195 ! 
    1196 !----------------------------------------------------------------------- 
    1197 &namobs        !  observation usage switch 
     1251&namobs        !  observation usage switch                              (default: NO) 
    11981252!----------------------------------------------------------------------- 
    11991253   ln_diaobs   = .false.             ! Logical switch for the observation operator 
     1254   ! 
    12001255   ln_t3d      = .false.             ! Logical switch for T profile observations 
    12011256   ln_s3d      = .false.             ! Logical switch for S profile observations 
     
    12461301   nn_2dint_sss = 0                  ! Horizontal interpolation method for SSS 
    12471302   nn_2dint_sic = 0                  ! Horizontal interpolation method for SIC 
    1248    nn_msshc    = 0                   ! MSSH correction scheme 
     1303   nn_msshc     = 0                  ! MSSH correction scheme 
    12491304   nn_profdavtypes = -1              ! Profile daily average types - array 
    12501305/ 
     
    12671322    nn_divdmp  = 0         !  Number of iterations of divergence damping operator 
    12681323/ 
    1269 !----------------------------------------------------------------------- 
    1270 &namdta_dyn    !   offline dynamics read in files                       ("key_offline") 
    1271 !----------------------------------------------------------------------- 
    1272 !          !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    1273 !          !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    1274    sn_tem  = 'dyna_grid_T' ,    120            , 'votemper' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1275    sn_sal  = 'dyna_grid_T' ,    120            , 'vosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1276    sn_mld  = 'dyna_grid_T' ,    120            , 'somixhgt' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1277    sn_emp  = 'dyna_grid_T' ,    120            , 'sowaflup' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1278    sn_fmf  = 'dyna_grid_T' ,    120            , 'iowaflup' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1279    sn_ice  = 'dyna_grid_T' ,    120            , 'soicecov' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1280    sn_qsr  = 'dyna_grid_T' ,    120            , 'soshfldo' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1281    sn_wnd  = 'dyna_grid_T' ,    120            , 'sowindsp' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1282    sn_uwd  = 'dyna_grid_U' ,    120            , 'uocetr_eff' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1283    sn_vwd  = 'dyna_grid_V' ,    120            , 'vocetr_eff' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1284    sn_wwd  = 'dyna_grid_W' ,    120            , 'wocetr_eff' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1285    sn_avt  = 'dyna_grid_W' ,    120            , 'voddmavs' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1286    sn_ubl  = 'dyna_grid_U' ,    120            , 'sobblcox' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1287    sn_vbl  = 'dyna_grid_V' ,    120            , 'sobblcoy' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    1288 ! 
    1289    cn_dir          = './'       !  root directory for the location of the dynamical files 
    1290    ln_dynrnf       =  .false.   !  runoffs option enabled (T) or not (F) 
    1291    ln_dynrnf_depth =  .false.   ! runoffs is spread in vertical (T) or not (F) 
    1292 !   fwbcorr      = 3.786e-06    ! annual global mean of empmr for ssh correction 
    1293 / 
Note: See TracChangeset for help on using the changeset viewer.