Changeset 12414


Ignore:
Timestamp:
2020-02-19T21:13:15+01:00 (10 months ago)
Author:
smueller
Message:

Reintegration of 2019 development branch /utils/tools_MERGE_2019 into the tools directory (/utils/tools)

Location:
utils/tools
Files:
48 deleted
5 edited
51 copied

Legend:

Unmodified
Added
Removed
  • utils/tools/DOMAINcfg/namelist_cfg

    r9051 r12414  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref 
     2!! NEMO/OCE :   Configuration namelist_cfg used to overwrite defaults value defined in namelist_ref 
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    4 ! 
     4!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namzgr, namzgr_sco ) 
     5!!              2 - diagnostics      (namnc4) 
     6!!              3 - miscellaneous    (nammpp, namctl) 
     7!! 
     8!! namelist skeleton : egrep -E '(^/ *$|^! *$|^ *$|&nam.*|!---.*|!! .*|!!==.*|!!>>>.*)' namelist_ref > namelist_skl 
     9!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    510!----------------------------------------------------------------------- 
    611&namrun        !   parameters of the run 
    712!----------------------------------------------------------------------- 
    8    nn_no       =       0   !  job number (no more used...) 
    9    cn_exp      =  "domaincfg"  !  experience name 
    10    nn_it000    =       1   !  first time step 
    11    nn_itend    =      75   !  last  time step (std 5475) 
     13/ 
     14!----------------------------------------------------------------------- 
     15&namdom        !   space and time domain (bathymetry, mesh, timestep) 
     16!----------------------------------------------------------------------- 
    1217/ 
    1318!----------------------------------------------------------------------- 
    1419&namcfg        !   parameters of the configuration 
    1520!----------------------------------------------------------------------- 
    16    ! 
    17    ln_e3_dep   = .true.    ! =T : e3=dk[depth] in discret sens.  
    18    !                       !      ===>>> will become the only possibility in v4.0 
    19    !                       ! =F : e3 analytical derivative of depth function 
    20    !                       !      only there for backward compatibility test with v3.6 
    21    !                       !       
    22    cp_cfg      =  "orca"   !  name of the configuration 
    23    jp_cfg      =       2   !  resolution of the configuration 
    24    jpidta      =     182   !  1st lateral dimension ( >= jpi ) 
    25    jpjdta      =     149   !  2nd    "         "    ( >= jpj ) 
    26    jpkdta      =      31   !  number of levels      ( >= jpk ) 
    27    jpiglo      =     182   !  1st dimension of global domain --> i =jpidta 
    28    jpjglo      =     149   !  2nd    -                  -    --> j  =jpjdta 
    29    jpizoom     =       1   !  left bottom (i,j) indices of the zoom 
    30    jpjzoom     =       1   !  in data domain indices 
    31    jperio      =       4   !  lateral cond. type (between 0 and 6) 
    3221/ 
    3322!----------------------------------------------------------------------- 
    34 &namzgr        !   vertical coordinate 
    35 !----------------------------------------------------------------------- 
    36    ln_zps      = .true.    !  z-coordinate - partial steps 
    37    ln_linssh   = .true.    !  linear free surface 
    38 / 
    39 !----------------------------------------------------------------------- 
    40 &namdom        !   space and time domain (bathymetry, mesh, timestep) 
    41 !----------------------------------------------------------------------- 
    42    jphgr_msh   =       0               !  type of horizontal mesh 
    43    ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
    44    ppgphi0     =  999999.0             ! latitude  of first raw and column T-point (jphgr_msh = 1) 
    45    ppe1_deg    =  999999.0             !  zonal      grid-spacing (degrees) 
    46    ppe2_deg    =  999999.0             !  meridional grid-spacing (degrees) 
    47    ppe1_m      =  999999.0             !  zonal      grid-spacing (degrees) 
    48    ppe2_m      =  999999.0             !  meridional grid-spacing (degrees) 
    49    ppsur       =   -4762.96143546300   !  ORCA r4, r2 and r05 coefficients 
    50    ppa0        =     255.58049070440   ! (default coefficients) 
    51    ppa1        =     245.58132232490   ! 
    52    ppkth       =      21.43336197938   ! 
    53    ppacr       =       3.0             ! 
    54    ppdzmin     =  999999.              !  Minimum vertical spacing 
    55    pphmax      =  999999.              !  Maximum depth 
    56    ldbletanh   =  .FALSE.              !  Use/do not use double tanf function for vertical coordinates 
    57    ppa2        =  999999.              !  Double tanh function parameters 
    58    ppkth2      =  999999.              ! 
    59    ppacr2      =  999999.              ! 
    60 / 
    61 !----------------------------------------------------------------------- 
    62 &namcrs        !   Grid coarsening for dynamics output and/or 
    63                !   passive tracer coarsened online simulations 
     23&namzgr        !   vertical coordinate                                  (default: NO selection) 
    6424!----------------------------------------------------------------------- 
    6525/ 
    6626!----------------------------------------------------------------------- 
    67 &namtsd    !   data : Temperature  & Salinity 
     27&namzgr_isf    !   isf cavity geometry definition 
    6828!----------------------------------------------------------------------- 
    6929/ 
    7030!----------------------------------------------------------------------- 
    71 &namsbc        !   Surface Boundary Condition (surface module) 
     31&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default F) 
    7232!----------------------------------------------------------------------- 
    7333/ 
    7434!----------------------------------------------------------------------- 
    75 &namsbc_core   !   namsbc_core  CORE bulk formulae 
     35&namclo ! (closed sea : need ln_domclo = .true. in namcfg) 
    7636!----------------------------------------------------------------------- 
    7737/ 
    7838!----------------------------------------------------------------------- 
    79 &namtra_qsr    !   penetrative solar radiation 
     39&namlbc        !   lateral momentum boundary condition                  (default: NO selection) 
    8040!----------------------------------------------------------------------- 
    8141/ 
    8242!----------------------------------------------------------------------- 
    83 &namsbc_rnf    !   runoffs namelist surface boundary condition 
     43&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    8444!----------------------------------------------------------------------- 
    8545/ 
    8646!----------------------------------------------------------------------- 
    87 &namsbc_ssr    !   surface boundary condition : sea surface restoring 
     47&nambdy        !  unstructured open boundaries                          (default: OFF) 
    8848!----------------------------------------------------------------------- 
    8949/ 
    9050!----------------------------------------------------------------------- 
    91 &namsbc_alb    !   albedo parameters 
     51&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
    9252!----------------------------------------------------------------------- 
    9353/ 
    9454!----------------------------------------------------------------------- 
    95 &namberg       !   iceberg parameters 
     55&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    9656!----------------------------------------------------------------------- 
    9757/ 
    9858!----------------------------------------------------------------------- 
    99 &namlbc        !   lateral momentum boundary condition 
     59&namctl        !   Control prints                                       (default: OFF) 
    10060!----------------------------------------------------------------------- 
    10161/ 
    102 !----------------------------------------------------------------------- 
    103 &nambfr        !   bottom friction 
    104 !----------------------------------------------------------------------- 
    105 / 
    106 !----------------------------------------------------------------------- 
    107 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    108 !----------------------------------------------------------------------- 
    109    ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    110 / 
    111 !----------------------------------------------------------------------- 
    112 &nambbl        !   bottom boundary layer scheme 
    113 !----------------------------------------------------------------------- 
    114 / 
    115 !----------------------------------------------------------------------- 
    116 &nameos        !   ocean physical parameters 
    117 !----------------------------------------------------------------------- 
    118    ln_teos10    = .true.         !  = Use TEOS-10 equation of state 
    119 / 
    120 !----------------------------------------------------------------------- 
    121 &namtra_adv    !   advection scheme for tracer 
    122 !----------------------------------------------------------------------- 
    123    ln_traadv_fct =  .true.    !  FCT scheme 
    124       nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
    125       nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
    126       nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
    127       !                             !        (number of sub-timestep = nn_fct_zts) 
    128 / 
    129 !----------------------------------------------------------------------- 
    130 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    131 !----------------------------------------------------------------------- 
    132 / 
    133 !---------------------------------------------------------------------------------- 
    134 &namtra_ldf    !   lateral diffusion scheme for tracers 
    135 !---------------------------------------------------------------------------------- 
    136    !                       !  Operator type: 
    137    ln_traldf_lap   =  .true.   !    laplacian operator 
    138    ln_traldf_blp   =  .false.  !  bilaplacian operator 
    139    !                       !  Direction of action: 
    140    ln_traldf_lev   =  .false.  !  iso-level 
    141    ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    142    ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
    143    ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
    144    ! 
    145    !                       !  iso-neutral options:         
    146    ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
    147    rn_slpmax       =   0.01    !  slope limit                      (both operators) 
    148    ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
    149    rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
    150    ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
    151    ! 
    152    !                       !  Coefficients: 
    153    nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
    154    !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    155    !                                !   =  0           constant 
    156    !                                !   = 10 F(k)      =ldf_c1d 
    157    !                                !   = 20 F(i,j)    =ldf_c2d 
    158    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    159    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    160    !                                !   = 31 F(i,j,k,t)=F(local velocity) 
    161    rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
    162    rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    163 / 
    164 !---------------------------------------------------------------------------------- 
    165 &namtra_ldfeiv !   eddy induced velocity param. 
    166 !---------------------------------------------------------------------------------- 
    167    ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
    168    ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
    169    rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
    170    nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
    171    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    172    !                                !   =  0           constant 
    173    !                                !   = 10 F(k)      =ldf_c1d 
    174    !                                !   = 20 F(i,j)    =ldf_c2d 
    175    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    176    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    177 / 
    178 !----------------------------------------------------------------------- 
    179 &namtra_dmp    !   tracer: T & S newtonian damping 
    180 !----------------------------------------------------------------------- 
    181 / 
    182 !----------------------------------------------------------------------- 
    183 &namdyn_adv    !   formulation of the momentum advection 
    184 !----------------------------------------------------------------------- 
    185 / 
    186 !----------------------------------------------------------------------- 
    187 &namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    188 !----------------------------------------------------------------------- 
    189    ln_dynvor_ene = .false. !  enstrophy conserving scheme 
    190    ln_dynvor_ens = .false. !  energy conserving scheme 
    191    ln_dynvor_mix = .false. !  mixed scheme 
    192    ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    193       nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    194 / 
    195 !----------------------------------------------------------------------- 
    196 &namdyn_hpg    !   Hydrostatic pressure gradient option 
    197 !----------------------------------------------------------------------- 
    198 / 
    199 !----------------------------------------------------------------------- 
    200 &namdyn_spg    !   Surface pressure gradient 
    201 !----------------------------------------------------------------------- 
    202    ln_dynspg_ts = .true.   !  Split-explicit free surface 
    203 / 
    204 !----------------------------------------------------------------------- 
    205 &namdyn_ldf    !   lateral diffusion on momentum 
    206 !----------------------------------------------------------------------- 
    207    !                       !  Type of the operator : 
    208    !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
    209    ln_dynldf_lap =  .true.     !    laplacian operator 
    210    ln_dynldf_blp =  .false.    !  bilaplacian operator 
    211    !                       !  Direction of action  : 
    212    ln_dynldf_lev =  .true.     !  iso-level 
    213    ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
    214    ln_dynldf_iso =  .false.    !  iso-neutral 
    215    !                       !  Coefficient 
    216    nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
    217    !                                !  =-30  read in eddy_viscosity_3D.nc file 
    218    !                                !  =-20  read in eddy_viscosity_2D.nc file 
    219    !                                !  =  0  constant  
    220    !                                !  = 10  F(k)=c1d 
    221    !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
    222    !                                !  = 30  F(i,j,k)=c2d*c1d 
    223    !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
    224    rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
    225    rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
    226    rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    227 / 
    228 !----------------------------------------------------------------------- 
    229 &namzdf        !   vertical physics 
    230 !----------------------------------------------------------------------- 
    231 / 
    232 !----------------------------------------------------------------------- 
    233 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    234 !----------------------------------------------------------------------- 
    235 / 
    236 !----------------------------------------------------------------------- 
    237 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    238 !----------------------------------------------------------------------- 
    239 / 
    240 !----------------------------------------------------------------------- 
    241 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    242 !----------------------------------------------------------------------- 
    243 / 
    244 !----------------------------------------------------------------------- 
    245 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    246 !----------------------------------------------------------------------- 
    247 / 
    248 !----------------------------------------------------------------------- 
    249 &namctl        !   Control prints & Benchmark 
    250 !----------------------------------------------------------------------- 
    251 / 
    252 !----------------------------------------------------------------------- 
    253 &namptr       !   Poleward Transport Diagnostic 
    254 !----------------------------------------------------------------------- 
    255 / 
    256 !----------------------------------------------------------------------- 
    257 &namhsb       !  Heat and salt budgets 
    258 !----------------------------------------------------------------------- 
    259 / 
    260 !----------------------------------------------------------------------- 
    261 &namobs       !  observation usage 
    262 !----------------------------------------------------------------------- 
    263 / 
    264 !----------------------------------------------------------------------- 
    265 &nam_asminc   !   assimilation increments                               ('key_asminc') 
    266 !----------------------------------------------------------------------- 
    267 / 
  • utils/tools/DOMAINcfg/namelist_ref

    r12207 r12414  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !!                            namelist_ref 
     2!! NEMO/OCE :   Reference namelist_ref                                !! 
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    4 !! NEMO/OPA  :  1 - run manager      (namrun) 
    5 !! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
    6 !!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    7 !!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    8 !!                                    namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave) 
    9 !!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    10 !!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    11 !!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
    12 !!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    13 !!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 
    14 !!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    15 !!             10 - miscellaneous    (nammpp, namctl) 
    16 !!             11 - Obs & Assim      (namobs, nam_asminc) 
     4!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namzgr, namzgr_sco ) 
     5!!              2 - diagnostics      (namnc4) 
     6!!              3 - miscellaneous    (nammpp, namctl) 
    177!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    18  
    19 !!====================================================================== 
    20 !!                   ***  Run management namelists  *** 
    21 !!====================================================================== 
    22 !!   namrun       parameters of the run 
    23 !!====================================================================== 
    24 ! 
    258!----------------------------------------------------------------------- 
    269&namrun        !   parameters of the run 
    2710!----------------------------------------------------------------------- 
    28    nn_no       =       0   !  job number (no more used...) 
     11   nn_no       =       0   !  Assimilation cycle index 
    2912   cn_exp      =  "ORCA2"  !  experience name 
    3013   nn_it000    =       1   !  first time step 
    31    nn_itend    =    5475   !  last  time step (std 5475) 
     14   nn_itend    =    5840   !  last  time step (std 5840) 
    3215   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    3316   nn_time0    =       0   !  initial time of day in hhmm 
    3417   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    3518   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    36       nn_euler    =    1            !  = 0 : start with forward time step if ln_rstart=T 
    37       nn_rstctl   =    0            !  restart control ==> activated only if ln_rstart=T 
    38       !                             !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    39       !                             !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
    40       !                             !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
     19      nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
     20      nn_rstctl   =    0      !  restart control ==> activated only if ln_rstart=T 
     21      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
     22      !                          !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     23      !                          !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
    4124      cn_ocerst_in    = "restart"   !  suffix of ocean restart name (input) 
    4225      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts 
    4326      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output) 
    44       cn_ocerst_outdir= "."         !  directory in which to write output ocean restarts 
     27      cn_ocerst_outdir = "."         !  directory in which to write output ocean restarts 
    4528   ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model 
    4629   nn_istate   =       0   !  output the initial state (1) or not (0) 
    4730   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
    48    nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1) 
     31   nn_stock    =    5840   !  frequency of creation of a restart file (modulo referenced to 1) 
    4932   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 
    5033   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000) 
     
    5437   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    5538/ 
    56 ! 
    57 !!====================================================================== 
    58 !!                      ***  Domain namelists  *** 
    59 !!====================================================================== 
    60 !!   namcfg       parameters of the configuration 
    61 !!   namzgr       vertical coordinate                                   (default: NO selection) 
    62 !!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    63 !!   namdom       space and time domain (bathymetry, mesh, timestep) 
    64 !!   namwad       Wetting and drying                                    (default F) 
    65 !!   namtsd       data: temperature & salinity 
    66 !!   namcrs       coarsened grid (for outputs and/or TOP)               ("key_crs") 
    67 !!   namc1d       1D configuration options                              ("key_c1d") 
    68 !!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
    69 !!   namc1d_uvd   1D data (currents)                                    ("key_c1d") 
    70 !!====================================================================== 
    71 ! 
    72 !----------------------------------------------------------------------- 
    73 &namcfg        !   parameters of the configuration 
    74 !----------------------------------------------------------------------- 
    75    ! 
    76    ln_e3_dep   = .true.    ! =T : e3=dk[depth] in discret sens.  
    77    !                       !      ===>>> will become the only possibility in v4.0 
    78    !                       ! =F : e3 analytical derivative of depth function 
    79    !                       !      only there for backward compatibility test with v3.6 
    80    ! 
    81    cp_cfg      = "default" !  name of the configuration 
    82    cp_cfz      = "no zoom" !  name of the zoom of configuration 
    83    jp_cfg      =      0    !  resolution of the configuration 
    84    jpidta      =     10    !  1st lateral dimension ( >= jpi ) 
    85    jpjdta      =     12    !  2nd    "         "    ( >= jpj ) 
    86    jpkdta      =     31    !  number of levels      ( >= jpk ) 
    87    jpiglo      =     10    !  1st dimension of global domain --> i =jpidta 
    88    jpjglo      =     12    !  2nd    -                  -    --> j =jpjdta 
    89    jpizoom     =      1    !  left bottom (i,j) indices of the zoom 
    90    jpjzoom     =      1    !  in data domain indices 
    91    jperio      =      0    !  lateral cond. type (between 0 and 6) 
    92                                  !  = 0 closed                 ;   = 1 cyclic East-West 
    93                                  !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
    94                                  !  = 4 cyclic East-West AND North fold T-point pivot 
    95                                  !  = 5 North fold F-point pivot 
    96                                  !  = 6 cyclic East-West AND North fold F-point pivot 
    97    ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present 
    98                            !  in netcdf input files, as the start j-row for reading 
    99 / 
    100 !----------------------------------------------------------------------- 
    101 &namzgr        !   vertical coordinate                                  (default: NO selection) 
    102 !----------------------------------------------------------------------- 
    103    ln_zco      = .false.   !  z-coordinate - full    steps 
    104    ln_zps      = .false.   !  z-coordinate - partial steps 
    105    ln_sco      = .false.   !  s- or hybrid z-s-coordinate 
    106    ln_isfcav   = .false.   !  ice shelf cavity 
    107    ln_linssh   = .false.   !  linear free surface 
    108 / 
    109 !----------------------------------------------------------------------- 
    110 &namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default F) 
    111 !----------------------------------------------------------------------- 
    112    ln_s_sh94   = .false.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
    113    ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
    114    ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
    115                            !  stretching coefficients for all functions 
    116    rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m) 
    117    rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
    118    rn_hc       =  150.0    !  critical depth for transition to stretched coordinates 
    119                         !!!!!!!  Envelop bathymetry 
    120    rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1) 
    121                         !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
    122    rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
    123    rn_bb       =    0.8    !  stretching with SH94 s-sigma 
    124                         !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
    125    rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
    126    rn_efold    =    0.0    !  efold length scale for transition to stretched coord 
    127    rn_zs       =    1.0    !  depth of surface grid box 
    128                            !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 
    129    rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb 
    130    rn_zb_b     =   -0.2    !  offset for calculating Zb 
    131                         !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
    132    rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1) 
    133 / 
    13439!----------------------------------------------------------------------- 
    13540&namdom        !   space and time domain (bathymetry, mesh, timestep) 
     
    13742   nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file 
    13843   rn_bathy    =    0.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
    139    nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    140    nn_msh      =    1      !  create (=1) a mesh file or not (=0) 
     44   nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
    14145   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
    142    rn_isfhmin  =    1.00   !  treshold (m) to discriminate grounding ice to floating ice 
    14346   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
    14447   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
     
    17275/ 
    17376!----------------------------------------------------------------------- 
    174 &namwad        !   Wetting and drying                                   (default F) 
    175 !----------------------------------------------------------------------- 
    176    ln_wd       = .false.   !  T/F activation of wetting and drying 
    177    rn_wdmin1   =  0.1      !  Minimum wet depth on dried cells 
    178    rn_wdmin2   =  0.01     !  Tolerance of min wet depth on dried cells 
    179    rn_wdld     =  20.0     !  Land elevation below which wetting/drying is allowed 
    180    nn_wdit     =  10       !  Max iterations for W/D limiter 
    181 / 
    182 !----------------------------------------------------------------------- 
    183 &namtsd        !   data : Temperature  & Salinity 
    184 !----------------------------------------------------------------------- 
    185 !              !  file name                 ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    186 !              !                            !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    187    sn_tem = 'data_1m_potential_temperature_nomask',     -1      ,'votemper',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    188    sn_sal = 'data_1m_salinity_nomask'             ,     -1      ,'vosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    189    ! 
    190    cn_dir      = './'      !  root directory for the location of the runoff files 
    191    ln_tsd_init = .true.    !  Initialisation of ocean T & S with T & S input data (T) or not (F) 
    192    ln_tsd_tradmp = .true.  !  damping of ocean T & S toward T & S input data (T) or not (F) 
    193 / 
    194 !----------------------------------------------------------------------- 
    195 &namcrs        !   coarsened grid (for outputs and/or TOP)              ("key_crs") 
    196 !----------------------------------------------------------------------- 
    197    nn_factx    = 3         !  Reduction factor of x-direction 
    198    nn_facty    = 3         !  Reduction factor of y-direction 
    199    nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT 
    200                            !  0, coarse grid is binned with preferential treatment of the north fold 
    201                            !  1, coarse grid is binned with centering at the equator 
    202                            !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty. 
    203    nn_msh_crs  = 1         !  create (=1) a mesh file or not (=0) 
    204    nn_crs_kz   = 0         ! 0, MEAN of volume boxes 
    205                            ! 1, MAX of boxes 
    206                            ! 2, MIN of boxes 
    207    ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F ) 
    208 / 
    209 !----------------------------------------------------------------------- 
    210 &namc1d        !   1D configuration options                             ("key_c1d") 
    211 !----------------------------------------------------------------------- 
    212    rn_lat1d    =      50   !  Column latitude (default at PAPA station) 
    213    rn_lon1d    =    -145   !  Column longitude (default at PAPA station) 
    214    ln_c1d_locpt=  .true.   ! Localization of 1D config in a grid (T) or independant point (F) 
    215 / 
    216 !----------------------------------------------------------------------- 
    217 &namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d") 
    218 !----------------------------------------------------------------------- 
    219    ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    220 / 
    221 !----------------------------------------------------------------------- 
    222 &namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
    223 !----------------------------------------------------------------------- 
    224 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    225 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    226    sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , '' 
    227    sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , '' 
    228 ! 
    229    cn_dir        = './'    !  root directory for the location of the files 
    230    ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F) 
    231    ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F) 
    232 / 
    233  
    234 !!====================================================================== 
    235 !!            ***  Surface Boundary Condition namelists  *** 
    236 !!====================================================================== 
    237 !!   namsbc          surface boundary condition 
    238 !!   namsbc_ana      analytical         formulation                     (ln_ana     =T) 
    239 !!   namsbc_flx      flux               formulation                     (ln_flx     =T) 
    240 !!   namsbc_clio     CLIO bulk formulae formulation                     (ln_blk_clio=T) 
    241 !!   namsbc_core     CORE bulk formulae formulation                     (ln_blk_core=T) 
    242 !!   namsbc_mfs      MFS  bulk formulae formulation                     (ln_blk_mfs =T) 
    243 !!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" ) 
    244 !!   namsbc_sas      StAndalone Surface module 
    245 !!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T) 
    246 !!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    247 !!   namsbc_isf      ice shelf melting/freezing                         (nn_isf     >0) 
    248 !!   namsbc_iscpl    coupling option between land ice model and ocean 
    249 !!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    250 !!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
    251 !!   namsbc_alb      albedo parameters 
    252 !!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    253 !!   namberg         iceberg floats                                     (ln_icebergs=T) 
    254 !!====================================================================== 
    255 ! 
    256 !----------------------------------------------------------------------- 
    257 &namsbc        !   Surface Boundary Condition (surface module) 
    258 !----------------------------------------------------------------------- 
    259    nn_fsbc     = 5         !  frequency of surface boundary condition computation 
    260                            !     (also = the frequency of sea-ice & iceberg model call) 
    261                      ! Type of air-sea fluxes  
    262    ln_ana      = .false.   !  analytical formulation                    (T => fill namsbc_ana ) 
    263    ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    264    ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio) 
    265    ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    266    ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
    267                      ! Type of coupling (Ocean/Ice/Atmosphere) : 
    268    ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
    269    ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 ) 
    270    nn_components = 0       !  configuration of the opa-sas OASIS coupling 
    271                            !  =0 no opa-sas OASIS coupling: default single executable configuration 
    272                            !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component 
    273                            !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component  
    274    nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
    275                            !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
    276                            !  = 0  Average per-category fluxes (forced and coupled mode) 
    277                            !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
    278                            !  = 2  Redistribute a single flux over categories (coupled mode only) 
    279                      ! Sea-ice : 
    280    nn_ice      = 2         !  =0 no ice boundary condition   , 
    281                            !  =1 use observed ice-cover      , 
    282                            !  =2 ice-model used                         ("key_lim3", "key_lim2", "key_cice") 
    283    nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    284                            !  =1 levitating ice with mass and salt exchange but no presure effect 
    285                            !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    286                      ! Misc. options of sbc :  
    287    ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr ) 
    288    ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    289    ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
    290    ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    291    nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
    292                            !     =1 global mean of e-p-r set to zero at each time step 
    293                            !     =2 annual global mean of e-p-r set to zero 
    294    ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    295    ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf) 
    296    ln_wave     = .false.   !  coupling with surface wave                (T => fill namsbc_wave) 
    297    nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
    298                            !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
    299 / 
    300 !----------------------------------------------------------------------- 
    301 &namsbc_ana    !   analytical surface boundary condition 
    302 !----------------------------------------------------------------------- 
    303    nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps 
    304    rn_utau0    =   0.5     !  uniform value for the i-stress 
    305    rn_vtau0    =   0.e0    !  uniform value for the j-stress 
    306    rn_qns0     =   0.e0    !  uniform value for the total heat flux 
    307    rn_qsr0     =   0.e0    !  uniform value for the solar radiation 
    308    rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P) 
    309 / 
    310 !----------------------------------------------------------------------- 
    311 &namsbc_flx    !   surface boundary condition : flux formulation 
    312 !----------------------------------------------------------------------- 
    313 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    314 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    315    sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    316    sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    317    sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    318    sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    319    sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    320  
    321    cn_dir      = './'      !  root directory for the location of the flux files 
    322 / 
    323 !----------------------------------------------------------------------- 
    324 &namsbc_clio   !   namsbc_clio  CLIO bulk formulae 
    325 !----------------------------------------------------------------------- 
    326 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    327 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    328    sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    329    sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    330    sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    331    sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    332    sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    333    sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    334    sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    335  
    336    cn_dir      = './'      !  root directory for the location of the bulk files are 
    337 / 
    338 !----------------------------------------------------------------------- 
    339 &namsbc_core   !   namsbc_core  CORE bulk formulae 
    340 !----------------------------------------------------------------------- 
    341 !              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
    342 !              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
    343    sn_wndi     = 'u_10.15JUNE2009_fill'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , '' 
    344    sn_wndj     = 'v_10.15JUNE2009_fill'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , '' 
    345    sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    346    sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    347    sn_tair     = 't_10.15JUNE2009_fill'       ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    348    sn_humi     = 'q_10.15JUNE2009_fill'       ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    349    sn_prec     = 'ncar_precip.15JUNE2009_fill',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    350    sn_snow     = 'ncar_precip.15JUNE2009_fill',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    351    sn_tdif     = 'taudif_core'                ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    352  
    353    cn_dir      = './'      !  root directory for the location of the bulk files 
    354    ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    355    rn_zqt      = 10.       !  Air temperature and humidity reference height (m) 
    356    rn_zu       = 10.       !  Wind vector reference height (m) 
    357    rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    358    rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    359    rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity 
    360                            !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    361 / 
    362 !----------------------------------------------------------------------- 
    363 &namsbc_mfs   !   namsbc_mfs  MFS bulk formulae 
    364 !----------------------------------------------------------------------- 
    365 !              !  file name  ! frequency (hours) ! variable ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask ! 
    366 !              !             !  (if <0  months)  !   name   !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      ! 
    367    sn_wndi     =   'ecmwf'   ,        6          , 'u10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    368    sn_wndj     =   'ecmwf'   ,        6          , 'v10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    369    sn_clc      =   'ecmwf'   ,        6          , 'clc'    ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,   '' 
    370    sn_msl      =   'ecmwf'   ,        6          , 'msl'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    371    sn_tair     =   'ecmwf'   ,        6          , 't2'     ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    372    sn_rhm      =   'ecmwf'   ,        6          , 'rh'     ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,   '' 
    373    sn_prec     =   'ecmwf'   ,        6          , 'precip' ,    .true.    , .true. , 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    374  
    375    cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files 
    376 / 
    377 !----------------------------------------------------------------------- 
    378 &namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
    379 !----------------------------------------------------------------------- 
    380 !                    !     description      !  multiple  !    vector   !      vector          ! vector ! 
    381 !                    !                      ! categories !  reference  !    orientation       ! grids  ! 
    382 ! send 
    383    sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    384    sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    385    sn_snd_thick  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    386    sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    387    sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    388 ! receive 
    389    sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    390    sn_rcv_taumod =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    391    sn_rcv_tau    =   'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
    392    sn_rcv_dqnsdt =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    393    sn_rcv_qsr    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    394    sn_rcv_qns    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    395    sn_rcv_emp    =   'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    396    sn_rcv_rnf    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    397    sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    398    sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    399 ! 
    400    nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
    401    ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models 
    402    !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    403 / 
    404 !----------------------------------------------------------------------- 
    405 &namsbc_sas    !   analytical surface boundary condition 
    406 !----------------------------------------------------------------------- 
    407 !              !  file name  ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    408 !              !             !  (if <0  months)  !   name    !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    409    sn_usp      = 'sas_grid_U',     120           , 'vozocrtx',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    410    sn_vsp      = 'sas_grid_V',     120           , 'vomecrty',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    411    sn_tem      = 'sas_grid_T',     120           , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    412    sn_sal      = 'sas_grid_T',     120           , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    413    sn_ssh      = 'sas_grid_T',     120           , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    414    sn_e3t      = 'sas_grid_T',     120           , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    415    sn_frq      = 'sas_grid_T',     120           , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    416  
    417    ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field 
    418    ln_read_frq = .false.   !  specify whether we must read frq or not 
    419    cn_dir      = './'      !  root directory for the location of the bulk files are 
    420 / 
    421 !----------------------------------------------------------------------- 
    422 &namtra_qsr    !   penetrative solar radiation                          (ln_traqsr=T) 
    423 !----------------------------------------------------------------------- 
    424 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    425 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    426    sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    427  
    428    cn_dir      = './'      !  root directory for the location of the runoff files 
    429    ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration 
    430    ln_qsr_2bd  = .false.   !  2 bands              light penetration 
    431    ln_qsr_bio  = .false.   !  bio-model light penetration 
    432    nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0) 
    433    rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1) 
    434    rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction 
    435    rn_si1      =   23.0    !  2 bands: longest depth of extinction 
    436    ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3 
    437 / 
    438 !----------------------------------------------------------------------- 
    439 &namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf=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_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    444    sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    445    sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    446    sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    447    sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    448  
    449    cn_dir      = './'      !  root directory for the location of the runoff files 
    450    ln_rnf_mouth= .true.    !  specific treatment at rivers mouths 
    451       rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T) 
    452       rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) 
    453    rn_rfact    =   1.e0    !  multiplicative factor for runoff 
    454    ln_rnf_depth= .false.   !  read in depth information for runoff 
    455    ln_rnf_tem  = .false.   !  read in temperature information for runoff 
    456    ln_rnf_sal  = .false.   !  read in salinity information for runoff 
    457    ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file 
    458       rn_rnf_max  = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) 
    459       rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) 
    460       nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0) 
    461 / 
    462 !----------------------------------------------------------------------- 
    463 &namsbc_isf    !  Top boundary layer (ISF)                              (nn_isf >0) 
    464 !----------------------------------------------------------------------- 
    465 !              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    466 !              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    467 ! nn_isf == 4 
    468    sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    469 ! nn_isf == 3 
    470    sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    471 ! nn_isf == 2 and 3 
    472    sn_depmax_isf='rnfisf'  ,         -12       ,'sozisfmax',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    473    sn_depmin_isf='rnfisf'  ,         -12       ,'sozisfmin',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    474 ! nn_isf == 2 
    475    sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    '' 
    476 ! 
    477 ! for all case 
    478    nn_isf      = 1         !  ice shelf melting/freezing 
    479                            !  1 = presence of ISF    2 = bg03 parametrisation  
    480                            !  3 = rnf file for isf   4 = ISF fwf specified 
    481                            !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    482 ! only for nn_isf = 1 or 2 
    483    rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    484    rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    485 ! only for nn_isf = 1 or 4 
    486    rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    487    !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    488 ! only for nn_isf = 1 
    489    nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    490    !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    491    nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    492    !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    493    !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    494 / 
    495 !----------------------------------------------------------------------- 
    496 &namsbc_iscpl  !   land ice / ocean coupling option                      
    497 !----------------------------------------------------------------------- 
    498    nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    499    ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
    500    nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
    501 / 
    502 !----------------------------------------------------------------------- 
    503 &namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk 
    504 !----------------------------------------------------------------------- 
    505 !              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    506 !              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    507    sn_apr      = 'patm'    ,         -1        ,'somslpre',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,      '' 
    508  
    509    cn_dir      = './'      !  root directory for the location of the bulk files 
    510    rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/ 
    511    ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F) 
    512    ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data 
    513 / 
    514 !----------------------------------------------------------------------- 
    515 &namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr=T) 
    516 !----------------------------------------------------------------------- 
    517 !              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    518 !              !           !  (if <0  months)  !   name   !   (logical) !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    519    sn_sst      = 'sst_data',        24         ,  'sst'   ,    .false.  , .false., 'yearly'  ,    ''    ,    ''    ,     '' 
    520    sn_sss      = 'sss_data',        -1         ,  'sss'   ,    .true.   , .true. , 'yearly'  ,    ''    ,    ''    ,     '' 
    521  
    522    cn_dir      = './'      !  root directory for the location of the runoff files 
    523    nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0) 
    524    nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2) 
    525                            !  or to SSS only (=1) or no damping term (=0) 
    526    rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K] 
    527    rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day] 
    528    ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2) 
    529    rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
    530 / 
    531 !----------------------------------------------------------------------- 
    532 &namsbc_alb    !   albedo parameters 
    533 !----------------------------------------------------------------------- 
    534    nn_ice_alb  =    0      !  parameterization of ice/snow albedo 
    535                            !     0: Shine & Henderson-Sellers (JGR 1985) 
    536                            !     1: "home made" based on Brandt et al. (J. Climate 2005) 
    537                            !                         and Grenfell & Perovich (JGR 2004) 
    538    rn_albice   =  0.53     !  albedo of bare puddled ice (values from 0.49 to 0.58) 
    539                            !     0.53 (default) => if nn_ice_alb=0 
    540                            !     0.50 (default) => if nn_ice_alb=1 
    541 / 
    542 !----------------------------------------------------------------------- 
    543 &namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
    544 !----------------------------------------------------------------------- 
    545 !              ! file name ! frequency (hours) ! variable    ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    546 !              !           !  (if <0  months)  !   name      !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    547    sn_cdg      = 'cdg_wave',        1          , 'drag_coeff',   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    548    sn_usd      = 'sdw_wave',        1          , 'u_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    549    sn_vsd      = 'sdw_wave',        1          , 'v_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    550    sn_wn       = 'sdw_wave',        1          , 'wave_num'  ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    551 ! 
    552    cn_dir_cdg  = './'      !  root directory for the location of drag coefficient files 
    553    ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model 
    554    ln_sdw      = .false.   !  Computation of 3D stokes drift                
    555 / 
    556 !----------------------------------------------------------------------- 
    557 &namberg       !   iceberg parameters                                   (default: No iceberg) 
    558 !----------------------------------------------------------------------- 
    559       ln_icebergs              = .false.              ! iceberg floats or not 
    560       ln_bergdia               = .true.               ! Calculate budgets 
    561       nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
    562       nn_verbose_write         = 15                   ! Timesteps between verbose messages 
    563       nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
    564                                                       ! Initial mass required for an iceberg of each class 
    565       rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
    566                                                       ! Proportion of calving mass to apportion to each class 
    567       rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    568                                                       ! Ratio between effective and real iceberg mass (non-dim) 
    569                                                       ! i.e. number of icebergs represented at a point 
    570       rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    571                                                       ! thickness of newly calved bergs (m) 
    572       rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
    573       rn_rho_bergs             = 850.                 ! Density of icebergs 
    574       rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
    575       ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
    576       rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
    577       rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
    578       ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
    579       nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
    580                                                       ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
    581       rn_test_box              = 108.0,  116.0, -66.0, -58.0 
    582       rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
    583  
    584 !            ! file name ! frequency (hours) !   variable   ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    585 !            !           !  (if <0  months)  !     name     !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    586       sn_icb =  'calving',       -1          , 'calvingmask',   .true.    , .true.  , 'yearly'  ,    ''    ,    ''    ,     '' 
    587  
    588       cn_dir = './' 
    589 / 
    590  
    591 !!====================================================================== 
    592 !!               ***  Lateral boundary condition  *** 
    593 !!====================================================================== 
    594 !!   namlbc        lateral momentum boundary condition 
    595 !!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
    596 !!   nam_tide      Tidal forcing  
    597 !!   nambdy        Unstructured open boundaries                         ("key_bdy") 
    598 !!   nambdy_dta    Unstructured open boundaries - external data         ("key_bdy") 
    599 !!   nambdy_tide   tidal forcing at open boundaries                     ("key_bdy_tides") 
    600 !!====================================================================== 
    601 ! 
    602 !----------------------------------------------------------------------- 
    603 &namlbc        !   lateral momentum boundary condition 
     77&namcfg        !   parameters of the configuration 
     78!----------------------------------------------------------------------- 
     79   ! 
     80   ln_e3_dep   = .true.    ! =T : e3=dk[depth] in discret sens. 
     81   !                       !      ===>>> will become the only possibility in v4.0 
     82   !                       ! =F : e3 analytical derivative of depth function 
     83   !                       !      only there for backward compatibility test with v3.6 
     84   ! 
     85   cp_cfg      = "default" !  name of the configuration 
     86   cp_cfz      = "no zoom" !  name of the zoom of configuration 
     87   jp_cfg      =      0    !  resolution of the configuration 
     88   jpkdta      =     31    !  number of levels      ( >= jpk ) 
     89   jpiglo      =     10    !  1st dimension of global domain --> i =jpidta 
     90   jpjglo      =     12    !  2nd    -                  -    --> j =jpjdta 
     91   jperio      =      0    !  lateral cond. type (between 0 and 6) 
     92                                 !  = 0 closed                 ;   = 1 cyclic East-West 
     93                                 !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
     94                                 !  = 4 cyclic East-West AND North fold T-point pivot 
     95                                 !  = 5 North fold F-point pivot 
     96                                 !  = 6 cyclic East-West AND North fold F-point pivot 
     97   ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present 
     98                           !  in netcdf input files, as the start j-row for reading 
     99   ln_domclo = .false.     ! computation of closed sea masks (see namclo) 
     100/ 
     101!----------------------------------------------------------------------- 
     102&namzgr        !   vertical coordinate                                  (default: NO selection) 
     103!----------------------------------------------------------------------- 
     104   ln_zco      = .false.   !  z-coordinate - full    steps 
     105   ln_zps      = .false.   !  z-coordinate - partial steps 
     106   ln_sco      = .false.   !  s- or hybrid z-s-coordinate 
     107   ln_isfcav   = .false.   !  ice shelf cavity             (T: see namzgr_isf) 
     108   ln_linssh   = .false.   !  linear free surface 
     109/ 
     110!----------------------------------------------------------------------- 
     111&namzgr_isf    !   isf cavity geometry definition                       (default: OFF) 
     112!----------------------------------------------------------------------- 
     113   rn_isfdep_min    = 10.         ! minimum isf draft tickness (if lower, isf draft set to this value) 
     114   rn_glhw_min      = 1.e-3       ! minimum water column thickness to define the grounding line 
     115   rn_isfhw_min     = 10          ! minimum water column thickness in the cavity once the grounding line defined. 
     116   ln_isfchannel    = .false.     ! remove channel (based on 2d mask build from isfdraft-bathy) 
     117   ln_isfconnect    = .false.     ! force connection under the ice shelf (based on 2d mask build from isfdraft-bathy) 
     118      nn_kisfmax       = 999         ! limiter in level on the previous condition. (if change larger than this number, get back to value before we enforce the connection) 
     119      rn_zisfmax       = 7000.       ! limiter in m     on the previous condition. (if change larger than this number, get back to value before we enforce the connection) 
     120   ln_isfcheminey   = .false.     ! close cheminey 
     121   ln_isfsubgl      = .false.     ! remove subglacial lake created by the remapping process 
     122      rn_isfsubgllon   =    0.0      !  longitude of the seed to determine the open ocean 
     123      rn_isfsubgllat   =    0.0      !  latitude  of the seed to determine the open ocean 
     124/ 
     125!----------------------------------------------------------------------- 
     126&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default: OFF) 
     127!----------------------------------------------------------------------- 
     128   ln_s_sh94   = .false.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
     129   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
     130   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
     131                           !  stretching coefficients for all functions 
     132   rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m) 
     133   rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
     134   rn_hc       =  150.0    !  critical depth for transition to stretched coordinates 
     135                        !!!!!!!  Envelop bathymetry 
     136   rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1) 
     137                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
     138   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
     139   rn_bb       =    0.8    !  stretching with SH94 s-sigma 
     140                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
     141   rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
     142   rn_efold    =    0.0    !  efold length scale for transition to stretched coord 
     143   rn_zs       =    1.0    !  depth of surface grid box 
     144                           !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 
     145   rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb 
     146   rn_zb_b     =   -0.2    !  offset for calculating Zb 
     147                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
     148   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1) 
     149/ 
     150!----------------------------------------------------------------------- 
     151&namclo ! (closed sea : need ln_domclo = .true. in namcfg)              (default: OFF) 
     152!----------------------------------------------------------------------- 
     153   rn_lon_opnsea =  0.0    ! longitude seed of open ocean 
     154   rn_lat_opnsea =  0.0    ! latitude  seed of open ocean 
     155   nn_closea = 8           ! number of closed seas ( = 0; only the open_sea mask will be computed) 
     156   ! 
     157   !                name   ! lon_src ! lat_src ! lon_trg ! lat_trg ! river mouth area   ! correction scheme ! radius trg   ! id trg 
     158   !                       ! (degree)! (degree)! (degree)! (degree)! local/coast/global ! (glo/rnf/emp)     !     (m)      ! 
     159   ! North American lakes 
     160   sn_lake(1) = 'superior' ,  -86.57 ,  47.30  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0   , 2     
     161   sn_lake(2) = 'michigan' ,  -87.06 ,  42.74  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0   , 2     
     162   sn_lake(3) = 'huron'    ,  -82.51 ,  44.74  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0   , 2     
     163   sn_lake(4) = 'erie'     ,  -81.13 ,  42.25  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0   , 2     
     164   sn_lake(5) = 'ontario'  ,  -77.72 ,  43.62  , -66.49  , 50.45   , 'local'            , 'rnf'             ,   550000.0   , 2     
     165   ! African Lake 
     166   sn_lake(6) = 'victoria' ,   32.93 ,  -1.08  ,  30.44  , 31.37   , 'coast'            , 'emp'             ,   100000.0   , 3     
     167   ! Asian Lakes 
     168   sn_lake(7) = 'caspian'  ,   50.0  ,  44.0   ,   0.0   ,  0.0    , 'global'           , 'glo'             ,        0.0   , 1      
     169   sn_lake(8) = 'aral'     ,   60.0  ,  45.0   ,   0.0   ,  0.0    , 'global'           , 'glo'             ,        0.0   , 1     
     170/ 
     171!----------------------------------------------------------------------- 
     172&namlbc        !   lateral momentum boundary condition                  (default: NO selection) 
    604173!----------------------------------------------------------------------- 
    605174   !                       !  free slip  !   partial slip  !   no slip   ! strong slip 
    606    rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
     175   rn_shlat    =  0        !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    607176   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs. 
    608177/ 
     
    610179&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    611180!----------------------------------------------------------------------- 
    612    nn_cln_update =    3    !  baroclinic update frequency 
    613181   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    614182   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    615183   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
    616    ln_chk_bathy  = .FALSE. ! 
    617 / 
    618 !----------------------------------------------------------------------- 
    619 &nam_tide      !   tide parameters                                      ("key_tide") 
    620 !----------------------------------------------------------------------- 
    621    ln_tide_pot = .true.    !  use tidal potential forcing 
    622    ln_tide_ramp= .false.   ! 
    623    rdttideramp =    0.     ! 
    624    clname(1)   = 'DUMMY'   !  name of constituent - all tidal components must be set in namelist_cfg 
    625 / 
    626 !----------------------------------------------------------------------- 
    627 &nambdy        !  unstructured open boundaries                          ("key_bdy") 
    628 !----------------------------------------------------------------------- 
    629     nb_bdy         = 0                    !  number of open boundary sets 
    630     ln_coords_file = .true.               !  =T : read bdy coordinates from file 
    631     cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files 
    632     ln_mask_file   = .false.              !  =T : read mask from file 
    633     cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.) 
    634     cn_dyn2d       = 'none'               ! 
    635     nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state 
    636                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    637                                           !  = 2, use tidal harmonic forcing data from files 
    638                                           !  = 3, use external data AND tidal harmonic forcing 
    639     cn_dyn3d      =  'none'               ! 
    640     nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state 
    641                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    642     cn_tra        =  'none'               ! 
    643     nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state 
    644                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    645     cn_ice_lim      =  'none'             ! 
    646     nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state 
    647                                           !  = 1, bdy data are read in 'bdydata   .nc' files 
    648     rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice 
    649     rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           -- 
    650     rn_ice_age      = 30.                 !  lim3 only:      --   age                -- 
    651  
    652     ln_tra_dmp    =.false.                !  open boudaries conditions for tracers 
    653     ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities 
    654     rn_time_dmp   =  1.                   ! Damping time scale in days 
    655     rn_time_dmp_out =  1.                 ! Outflow damping time scale 
    656     nn_rimwidth   = 10                    !  width of the relaxation zone 
    657     ln_vol        = .false.               !  total volume correction (see nn_volctl parameter) 
    658     nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero 
    659 / 
    660 !----------------------------------------------------------------------- 
    661 &nambdy_dta    !  open boundaries - external data                       ("key_bdy") 
    662 !----------------------------------------------------------------------- 
    663 !              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    664 !              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    665    bn_ssh      = 'amm12_bdyT_u2d',         24        , 'sossheig',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    666    bn_u2d      = 'amm12_bdyU_u2d',         24        , 'vobtcrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    667    bn_v2d      = 'amm12_bdyV_u2d',         24        , 'vobtcrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    668    bn_u3d      = 'amm12_bdyU_u3d',         24        , 'vozocrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    669    bn_v3d      = 'amm12_bdyV_u3d',         24        , 'vomecrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    670    bn_tem      = 'amm12_bdyT_tra',         24        , 'votemper',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    671    bn_sal      = 'amm12_bdyT_tra',         24        , 'vosaline',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    672 ! for lim2 
    673 !   bn_frld    = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    674 !   bn_hicif   = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    675 !   bn_hsnif   = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    676 ! for lim3 
    677 !   bn_a_i     = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    678 !   bn_ht_i    = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    679 !   bn_ht_s    = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     '' 
    680  
    681    cn_dir      = 'bdydta/' !  root directory for the location of the bulk files 
    682    ln_full_vel = .false.   !   
    683 / 
    684 !----------------------------------------------------------------------- 
    685 &nambdy_tide   !  tidal forcing at open boundaries 
    686 !----------------------------------------------------------------------- 
    687    filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
    688    ln_bdytide_2ddta = .false.                   ! 
    689    ln_bdytide_conj  = .false.                   !  
    690 / 
    691  
    692 !!====================================================================== 
    693 !!                 ***  Bottom boundary condition  *** 
    694 !!====================================================================== 
    695 !!   nambfr        bottom friction 
    696 !!   nambbc        bottom temperature boundary condition 
    697 !!   nambbl        bottom boundary layer scheme                         ("key_trabbl") 
    698 !!====================================================================== 
    699 ! 
    700 !----------------------------------------------------------------------- 
    701 &nambfr        !   bottom friction                                      (default: linear) 
    702 !----------------------------------------------------------------------- 
    703    nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    704                            !                              = 2 : nonlinear friction 
    705    rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case) 
    706    rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
    707    rn_bfri2_max=    1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T) 
    708    rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    709    rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T 
    710    ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    711    rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
    712    rn_tfri1    =    4.e-4  !  top drag coefficient (linear case) 
    713    rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
    714    rn_tfri2_max=    1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T) 
    715    rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2) 
    716    rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T 
    717    ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file ) 
    718    rn_tfrien   =   50.     !  local multiplying factor of tfr (ln_tfr2d=T) 
    719  
    720    ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true) 
    721    ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    722 / 
    723 !----------------------------------------------------------------------- 
    724 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    725 !----------------------------------------------------------------------- 
    726 !              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    727 !              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    728    sn_qgh      ='geothermal_heating.nc',  -12.       , 'heatflow',   .false.   , .true.  , 'yearly'  ,   ''     ,   ''     ,   '' 
    729    ! 
    730    ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    731    nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux 
    732                            !     = 1 constant flux 
    733                            !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    734    rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    735    cn_dir      = './'      !  root directory for the location of the runoff files  
    736 / 
    737 !----------------------------------------------------------------------- 
    738 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    739 !----------------------------------------------------------------------- 
    740    nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    741    nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
    742    rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
    743    rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
    744 / 
    745  
    746 !!====================================================================== 
    747 !!                        Tracer (T & S ) namelists 
    748 !!====================================================================== 
    749 !!   nameos           equation of state 
    750 !!   namtra_adv       advection scheme 
    751 !!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    752 !!   namtra_ldf       lateral diffusion scheme 
    753 !!   namtra_ldfeiv    eddy induced velocity param. 
    754 !!   namtra_dmp       T & S newtonian damping 
    755 !!====================================================================== 
    756 ! 
    757 !----------------------------------------------------------------------- 
    758 &nameos        !   ocean physical parameters 
    759 !----------------------------------------------------------------------- 
    760    ln_teos10   = .false.         !  = Use TEOS-10 equation of state 
    761    ln_eos80    = .false.         !  = Use EOS80 equation of state 
    762    ln_seos     = .false.         !  = Use simplified equation of state (S-EOS) 
    763                                  ! 
    764    !                     ! S-EOS coefficients (ln_seos=T): 
    765    !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    766    rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
    767    rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
    768    rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
    769    rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
    770    rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos) 
    771    rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos) 
    772    rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos) 
    773 / 
    774 !----------------------------------------------------------------------- 
    775 &namtra_adv    !   advection scheme for tracer                          (default: NO advection) 
    776 !----------------------------------------------------------------------- 
    777    ln_traadv_cen = .false. !  2nd order centered scheme 
    778       nn_cen_h   =  4            !  =2/4, horizontal 2nd order CEN / 4th order CEN 
    779       nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
    780    ln_traadv_fct = .false. !  FCT scheme 
    781       nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order  
    782       nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order  
    783       nn_fct_zts =  0            !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
    784       !                          !        (number of sub-timestep = nn_fct_zts) 
    785    ln_traadv_mus = .false. !  MUSCL scheme 
    786       ln_mus_ups = .false.       !  use upstream scheme near river mouths 
    787    ln_traadv_ubs = .false. !  UBS scheme 
    788       nn_ubs_v   =  2            !  =2  , vertical 2nd order FCT / COMPACT 4th order 
    789    ln_traadv_qck = .false. !  QUICKEST scheme 
    790 / 
    791 !----------------------------------------------------------------------- 
    792 &namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) 
    793 !----------------------------------------------------------------------- 
    794    ln_mle      = .false.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
    795    rn_ce       = 0.06      ! magnitude of the MLE (typical value: 0.06 to 0.08) 
    796    nn_mle      = 1         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
    797    rn_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
    798    rn_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
    799    rn_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
    800    nn_mld_uv   = 0         ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
    801    nn_conv     = 0         ! =1 no MLE in case of convection ; =0 always MLE 
    802    rn_rho_c_mle= 0.01      ! delta rho criterion used to calculate MLD for FK 
    803 / 
    804 !----------------------------------------------------------------------- 
    805 &namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO diffusion) 
    806 !----------------------------------------------------------------------- 
    807    !                       !  Operator type: 
    808    !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
    809    ln_traldf_lap   =  .false.  !    laplacian operator 
    810    ln_traldf_blp   =  .false.  !  bilaplacian operator 
    811    ! 
    812    !                       !  Direction of action: 
    813    ln_traldf_lev   =  .false.  !  iso-level 
    814    ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    815    ln_traldf_iso   =  .false.  !  iso-neutral (standard operator) 
    816    ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
    817    ! 
    818    !                       !  iso-neutral options:         
    819    ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
    820    rn_slpmax       =   0.01    !  slope limit                      (both operators) 
    821    ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
    822    rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
    823    ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
    824    ! 
    825    !                       !  Coefficients: 
    826    nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
    827    !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    828    !                                !   =  0           constant  
    829    !                                !   = 10 F(k)      =ldf_c1d  
    830    !                                !   = 20 F(i,j)    =ldf_c2d  
    831    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    832    !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    833    !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    834    rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
    835    rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    836 / 
    837 !----------------------------------------------------------------------- 
    838 &namtra_ldfeiv !   eddy induced velocity param.                         (default: NO) 
    839 !----------------------------------------------------------------------- 
    840    ln_ldfeiv     =.false.  ! use eddy induced velocity parameterization 
    841    ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities 
    842    rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
    843    nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
    844    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    845    !                                !   =  0           constant  
    846    !                                !   = 10 F(k)      =ldf_c1d  
    847    !                                !   = 20 F(i,j)    =ldf_c2d  
    848    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    849    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    850 / 
    851 !----------------------------------------------------------------------- 
    852 &namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO) 
    853 !----------------------------------------------------------------------- 
    854    ln_tradmp   =  .true.   !  add a damping termn (T) or not (F) 
    855    nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column 
    856                            !                   =1 no damping in the mixing layer (kz  criteria) 
    857                            !                   =2 no damping in the mixed  layer (rho crieria) 
    858    cn_resto    ='resto.nc' !  Name of file containing restoration coeff. field (use dmp_tools to create this) 
    859 / 
    860  
    861 !!====================================================================== 
    862 !!                      ***  Dynamics namelists  *** 
    863 !!====================================================================== 
    864 !!   namdyn_adv    formulation of the momentum advection 
    865 !!   namdyn_vor    advection scheme 
    866 !!   namdyn_hpg    hydrostatic pressure gradient 
    867 !!   namdyn_spg    surface pressure gradient 
    868 !!   namdyn_ldf    lateral diffusion scheme 
    869 !!====================================================================== 
    870 ! 
    871 !----------------------------------------------------------------------- 
    872 &namdyn_adv    !   formulation of the momentum advection                (default: vector form) 
    873 !----------------------------------------------------------------------- 
    874    ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
    875    nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction 
    876    ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    877    ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
    878    ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection 
    879 / 
    880 !----------------------------------------------------------------------- 
    881 &nam_vvl    !   vertical coordinate options                             (default: zstar) 
    882 !----------------------------------------------------------------------- 
    883    ln_vvl_zstar  = .true.           !  zstar vertical coordinate 
    884    ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations 
    885    ln_vvl_layer  = .false.          !  full layer vertical coordinate 
    886    ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar 
    887    ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator 
    888    rn_ahe3       = 0.0e0            !  thickness diffusion coefficient 
    889    rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days] 
    890    rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days] 
    891    rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation 
    892    ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
    893 / 
    894 !----------------------------------------------------------------------- 
    895 &namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO) 
    896 !----------------------------------------------------------------------- 
    897    ln_dynvor_ene = .false. !  enstrophy conserving scheme 
    898    ln_dynvor_ens = .false. !  energy conserving scheme 
    899    ln_dynvor_mix = .false. !  mixed scheme 
    900    ln_dynvor_een = .false. !  energy & enstrophy scheme 
    901       nn_een_e3f = 1          ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    902    ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes)  ! PLEASE DO NOT ACTIVATE 
    903 / 
    904 !----------------------------------------------------------------------- 
    905 &namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: zps) 
    906 !----------------------------------------------------------------------- 
    907    ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    908    ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
    909    ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
    910    ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf 
    911    ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    912    ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    913 / 
    914 !----------------------------------------------------------------------- 
    915 &namdyn_spg    !   surface pressure gradient                            (default: NO) 
    916 !----------------------------------------------------------------------- 
    917    ln_dynspg_exp  = .false.   ! explicit free surface 
    918    ln_dynspg_ts   = .false.   ! split-explicit free surface 
    919       ln_bt_fw      = .true.     ! Forward integration of barotropic Eqs. 
    920       ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    921          nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    922          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    923          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    924       ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    925          rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed 
    926          nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds 
    927 / 
    928 !----------------------------------------------------------------------- 
    929 &namdyn_ldf    !   lateral diffusion on momentum                        (default: NO) 
    930 !----------------------------------------------------------------------- 
    931    !                       !  Type of the operator : 
    932    !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
    933    ln_dynldf_lap =  .false.    !    laplacian operator 
    934    ln_dynldf_blp =  .false.    !  bilaplacian operator 
    935    !                       !  Direction of action  : 
    936    ln_dynldf_lev =  .false.    !  iso-level 
    937    ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
    938    ln_dynldf_iso =  .false.    !  iso-neutral 
    939    !                       !  Coefficient 
    940    nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
    941    !                                !  =-30  read in eddy_viscosity_3D.nc file 
    942    !                                !  =-20  read in eddy_viscosity_2D.nc file 
    943    !                                !  =  0  constant  
    944    !                                !  = 10  F(k)=c1d 
    945    !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
    946    !                                !  = 30  F(i,j,k)=c2d*c1d 
    947    !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
    948    rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
    949    rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
    950    rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    951    ! 
    952    ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    953 / 
    954  
    955 !!====================================================================== 
    956 !!             Tracers & Dynamics vertical physics namelists 
    957 !!====================================================================== 
    958 !!    namzdf        vertical physics 
    959 !!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric") 
    960 !!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke") 
    961 !!    namzdf_gls    GLS vertical mixing                                 ("key_zdfgls") 
    962 !!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm") 
    963 !!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx") 
    964 !!====================================================================== 
    965 ! 
    966 !----------------------------------------------------------------------- 
    967 &namzdf        !   vertical physics 
    968 !----------------------------------------------------------------------- 
    969    rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    970    rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    971    nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0) 
    972    nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0) 
    973    ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F) 
    974       nn_evdm     =    0        ! evd apply on tracer (=0) or on tracer and momentum (=1) 
    975       rn_avevd    =  100.       !  evd mixing coefficient [m2/s] 
    976    ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F) 
    977       nn_npc      =    1        ! frequency of application of npc 
    978       nn_npcp     =  365        ! npc control print frequency 
    979    ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping 
    980       nn_zdfexp   =    3        ! number of sub-timestep for ln_zdfexp=T 
    981 / 
    982 !----------------------------------------------------------------------- 
    983 &namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
    984 !----------------------------------------------------------------------- 
    985    rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity 
    986    rn_alp      =    5.     !  coefficient of the parameterization 
    987    nn_ric      =    2      !  coefficient of the parameterization 
    988    rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation 
    989    rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m) 
    990    rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m) 
    991    rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer 
    992    rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer 
    993    ln_mldw     =  .true.   !  Flag to use or not the mixed layer depth param. 
    994 / 
    995 !----------------------------------------------------------------------- 
    996 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    997 !----------------------------------------------------------------------- 
    998    rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
    999    rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation 
    1000    rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) 
    1001    rn_emin     =   1.e-6   !  minimum value of tke [m2/s2] 
    1002    rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2] 
    1003    rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it) 
    1004    nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
    1005                            !                 = 1 bounded by the local vertical scale factor 
    1006                            !                 = 2 first vertical derivative of mixing length bounded by 1 
    1007                            !                 = 3 as =2 with distinct disspipative an mixing length scale 
    1008    nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
    1009    ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
    1010    rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
    1011    ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
    1012    rn_lc       =   0.15    !  coef. associated to Langmuir cells 
    1013    nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to near intertial waves 
    1014                            !        = 0 no penetration 
    1015                            !        = 1 add a tke source below the ML 
    1016                            !        = 2 add a tke source just at the base of the ML 
    1017                            !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T) 
    1018    rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) 
    1019    nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML 
    1020                            !        = 0  constant 10 m length scale 
    1021                            !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    1022 / 
    1023 !----------------------------------------------------------------------- 
    1024 &namzdf_gls    !   GLS vertical diffusion                               ("key_zdfgls") 
    1025 !----------------------------------------------------------------------- 
    1026    rn_emin       = 1.e-7   !  minimum value of e   [m2/s2] 
    1027    rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3] 
    1028    ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988) 
    1029    rn_clim_galp  = 0.267   !  galperin limit 
    1030    ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case 
    1031    rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux 
    1032    rn_charn      = 70000.  !  Charnock constant for wb induced roughness length 
    1033    rn_hsro       =  0.02   !  Minimum surface roughness 
    1034    rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2) 
    1035    nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2) 
    1036    nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum) 
    1037    nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum) 
    1038    nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) 
    1039    nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) 
    1040 / 
    1041 !----------------------------------------------------------------------- 
    1042 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    1043 !----------------------------------------------------------------------- 
    1044    rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity) 
    1045    rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio 
    1046 / 
    1047 !----------------------------------------------------------------------- 
    1048 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    1049 !----------------------------------------------------------------------- 
    1050    rn_htmx     = 500.      !  vertical decay scale for turbulence (meters) 
    1051    rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1) 
    1052    rn_tfe      = 0.333     !  tidal dissipation efficiency 
    1053    rn_me       = 0.2       !  mixing efficiency 
    1054    ln_tmx_itf  = .true.    !  ITF specific parameterisation 
    1055    rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency 
    1056 / 
    1057 !----------------------------------------------------------------------- 
    1058 &namzdf_tmx_new !   internal wave-driven mixing parameterization        ("key_zdftmx_new" & "key_zdfddm") 
    1059 !----------------------------------------------------------------------- 
    1060    nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) 
    1061    ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
    1062    ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    1063 / 
    1064  
    1065  
    1066 !!====================================================================== 
    1067 !!                  ***  Miscellaneous namelists  *** 
    1068 !!====================================================================== 
    1069 !!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    1070 !!   namctl            Control prints & Benchmark 
    1071 !!   namsto            Stochastic parametrization of EOS 
    1072 !!====================================================================== 
    1073 ! 
    1074 !----------------------------------------------------------------------- 
    1075 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     184   ln_chk_bathy  = .false. !  =T  check the parent bathymetry 
     185/ 
     186!----------------------------------------------------------------------- 
     187&nambdy        !  unstructured open boundaries                          (default: OFF) 
     188!----------------------------------------------------------------------- 
     189   ln_bdy         = .false.   !  Use unstructured open boundaries 
     190   nb_bdy         = 0         !  number of open boundary sets 
     191   ln_coords_file = .true.    !  =T : read bdy coordinates from file 
     192      cn_coords_file = 'coordinates.bdy.nc'  !  bdy coordinates files 
     193   ln_mask_file   = .false.   !  =T : read mask from file 
     194      cn_mask_file = ''        !  name of mask file (if ln_mask_file=.TRUE.) 
     195   cn_dyn2d    = 'none'       ! 
     196   nn_dyn2d_dta   =  0        !  = 0, bdy data are equal to the initial state 
     197      !                       !  = 1, bdy data are read in 'bdydata   .nc' files 
     198      !                       !  = 2, use tidal harmonic forcing data from files 
     199      !                       !  = 3, use external data AND tidal harmonic forcing 
     200   cn_dyn3d      =  'none'    ! 
     201   nn_dyn3d_dta  =  0         !  = 0, bdy data are equal to the initial state 
     202   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     203   cn_tra        =  'none'    ! 
     204   nn_tra_dta    =  0         !  = 0, bdy data are equal to the initial state 
     205   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     206   cn_ice        =  'none'    ! 
     207   nn_ice_dta    =  0         !  = 0, bdy data are equal to the initial state 
     208   !                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     209   rn_ice_tem    = 270.       !  si3 only: arbitrary temperature of incoming sea ice 
     210   rn_ice_sal    = 10.        !  si3 only:      --   salinity           -- 
     211   rn_ice_age    = 30.        !  si3 only:      --   age                -- 
     212   ! 
     213   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers 
     214   ln_dyn3d_dmp  =.false.     !  open boundary condition for baroclinic velocities 
     215   rn_time_dmp   =  1.        !  Damping time scale in days 
     216   rn_time_dmp_out = 1.        !  Outflow damping time scale 
     217   nn_rimwidth   = 10         !  width of the relaxation zone 
     218   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter) 
     219   nn_volctl     =  1         !  = 0, the total water flux across open boundaries is zero 
     220   nb_jpk_bdy    = -1         ! number of levels in the bdy data (set < 0 if consistent with planned run) 
     221/ 
     222!----------------------------------------------------------------------- 
     223&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
     224!----------------------------------------------------------------------- 
     225   nn_nchunks_i =   4       !  number of chunks in i-dimension 
     226   nn_nchunks_j =   4       !  number of chunks in j-dimension 
     227   nn_nchunks_k =   31      !  number of chunks in k-dimension 
     228   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
     229   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs 
     230   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
     231   !                       !  (F) ignore chunking information and produce netcdf3-compatible files 
     232/ 
     233!----------------------------------------------------------------------- 
     234&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    1076235!----------------------------------------------------------------------- 
    1077236   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send, 
    1078                            !  buffer blocking send or immediate non-blocking sends, resp. 
     237   !                       !  buffer blocking send or immediate non-blocking sends, resp. 
    1079238   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation 
    1080    ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold 
     239   ln_nnogather =  .true.  !  activate code to avoid mpi_allgather use at the northfold 
    1081240   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1) 
    1082241   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1) 
    1083    jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1) 
    1084 / 
    1085 !----------------------------------------------------------------------- 
    1086 &namctl        !   Control prints & Benchmark 
    1087 !----------------------------------------------------------------------- 
    1088    ln_ctl      = .false.   !  trends control print (expensive!) 
     242/ 
     243!----------------------------------------------------------------------- 
     244&namctl        !   Control prints                                       (default: OFF) 
     245!----------------------------------------------------------------------- 
     246   ln_ctl = .FALSE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 
     247     sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
     248       sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 
     249       sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
     250       sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
     251       sn_cfctl%l_layout  = .FALSE. ! 
     252       sn_cfctl%l_mppout  = .FALSE. ! 
     253       sn_cfctl%l_mpptop  = .FALSE. ! 
     254       sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
     255       sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
     256       sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1] 
     257       sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info 
    1089258   nn_print    =    0      !  level of print (0 no extra print) 
    1090259   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus 
     
    1094263   nn_isplt    =    1      !  number of processors in i-direction 
    1095264   nn_jsplt    =    1      !  number of processors in j-direction 
    1096    nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench 
    1097                            !     (no physical validity of the results) 
    1098    nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    1099    nn_diacfl   =    0      !  Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 
    1100 / 
    1101 !----------------------------------------------------------------------- 
    1102 &namsto        ! Stochastic parametrization of EOS                      (default: NO) 
    1103 !----------------------------------------------------------------------- 
    1104    ln_sto_eos  = .false.   ! stochastic equation of state 
    1105    nn_sto_eos  = 1         ! number of independent random walks 
    1106    rn_eos_stdxy= 1.4       ! random walk horz. standard deviation (in grid points) 
    1107    rn_eos_stdz = 0.7       ! random walk vert. standard deviation (in grid points) 
    1108    rn_eos_tcor = 1440.     ! random walk time correlation (in timesteps) 
    1109    nn_eos_ord  = 1         ! order of autoregressive processes 
    1110    nn_eos_flt  = 0         ! passes of Laplacian filter 
    1111    rn_eos_lim  = 2.0       ! limitation factor (default = 3.0) 
    1112    ln_rststo   = .false.   ! start from mean parameter (F) or from restart file (T) 
    1113    ln_rstseed  = .true.    ! read seed of RNG from restart file 
    1114    cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input) 
    1115    cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output) 
    1116 / 
    1117  
    1118 !!====================================================================== 
    1119 !!                  ***  Diagnostics namelists  *** 
    1120 !!====================================================================== 
    1121 !!   namtrd       dynamics and/or tracer trends                         (default F) 
    1122 !!   namptr       Poleward Transport Diagnostics                        (default F) 
    1123 !!   namhsb       Heat and salt budgets                                 (default F) 
    1124 !!   namdiu       Cool skin and warm layer models                       (default F) 
    1125 !!   namflo       float parameters                                      ("key_float") 
    1126 !!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
    1127 !!   namdct       transports through some sections                      ("key_diadct") 
    1128 !!   nam_dia25h   25h Mean Output                                       (default F) 
    1129 !!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    1130 !!====================================================================== 
    1131 ! 
    1132 !----------------------------------------------------------------------- 
    1133 &namtrd        !   trend diagnostics                                    (default F) 
    1134 !----------------------------------------------------------------------- 
    1135    ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
    1136    ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
    1137    ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
    1138    ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet) 
    1139    ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
    1140    ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
    1141    ln_tra_trd  = .false.   ! (T) 3D tracer trend output 
    1142    ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
    1143    nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
    1144 / 
    1145 !!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
    1146 !!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
    1147 !!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
    1148 !!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
    1149 !!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics 
    1150 !!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    1151 !!gm 
    1152 !----------------------------------------------------------------------- 
    1153 &namptr        !   Poleward Transport Diagnostic                         (default F) 
    1154 !----------------------------------------------------------------------- 
    1155    ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F) 
    1156    ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not 
    1157 / 
    1158 !----------------------------------------------------------------------- 
    1159 &namhsb        !  Heat and salt budgets                                  (default F) 
    1160 !----------------------------------------------------------------------- 
    1161    ln_diahsb   = .false.   !  check the heat and salt budgets (T) or not (F) 
    1162 / 
    1163 !----------------------------------------------------------------------- 
    1164 &namdiu        !   Cool skin and warm layer models                       (default F) 
    1165 !----------------------------------------------------------------------- 
    1166    ln_diurnal      = .false.   !  
    1167    ln_diurnal_only = .false.   ! 
    1168 / 
    1169 !----------------------------------------------------------------------- 
    1170 &namflo        !   float parameters                                      ("key_float") 
    1171 !----------------------------------------------------------------------- 
    1172    jpnfl       = 1         !  total number of floats during the run 
    1173    jpnnewflo   = 0         !  number of floats for the restart 
    1174    ln_rstflo   = .false.   !  float restart (T) or not (F) 
    1175    nn_writefl  =      75   !  frequency of writing in float output file 
    1176    nn_stockfl  =    5475   !  frequency of creation of the float restart file 
    1177    ln_argo     = .false.   !  Argo type floats (stay at the surface each 10 days) 
    1178    ln_flork4   = .false.   !  trajectories computed with a 4th order Runge-Kutta (T) 
    1179    !                       !  or computed with Blanke' scheme (F) 
    1180    ln_ariane   = .true.    !  Input with Ariane tool convention(T) 
    1181    ln_flo_ascii= .true.    !  Output with Ariane tool netcdf convention(F) or ascii file (T) 
    1182 / 
    1183 !----------------------------------------------------------------------- 
    1184 &nam_diaharm   !   Harmonic analysis of tidal constituents               ("key_diaharm") 
    1185 !----------------------------------------------------------------------- 
    1186     nit000_han = 1         ! First time step used for harmonic analysis 
    1187     nitend_han = 75        ! Last time step used for harmonic analysis 
    1188     nstep_han  = 15        ! Time step frequency for harmonic analysis 
    1189     tname(1)   = 'M2'      ! Name of tidal constituents 
    1190     tname(2)   = 'K1' 
    1191 / 
    1192 !----------------------------------------------------------------------- 
    1193 &namdct        ! transports through some sections                        ("key_diadct") 
    1194 !----------------------------------------------------------------------- 
    1195     nn_dct     = 15        !  time step frequency for transports computing 
    1196     nn_dctwri  = 15        !  time step frequency for transports writing 
    1197     nn_secdebug= 112       !      0 : no section to debug 
    1198     !                      !     -1 : debug all section 
    1199     !                      !  0 < n : debug section number n 
    1200 / 
    1201 !----------------------------------------------------------------------- 
    1202 &nam_dia25h    !  25h Mean Output                                        (default F) 
    1203 !----------------------------------------------------------------------- 
    1204    ln_dia25h   = .false.   ! Choose 25h mean output or not 
    1205 / 
    1206 !----------------------------------------------------------------------- 
    1207 &namnc4        !   netcdf4 chunking and compression settings             ("key_netcdf4") 
    1208 !----------------------------------------------------------------------- 
    1209    nn_nchunks_i=   4       !  number of chunks in i-dimension 
    1210    nn_nchunks_j=   4       !  number of chunks in j-dimension 
    1211    nn_nchunks_k=   31      !  number of chunks in k-dimension 
    1212    !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
    1213    !                       !  is optimal for postprocessing which works exclusively with horizontal slabs 
    1214    ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
    1215    !                       !  (F) ignore chunking information and produce netcdf3-compatible files 
    1216 / 
    1217  
    1218 !!====================================================================== 
    1219 !!               ***  Observation & Assimilation  *** 
    1220 !!====================================================================== 
    1221 !!   namobs       observation and model comparison 
    1222 !!   nam_asminc   assimilation increments                               ('key_asminc') 
    1223 !!====================================================================== 
    1224 ! 
    1225 !----------------------------------------------------------------------- 
    1226 &namobs        !  observation usage switch 
    1227 !----------------------------------------------------------------------- 
    1228    ln_diaobs   = .false.             ! Logical switch for the observation operator 
    1229    ln_t3d      = .false.             ! Logical switch for T profile observations 
    1230    ln_s3d      = .false.             ! Logical switch for S profile observations 
    1231    ln_sla      = .false.             ! Logical switch for SLA observations 
    1232    ln_sst      = .false.             ! Logical switch for SST observations 
    1233    ln_sic      = .false.             ! Logical switch for Sea Ice observations 
    1234    ln_vel3d    = .false.             ! Logical switch for velocity observations 
    1235    ln_altbias  = .false.             ! Logical switch for altimeter bias correction 
    1236    ln_nea      = .false.             ! Logical switch for rejection of observations near land 
    1237    ln_grid_global = .true.           ! Logical switch for global distribution of observations 
    1238    ln_grid_search_lookup = .false.   ! Logical switch for obs grid search w/lookup table 
    1239    ln_ignmis   = .true.              ! Logical switch for ignoring missing files 
    1240    ln_s_at_t   = .false.             ! Logical switch for computing model S at T obs if not there 
    1241    ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs 
    1242 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files 
    1243    cn_profbfiles = 'profiles_01.nc'  ! Profile feedback input observation file names 
    1244    cn_slafbfiles = 'sla_01.nc'       ! SLA feedback input observation file names 
    1245    cn_sstfbfiles = 'sst_01.nc'       ! SST feedback input observation file names 
    1246    cn_sicfbfiles = 'sic_01.nc'       ! SIC feedback input observation file names 
    1247    cn_velfbfiles = 'vel_01.nc'       ! Velocity feedback input observation file names 
    1248    cn_altbiasfile = 'altbias.nc'     ! Altimeter bias input file name 
    1249    cn_gridsearchfile='gridsearch.nc' ! Grid search file name 
    1250    rn_gridsearchres = 0.5            ! Grid search resolution 
    1251    rn_dobsini  = 00010101.000000     ! Initial date in window YYYYMMDD.HHMMSS 
    1252    rn_dobsend  = 00010102.000000     ! Final date in window YYYYMMDD.HHMMSS 
    1253    nn_1dint    = 0                   ! Type of vertical interpolation method 
    1254    nn_2dint    = 0                   ! Type of horizontal interpolation method 
    1255    nn_msshc    = 0                   ! MSSH correction scheme 
    1256    rn_mdtcorr  = 1.61                ! MDT  correction 
    1257    rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction 
    1258    nn_profdavtypes = -1              ! Profile daily average types - array 
    1259    ln_sstbias  = .false.             ! 
    1260    cn_sstbias_files = 'sstbias.nc'   ! 
    1261 / 
    1262 !----------------------------------------------------------------------- 
    1263 &nam_asminc    !   assimilation increments                              ('key_asminc') 
    1264 !----------------------------------------------------------------------- 
    1265     ln_bkgwri  = .false.   !  Logical switch for writing out background state 
    1266     ln_trainc  = .false.   !  Logical switch for applying tracer increments 
    1267     ln_dyninc  = .false.   !  Logical switch for applying velocity increments 
    1268     ln_sshinc  = .false.   !  Logical switch for applying SSH increments 
    1269     ln_asmdin  = .false.   !  Logical switch for Direct Initialization (DI) 
    1270     ln_asmiau  = .false.   !  Logical switch for Incremental Analysis Updating (IAU) 
    1271     nitbkg     = 0         !  Timestep of background in [0,nitend-nit000-1] 
    1272     nitdin     = 0         !  Timestep of background for DI in [0,nitend-nit000-1] 
    1273     nitiaustr  = 1         !  Timestep of start of IAU interval in [0,nitend-nit000-1] 
    1274     nitiaufin  = 15        !  Timestep of end of IAU interval in [0,nitend-nit000-1] 
    1275     niaufn     = 0         !  Type of IAU weighting function 
    1276     ln_salfix  = .false.   !  Logical switch for ensuring that the sa > salfixmin 
    1277     salfixmin  = -9999     !  Minimum salinity after applying the increments 
    1278     nn_divdmp  = 0         !  Number of iterations of divergence damping operator 
    1279 / 
     265   ln_timing   = .false.   !  timing by routine write out in timing.output file 
     266   ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii 
     267/ 
  • utils/tools/DOMAINcfg/src/domcfg.f90

    r9598 r12414  
    1515   USE in_out_manager  ! I/O manager 
    1616   USE lib_mpp         ! distributed memory computing library 
    17    USE timing          ! Timing 
    1817 
    1918   IMPLICIT NONE 
     
    3635      !! 
    3736      !!---------------------------------------------------------------------- 
    38       ! 
    39       IF( nn_timing == 1 )  CALL timing_start('dom_cfg') 
    4037      ! 
    4138      IF(lwp) THEN                   ! Control print 
     
    6057      CALL dom_glo                   ! global domain versus zoom and/or local domain 
    6158      ! 
    62       IF( nn_timing == 1 )  CALL timing_stop('dom_cfg') 
    63       ! 
    6459   END SUBROUTINE dom_cfg 
    65  
    6660 
    6761   SUBROUTINE dom_glo 
     
    6963      !!                     ***  ROUTINE dom_glo  *** 
    7064      !! 
    71       !! ** Purpose :   initialization for global domain, zoom and local domain 
     65      !! ** Purpose :   initialization of global domain <--> local domain indices 
    7266      !! 
    7367      !! ** Method  :    
    7468      !! 
    75       !! ** Action  : - mig  , mjg :  
    76       !!              - mi0  , mi1   : 
    77       !!              - mj0, , mj1   : 
     69      !! ** Action  : - mig , mjg : local  domain indices ==> global domain indices 
     70      !!              - mi0 , mi1 : global domain indices ==> local  domain indices 
     71      !!              - mj0,, mj1   (global point not in the local domain ==> mi0>mi1 and/or mj0>mj1) 
    7872      !!---------------------------------------------------------------------- 
    7973      INTEGER ::   ji, jj   ! dummy loop argument 
    8074      !!---------------------------------------------------------------------- 
    81       !                              ! recalculate jpizoom/jpjzoom given lat/lon 
    8275      ! 
    83       !                        ! ============== ! 
    84       !                        !  Local domain  !  
    85       !                        ! ============== ! 
    86       DO ji = 1, jpi                 ! local domain indices ==> data domain indices 
    87         mig(ji) = ji + jpizoom - 1 + nimpp - 1 
     76      DO ji = 1, jpi                 ! local domain indices ==> global domain indices 
     77        mig(ji) = ji + nimpp - 1 
    8878      END DO 
    8979      DO jj = 1, jpj 
    90         mjg(jj) = jj + jpjzoom - 1 + njmpp - 1 
     80        mjg(jj) = jj + njmpp - 1 
    9181      END DO 
    92       ! 
    93       !                              ! data domain indices ==> local domain indices 
     82      !                              ! global domain indices ==> local domain indices 
    9483      !                                   ! (return (m.0,m.1)=(1,0) if data domain gridpoint is to the west/south of the  
    95       !                                   !local domain, or (m.0,m.1)=(jp.+1,jp.) to the east/north of local domain.  
    96       DO ji = 1, jpidta 
    97         mi0(ji) = MAX( 1, MIN( ji - jpizoom + 1 - nimpp + 1, jpi+1 ) ) 
    98         mi1(ji) = MAX( 0, MIN( ji - jpizoom + 1 - nimpp + 1, jpi   ) ) 
     84      !                                   ! local domain, or (m.0,m.1)=(jp.+1,jp.) to the east/north of local domain.  
     85      DO ji = 1, jpiglo 
     86        mi0(ji) = MAX( 1 , MIN( ji - nimpp + 1, jpi+1 ) ) 
     87        mi1(ji) = MAX( 0 , MIN( ji - nimpp + 1, jpi   ) ) 
    9988      END DO 
    100       DO jj = 1, jpjdta 
    101         mj0(jj) = MAX( 1, MIN( jj - jpjzoom + 1 - njmpp + 1, jpj+1 ) ) 
    102         mj1(jj) = MAX( 0, MIN( jj - jpjzoom + 1 - njmpp + 1, jpj   ) ) 
     89      DO jj = 1, jpjglo 
     90        mj0(jj) = MAX( 1 , MIN( jj - njmpp + 1, jpj+1 ) ) 
     91        mj1(jj) = MAX( 0 , MIN( jj - njmpp + 1, jpj   ) ) 
    10392      END DO 
    10493      IF(lwp) THEN                   ! control print 
    10594         WRITE(numout,*) 
    106          WRITE(numout,*) 'dom_glo : domain: data / local ' 
     95         WRITE(numout,*) 'dom_glo : domain: global <<==>> local ' 
    10796         WRITE(numout,*) '~~~~~~~ ' 
    108          WRITE(numout,*) '          data input domain    : jpidta = ', jpidta,   & 
    109             &                                            ' jpjdta = ', jpjdta, ' jpkdta = ', jpkdta 
    110          WRITE(numout,*) '          global or zoom domain: jpiglo = ', jpiglo,   & 
    111             &                                            ' jpjglo = ', jpjglo, ' jpk    = ', jpk 
    112          WRITE(numout,*) '          local domain         : jpi    = ', jpi   ,   & 
    113             &                                            ' jpj    = ', jpj   , ' jpk    = ', jpk 
     97         WRITE(numout,*) '   global domain:   jpiglo = ', jpiglo, ' jpjglo = ', jpjglo, ' jpkglo = ', jpkglo 
     98         WRITE(numout,*) '   local  domain:   jpi    = ', jpi   , ' jpj    = ', jpj   , ' jpk    = ', jpk 
    11499         WRITE(numout,*) 
    115          WRITE(numout,*) '          south-west indices    jpizoom = ', jpizoom,   & 
    116             &                                           ' jpjzoom = ', jpjzoom 
     100         WRITE(numout,*) '   conversion from local to global domain indices (and vise versa) done' 
    117101         IF( nn_print >= 1 ) THEN 
    118102            WRITE(numout,*) 
    119             WRITE(numout,*) '          conversion local  ==> data i-index domain' 
     103            WRITE(numout,*) '          conversion local  ==> global i-index domain (mig)' 
    120104            WRITE(numout,25)              (mig(ji),ji = 1,jpi) 
    121105            WRITE(numout,*) 
    122             WRITE(numout,*) '          conversion data  ==> local  i-index domain' 
    123             WRITE(numout,*) '             starting index' 
    124             WRITE(numout,25)              (mi0(ji),ji = 1,jpidta) 
    125             WRITE(numout,*) '             ending index' 
    126             WRITE(numout,25)              (mi1(ji),ji = 1,jpidta) 
     106            WRITE(numout,*) '          conversion global ==> local  i-index domain' 
     107            WRITE(numout,*) '             starting index (mi0)' 
     108            WRITE(numout,25)              (mi0(ji),ji = 1,jpiglo) 
     109            WRITE(numout,*) '             ending index (mi1)' 
     110            WRITE(numout,25)              (mi1(ji),ji = 1,jpiglo) 
    127111            WRITE(numout,*) 
    128             WRITE(numout,*) '          conversion local  ==> data j-index domain' 
     112            WRITE(numout,*) '          conversion local  ==> global j-index domain (mjg)' 
    129113            WRITE(numout,25)              (mjg(jj),jj = 1,jpj) 
    130114            WRITE(numout,*) 
    131             WRITE(numout,*) '          conversion data  ==> local j-index domain' 
    132             WRITE(numout,*) '             starting index' 
    133             WRITE(numout,25)              (mj0(jj),jj = 1,jpjdta) 
    134             WRITE(numout,*) '             ending index' 
    135             WRITE(numout,25)              (mj1(jj),jj = 1,jpjdta) 
     115            WRITE(numout,*) '          conversion global ==> local j-index domain' 
     116            WRITE(numout,*) '             starting index (mj0)' 
     117            WRITE(numout,25)              (mj0(jj),jj = 1,jpjglo) 
     118            WRITE(numout,*) '             ending index (mj1)' 
     119            WRITE(numout,25)              (mj1(jj),jj = 1,jpjglo) 
    136120         ENDIF 
    137121      ENDIF 
    138122 25   FORMAT( 100(10x,19i4,/) ) 
    139  
    140       !                        ! ============== ! 
    141       !                        !  Zoom domain   ! 
    142       !                        ! ============== ! 
    143       !                              ! zoom control 
    144       IF( jpiglo + jpizoom - 1  >  jpidta .OR.   & 
    145           jpjglo + jpjzoom - 1  >  jpjdta      ) & 
    146           &   CALL ctl_stop( ' global or zoom domain exceed the data domain ! ' ) 
    147  
    148       !                              ! set zoom flag 
    149       IF( jpiglo < jpidta .OR. jpjglo < jpjdta )   lzoom = .TRUE. 
    150  
    151       !                              ! set zoom type flags 
    152       IF( lzoom .AND. jpizoom /= 1 )   lzoom_w = .TRUE.                     !  
    153       IF( lzoom .AND. jpjzoom /= 1 )   lzoom_s = .TRUE. 
    154       IF( lzoom .AND. jpiglo + jpizoom -1 /= jpidta )   lzoom_e = .TRUE. 
    155       IF( lzoom .AND. jpjglo + jpjzoom -1 /= jpjdta )   lzoom_n = .TRUE. 
    156       IF(lwp) THEN 
    157          WRITE(numout,*) 
    158          WRITE(numout,*) '          zoom flags : ' 
    159          WRITE(numout,*) '             lzoom   = ', lzoom  , ' (T = zoom, F = global )' 
    160          WRITE(numout,*) '             lzoom_e = ', lzoom_e, ' (T = forced closed east  boundary)' 
    161          WRITE(numout,*) '             lzoom_w = ', lzoom_w, ' (T = forced closed west  boundary)' 
    162          WRITE(numout,*) '             lzoom_s = ', lzoom_s, ' (T = forced closed South boundary)' 
    163          WRITE(numout,*) '             lzoom_n = ', lzoom_n, ' (T = forced closed North boundary)' 
    164       ENDIF 
    165       IF(  ( lzoom_e .OR. lzoom_w )  .AND.  ( jperio == 1 .OR. jperio == 4 .OR. jperio == 6 )  )   & 
    166            &   CALL ctl_stop( ' Your zoom choice is inconsistent with east-west cyclic boundary condition' ) 
    167       IF(  lzoom_n  .AND.  ( 3 <= jperio .AND. jperio <= 6 )  )   & 
    168            &   CALL ctl_stop( ' Your zoom choice is inconsistent with North fold boundary condition' ) 
    169  
    170       !                              ! Pre-defined arctic/antarctic zoom of ORCA configuration flag 
    171       IF( cp_cfg == "orca" ) THEN 
    172          SELECT CASE ( jp_cfg ) 
    173          CASE ( 2 )                               !  ORCA_R2 configuration 
    174             IF(  cp_cfz == "arctic"    .AND. jpiglo  == 142    .AND. jpjglo  ==  53 .AND.   & 
    175                & jpizoom ==  21    .AND. jpjzoom ==  97         )   THEN 
    176               IF(lwp) WRITE(numout,*) '          ORCA configuration: arctic zoom ' 
    177             ENDIF 
    178             IF(  cp_cfz == "antarctic" .AND. jpiglo  == jpidta .AND. jpjglo  ==  50 .AND.   & 
    179                & jpizoom ==   1    .AND. jpjzoom ==   1         )   THEN 
    180               IF(lwp) WRITE(numout,*) '          ORCA configuration: antarctic zoom ' 
    181             ENDIF 
    182             !                              
    183          CASE ( 05 )                              !  ORCA_R05 configuration 
    184             IF(    cp_cfz == "arctic"    .AND. jpiglo  == 562    .AND. jpjglo  == 202 .AND.   & 
    185                & jpizoom ==  81    .AND. jpjzoom == 301         )   THEN 
    186               IF(lwp) WRITE(numout,*) '          ORCA configuration: arctic zoom ' 
    187             ENDIF 
    188             IF(    cp_cfz == "antarctic" .AND. jpiglo  == jpidta .AND. jpjglo  == 187 .AND.   & 
    189                & jpizoom ==   1    .AND. jpjzoom ==   1         )   THEN 
    190               IF(lwp) WRITE(numout,*) '          ORCA configuration: antarctic zoom ' 
    191             ENDIF 
    192          END SELECT 
    193          ! 
    194       ENDIF 
    195123      ! 
    196124   END SUBROUTINE dom_glo 
    197  
    198125   !!====================================================================== 
    199126END MODULE domcfg 
  • utils/tools/DOMAINcfg/src/ioipsl.f90

    r6951 r12414  
    66! See IOIPSL/IOIPSL_License_CeCILL.txt 
    77! 
    8   USE errioipsl  
     8  USE errioipsl    
     9  USE calendar    
    910  USE stringop 
    10   USE mathelp     
    11   USE getincom 
    12   USE calendar    
    1311  USE fliocom     
    14   USE flincom     
    15   USE histcom     
    16   USE restcom 
     12 
    1713END MODULE ioipsl 
  • utils/tools/DOMAINcfg/src/par_oce.f90

    r9598 r12414  
    1313   PUBLIC 
    1414 
     15   CHARACTER(lc) ::   cp_cfg           !: name of the configuration 
     16   CHARACTER(lc) ::   cp_cfz           !: name of the zoom of configuration 
     17   INTEGER       ::   jp_cfg           !: resolution of the configuration 
     18 
     19   ! data size                                       !!! * size of all input files * 
     20   INTEGER       ::   jpidta           !: 1st lateral dimension ( >= jpi ) 
     21   INTEGER       ::   jpjdta           !: 2nd    "         "    ( >= jpj ) 
     22   INTEGER       ::   jpkdta           !: number of levels      ( >= jpk ) 
     23   LOGICAL       ::   ln_e3_dep        ! e3. definition flag 
     24   REAL(wp)      ::   pp_not_used       = 999999._wp   !: vertical grid parameter 
     25   REAL(wp)      ::   pp_to_be_computed = 999999._wp   !:    -      -       - 
     26   !!---------------------------------------------------------------------- 
     27   !!                   namcfg namelist parameters 
     28   !!---------------------------------------------------------------------- 
     29   LOGICAL       ::   ln_read_cfg      !: (=T) read the domain configuration file or (=F) not 
     30   CHARACTER(lc) ::      cn_domcfg        !: filename the configuration file to be read 
     31   LOGICAL       ::   ln_write_cfg     !: (=T) create the domain configuration file 
     32   CHARACTER(lc) ::      cn_domcfg_out    !: filename the configuration file to be read 
     33   ! 
     34   LOGICAL       ::   ln_use_jattr     !: input file read offset 
     35   !                                   !  Use file global attribute: open_ocean_jstart to determine start j-row  
     36   !                                   !  when reading input from those netcdf files that have the  
     37   !                                   !  attribute defined. This is designed to enable input files associated  
     38   !                                   !  with the extended grids used in the under ice shelf configurations to  
     39   !                                   !  be used without redundant rows when the ice shelves are not in use. 
     40   !  
     41 
     42   !!--------------------------------------------------------------------- 
     43   !! Domain Matrix size  
     44   !!--------------------------------------------------------------------- 
     45   ! configuration name & resolution   (required only in ORCA family case) 
     46   CHARACTER(lc) ::   cn_cfg           !: name of the configuration 
     47   INTEGER       ::   nn_cfg           !: resolution of the configuration  
     48 
     49   ! global domain size               !!! * total computational domain * 
     50   INTEGER       ::   jpiglo           !: 1st dimension of global domain --> i-direction 
     51   INTEGER       ::   jpjglo           !: 2nd    -                  -    --> j-direction 
     52   INTEGER       ::   jpkglo           !: 3nd    -                  -    --> k levels 
     53 
     54   ! global domain size for AGRIF     !!! * total AGRIF computational domain * 
     55   INTEGER, PUBLIC            ::   nbug_in_agrif_conv_do_not_remove_or_modify = 1 - 1 
     56   INTEGER, PUBLIC, PARAMETER ::   nbghostcells = 3                             !: number of ghost cells 
     57   INTEGER, PUBLIC            ::   nbcellsx   ! = jpiglo - 2 - 2*nbghostcells   !: number of cells in i-direction 
     58   INTEGER, PUBLIC            ::   nbcellsy   ! = jpjglo - 2 - 2*nbghostcells   !: number of cells in j-direction 
     59 
     60   ! local domain size                !!! * local computational domain * 
     61   INTEGER, PUBLIC ::   jpi   !                                                    !: first  dimension 
     62   INTEGER, PUBLIC ::   jpj   !                                                    !: second dimension 
     63   INTEGER, PUBLIC ::   jpk   ! = jpkglo                                           !: third  dimension 
     64   INTEGER, PUBLIC ::   jpim1 ! = jpi-1                                            !: inner domain indices 
     65   INTEGER, PUBLIC ::   jpjm1 ! = jpj-1                                            !:   -     -      - 
     66   INTEGER, PUBLIC ::   jpkm1 ! = jpk-1                                            !:   -     -      - 
     67   INTEGER, PUBLIC ::   jpij  ! = jpi*jpj                                          !:  jpi x jpj 
     68   INTEGER, PUBLIC ::   jpimax! = ( jpiglo-2*nn_hls + (jpni-1) ) / jpni + 2*nn_hls !: maximum jpi 
     69   INTEGER, PUBLIC ::   jpjmax! = ( jpjglo-2*nn_hls + (jpnj-1) ) / jpnj + 2*nn_hls !: maximum jpj 
     70 
     71   !!--------------------------------------------------------------------- 
     72   !! Active tracer parameters 
     73   !!--------------------------------------------------------------------- 
     74   INTEGER, PUBLIC, PARAMETER ::   jpts   = 2    !: Number of active tracers (=2, i.e. T & S ) 
     75   INTEGER, PUBLIC, PARAMETER ::   jp_tem = 1    !: indice for temperature 
     76   INTEGER, PUBLIC, PARAMETER ::   jp_sal = 2    !: indice for salinity 
     77 
    1578   !!---------------------------------------------------------------------- 
    1679   !!   Domain decomposition 
     
    2285   INTEGER, PUBLIC, PARAMETER ::   jpr2di = 0   !: number of columns for extra outer halo  
    2386   INTEGER, PUBLIC, PARAMETER ::   jpr2dj = 0   !: number of rows    for extra outer halo  
    24    INTEGER, PUBLIC, PARAMETER ::   jpreci = 1   !: number of columns for overlap  
    25    INTEGER, PUBLIC, PARAMETER ::   jprecj = 1   !: number of rows    for overlap  
    26  
    27    !!---------------------------------------------------------------------- 
    28    !!                   namcfg namelist parameters 
    29    !!---------------------------------------------------------------------- 
    30    ! 
    31    LOGICAL       ::   ln_e3_dep        ! e3. definition flag 
    32    ! 
    33    CHARACTER(lc) ::   cp_cfg           !: name of the configuration 
    34    CHARACTER(lc) ::   cp_cfz           !: name of the zoom of configuration 
    35    INTEGER       ::   jp_cfg           !: resolution of the configuration 
    36  
    37    ! data size                                       !!! * size of all input files * 
    38    INTEGER       ::   jpidta           !: 1st lateral dimension ( >= jpi ) 
    39    INTEGER       ::   jpjdta           !: 2nd    "         "    ( >= jpj ) 
    40    INTEGER       ::   jpkdta           !: number of levels      ( >= jpk ) 
    41  
    42    ! global or zoom domain size                      !!! * computational domain * 
    43    INTEGER       ::   jpiglo           !: 1st dimension of global domain --> i 
    44    INTEGER       ::   jpjglo           !: 2nd    -                  -    --> j 
    45  
    46    ! zoom starting position  
    47    INTEGER       ::   jpizoom          !: left bottom (i,j) indices of the zoom 
    48    INTEGER       ::   jpjzoom          !: in data domain indices 
    49  
    50    ! Domain characteristics 
    51    INTEGER       ::   jperio           !: lateral cond. type (between 0 and 6) 
    52    !                                       !  = 0 closed                 ;   = 1 cyclic East-West 
    53    !                                       !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
    54    !                                       !  = 4 cyclic East-West AND North fold T-point pivot 
    55    !                                       !  = 5 North fold F-point pivot 
    56    !                                       !  = 6 cyclic East-West AND North fold F-point pivot 
    57  
    58    ! Input file read offset 
    59    LOGICAL       ::   ln_use_jattr     !: Use file global attribute: open_ocean_jstart to determine start j-row  
    60                                            ! when reading input from those netcdf files that have the  
    61                                            ! attribute defined. This is designed to enable input files associated  
    62                                            ! with the extended grids used in the under ice shelf configurations to  
    63                                            ! be used without redundant rows when the ice shelves are not in use. 
    64  
    65    !!  Values set to pp_not_used indicates that this parameter is not used in THIS config. 
    66    !!  Values set to pp_to_be_computed  indicates that variables will be computed in domzgr 
    67    REAL(wp)      ::   pp_not_used       = 999999._wp   !: vertical grid parameter 
    68    REAL(wp)      ::   pp_to_be_computed = 999999._wp   !:    -      -       - 
    69  
    70  
    71  
    72  
    73    !!--------------------------------------------------------------------- 
    74    !! Active tracer parameters 
    75    !!--------------------------------------------------------------------- 
    76    INTEGER, PUBLIC, PARAMETER ::   jpts   = 2    !: Number of active tracers (=2, i.e. T & S ) 
    77    INTEGER, PUBLIC, PARAMETER ::   jp_tem = 1    !: indice for temperature 
    78    INTEGER, PUBLIC, PARAMETER ::   jp_sal = 2    !: indice for salinity 
    79  
    80    !!--------------------------------------------------------------------- 
    81    !! Domain Matrix size  (if AGRIF, they are not all parameters) 
    82    !!--------------------------------------------------------------------- 
    83  
    84  
    85  
    86  
    87  
    88  
    89    INTEGER, PUBLIC  ::   jpi   ! = ( jpiglo-2*jpreci + (jpni-1) ) / jpni + 2*jpreci   !: first  dimension 
    90    INTEGER, PUBLIC  ::   jpj   ! = ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj   !: second dimension 
    91    INTEGER, PUBLIC  ::   jpk   ! = jpkdta 
    92    INTEGER, PUBLIC  ::   jpim1 ! = jpi-1                                            !: inner domain indices 
    93    INTEGER, PUBLIC  ::   jpjm1 ! = jpj-1                                            !:   -     -      - 
    94    INTEGER, PUBLIC  ::   jpkm1 ! = jpk-1                                            !:   -     -      - 
    95    INTEGER, PUBLIC  ::   jpij  ! = jpi*jpj                                          !:  jpi x jpj 
     87   INTEGER, PUBLIC, PARAMETER ::   nn_hls = 1   !: halo width (applies to both rows and columns) 
    9688 
    9789   !!---------------------------------------------------------------------- 
    9890   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    99    !! $Id: par_oce.F90 5836 2015-10-26 14:49:40Z cetlod $  
    100    !! Software governed by the CeCILL licence (./LICENSE) 
     91   !! $Id: par_oce.F90 10068 2018-08-28 14:09:04Z nicolasmartin $  
     92   !! Software governed by the CeCILL license (see ./LICENSE) 
    10193   !!====================================================================== 
    10294END MODULE par_oce 
Note: See TracChangeset for help on using the changeset viewer.