Changeset 907


Ignore:
Timestamp:
2008-04-25T09:35:26+02:00 (13 years ago)
Author:
ctlod
Message:

Re-organisation of ocean namelists to improve readability, see ticket: #125

Location:
trunk/CONFIG
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/CONFIG/GYRE/EXP00/namelist

    r889 r907  
    1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! OPA namelist :  model option and parameter input 
    3 ! ------------- 
    4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    5 ! 
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/OPA  :       1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 
     3!! namelists         2 - Domain           (nam_zgr, nam_zgr_sco, namdom) 
     4!!                   3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
     5!!                                         namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 
     6!!                   4 - lateral boundary (namlbc, namcla, namobc) 
     7!!                   5 - bottom  boundary (nambfr, nambbc, nambbl) 
     8!!                   6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp) 
     9!!                   7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) 
     10!!                   8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm) 
     11!!                   9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr) 
     12!!                   9 - miscellaneous    (namsol) 
     13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     14!  CAUTION: some scripts does not support CAPITALs for logical definition  
     15!  *******       use .true. or .false. and NOT .TRUE. or .FALSE. 
     16 
     17!!====================================================================== 
     18!!                   ***  Run management namelists  *** 
     19!!====================================================================== 
     20!!   namrun            parameters of the run 
     21!!   nam_mpp           Massively Parallel Processing 
     22!!   nam_mpp_dyndist    
     23!!   nam_ctl           Control prints & Benchmark 
     24!!====================================================================== 
     25 
    626!----------------------------------------------------------------------- 
    727!       namrun   parameters of the run 
     
    1030!  cexper     experience name for vairmer format 
    1131!  ln_rstart  boolean term for restart (true or false) 
    12 !  nrstdt     control of the restart timestep: 
    13 !                = 0 restart, do not control nit000 in the restart file. 
    14 !                = 1 restart, control nit000 in the restart file. Do not 
    15 !                    use the date in the restart file (use ndate0 in namelist) 
    16 !                = 2 restart, control nit000 in the restart file, use the date 
    17 !                    in the restart file. ndate0 in the namelist is ignored. 
     32!  nrstdt     restart control = 0 restart, do not control nit000 in the restart file. 
     33!             !               = 1 restart, control nit000 in the restart file. Do not 
     34!             !                   use the date in the restart file (use ndate0 in namelist) 
     35!             !               = 2 restart, control nit000 in the restart file, use the date 
     36!             !                   in the restart file. ndate0 in the namelist is ignored. 
    1837!  nit000     number of the first time step 
    1938!  nitend     number of the last time step 
     
    2443!  nwrite     frequency of OUTPUT file 
    2544!  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 
    26 ! 
    27 !  CAUTION: for usual run scripts, logical value of 
    28 !  *******  ln_rstart must be .true. or .false. 
    29 !                     and NOT .TRUE. or .FALSE. 
    3045&namrun 
    3146   no         =       0 
     
    4055   nstock     =    4320  
    4156   nwrite     =      60  
    42    ln_dimgnnn = .FALSE. 
     57   ln_dimgnnn = .false. 
     58/ 
     59!----------------------------------------------------------------------- 
     60!       nam_mpp      Massively Parallel Processing 
     61!----------------------------------------------------------------------- 
     62!  c_mpi_send         mpi send/recieve type 
     63!                      = 'S'  : standard blocking send 
     64!                      = 'B'  : buffer blocking send 
     65!                      = 'I'  : immediate non-blocking send 
     66&nam_mpp 
     67   c_mpi_send =  'S' 
     68/ 
     69!----------------------------------------------------------------------- 
     70!       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist") 
     71!----------------------------------------------------------------------- 
     72!  jpni   number of processors following i 
     73!  jpnj   number of processors following j 
     74!  jpnij  number of local domains 
     75&nam_mpp_dyndist 
     76   jpni  = 1 
     77   jpnj  = 1 
     78   jpnij = 1 
    4379/ 
    4480!----------------------------------------------------------------------- 
     
    69105   nbit_cmp =     0 
    70106/ 
    71 !----------------------------------------------------------------------- 
    72 !       nam_mpp      Massively Parallel Processing 
    73 !----------------------------------------------------------------------- 
    74 !  c_mpi_send         mpi send/recieve type 
    75 !                      = 'S'  : standard blocking send 
    76 !                      = 'B'  : buffer blocking send 
    77 !                      = 'I'  : immediate non-blocking send 
    78 &nam_mpp 
    79    c_mpi_send =  'S' 
    80 / 
    81 !----------------------------------------------------------------------- 
    82 !       nam_mpp_dyndist      Massively Parallel Distribution 
    83 !              #ifdef ( key_agrif && key_mpp_dyndist ) 
    84 !----------------------------------------------------------------------- 
    85 !  jpni   number of processors following i 
    86 !  jpnj   number of processors following j 
    87 !  jpnij  number of local domains 
    88 &nam_mpp_dyndist 
    89    jpni  = 1 
    90    jpnj  = 1 
    91    jpnij = 1 
    92 / 
     107 
     108!!====================================================================== 
     109!!                      ***  Domain namelists  *** 
     110!!====================================================================== 
     111!!   nam_zgr       vertical coordinate 
     112!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate 
     113!!   namdom        space and time domain (bathymetry, mesh, timestep) 
     114!!====================================================================== 
     115 
    93116!----------------------------------------------------------------------- 
    94117!       nam_zgr       vertical coordinate 
    95118!----------------------------------------------------------------------- 
    96 !  ln_zco     z-coordinate - full steps      (T/F) 
     119!  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined) 
    97120!  ln_zps     z-coordinate - partial steps   (T/F) 
    98121!  ln_sco     s- or hybrid z-s-coordinate    (T/F) 
     
    118141/ 
    119142!----------------------------------------------------------------------- 
    120 !       nam_traadv   advection scheme for tracer (option not control by CPP keys) 
    121 !----------------------------------------------------------------------- 
    122 !  ln_traadv_cen2     2nd order centered scheme    (default T) 
    123 !  ln_traadv_tvd      TVD scheme                   (default F) 
    124 !  ln_traadv_muscl    MUSCL scheme                 (default F) 
    125 !  ln_traadv_muscl2   MUSCL2 scheme                (default F) 
    126 !  ln_traadv_ubs      UBS scheme                   (default F) 
    127 !  ln_traadv_qck      QUICKEST scheme              (default F) 
    128 &nam_traadv 
    129    ln_traadv_cen2   =  .false. 
    130    ln_traadv_tvd    =  .true. 
    131    ln_traadv_muscl  =  .false. 
    132    ln_traadv_muscl2 =  .false. 
    133    ln_traadv_ubs    =  .false. 
    134    ln_traadv_qck    =  .false. 
    135 / 
    136 !----------------------------------------------------------------------- 
    137 !       nam_traldf   lateral diffusion scheme for tracer (option not control by CPP keys) 
    138 !----------------------------------------------------------------------- 
    139 !  Type of the operator :  
    140 !     ln_traldf_lap    laplacian operator          (default T) 
    141 !     ln_traldf_bilap  bilaplacian operator        (default F) 
    142 !  Direction of action  :  
    143 !     ln_traldf_level  iso-level                   (default F) 
    144 !     ln_traldf_hor    horizontal (geopotential)   (default F)^** 
    145 !     ln_traldf_iso    iso-neutral                 (default T)^* 
    146 !  Coefficient 
    147 !     aht0    horizontal eddy diffusivity for tracers (m2/s) 
    148 !     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
    149 !     aeiv0   eddy induced velocity coefficient (m2/s) 
    150 ! ^* require key_ldfslp to compute the direction of the lateral diffusion 
    151 ! ^** require key_ldfslp in s-coordinate 
    152 &nam_traldf 
    153    ln_traldf_lap    =  .true. 
    154    ln_traldf_bilap  =  .false. 
    155    ln_traldf_level  =  .false. 
    156    ln_traldf_hor    =  .false. 
    157    ln_traldf_iso    =  .true. 
    158    aht0    =  1000. 
    159    ahtb0   =     0. 
    160    aeiv0   =  1000. 
    161 / 
    162 !----------------------------------------------------------------------- 
    163 !       nam_dynldf   lateral diffusion on momentum 
    164 !----------------------------------------------------------------------- 
    165 !  Type of the operator :  
    166 !     ln_dynldf_lap    laplacian operator          (default T) 
    167 !     ln_dynldf_bilap  bilaplacian operator        (default F) 
    168 !  Direction of action  :  
    169 !     ln_dynldf_level  iso-level                   (default F) 
    170 !     ln_dynldf_hor    horizontal (geopotential)   (default F)^** 
    171 !     ln_dynldf_iso    iso-neutral                 (default T)^* 
    172 !  Coefficient 
    173 !  ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
    174 !  ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
    175 &nam_dynldf 
    176    ln_dynldf_lap    =  .true. 
    177    ln_dynldf_bilap  =  .false. 
    178    ln_dynldf_level  =  .false. 
    179    ln_dynldf_hor    =  .true. 
    180    ln_dynldf_iso    =  .false. 
    181    ahm0    = 100000. 
    182    ahmb0   =      0. 
    183 / 
    184 !----------------------------------------------------------------------- 
    185 !       namflg   algorithm flags (algorithm not control by CPP keys) 
    186 !----------------------------------------------------------------------- 
    187 !  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
    188 !                                                  centered      time scheme  (F) 
    189 !   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
    190 &namflg 
    191    ln_dynhpg_imp   =  .false. 
    192    nn_dynhpg_rst   =  0 
    193 / 
    194 !----------------------------------------------------------------------- 
    195 !       nam_dynhpg   Hydrostatic pressure gradient option 
    196 !----------------------------------------------------------------------- 
    197 !  type of pressure gradient scheme (choose one only!) 
    198 !     ln_hpg_zco    z-coordinate - full steps                   (default T) 
    199 !     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
    200 !     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
    201 !     ln_hpg_hel    s-coordinate (helsinki modification) 
    202 !     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
    203 !     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
    204 !     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
    205 !  parameters 
    206 !    gamm          weighting coefficient (wdj scheme) 
    207 &nam_dynhpg 
    208    ln_hpg_zco = .true. 
    209    ln_hpg_zps = .false. 
    210    ln_hpg_sco = .false. 
    211    ln_hpg_hel = .false. 
    212    ln_hpg_wdj = .false. 
    213    ln_hpg_djc = .false. 
    214    ln_hpg_rot = .false. 
    215    gamm       = 0.e0 
    216 / 
    217 !----------------------------------------------------------------------- 
    218 !       nam_dynadv   option of physics/algorithm (not control by CPP keys) 
    219 !----------------------------------------------------------------------- 
    220 !  ln_dynadv_vec      vector form flag 
    221 !  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
    222 !  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
    223 &nam_dynadv    
    224    ln_dynadv_vec   = .TRUE. 
    225    ln_dynadv_cen2  = .FALSE. 
    226    ln_dynadv_ubs   = .FALSE. 
    227  
    228 !----------------------------------------------------------------------- 
    229 !       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
    230 !----------------------------------------------------------------------- 
    231 !  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
    232 !  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
    233 !  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
    234 !  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
    235 &nam_dynvor 
    236    ln_dynvor_ene = .TRUE. 
    237    ln_dynvor_ens = .FALSE. 
    238    ln_dynvor_mix = .FALSE. 
    239    ln_dynvor_een = .FALSE. 
    240 / 
    241 !----------------------------------------------------------------------- 
    242 !       namalb   albedo parameters 
    243 !----------------------------------------------------------------------- 
    244 !  cgren    correction of the snow or ice albedo to take into account 
    245 !  albice   albedo of melting ice in the arctic and antarctic 
    246 !  alphd    coefficients for linear interpolation used to compute albedo 
    247 !           between two extremes values (Pyane, 1972) 
    248 !  alphc     "                                         " 
    249 !  alphdi    "                                         " 
    250 &namalb 
    251    cgren    =      0.06 
    252    albice   =      0.5 
    253    alphd    =      0.80 
    254    alphc    =      0.65 
    255    alphdi   =      0.72 
    256 / 
    257 !----------------------------------------------------------------------- 
    258143!       namdom   space and time domain (bathymetry, mesh, timestep) 
    259144!----------------------------------------------------------------------- 
    260 !  ntopo      = 0/1 ,compute/read the bathymetry file 
    261 !               (mbathy, nb of T-ocean levels) 
     145!  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 
    262146!  e3zps_min  the thickness of the partial step is set larger than the 
    263 !  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t 
    264 !                (N.B. 0<e3zps_rat<1) 
     147!  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 
    265148!  nmsh       =1 create a mesh file (coordinates, scale factors, masks) 
    266 !  nacc       the acceleration of convergence method 
    267 !             = 0, no acceleration, rdt = rdttra 
    268 !             = 1, acceleration used, rdt < rdttra(k) 
     149!  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra 
     150!                                                  = 1,     used, rdt < rdttra(k) 
    269151!  atfp       asselin time filter parameter 
    270152!  rdt        time step for the dynamics (and tracer if nacc=0) 
     
    273155!  rdth       depth variation of tracer time step 
    274156!  rdtbt      barotropic time step (for the time splitting algorithm) 
    275 !  nclosea    0 no closed sea  
    276 !             1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
     157!  nclosea    =0 remove the closed sea from the global domain (orca configuration) 
     158!             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
    277159&namdom 
    278160   ntopo     =     0 
     
    289171   nclosea   =     0 
    290172/ 
    291 !----------------------------------------------------------------------- 
    292 !       namfwb   freshwater budget correction 
    293 !----------------------------------------------------------------------- 
    294 !  ln_fwb     logical flag for freshwater budget correction (0 annual mean) 
    295 &namfwb 
    296    ln_fwb    = .false. 
    297 / 
    298 !----------------------------------------------------------------------- 
    299 !       namptr   Poleward Transport Diagnostic 
    300 !----------------------------------------------------------------------- 
    301 !  ln_diaptr  logical flag for Poleward transport computation 
    302 !  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
    303 !             need input basins mask file named "subbasins.nc" 
    304 !  nf_ptr     Frequency of computation 
    305 &namptr 
    306    ln_diaptr = .false. 
    307    ln_subbas = .false. 
    308    nf_ptr    =  15 
    309 / 
    310 !----------------------------------------------------------------------- 
    311 !       namcro   cross land advection 
    312 !----------------------------------------------------------------------- 
    313 !  n_cla   advection between 2 ocean pts separates by land  
    314 &namcla 
    315    n_cla   = 0 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 !       namzdf   vertical physics 
    319 !----------------------------------------------------------------------- 
    320 !  ln_zdfevd  enhanced vertical diffusion         (default T) 
    321 !  ln_zdfnpc  Non-Penetrative Convection          (default T) 
    322 !  avm0       vertical eddy viscosity for the dynamic (m2/s) 
    323 !  avt0       vertical eddy diffusivity for tracers (m2/s) 
    324 !  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
    325 !  nevdm      = 0  apply enhanced mixing on tracer only 
    326 !             = 1  apply enhanced mixing on both tracer and momentum 
    327 !  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
    328 !                               (=F)  euler backward (F) 
    329 !  n_zdfexp   number of sub-timestep for time splitting scheme 
    330 &namzdf 
    331    ln_zdfevd = .true. 
    332    ln_zdfnpc = .false. 
    333    avm0      = 1.2e-4 
    334    avt0      = 1.2e-5 
    335    avevd     =   100. 
    336    n_evdm    =     1 
    337    ln_zdfexp =  .false. 
    338    n_zdfexp  =     3 
    339 / 
    340 !----------------------------------------------------------------------- 
    341 !       namnpc   vnon penetrative convection 
    342 !----------------------------------------------------------------------- 
    343 !  nnpc1   non penetrative convective scheme frequency 
    344 !  nnpc2   non penetrative convective scheme print frequency 
    345 &namnpc 
    346    nnpc1  =      1 
    347    nnpc2  =    365 
    348 / 
    349 !----------------------------------------------------------------------- 
    350 !       nambbl   bottom boundary layer scheme 
    351 !----------------------------------------------------------------------- 
    352 !  atrbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s) 
    353 &nambbl 
    354    atrbbl = 10000. 
    355 / 
    356 !----------------------------------------------------------------------- 
    357 !       namric   richardson number dependent vertical diffusion 
    358 !                ( #ifdef "key_zdfrichardson" ) 
    359 !----------------------------------------------------------------------- 
    360 !  avmri   maximum value of the vertical viscosity 
    361 !  alp     coefficient of the parameterization 
    362 !  nric    coefficient of the parameterization 
    363 &namric 
    364    avmri = 100.e-4 
    365    alp   =      5. 
    366    nric  =       2 
    367 / 
    368 !----------------------------------------------------------------------- 
    369 !       namtke   turbulent eddy kinetic dependent vertical diffusion 
    370 !                ( #ifdef "key_zdftke" ) 
    371 !----------------------------------------------------------------------- 
    372 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    373 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    374 !  ediss    coef. of the Kolmogoroff dissipation   
    375 !  ebb      coef. of the surface input of tke 
    376 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    377 !  emin     minimum value of tke (m^2/s^2) 
    378 !  emin0    surface minimum value of tke (m^2/s^2) 
    379 !  nitke    number of restart iterative loops 
    380 !  ri_c     critic richardson number 
    381 !  nmxl     flag on mixing length used 
    382 !           = 0 bounded by the distance to surface and bottom 
    383 !           = 1 bounded by the local vertical scale factor 
    384 !           = 2 first vertical derivative of mixing length bounded by 1 
    385 !  npdl     flag on prandtl number 
    386 !           = 0 no vertical prandtl number (avt=avm) 
    387 !           = 1 prandtl number function of richarson number (avt=pdl*avm) 
    388 !           = 2 same as = 1 but a shapiro filter is applied on pdl 
    389 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    390 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    391 &namtke 
    392    ln_rstke = .false. 
    393    ediff =       0.1 
    394    ediss =       0.7 
    395    ebb   =      3.75 
    396    efave =        1. 
    397    emin  =     1.e-5  
    398    emin0 =     1.e-4 
    399    nitke =        50 
    400    nmxl  =         2 
    401    npdl  =         1 
    402    navb  =         0 
    403 / 
    404 !----------------------------------------------------------------------- 
    405 !       namkpp   K-Profile Parameterization dependent vertical diffusion 
    406 !                ( #ifdef "key_zdfkpp" ) 
    407 !----------------------------------------------------------------------- 
    408 !  ln_kpprimix   shear instability mixing  (default T) 
    409 !  difmiw        constant internal wave viscosity (m2/s) 
    410 !  difsiw        constant internal wave diffusivity (m2/s) 
    411 !  Riinfty       local Richardson Number limit for shear instability 
    412 !  difri         maximum shear mixing at Rig = 0    (m2/s) 
    413 !  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
    414 !  difcon        maximum mixing in interior convection (m2/s) 
    415 !  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
    416 !  navb          = 0/1 flag for constant or profile background avt 
    417 &namkpp 
    418    ln_kpprimix  = .true. 
    419    difmiw       =  1.e-04 
    420    difsiw       =  0.1e-04 
    421    Riinfty      =  0.8 
    422    difri        =  0.0050 
    423    bvsqcon      = -0.01e-07 
    424    difcon       =  1. 
    425    navb         =  0 
    426    nave         =  1 
    427 / 
    428 !----------------------------------------------------------------------- 
    429 !       namddm   double diffusive mixing parameterization 
    430 !----------------------------------------------------------------------- 
    431 !   avts    maximum avs for dd mixing  
    432 !   hsbfr   heat/salt buoyancy flux ratio 
    433 &namddm 
    434       avts  = 1.e-4 
    435       hsbfr = 1.6 
    436 / 
    437 !----------------------------------------------------------------------- 
    438 !       namlbc   lateral momentum boundary condition 
    439 !----------------------------------------------------------------------- 
    440 !  shlat   lateral boundary condition on velocity 
    441 !                   shlat = 0 , free slip 
    442 !               0 < shlat < 2 , partial slip 
    443 !                   shlat = 2 , no slip 
    444 !               2 < shlat     , strong slip 
    445 &namlbc 
    446    shlat  =      0. 
    447 / 
    448 !----------------------------------------------------------------------- 
    449 !       nambfr   bottom friction 
    450 !----------------------------------------------------------------------- 
    451 !  nbotfr  type of bottom friction  
    452 !                  nbotfr = 0 , no slip  
    453 !                  nbotfr = 1 , linear friction 
    454 !                  nbotfr = 2 , nonlinear friction 
    455 !                  nbotfr = 3 , free slip 
    456 !  bfri1   bottom drag coefficient (linear case) 
    457 !  bfri2   bottom drag coefficient (non linear case) 
    458 !  bfeb2   bottom turbulent kinetic energy  (m^2/s^2) 
    459 &nambfr 
    460    nbotfr =       2 
    461    bfri1  =   4.e-4 
    462    bfri2  =   1.e-3 
    463    bfeb2  =  2.5e-3 
    464 / 
    465 !----------------------------------------------------------------------- 
    466 !       nambbc   bottom temperature boundary condition 
    467 !----------------------------------------------------------------------- 
    468 !  ngeo_flux  = 0 no geothermal heat flux 
    469 !             = 1 constant geothermal heat flux 
    470 !             = 2 variable geothermal heat flux (read in geothermal_heating.nc)  
    471 !                 ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 
    472 !  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
    473 &nambbc 
    474    ngeo_flux =  0 
    475    ngeo_flux_const = 86.4e-3 
     173 
     174!!====================================================================== 
     175!!            ***  Surface Boundary Condition namelists  *** 
     176!!====================================================================== 
     177!!   namsbc        surface boundary condition 
     178!!   namsbc_ana    analytical         formulation 
     179!!   namsbc_flx    flux               formulation 
     180!!   namsbc_clio   CLIO bulk formulea formulation 
     181!!   namsbc_core   CORE bulk formulea formulation 
     182!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled") 
     183!!   namqsr        penetrative solar radiation 
     184!!   namsbc_rnf    river runoffs 
     185!!   namsbc_ssr    sea surface restoring term (for T and/or S) 
     186!!   namalb        albedo parameters 
     187!!====================================================================== 
     188 
     189!----------------------------------------------------------------------- 
     190!       namsbc       Surface Boundary Condition (surface module) 
     191!----------------------------------------------------------------------- 
     192!  nn_fsbc       frequency of surface boundary condition computation  
     193!                (= the frequency of sea-ice model call) 
     194!  ln_ana        analytical formulation (fill namsbc_ana )  
     195!  ln_flx        flux formulation       (fill namsbc_flx )  
     196!  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core)  
     197!  ln_blk_core   CORE bulk formulation  (fill namsbc_clio)  
     198!  ln_cpl        Coupled formulation    (fill namsbc_cpl )  
     199!  nn_ice        =0 no ice boundary condition   , 
     200!                =1 observed ice-cover          , 
     201!                =2 ice-model used                                      ("key_lim3" or "key_lim2) 
     202!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
     203!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr) 
     204!  nn_fwb       FreshWater Budget: =0 no control                              ,  
     205!                                  =1 annual global mean of e-p-r set to zero   , 
     206!                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step 
     207&namsbc 
     208   nn_fsbc     = 1 
     209   ln_ana      = .true. 
     210   ln_flx      = .false. 
     211   ln_blk_clio = .false. 
     212   ln_blk_core = .false. 
     213   ln_cpl      = .false. 
     214   nn_ice      = 0 
     215   ln_dm2dc    = .false. 
     216   ln_rnf      = .false. 
     217   ln_ssr      = .false. 
     218   nn_fwb      = 0 
     219/ 
     220!----------------------------------------------------------------------- 
     221!       namsbc_ana   analytical surface boundary condition 
     222!----------------------------------------------------------------------- 
     223!  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
     224!  rn_utau0    uniform value for the i-stress 
     225!  rn_vtau0    uniform value for the j-stress 
     226!  rn_q0       uniform value for the total heat flux 
     227!  rn_qsr0     uniform value for the solar radiation 
     228!  rn_emp0     uniform value for the freswater budget (E-P) 
     229&namsbc_ana 
     230   nn_tau000  =   100 
     231   rn_utau0   =   0.1e0 
     232   rn_vtau0   =   0.e0 
     233   rn_q0      =   0.e0 
     234   rn_qsr0    =   0.e0 
     235   rn_emp0    =   0.e0 
     236/ 
     237!----------------------------------------------------------------------- 
     238!       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
     239!----------------------------------------------------------------------- 
     240!  cn_dir   directory in which the model is executed 
     241! 
     242!  THE ORDER OF THE FILES MATTER: 
     243!               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     244!               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     245&namsbc_flx 
     246   cn_dir       = './' 
     247   sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0 
     248   sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0 
     249   sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0 
     250   sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0 
     251   sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0 
     252/       
     253!----------------------------------------------------------------------- 
     254!       namsbc_clio  CLIO bulk formulea 
     255!----------------------------------------------------------------------- 
     256!  cn_dir   directory in which the model is executed 
     257! 
     258!  THE ORDER OF THE FILES MATTER: 
     259!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     260!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     261&namsbc_clio 
     262   cn_dir       = './' 
     263   sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0     
     264   sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0     
     265   sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0       
     266   sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0      
     267   sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0     
     268   sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0       
     269   sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0        
     270/ 
     271!----------------------------------------------------------------------- 
     272!       namsbc_core  CORE bulk formulea 
     273!----------------------------------------------------------------------- 
     274!  cn_dir         directory in which the model is executed 
     275!  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
     276!  alpha_precip   multiplication factor for precipitation (total & snow) 
     277! 
     278!  THE ORDER OF THE FILES MATTER: 
     279!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     280!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     281&namsbc_core 
     282   cn_dir       = './' 
     283   ln_2m        = .false. 
     284   alpha_precip = 1. 
     285   sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0     
     286   sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0     
     287   sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0       
     288   sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0       
     289   sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0      
     290   sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0     
     291   sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0       
     292   sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0        
     293/ 
     294!----------------------------------------------------------------------- 
     295!       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled") 
     296!----------------------------------------------------------------------- 
     297&namsbc_cpl 
    476298/ 
    477299!----------------------------------------------------------------------- 
     
    489311/ 
    490312!----------------------------------------------------------------------- 
    491 !       namtdp   tracer newtonian damping ('key_tradmp') 
    492 !----------------------------------------------------------------------- 
    493 !  ndmp    type of damping in temperature and salinity  
    494 !          (='latitude', damping poleward of 'ndmp' degrees and function  
    495 !             of the distance-to-coast. Red and Med Seas as ndmp=-1) 
    496 !          (=-1 damping only in Med and Red Seas) 
    497 !  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
    498 !  nmldmp  type of damping in the mixed layer 
    499 !          (=0 damping throughout the water column) 
    500 !     (=1 no damping in the mixed layer defined by avt >5cm2/s ) 
    501 !     (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
    502 !  sdmp    surface time scale for internal damping (days) 
    503 !  bdmp    bottom time scale for internal damping (days) 
    504 !  hdmp    depth of transition between sdmp and bdmp (meters) 
    505 &namtdp 
    506    ndmp   =   -1 
    507    ndmpf  =    1 
    508    nmldmp =    1 
    509    sdmp   =  50. 
    510    bdmp   = 360. 
    511    hdmp   = 800. 
    512 / 
    513 !----------------------------------------------------------------------- 
    514 !       nameos   ocean physical parameters 
    515 !----------------------------------------------------------------------- 
    516 !  neos    type of equation of state and Brunt-Vaisala frequency 
    517 !          = 0, UNESCO (formulation of Jackett and McDougall (1994) 
    518 !                                         and of McDougall (1987) ) 
    519 !          = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    520 !          = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    521 !                               with rau0=1020 set in parcst routine 
    522 !  ralpha  thermal expension coefficient (linear equation of state) 
    523 !  rbeta   saline  expension coefficient (linear equation of state) 
    524 &nameos 
    525    neos   =      2 
    526    ralpha =  2.e-4 
    527    rbeta  =  7.7e-4 
    528 / 
    529 !----------------------------------------------------------------------- 
    530 !       namsol   elliptic solver / island / free surface  
    531 !----------------------------------------------------------------------- 
    532 !  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
    533 !                            (=2 successive-over-relaxation: sor) 
    534 !                            (=3 FETI: fet, all require "key_feti" defined) 
    535 !  nsol_arp  absolute/relative (0/1) precision convergence test 
    536 !  nmin      minimum of iterations for the SOR solver 
    537 !  nmax      maximum of iterations for the SOR solver 
    538 !  nmod      frequency of test for the SOR solver 
    539 !  eps       absolute precision of the solver 
    540 !  resmax    absolute precision for the SOR solver 
    541 !  sor       optimal coefficient for SOR solver 
    542 !  epsisl    absolute precision on stream function solver 
    543 !  nmisl     maximum pcg iterations for island 
    544 !  rnu       strength of the additional force used in free surface b.c. 
    545 &namsol 
    546    nsolv     =      2 
    547    nsol_arp  =      0 
    548    nmin      =    210 
    549    nmax      =    800 
    550    nmod      =     10 
    551    eps       =  1.E-6 
    552    resmax    = 1.E-10 
    553    sor       =   1.96 
    554    epsisl    = 1.e-10 
    555    nmisl     =   4000 
    556    rnu       =     1. 
    557 / 
    558 !======================================================================= 
    559 !   Diagnostics namelists 
    560 !       namtrd    dynamics and/or tracer trends 
    561 !       namgap    level mean model-data gap 
    562 !       namznl    zonal mean heat & freshwater fluxes computation 
    563 !       namspr    surface pressure in rigid-lid 
    564 !======================================================================= 
    565 !----------------------------------------------------------------------- 
    566 !       namtrd    diagnostics on dynamics and/or tracer trends 
    567 !                         ('key_trdyn' and/or 'key_trdtra') 
    568 !                 or mixed-layer trends ('key_trdmld') 
    569 !----------------------------------------------------------------------- 
    570 !  ntrd              time step frequency dynamics and tracers trends 
    571 !  nctls             control surface type in mixed-layer trends (0,1 or n<jpk) 
    572 !  ln_trdmld_restart restart for ML diagnostics 
    573 !  ucf               unit conversion factor (=1 -> /seconds | =86400. -> /day) 
    574 !  ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 
    575 &namtrd 
    576    ntrd  = 365 
    577    nctls =   0 
    578    ln_trdmld_restart = .false. 
    579    ucf   =  1. 
    580    ln_trdmld_instant = .false. 
    581 / 
    582 !----------------------------------------------------------------------- 
    583 !       namgap    level mean model-data gap ('key_diagap') 
    584 !----------------------------------------------------------------------- 
    585 !  ngap    time-step frequency of model-data gap computation 
    586 !  nprg    time-step frequency of gap print in model output 
    587 &namgap 
    588    ngap =  15 
    589    nprg =  10 
    590 / 
    591 !----------------------------------------------------------------------- 
    592 !       namznl    zonal mean heat & freshwater fluxes computation 
    593 !                 (#ifdef "key_diaznl") 
    594 !----------------------------------------------------------------------- 
    595 !  nfznl   time-step frequency of zonal mean fluxes computation 
    596 &namznl 
    597    nfznl =  15 
    598 / 
    599 !----------------------------------------------------------------------- 
    600 !       namspr  surface pressure diagnostic 
    601 !----------------------------------------------------------------------- 
    602 !  nmaxp   maximum of iterations for the solver 
    603 !  epsp    absolute precision of the solver 
    604 !  niterp  number of iteration done by the solver 
    605 &namspr 
    606    nmaxp   =   1000 
    607    epsp    =  1.e-3 
    608    niterp  =    400 
     313!       namsbc_rnf   runoffs namelist surface boundary condition 
     314!----------------------------------------------------------------------- 
     315!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
     316!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
     317!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
     318!cn_dir      Root directory for location of ssr files 
     319!sn_rnf      informations about the runoff file to be read 
     320!sn_cnf      informations about the runoff mouth file to be read 
     321!  THE ORDER OF THE FILES MATTER: 
     322!               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
     323!               !                    !   (if <0  months) !   name     !     (logical)      ! 
     324&namsbc_rnf 
     325   nn_runoff  = 0 
     326   rn_hrnf    = 0.e0 
     327   rn_avt_rnf = 1.e-3 
     328   cn_dir     = './' 
     329   sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0      
     330   sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0               
     331/ 
     332!----------------------------------------------------------------------- 
     333!       namsbc_ssr   surface boundary condition : sea surface restoring 
     334!----------------------------------------------------------------------- 
     335!  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
     336!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
     337!             (=1, salt flux, concentration/dillution emps only) 
     338!             (=2, volume flux, both emp and emps are updated)  
     339!  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
     340!  deds       magnitude of the damping on salinity [mm/day] 
     341!  THE ORDER OF THE FILES MATTER: 
     342!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     343!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     344&namsbc_ssr 
     345   nn_sstr      = 0 
     346   nn_sssr      = 0 
     347   dqdt         = -40. 
     348   deds         =  27.7 
     349   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0 
     350   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0 
     351/       
     352!----------------------------------------------------------------------- 
     353!       namalb   albedo parameters 
     354!----------------------------------------------------------------------- 
     355!  cgren    correction of the snow or ice albedo to take into account 
     356!  albice   albedo of melting ice in the arctic and antarctic 
     357!  alphd    coefficients for linear interpolation used to compute albedo 
     358!           between two extremes values (Pyane, 1972) 
     359!  alphc     "                                         " 
     360!  alphdi    "                                         " 
     361&namalb 
     362   cgren    =      0.06 
     363   albice   =      0.5 
     364   alphd    =      0.80 
     365   alphc    =      0.65 
     366   alphdi   =      0.72 
     367/ 
     368 
     369!!====================================================================== 
     370!!               ***  Lateral boundary condition  *** 
     371!!====================================================================== 
     372!!   namlbc        lateral momentum boundary condition 
     373!!   namcla        cross land advection 
     374!!   namobc        open boundaries parameters                           ("key_obc") 
     375!!====================================================================== 
     376 
     377!----------------------------------------------------------------------- 
     378!       namlbc   lateral momentum boundary condition 
     379!----------------------------------------------------------------------- 
     380!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip 
     381!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip 
     382&namlbc 
     383   shlat  =      0. 
     384/ 
     385!----------------------------------------------------------------------- 
     386!       namcla   cross land advection 
     387!----------------------------------------------------------------------- 
     388!  n_cla   advection between 2 ocean pts separates by land  
     389&namcla 
     390   n_cla   = 0 
    609391/ 
    610392!----------------------------------------------------------------------- 
     
    638420    ln_vol_cst  = .false. 
    639421/ 
    640 !----------------------------------------------------------------------- 
    641 !       namflo    float parameters (#ifdef key_float) 
     422 
     423!!====================================================================== 
     424!!                 ***  Bottom boundary condition  *** 
     425!!====================================================================== 
     426!!   nambfr        bottom friction 
     427!!   nambbc        bottom temperature boundary condition                ("key_trabbc") 
     428!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv") 
     429!!====================================================================== 
     430!----------------------------------------------------------------------- 
     431!       nambfr   bottom friction 
     432!----------------------------------------------------------------------- 
     433!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction 
     434!                                     = 3 : free slip  ;  = 2 : non linear friction 
     435!  bfri1    bottom drag coefficient (linear case) 
     436!  bfri2    bottom drag coefficient (non linear case) 
     437!  bfeb2    bottom turbulent kinetic energy  (m2/s2) 
     438&nambfr 
     439   nbotfr =       2 
     440   bfri1  =   4.e-4 
     441   bfri2  =   1.e-3 
     442   bfeb2  =  2.5e-3 
     443/ 
     444!----------------------------------------------------------------------- 
     445!       nambbc   bottom temperature boundary condition 
     446!----------------------------------------------------------------------- 
     447!  ngeo_flux  = 0 no geothermal heat flux 
     448!             = 1 constant geothermal heat flux 
     449!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 
     450!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
     451&nambbc 
     452   ngeo_flux =  0 
     453   ngeo_flux_const = 86.4e-3 
     454/ 
     455!----------------------------------------------------------------------- 
     456!       nambbl   bottom boundary layer scheme 
     457!----------------------------------------------------------------------- 
     458!                                !  diffusive bbl                       ("key_trabbl") 
     459!                                !  advective bbl                       ("key_trabbl_adv") 
     460!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s) 
     461&nambbl 
     462   atrbbl = 10000. 
     463/ 
     464 
     465!!====================================================================== 
     466!!                        Tracer (T & S ) namelists 
     467!!====================================================================== 
     468!!   nameos        equation of state 
     469!!   nam_traadv    advection scheme 
     470!!   nam_traldf    lateral diffusion scheme 
     471!!   namtdp        tracer newtonian damping                             ("key_tradmp") 
     472!!====================================================================== 
     473 
     474!----------------------------------------------------------------------- 
     475!       nameos   ocean physical parameters 
     476!----------------------------------------------------------------------- 
     477!  neos    type of equation of state and Brunt-Vaisala frequency 
     478!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
     479!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
     480!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
     481!  ralpha  thermal expension coefficient (neos= 1 or 2) 
     482!  rbeta   saline  expension coefficient (neos= 2) 
     483&nameos 
     484   neos   =      2 
     485   ralpha =  2.e-4 
     486   rbeta  =  7.7e-4 
     487/ 
     488!----------------------------------------------------------------------- 
     489!       nam_traadv   advection scheme for tracer 
     490!----------------------------------------------------------------------- 
     491!  ln_traadv_cen2     2nd order centered scheme          (default T) 
     492!  ln_traadv_tvd      TVD scheme                         (default F) 
     493!  ln_traadv_muscl    MUSCL scheme                       (default F) 
     494!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F) 
     495!  ln_traadv_ubs      UBS scheme                         (default F) 
     496!  ln_traadv_qck      QUICKEST scheme                    (default F) 
     497&nam_traadv 
     498   ln_traadv_cen2   =  .false. 
     499   ln_traadv_tvd    =  .true. 
     500   ln_traadv_muscl  =  .false. 
     501   ln_traadv_muscl2 =  .false. 
     502   ln_traadv_ubs    =  .false. 
     503   ln_traadv_qck    =  .false. 
     504/ 
     505!----------------------------------------------------------------------- 
     506!       nam_traldf   lateral diffusion scheme for tracer 
     507!----------------------------------------------------------------------- 
     508!  Type of the operator :  
     509!     ln_traldf_lap    laplacian operator          (default T) 
     510!     ln_traldf_bilap  bilaplacian operator        (default F) 
     511!  Direction of action  :  
     512!     ln_traldf_level  iso-level                   (default F) 
     513!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T) 
     514!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
     515!  Coefficient 
     516!     aht0    horizontal eddy diffusivity for tracers (m2/s) 
     517!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
     518!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv") 
     519&nam_traldf 
     520   ln_traldf_lap    =  .true. 
     521   ln_traldf_bilap  =  .false. 
     522   ln_traldf_level  =  .false. 
     523   ln_traldf_hor    =  .false. 
     524   ln_traldf_iso    =  .true. 
     525   aht0    =  1000. 
     526   ahtb0   =     0. 
     527   aeiv0   =  1000. 
     528/ 
     529!----------------------------------------------------------------------- 
     530!       namtdp   tracer newtonian damping ('key_tradmp') 
     531!----------------------------------------------------------------------- 
     532!  ndmp    type of damping in temperature and salinity  
     533!          !   ='latitude', damping poleward of 'ndmp' degrees and function  
     534!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1 
     535!          !   =-1          damping only in Med and Red Seas 
     536!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
     537!  nmldmp  type of damping in the mixed layer 
     538!          !   =0 damping throughout the water column 
     539!          !   =1 no damping in the mixed layer defined by avt >5cm2/s  
     540!          !   =2 no damping in the mixed layer defined rho<rho(surf)+.01  
     541!  sdmp    surface time scale for internal damping (days) 
     542!  bdmp    bottom time scale for internal damping (days) 
     543!  hdmp    depth of transition between sdmp and bdmp (meters) 
     544&namtdp 
     545   ndmp   =   -1 
     546   ndmpf  =    1 
     547   nmldmp =    1 
     548   sdmp   =  50. 
     549   bdmp   = 360. 
     550   hdmp   = 800. 
     551/ 
     552 
     553!!====================================================================== 
     554!!                      ***  Dynamics namelists  *** 
     555!!====================================================================== 
     556!!   nam_dynadv    formulation of the momentum advection 
     557!!   nam_dynvor    advection scheme 
     558!!   nam_dynhpg    hydrostatic pressure gradient 
     559!!   namflg        hydrostatic pressure gradient time stepping 
     560!!   nam_dynspg    surface pressure gradient                            (CPP key only) 
     561!!   nam_dynldf    lateral diffusion scheme 
     562!!====================================================================== 
     563 
     564!----------------------------------------------------------------------- 
     565!       nam_dynadv   formulation of the momentum advection 
     566!----------------------------------------------------------------------- 
     567!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T) 
     568!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
     569!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
     570&nam_dynadv    
     571   ln_dynadv_vec   = .true. 
     572   ln_dynadv_cen2  = .false. 
     573   ln_dynadv_ubs   = .false. 
     574 
     575!----------------------------------------------------------------------- 
     576!       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
     577!----------------------------------------------------------------------- 
     578!  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
     579!  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
     580!  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
     581!  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
     582&nam_dynvor 
     583   ln_dynvor_ene = .true. 
     584   ln_dynvor_ens = .false. 
     585   ln_dynvor_mix = .false. 
     586   ln_dynvor_een = .false. 
     587/ 
     588!----------------------------------------------------------------------- 
     589!       nam_dynhpg   Hydrostatic pressure gradient option 
     590!----------------------------------------------------------------------- 
     591!  type of pressure gradient scheme (choose one only!) 
     592!     ln_hpg_zco    z-coordinate - full steps                   (default T) 
     593!     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
     594!     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
     595!     ln_hpg_hel    s-coordinate (helsinki modification) 
     596!     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
     597!     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
     598!     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
     599!  parameters 
     600!     gamm          weighting coefficient (wdj scheme) 
     601&nam_dynhpg 
     602   ln_hpg_zco = .true. 
     603   ln_hpg_zps = .false. 
     604   ln_hpg_sco = .false. 
     605   ln_hpg_hel = .false. 
     606   ln_hpg_wdj = .false. 
     607   ln_hpg_djc = .false. 
     608   ln_hpg_rot = .false. 
     609   gamm       = 0.e0 
     610/ 
     611!----------------------------------------------------------------------- 
     612!       namflg   algorithm flags (algorithm not control by CPP keys) 
     613!----------------------------------------------------------------------- 
     614!  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
     615!                                                  centered     time scheme  (F) 
     616!   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
     617&namflg 
     618   ln_dynhpg_imp   =  .false. 
     619   nn_dynhpg_rst   =  0 
     620/ 
     621!----------------------------------------------------------------------- 
     622!       nam_dynspg   surface pressure gradient   (CPP key only) 
     623!----------------------------------------------------------------------- 
     624!                                !  explicit free surface               ("key_dynspg_exp") 
     625!                                !  filtered free surface               ("key_dynspg_flt") 
     626!                                !  split-explicit free surface         ("key_dynspg_ts") 
     627!                                !  rigid-lid                           ("key_dynspg_rl") 
     628 
     629!----------------------------------------------------------------------- 
     630!       nam_dynldf   lateral diffusion on momentum 
     631!----------------------------------------------------------------------- 
     632!  Type of the operator :  
     633!     ln_dynldf_lap    laplacian operator          (default T) 
     634!     ln_dynldf_bilap  bilaplacian operator        (default F) 
     635!  Direction of action  :  
     636!     ln_dynldf_level  iso-level                   (default F) 
     637!     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T) 
     638!     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
     639!  Coefficient 
     640!     ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
     641!     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
     642&nam_dynldf 
     643   ln_dynldf_lap    =  .true. 
     644   ln_dynldf_bilap  =  .false. 
     645   ln_dynldf_level  =  .false. 
     646   ln_dynldf_hor    =  .true. 
     647   ln_dynldf_iso    =  .false. 
     648   ahm0    = 100000. 
     649   ahmb0   =      0. 
     650/ 
     651 
     652!!====================================================================== 
     653!!             Tracers & Dynamics vertical physics namelists 
     654!!====================================================================== 
     655!!       namzdf        vertical physics 
     656!!       namnpc        non penetrative convection                        
     657!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      ) 
     658!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      ) 
     659!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      ) 
     660!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      ) 
     661!!====================================================================== 
     662 
     663!----------------------------------------------------------------------- 
     664!       namzdf   vertical physics 
     665!----------------------------------------------------------------------- 
     666!  ln_zdfevd  enhanced vertical diffusion         (default T) 
     667!  ln_zdfnpc  Non-Penetrative Convection          (default T) 
     668!  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst") 
     669!  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst") 
     670!  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
     671!  n_evdm     = 0  apply enhanced mixing on tracer only 
     672!             = 1  apply enhanced mixing on both tracer and momentum 
     673!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
     674!                               (=F)  euler backward (F) 
     675!  n_zdfexp   number of sub-timestep for time splitting scheme 
     676&namzdf 
     677   ln_zdfevd = .true. 
     678   ln_zdfnpc = .false. 
     679   avm0      = 1.2e-4 
     680   avt0      = 1.2e-5 
     681   avevd     =   100. 
     682   n_evdm    =     1 
     683   ln_zdfexp =  .false. 
     684   n_zdfexp  =     3 
     685/ 
     686!----------------------------------------------------------------------- 
     687!       namnpc   non penetrative convection 
     688!----------------------------------------------------------------------- 
     689!  nnpc1   non penetrative convective scheme frequency    
     690!  nnpc2   non penetrative convective scheme print frequency 
     691&namnpc 
     692   nnpc1  =      1 
     693   nnpc2  =    365 
     694/ 
     695!----------------------------------------------------------------------- 
     696!       namric   richardson number dependent vertical diffusion         ("key_zdfric" ) 
     697!----------------------------------------------------------------------- 
     698!  avmri   maximum value of the vertical viscosity 
     699!  alp     coefficient of the parameterization 
     700!  nric    coefficient of the parameterization 
     701&namric 
     702   avmri = 100.e-4 
     703   alp   =      5. 
     704   nric  =       2 
     705/ 
     706!----------------------------------------------------------------------- 
     707!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
     708!----------------------------------------------------------------------- 
     709!  ln_rstke flag to restart with tke from a run without tke (default F) 
     710!  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
     711!  ediss    coef. of the Kolmogoroff dissipation   
     712!  ebb      coef. of the surface input of tke 
     713!  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
     714!  emin     minimum value of tke (m^2/s^2) 
     715!  emin0    surface minimum value of tke (m^2/s^2) 
     716!  nitke    number of restart iterative loops 
     717!  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
     718!  nmxl     length used   = 0 bounded by the distance to surface and bottom 
     719!           !             = 1 bounded by the local vertical scale factor 
     720!           !             = 2 first vertical derivative of mixing length bounded by 1 
     721!  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
     722!           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
     723!           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
     724!  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
     725!  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
     726&namtke 
     727   ln_rstke = .false. 
     728   ediff =       0.1 
     729   ediss =       0.7 
     730   ebb   =      3.75 
     731   efave =        1. 
     732   emin  =     1.e-5 
     733   emin0 =     1.e-4 
     734   nitke =        50 
     735   nmxl  =         2 
     736   npdl  =         1 
     737   navb  =         0 
     738/ 
     739!----------------------------------------------------------------------- 
     740!       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" ) 
     741!----------------------------------------------------------------------- 
     742!  ln_kpprimix   shear instability mixing  (default T) 
     743!  difmiw        constant internal wave viscosity (m2/s) 
     744!  difsiw        constant internal wave diffusivity (m2/s) 
     745!  Riinfty       local Richardson Number limit for shear instability 
     746!  difri         maximum shear mixing at Rig = 0    (m2/s) 
     747!  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
     748!  difcon        maximum mixing in interior convection (m2/s) 
     749!  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
     750!  navb          = 0/1 flag for constant or profile background avt 
     751&namkpp 
     752   ln_kpprimix  = .true. 
     753   difmiw       =  1.e-04 
     754   difsiw       =  0.1e-04 
     755   Riinfty      =  0.8 
     756   difri        =  0.0050 
     757   bvsqcon      = -0.01e-07 
     758   difcon       =  1. 
     759   navb         =  0 
     760   nave         =  1 
     761/ 
     762!----------------------------------------------------------------------- 
     763!       namddm   double diffusive mixing parameterization               ("key_zdfddm") 
     764!----------------------------------------------------------------------- 
     765!   avts    maximum avs for dd mixing  
     766!   hsbfr   heat/salt buoyancy flux ratio 
     767&namddm 
     768      avts  = 1.e-4 
     769      hsbfr = 1.6 
     770/ 
     771 
     772!!====================================================================== 
     773!!                  ***  Miscelaneous namelists  *** 
     774!!====================================================================== 
     775!!   namsol      elliptic solver / island / free surface  
     776!!====================================================================== 
     777 
     778!----------------------------------------------------------------------- 
     779!       namsol   elliptic solver / island / free surface  
     780!----------------------------------------------------------------------- 
     781!  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
     782!                            (=2 successive-over-relaxation: sor) 
     783!                            (=3 FETI currently it does not work!       ("key_feti") 
     784!  nsol_arp  absolute/relative (0/1) precision convergence test 
     785!  nmin      minimum of iterations for the SOR solver 
     786!  nmax      maximum of iterations for the SOR solver 
     787!  nmod      frequency of test for the SOR solver 
     788!  eps       absolute precision of the solver 
     789!  resmax    absolute precision for the SOR solver 
     790!  sor       optimal coefficient for SOR solver 
     791!  epsisl    absolute precision on stream function solver 
     792!  nmisl     maximum pcg iterations for island 
     793!  rnu       strength of the additional force used in free surface b.c. 
     794&namsol 
     795   nsolv     =      2 
     796   nsol_arp  =      0 
     797   nmin      =    210 
     798   nmax      =    800 
     799   nmod      =     10 
     800   eps       =  1.E-6 
     801   resmax    = 1.E-10 
     802   sor       =   1.96 
     803   epsisl    = 1.e-10 
     804   nmisl     =   4000 
     805   rnu       =     1. 
     806/ 
     807!!====================================================================== 
     808!!                       ***  Diagnostics namelists  *** 
     809!!====================================================================== 
     810!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld") 
     811!!   namgap       level mean model-data gap                             ("key_diagap") 
     812!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr") 
     813!!   namflo       float parameters                                      ("key_float") 
     814!!   namptr       Poleward Transport Diagnostics 
     815!!====================================================================== 
     816 
     817!----------------------------------------------------------------------- 
     818!       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra") 
     819!                          or mixed-layer trends                        ('key_trdmld') 
     820!                          or barotropic vorticity                      ("key_trdvor") 
     821!----------------------------------------------------------------------- 
     822!  ntrd                time step frequency dynamics and tracers trends 
     823!  nctls               control surface type in mixed-layer trends (0,1 or n<jpk) 
     824!  ln_trdmld_restart   restart for ML diagnostics 
     825!  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day) 
     826!  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S 
     827&namtrd 
     828   ntrd              = 365 
     829   nctls             =   0 
     830   ln_trdmld_restart = .false. 
     831   ucf               =  1. 
     832   ln_trdmld_instant = .false. 
     833/ 
     834!----------------------------------------------------------------------- 
     835!       namgap    level mean model-data gap ('key_diagap') 
     836!----------------------------------------------------------------------- 
     837!  ngap    time-step frequency of model-data gap computation 
     838!  nprg    time-step frequency of gap print in model output 
     839&namgap 
     840   ngap =  15 
     841   nprg =  10 
     842/ 
     843!----------------------------------------------------------------------- 
     844!       namspr  surface pressure diagnostic 
     845!----------------------------------------------------------------------- 
     846!  nmaxp   maximum of iterations for the solver 
     847!  epsp    absolute precision of the solver 
     848!  niterp  number of iteration done by the solver 
     849&namspr 
     850   nmaxp   =   1000 
     851   epsp    =  1.e-3 
     852   niterp  =    400 
     853/ 
     854!----------------------------------------------------------------------- 
     855!       namflo    float parameters                                      ("key_float") 
    642856!----------------------------------------------------------------------- 
    643857!  ln_rstflo   boolean term for float restart (true or false) 
     
    654868    ln_flork4 = .false. 
    655869/ 
    656 !!====================================================================== 
    657 !!              Surface Boundary Condition namelists 
    658 !! 
    659 !!       namsbc    surface boundary condition 
    660 !!       namsbc_ana    analytical         formulation of sbc 
    661 !!       namsbc_flx    flux               formulation of sbc 
    662 !!       namsbc_core   CORE bulk formulea formulation of sbc 
    663 !!       namsbc_clio   CLIO bulk formulea formulation of sbc 
    664 !!       namsbc_cpl    CouPLed            formulation of sbc 
    665 !!       namsbc_ssr    sea surface restoring term (for T and/or S) 
    666 !!====================================================================== 
    667 !----------------------------------------------------------------------- 
    668 !       namsbc       Surface Boundary Condition (surface module) 
    669 !----------------------------------------------------------------------- 
    670 !  nn_fsbc          frequency of surface boundary condition computation  
    671 !                   (also correspond to the frequency of sea-ice model call) 
    672 ! 
    673 !  ln_ana        analytical formulation (fill namsbc_ana )  
    674 !  ln_flx        flux formulation       (fill namsbc_flx ) (overwritten by key_sbc_flux     ) 
    675 !  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core) (overwritten by key_sbc_bulk_clio) 
    676 !  ln_blk_core   CORE bulk formulation  (fill namsbc_clio) (overwritten by key_sbc_bulk_core) 
    677 !  ln_cpl        Coupled formulation    (fill namsbc_cpl ) (overwritten by key_sbc_coupled  ) 
    678 ! 
    679 !  nn_ice        =0 no ice boundary condition   , 
    680 !                =1 observed ice-cover          , 
    681 !                =2 ice-model ice-cover (overwritten by key_ice_lim) 
    682 ! 
    683 !  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
    684 ! 
    685 !  ln_ssr       Sea Surface restoring on SST and/or SSS (fill namsbc_ssr) 
    686 !               (overwritten by key_dta_sst and/or key_dtasss) 
    687 ! 
    688 !  nn_fwb       FreshWater Budget: =0 no control                              ,  
    689 !                                  =1 annual global mean of emp set to zero   , 
    690 !                                  =2 global mean of emp set to zero at each nn_fsbc time step 
    691 &namsbc 
    692    nn_fsbc     = 0 
    693    ln_ana      = .true. 
    694    ln_flx      = .false. 
    695    ln_blk_clio = .false. 
    696    ln_blk_core = .false. 
    697    ln_cpl      = .false. 
    698    nn_ice      = 0 
    699    ln_dm2dc    = .false. 
    700    ln_rnf      = .false. 
    701    ln_ssr      = .false. 
    702    nn_fwb      = 0 
    703 / 
    704 !----------------------------------------------------------------------- 
    705 !       namsbc_rnf   runoffs namelist surface boundary condition 
    706 !----------------------------------------------------------------------- 
    707 !nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
    708 !rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
    709 !rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
    710 !cn_dir      Root directory for location of ssr files 
    711 !sn_rnf      informations about the runoff file to be read 
    712 !sn_cnf      informations about the runoff mouth file to be read 
    713 !  THE ORDER OF THE FILES MATTER: 
    714 !               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
    715 !               !                    !   (if <0  months) !   name     !     (logical)      ! 
    716 &namsbc_rnf 
    717    nn_runoff  = 0 
    718    rn_hrnf    = 0.e0 
    719    rn_avt_rnf = 1.e-3 
    720    cn_dir     = './' 
    721    sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .TRUE.   ,  1  ,  0      
    722    sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .FALSE.  ,  1  ,  0               
    723 / 
    724 !----------------------------------------------------------------------- 
    725 !       namsbc_ana   analytical surface boundary condition 
    726 !----------------------------------------------------------------------- 
    727 !  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
    728 !  rn_utau0    uniform value used as default wind stress -u 
    729 !  rn_vtau0    uniform value used as default wind stress -v 
    730 !  rn_q0       uniform value used as default surface heat flux 
    731 !  rn_qsr0     uniform value used as default solar radiation flux 
    732 !  rn_emp0     uniform value used as default surface freswater budget (E-P) 
    733 &namsbc_ana 
    734    nn_tau000  =   100 
    735    rn_utau0   =   0.1e0 
    736    rn_vtau0   =   0.e0 
    737    rn_q0      =   0.e0 
    738    rn_qsr0    =   0.e0 
    739    rn_emp0    =   0.e0 
    740 / 
    741 !----------------------------------------------------------------------- 
    742 !       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
    743 !----------------------------------------------------------------------- 
    744 !  cn_dir   directory in which the model is executed 
    745 ! 
    746 !  THE ORDER OF THE FILES MATTER: 
    747 !               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    748 !               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    749 &namsbc_flx 
    750    cn_dir       = './' 
    751    sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .FALSE.       ,    0   ,     0 
    752    sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .FALSE.       ,    0   ,     0 
    753    sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .FALSE.       ,    0   ,     0 
    754    sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .FALSE.       ,    0   ,     0 
    755    sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .FALSE.       ,    0   ,     0 
    756  
    757 /       
    758 !----------------------------------------------------------------------- 
    759 !       namsbc_clio  CLIO bulk formulea 
    760 !----------------------------------------------------------------------- 
    761 !  cn_dir   directory in which the model is executed 
    762 ! 
    763 !  THE ORDER OF THE FILES MATTER: 
    764 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    765 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    766 &namsbc_clio 
    767    cn_dir       = './' 
    768    sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .FALSE.        ,    1   ,    0     
    769    sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .FALSE.        ,    1   ,    0     
    770    sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .FALSE.        ,    1   ,    0       
    771    sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .FALSE.        ,    1   ,    0      
    772    sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .FALSE.        ,    1   ,    0     
    773    sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .FALSE.        ,    1   ,    0       
    774    sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .FALSE.        ,    1   ,    0        
    775 / 
    776 !----------------------------------------------------------------------- 
    777 !       namsbc_core  CORE bulk formulea 
    778 !----------------------------------------------------------------------- 
    779 !  cn_dir         directory in which the model is executed 
    780 !  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
    781 !  alpha_precip   multiplication factor for precipitation (total & snow) 
    782 ! 
    783 !  THE ORDER OF THE FILES MATTER: 
    784 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    785 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    786 &namsbc_core 
    787    cn_dir       = './' 
    788    ln_2m        = .FALSE. 
    789    alpha_precip = 1. 
    790    sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .FALSE.        ,    1   ,    0     
    791    sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .FALSE.        ,    1   ,    0     
    792    sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .FALSE.        ,    1   ,    0       
    793    sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .FALSE.        ,    1   ,    0       
    794    sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .FALSE.        ,    1   ,    0      
    795    sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .FALSE.        ,    1   ,    0     
    796    sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .FALSE.        ,    1   ,    0       
    797    sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .FALSE.        ,    1   ,    0        
    798 / 
    799 !----------------------------------------------------------------------- 
    800 !       namsbc_cpl    coupled ocean/atmosphere model  (#ifdef "key_coupled") 
    801 !----------------------------------------------------------------------- 
    802 &namsbc_cpl 
    803 / 
    804 !----------------------------------------------------------------------- 
    805 !       namsbc_ssr   surface boundary condition : sea surface restoring 
    806 !----------------------------------------------------------------------- 
    807 !  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
    808 !  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
    809 !             (=1, salt flux, concentration/dillution emps only) 
    810 !             (=2, volume flux, both emp and emps are updated)  
    811 !  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
    812 !  deds       magnitude of the damping on salinity [mm/day] 
    813 !  THE ORDER OF THE FILES MATTER: 
    814 !               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    815 !               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    816 &namsbc_ssr 
    817    nn_sstr      = 0 
    818    nn_sssr      = 0 
    819    dqdt         = -40. 
    820    deds         =  27.7 
    821    sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .FALSE.       ,    0   ,     0 
    822    sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .TRUE.        ,    0   ,     0 
    823 /       
     870!----------------------------------------------------------------------- 
     871!       namptr   Poleward Transport Diagnostic 
     872!----------------------------------------------------------------------- 
     873!  ln_diaptr  logical flag for Poleward transport computation 
     874!  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
     875!             need input basins mask file named "subbasins.nc" 
     876!  nf_ptr     Frequency of computation 
     877&namptr 
     878   ln_diaptr = .false. 
     879   ln_subbas = .false. 
     880   nf_ptr    =  15 
     881/ 
  • trunk/CONFIG/GYRE_LOBSTER/EXP00/namelist

    r889 r907  
    1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! OPA namelist :  model option and parameter input 
    3 ! ------------- 
    4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    5 ! 
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/OPA  :       1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 
     3!! namelists         2 - Domain           (nam_zgr, nam_zgr_sco, namdom) 
     4!!                   3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
     5!!                                         namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 
     6!!                   4 - lateral boundary (namlbc, namcla, namobc) 
     7!!                   5 - bottom  boundary (nambfr, nambbc, nambbl) 
     8!!                   6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp) 
     9!!                   7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) 
     10!!                   8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm) 
     11!!                   9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr) 
     12!!                   9 - miscellaneous    (namsol) 
     13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     14!  CAUTION: some scripts does not support CAPITALs for logical definition  
     15!  *******       use .true. or .false. and NOT .TRUE. or .FALSE. 
     16 
     17!!====================================================================== 
     18!!                   ***  Run management namelists  *** 
     19!!====================================================================== 
     20!!   namrun            parameters of the run 
     21!!   nam_mpp           Massively Parallel Processing 
     22!!   nam_mpp_dyndist    
     23!!   nam_ctl           Control prints & Benchmark 
     24!!====================================================================== 
     25 
    626!----------------------------------------------------------------------- 
    727!       namrun   parameters of the run 
     
    1030!  cexper     experience name for vairmer format 
    1131!  ln_rstart  boolean term for restart (true or false) 
    12 !  nrstdt     control of the restart timestep: 
    13 !                = 0 restart, do not control nit000 in the restart file. 
    14 !                = 1 restart, control nit000 in the restart file. Do not 
    15 !                    use the date in the restart file (use ndate0 in namelist) 
    16 !                = 2 restart, control nit000 in the restart file, use the date 
    17 !                    in the restart file. ndate0 in the namelist is ignored. 
     32!  nrstdt     restart control = 0 restart, do not control nit000 in the restart file. 
     33!             !               = 1 restart, control nit000 in the restart file. Do not 
     34!             !                   use the date in the restart file (use ndate0 in namelist) 
     35!             !               = 2 restart, control nit000 in the restart file, use the date 
     36!             !                   in the restart file. ndate0 in the namelist is ignored. 
    1837!  nit000     number of the first time step 
    1938!  nitend     number of the last time step 
     
    2443!  nwrite     frequency of OUTPUT file 
    2544!  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 
    26 ! 
    27 !  CAUTION: for usual run scripts, logical value of 
    28 !  *******  ln_rstart must be .true. or .false. 
    29 !                     and NOT .TRUE. or .FALSE. 
    3045&namrun 
    3146   no         =       0 
     
    3954   ninist     =       0 
    4055   nstock     =    4320  
    41    nwrite     =      360  
    42    ln_dimgnnn = .FALSE. 
     56   nwrite     =     360  
     57   ln_dimgnnn = .false. 
     58/ 
     59!----------------------------------------------------------------------- 
     60!       nam_mpp      Massively Parallel Processing 
     61!----------------------------------------------------------------------- 
     62!  c_mpi_send         mpi send/recieve type 
     63!                      = 'S'  : standard blocking send 
     64!                      = 'B'  : buffer blocking send 
     65!                      = 'I'  : immediate non-blocking send 
     66&nam_mpp 
     67   c_mpi_send =  'S' 
     68/ 
     69!----------------------------------------------------------------------- 
     70!       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist") 
     71!----------------------------------------------------------------------- 
     72!  jpni   number of processors following i 
     73!  jpnj   number of processors following j 
     74!  jpnij  number of local domains 
     75&nam_mpp_dyndist 
     76   jpni  = 1 
     77   jpnj  = 1 
     78   jpnij = 1 
    4379/ 
    4480!----------------------------------------------------------------------- 
     
    69105   nbit_cmp =     0 
    70106/ 
    71 !----------------------------------------------------------------------- 
    72 !       nam_mpp      Massively Parallel Processing 
    73 !----------------------------------------------------------------------- 
    74 !  c_mpi_send         mpi send/recieve type 
    75 !                      = 'S'  : standard blocking send 
    76 !                      = 'B'  : buffer blocking send 
    77 !                      = 'I'  : immediate non-blocking send 
    78 &nam_mpp 
    79    c_mpi_send =  'S' 
    80 / 
    81 !----------------------------------------------------------------------- 
    82 !       nam_mpp_dyndist      Massively Parallel Distribution 
    83 !              #ifdef ( key_agrif && key_mpp_dyndist ) 
    84 !----------------------------------------------------------------------- 
    85 !  jpni   number of processors following i 
    86 !  jpnj   number of processors following j 
    87 !  jpnij  number of local domains 
    88 &nam_mpp_dyndist 
    89    jpni  = 1 
    90    jpnj  = 1 
    91    jpnij = 1 
    92 / 
     107 
     108!!====================================================================== 
     109!!                      ***  Domain namelists  *** 
     110!!====================================================================== 
     111!!   nam_zgr       vertical coordinate 
     112!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate 
     113!!   namdom        space and time domain (bathymetry, mesh, timestep) 
     114!!====================================================================== 
     115 
    93116!----------------------------------------------------------------------- 
    94117!       nam_zgr       vertical coordinate 
    95118!----------------------------------------------------------------------- 
    96 !  ln_zco     z-coordinate - full steps      (T/F) 
     119!  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined) 
    97120!  ln_zps     z-coordinate - partial steps   (T/F) 
    98121!  ln_sco     s- or hybrid z-s-coordinate    (T/F) 
     
    118141/ 
    119142!----------------------------------------------------------------------- 
    120 !       nam_traadv   advection scheme for tracer (option not control by CPP keys) 
    121 !----------------------------------------------------------------------- 
    122 !  ln_traadv_cen2     2nd order centered scheme    (default T) 
    123 !  ln_traadv_tvd      TVD scheme                   (default F) 
    124 !  ln_traadv_muscl    MUSCL scheme                 (default F) 
    125 !  ln_traadv_muscl2   MUSCL2 scheme                (default F) 
    126 !  ln_traadv_ubs      UBS scheme                   (default F) 
    127 !  ln_traadv_qck      QUICKEST scheme              (default F) 
    128 &nam_traadv 
    129    ln_traadv_cen2   =  .false. 
    130    ln_traadv_tvd    =  .true. 
    131    ln_traadv_muscl  =  .false. 
    132    ln_traadv_muscl2 =  .false. 
    133    ln_traadv_ubs    =  .false. 
    134    ln_traadv_qck    =  .false. 
    135 / 
    136 !----------------------------------------------------------------------- 
    137 !       nam_traldf   lateral diffusion scheme for tracer (option not control by CPP keys) 
    138 !----------------------------------------------------------------------- 
    139 !  Type of the operator :  
    140 !     ln_traldf_lap    laplacian operator          (default T) 
    141 !     ln_traldf_bilap  bilaplacian operator        (default F) 
    142 !  Direction of action  :  
    143 !     ln_traldf_level  iso-level                   (default F) 
    144 !     ln_traldf_hor    horizontal (geopotential)   (default F)^** 
    145 !     ln_traldf_iso    iso-neutral                 (default T)^* 
    146 !  Coefficient 
    147 !     aht0    horizontal eddy diffusivity for tracers (m2/s) 
    148 !     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
    149 !     aeiv0   eddy induced velocity coefficient (m2/s) 
    150 ! ^* require key_ldfslp to compute the direction of the lateral diffusion 
    151 ! ^** require key_ldfslp in s-coordinate 
    152 &nam_traldf 
    153    ln_traldf_lap    =  .true. 
    154    ln_traldf_bilap  =  .false. 
    155    ln_traldf_level  =  .false. 
    156    ln_traldf_hor    =  .false. 
    157    ln_traldf_iso    =  .true. 
    158    aht0    =  1000. 
    159    ahtb0   =     0. 
    160    aeiv0   =     0. 
    161 / 
    162 !----------------------------------------------------------------------- 
    163 !       nam_dynldf   lateral diffusion on momentum 
    164 !----------------------------------------------------------------------- 
    165 !  Type of the operator :  
    166 !     ln_dynldf_lap    laplacian operator          (default T) 
    167 !     ln_dynldf_bilap  bilaplacian operator        (default F) 
    168 !  Direction of action  :  
    169 !     ln_dynldf_level  iso-level                   (default F) 
    170 !     ln_dynldf_hor    horizontal (geopotential)   (default F)^** 
    171 !     ln_dynldf_iso    iso-neutral                 (default T)^* 
    172 !  Coefficient 
    173 !  ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
    174 !  ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
    175 &nam_dynldf 
    176    ln_dynldf_lap    =  .true. 
    177    ln_dynldf_bilap  =  .false. 
    178    ln_dynldf_level  =  .false. 
    179    ln_dynldf_hor    =  .true. 
    180    ln_dynldf_iso    =  .false. 
    181    ahm0    = 100000. 
    182    ahmb0   =      0. 
    183 / 
    184 !----------------------------------------------------------------------- 
    185 !       namflg   algorithm flags (algorithm not control by CPP keys) 
    186 !----------------------------------------------------------------------- 
    187 !  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
    188 !                                                  centered      time scheme  (F) 
    189 !   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
    190 &namflg 
    191    ln_dynhpg_imp   =  .false. 
    192    nn_dynhpg_rst   =  0 
    193 / 
    194 !----------------------------------------------------------------------- 
    195 !       nam_dynhpg   Hydrostatic pressure gradient option 
    196 !----------------------------------------------------------------------- 
    197 !  type of pressure gradient scheme (choose one only!) 
    198 !     ln_hpg_zco    z-coordinate - full steps                   (default T) 
    199 !     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
    200 !     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
    201 !     ln_hpg_hel    s-coordinate (helsinki modification) 
    202 !     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
    203 !     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
    204 !     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
    205 !  parameters 
    206 !    gamm          weighting coefficient (wdj scheme) 
    207 &nam_dynhpg 
    208    ln_hpg_zco = .true. 
    209    ln_hpg_zps = .false. 
    210    ln_hpg_sco = .false. 
    211    ln_hpg_hel = .false. 
    212    ln_hpg_wdj = .false. 
    213    ln_hpg_djc = .false. 
    214    ln_hpg_rot = .false. 
    215    gamm       = 0.e0 
    216 / 
    217 !----------------------------------------------------------------------- 
    218 !       nam_dynadv   option of physics/algorithm (not control by CPP keys) 
    219 !----------------------------------------------------------------------- 
    220 !  ln_dynadv_vec      vector form flag 
    221 !  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
    222 !  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
    223 &nam_dynadv    
    224    ln_dynadv_vec   = .TRUE. 
    225    ln_dynadv_cen2  = .FALSE. 
    226    ln_dynadv_ubs   = .FALSE. 
    227  
    228 !----------------------------------------------------------------------- 
    229 !       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
    230 !----------------------------------------------------------------------- 
    231 !  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
    232 !  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
    233 !  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
    234 !  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
    235 &nam_dynvor 
    236    ln_dynvor_ene = .TRUE. 
    237    ln_dynvor_ens = .FALSE. 
    238    ln_dynvor_mix = .FALSE. 
    239    ln_dynvor_een = .FALSE. 
    240 / 
    241 !----------------------------------------------------------------------- 
    242 !       namalb   albedo parameters 
    243 !----------------------------------------------------------------------- 
    244 !  cgren    correction of the snow or ice albedo to take into account 
    245 !  albice   albedo of melting ice in the arctic and antarctic 
    246 !  alphd    coefficients for linear interpolation used to compute albedo 
    247 !           between two extremes values (Pyane, 1972) 
    248 !  alphc     "                                         " 
    249 !  alphdi    "                                         " 
    250 &namalb 
    251    cgren    =      0.06 
    252    albice   =      0.5 
    253    alphd    =      0.80 
    254    alphc    =      0.65 
    255    alphdi   =      0.72 
    256 / 
    257 !----------------------------------------------------------------------- 
    258143!       namdom   space and time domain (bathymetry, mesh, timestep) 
    259144!----------------------------------------------------------------------- 
    260 !  ntopo      = 0/1 ,compute/read the bathymetry file 
    261 !               (mbathy, nb of T-ocean levels) 
     145!  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 
    262146!  e3zps_min  the thickness of the partial step is set larger than the 
    263 !  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t 
    264 !                (N.B. 0<e3zps_rat<1) 
     147!  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 
    265148!  nmsh       =1 create a mesh file (coordinates, scale factors, masks) 
    266 !  nacc       the acceleration of convergence method 
    267 !             = 0, no acceleration, rdt = rdttra 
    268 !             = 1, acceleration used, rdt < rdttra(k) 
     149!  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra 
     150!                                                  = 1,     used, rdt < rdttra(k) 
    269151!  atfp       asselin time filter parameter 
    270152!  rdt        time step for the dynamics (and tracer if nacc=0) 
     
    273155!  rdth       depth variation of tracer time step 
    274156!  rdtbt      barotropic time step (for the time splitting algorithm) 
    275 !  nclosea    0 no closed sea  
    276 !             1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
     157!  nclosea    =0 remove the closed sea from the global domain (orca configuration) 
     158!             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
    277159&namdom 
    278160   ntopo     =     0 
     
    289171   nclosea   =     0 
    290172/ 
    291 !----------------------------------------------------------------------- 
    292 !       namfwb   freshwater budget correction 
    293 !----------------------------------------------------------------------- 
    294 !  ln_fwb     logical flag for freshwater budget correction (0 annual mean) 
    295 &namfwb 
    296    ln_fwb    = .false. 
    297 / 
    298 !----------------------------------------------------------------------- 
    299 !       namptr   Poleward Transport Diagnostic 
    300 !----------------------------------------------------------------------- 
    301 !  ln_diaptr  logical flag for Poleward transport computation 
    302 !  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
    303 !             need input basins mask file named "subbasins.nc" 
    304 !  nf_ptr     Frequency of computation 
    305 &namptr 
    306    ln_diaptr = .false. 
    307    ln_subbas = .false. 
    308    nf_ptr    =  15 
    309 / 
    310 !----------------------------------------------------------------------- 
    311 !       namcro   cross land advection 
    312 !----------------------------------------------------------------------- 
    313 !  n_cla   advection between 2 ocean pts separates by land  
    314 &namcla 
    315    n_cla   = 0 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 !       namzdf   vertical physics 
    319 !----------------------------------------------------------------------- 
    320 !  ln_zdfevd  enhanced vertical diffusion         (default T) 
    321 !  ln_zdfnpc  Non-Penetrative Convection          (default T) 
    322 !  avm0       vertical eddy viscosity for the dynamic (m2/s) 
    323 !  avt0       vertical eddy diffusivity for tracers (m2/s) 
    324 !  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
    325 !  n_evdm     = 0  apply enhanced mixing on tracer only 
    326 !             = 1  apply enhanced mixing on both tracer and momentum 
    327 !  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
    328 !                               (=F)  euler backward (F) 
    329 !  n_zdfexp   number of sub-timestep for time splitting scheme 
    330 &namzdf 
    331    ln_zdfevd = .true. 
    332    ln_zdfnpc = .false. 
    333    avm0      = 1.2e-4 
    334    avt0      = 1.2e-5 
    335    avevd     =   100. 
    336    n_evdm    =     1 
    337    ln_zdfexp =  .false. 
    338    n_zdfexp  =     3 
    339 / 
    340 !----------------------------------------------------------------------- 
    341 !       namnpc   vnon penetrative convection 
    342 !----------------------------------------------------------------------- 
    343 !  nnpc1   non penetrative convective scheme frequency 
    344 !  nnpc2   non penetrative convective scheme print frequency 
    345 &namnpc 
    346    nnpc1  =      1 
    347    nnpc2  =    365 
    348 / 
    349 !----------------------------------------------------------------------- 
    350 !       nambbl   bottom boundary layer scheme 
    351 !----------------------------------------------------------------------- 
    352 !  atrbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s) 
    353 &nambbl 
    354    atrbbl = 10000. 
    355 / 
    356 !----------------------------------------------------------------------- 
    357 !       namric   richardson number dependent vertical diffusion 
    358 !                ( #ifdef "key_zdfrichardson" ) 
    359 !----------------------------------------------------------------------- 
    360 !  avmri   maximum value of the vertical viscosity 
    361 !  alp     coefficient of the parameterization 
    362 !  nric    coefficient of the parameterization 
    363 &namric 
    364    avmri = 100.e-4 
    365    alp   =      5. 
    366    nric  =       2 
    367 / 
    368 !----------------------------------------------------------------------- 
    369 !       namtke   turbulent eddy kinetic dependent vertical diffusion 
    370 !                ( #ifdef "key_zdftke" ) 
    371 !----------------------------------------------------------------------- 
    372 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    373 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    374 !  ediss    coef. of the Kolmogoroff dissipation   
    375 !  ebb      coef. of the surface input of tke 
    376 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    377 !  emin     minimum value of tke (m^2/s^2) 
    378 !  emin0    surface minimum value of tke (m^2/s^2) 
    379 !  nitke    number of restart iterative loops 
    380 !  ri_c     critic richardson number 
    381 !  nmxl     flag on mixing length used 
    382 !           = 0 bounded by the distance to surface and bottom 
    383 !           = 1 bounded by the local vertical scale factor 
    384 !           = 2 first vertical derivative of mixing length bounded by 1 
    385 !  npdl     flag on prandtl number 
    386 !           = 0 no vertical prandtl number (avt=avm) 
    387 !           = 1 prandtl number function of richarson number (avt=pdl*avm) 
    388 !           = 2 same as = 1 but a shapiro filter is applied on pdl 
    389 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    390 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    391 &namtke 
    392    ln_rstke = .false. 
    393    ediff =       0.1 
    394    ediss =       0.7 
    395    ebb   =      3.75 
    396    efave =        1. 
    397    emin  =     1.e-5  
    398    emin0 =     1.e-4 
    399    nitke =        50 
    400    nmxl  =         2 
    401    npdl  =         1 
    402    navb  =         0 
    403 / 
    404 !----------------------------------------------------------------------- 
    405 !       namkpp   K-Profile Parameterization dependent vertical diffusion 
    406 !                ( #ifdef "key_zdfkpp" ) 
    407 !----------------------------------------------------------------------- 
    408 !  ln_kpprimix   shear instability mixing  (default T) 
    409 !  difmiw        constant internal wave viscosity (m2/s) 
    410 !  difsiw        constant internal wave diffusivity (m2/s) 
    411 !  Riinfty       local Richardson Number limit for shear instability 
    412 !  difri         maximum shear mixing at Rig = 0    (m2/s) 
    413 !  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
    414 !  difcon        maximum mixing in interior convection (m2/s) 
    415 !  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
    416 !  navb          = 0/1 flag for constant or profile background avt 
    417 &namkpp 
    418    ln_kpprimix  = .true. 
    419    difmiw       =  1.e-04 
    420    difsiw       =  0.1e-04 
    421    Riinfty      =  0.8 
    422    difri        =  0.0050 
    423    bvsqcon      = -0.01e-07 
    424    difcon       =  1. 
    425    navb         =  0 
    426    nave         =  1 
    427 / 
    428 !----------------------------------------------------------------------- 
    429 !       namddm   double diffusive mixing parameterization 
    430 !----------------------------------------------------------------------- 
    431 !   avts    maximum avs for dd mixing  
    432 !   hsbfr   heat/salt buoyancy flux ratio 
    433 &namddm 
    434       avts  = 1.e-4 
    435       hsbfr = 1.6 
    436 / 
    437 !----------------------------------------------------------------------- 
    438 !       namlbc   lateral momentum boundary condition 
    439 !----------------------------------------------------------------------- 
    440 !  shlat   lateral boundary condition on velocity 
    441 !                   shlat = 0 , free slip 
    442 !               0 < shlat < 2 , partial slip 
    443 !                   shlat = 2 , no slip 
    444 !               2 < shlat     , strong slip 
    445 &namlbc 
    446    shlat  =      0. 
    447 / 
    448 !----------------------------------------------------------------------- 
    449 !       nambfr   bottom friction 
    450 !----------------------------------------------------------------------- 
    451 !  nbotfr  type of bottom friction  
    452 !                  nbotfr = 0 , no slip  
    453 !                  nbotfr = 1 , linear friction 
    454 !                  nbotfr = 2 , nonlinear friction 
    455 !                  nbotfr = 3 , free slip 
    456 !  bfri1   bottom drag coefficient (linear case) 
    457 !  bfri2   bottom drag coefficient (non linear case) 
    458 !  bfeb2   bottom turbulent kinetic energy  (m^2/s^2) 
    459 &nambfr 
    460    nbotfr =       2 
    461    bfri1  =   4.e-4 
    462    bfri2  =   1.e-3 
    463    bfeb2  =  2.5e-3 
    464 / 
    465 !----------------------------------------------------------------------- 
    466 !       nambbc   bottom temperature boundary condition 
    467 !----------------------------------------------------------------------- 
    468 !  ngeo_flux  = 0 no geothermal heat flux 
    469 !             = 1 constant geothermal heat flux 
    470 !             = 2 variable geothermal heat flux (read in geothermal_heating.nc)  
    471 !                 ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 
    472 !  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
    473 &nambbc 
    474    ngeo_flux =  0 
    475    ngeo_flux_const = 86.4e-3 
     173 
     174!!====================================================================== 
     175!!            ***  Surface Boundary Condition namelists  *** 
     176!!====================================================================== 
     177!!   namsbc        surface boundary condition 
     178!!   namsbc_ana    analytical         formulation 
     179!!   namsbc_flx    flux               formulation 
     180!!   namsbc_clio   CLIO bulk formulea formulation 
     181!!   namsbc_core   CORE bulk formulea formulation 
     182!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled") 
     183!!   namqsr        penetrative solar radiation 
     184!!   namsbc_rnf    river runoffs 
     185!!   namsbc_ssr    sea surface restoring term (for T and/or S) 
     186!!   namalb        albedo parameters 
     187!!====================================================================== 
     188 
     189!----------------------------------------------------------------------- 
     190!       namsbc       Surface Boundary Condition (surface module) 
     191!----------------------------------------------------------------------- 
     192!  nn_fsbc       frequency of surface boundary condition computation  
     193!                (= the frequency of sea-ice model call) 
     194!  ln_ana        analytical formulation (fill namsbc_ana )  
     195!  ln_flx        flux formulation       (fill namsbc_flx )  
     196!  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core)  
     197!  ln_blk_core   CORE bulk formulation  (fill namsbc_clio)  
     198!  ln_cpl        Coupled formulation    (fill namsbc_cpl )  
     199!  nn_ice        =0 no ice boundary condition   , 
     200!                =1 observed ice-cover          , 
     201!                =2 ice-model used                                      ("key_lim3" or "key_lim2) 
     202!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
     203!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr) 
     204!  nn_fwb       FreshWater Budget: =0 no control                              ,  
     205!                                  =1 annual global mean of e-p-r set to zero   , 
     206!                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step 
     207&namsbc 
     208   nn_fsbc     = 1 
     209   ln_ana      = .true. 
     210   ln_flx      = .false. 
     211   ln_blk_clio = .false. 
     212   ln_blk_core = .false. 
     213   ln_cpl      = .false. 
     214   nn_ice      = 0 
     215   ln_dm2dc    = .false. 
     216   ln_rnf      = .false. 
     217   ln_ssr      = .false. 
     218   nn_fwb      = 0 
     219/ 
     220!----------------------------------------------------------------------- 
     221!       namsbc_ana   analytical surface boundary condition 
     222!----------------------------------------------------------------------- 
     223!  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
     224!  rn_utau0    uniform value for the i-stress 
     225!  rn_vtau0    uniform value for the j-stress 
     226!  rn_q0       uniform value for the total heat flux 
     227!  rn_qsr0     uniform value for the solar radiation 
     228!  rn_emp0     uniform value for the freswater budget (E-P) 
     229&namsbc_ana 
     230   nn_tau000  =   100 
     231   rn_utau0   =   0.1e0 
     232   rn_vtau0   =   0.e0 
     233   rn_q0      =   0.e0 
     234   rn_qsr0    =   0.e0 
     235   rn_emp0    =   0.e0 
     236/ 
     237!----------------------------------------------------------------------- 
     238!       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
     239!----------------------------------------------------------------------- 
     240!  cn_dir   directory in which the model is executed 
     241! 
     242!  THE ORDER OF THE FILES MATTER: 
     243!               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     244!               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     245&namsbc_flx 
     246   cn_dir       = './' 
     247   sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0 
     248   sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0 
     249   sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0 
     250   sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0 
     251   sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0 
     252/       
     253!----------------------------------------------------------------------- 
     254!       namsbc_clio  CLIO bulk formulea 
     255!----------------------------------------------------------------------- 
     256!  cn_dir   directory in which the model is executed 
     257! 
     258!  THE ORDER OF THE FILES MATTER: 
     259!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     260!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     261&namsbc_clio 
     262   cn_dir       = './' 
     263   sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0     
     264   sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0     
     265   sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0       
     266   sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0      
     267   sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0     
     268   sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0       
     269   sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0        
     270/ 
     271!----------------------------------------------------------------------- 
     272!       namsbc_core  CORE bulk formulea 
     273!----------------------------------------------------------------------- 
     274!  cn_dir         directory in which the model is executed 
     275!  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
     276!  alpha_precip   multiplication factor for precipitation (total & snow) 
     277! 
     278!  THE ORDER OF THE FILES MATTER: 
     279!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     280!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     281&namsbc_core 
     282   cn_dir       = './' 
     283   ln_2m        = .false. 
     284   alpha_precip = 1. 
     285   sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0     
     286   sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0     
     287   sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0       
     288   sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0       
     289   sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0      
     290   sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0     
     291   sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0       
     292   sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0        
     293/ 
     294!----------------------------------------------------------------------- 
     295!       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled") 
     296!----------------------------------------------------------------------- 
     297&namsbc_cpl 
    476298/ 
    477299!----------------------------------------------------------------------- 
     
    489311/ 
    490312!----------------------------------------------------------------------- 
    491 !       namtdp   tracer newtonian damping ('key_tradmp') 
    492 !----------------------------------------------------------------------- 
    493 !  ndmp    type of damping in temperature and salinity  
    494 !          (='latitude', damping poleward of 'ndmp' degrees and function  
    495 !             of the distance-to-coast. Red and Med Seas as ndmp=-1) 
    496 !          (=-1 damping only in Med and Red Seas) 
    497 !  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
    498 !  nmldmp  type of damping in the mixed layer 
    499 !          (=0 damping throughout the water column) 
    500 !     (=1 no damping in the mixed layer defined by avt >5cm2/s ) 
    501 !     (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
    502 !  sdmp    surface time scale for internal damping (days) 
    503 !  bdmp    bottom time scale for internal damping (days) 
    504 !  hdmp    depth of transition between sdmp and bdmp (meters) 
    505 &namtdp 
    506    ndmp   =   -1 
    507    ndmpf  =    1 
    508    nmldmp =    1 
    509    sdmp   =  50. 
    510    bdmp   = 360. 
    511    hdmp   = 800. 
    512 / 
    513 !----------------------------------------------------------------------- 
    514 !       nameos   ocean physical parameters 
    515 !----------------------------------------------------------------------- 
    516 !  neos    type of equation of state and Brunt-Vaisala frequency 
    517 !          = 0, UNESCO (formulation of Jackett and McDougall (1994) 
    518 !                                         and of McDougall (1987) ) 
    519 !          = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    520 !          = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    521 !                               with rau0=1020 set in parcst routine 
    522 !  ralpha  thermal expension coefficient (linear equation of state) 
    523 !  rbeta   saline  expension coefficient (linear equation of state) 
    524 &nameos 
    525    neos   =      2 
    526    ralpha =  2.e-4 
    527    rbeta  =  7.7e-4 
    528 / 
    529 !----------------------------------------------------------------------- 
    530 !       namsol   elliptic solver / island / free surface  
    531 !----------------------------------------------------------------------- 
    532 !  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
    533 !                            (=2 successive-over-relaxation: sor) 
    534 !                            (=3 FETI: fet, all require "key_feti" defined) 
    535 !                            (=4 sor with extra outer halo) 
    536 !  nsol_arp  absolute/relative (0/1) precision convergence test 
    537 !  nmin      minimum of iterations for the SOR solver 
    538 !  nmax      maximum of iterations for the SOR solver 
    539 !  nmod      frequency of test for the SOR solver 
    540 !  eps       absolute precision of the solver 
    541 !  resmax    absolute precision for the SOR solver 
    542 !  sor       optimal coefficient for SOR solver 
    543 !  epsisl    absolute precision on stream function solver 
    544 !  nmisl     maximum pcg iterations for island 
    545 !  rnu       strength of the additional force used in free surface b.c. 
    546 &namsol 
    547    nsolv     =      2 
    548    nsol_arp  =      0 
    549    nmin      =    100 
    550    nmax      =    800 
    551    nmod      =     10 
    552    eps       =  1.E-6 
    553    resmax    = 1.E-10 
    554    sor       =   1.96 
    555    epsisl    = 1.e-10 
    556    nmisl     =   4000 
    557    rnu       =     1. 
    558 / 
    559 !======================================================================= 
    560 !   Diagnostics namelists 
    561 !       namtrd    dynamics and/or tracer trends 
    562 !       namgap    level mean model-data gap 
    563 !       namznl    zonal mean heat & freshwater fluxes computation 
    564 !       namspr    surface pressure in rigid-lid 
    565 !======================================================================= 
    566 !----------------------------------------------------------------------- 
    567 !       namtrd    diagnostics on dynamics and/or tracer trends 
    568 !                         ('key_trdyn' and/or 'key_trdtra') 
    569 !                 or mixed-layer trends ('key_trdmld') 
    570 !----------------------------------------------------------------------- 
    571 !  ntrd              time step frequency dynamics and tracers trends 
    572 !  nctls             control surface type in mixed-layer trends (0,1 or n<jpk) 
    573 !  ln_trdmld_restart restart for ML diagnostics 
    574 !  ucf               unit conversion factor (=1 -> /seconds | =86400. -> /day) 
    575 !  ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 
    576 &namtrd 
    577    ntrd  = 365 
    578    nctls =   0 
    579    ln_trdmld_restart = .false. 
    580    ucf   =  1. 
    581    ln_trdmld_instant = .false. 
    582 / 
    583 !----------------------------------------------------------------------- 
    584 !       namgap    level mean model-data gap ('key_diagap') 
    585 !----------------------------------------------------------------------- 
    586 !  ngap    time-step frequency of model-data gap computation 
    587 !  nprg    time-step frequency of gap print in model output 
    588 &namgap 
    589    ngap =  15 
    590    nprg =  10 
    591 / 
    592 !----------------------------------------------------------------------- 
    593 !       namznl    zonal mean heat & freshwater fluxes computation 
    594 !                 (#ifdef "key_diaznl") 
    595 !----------------------------------------------------------------------- 
    596 !  nfznl   time-step frequency of zonal mean fluxes computation 
    597 &namznl 
    598    nfznl =  15 
    599 / 
    600 !----------------------------------------------------------------------- 
    601 !       namspr  surface pressure diagnostic 
    602 !----------------------------------------------------------------------- 
    603 !  nmaxp   maximum of iterations for the solver 
    604 !  epsp    absolute precision of the solver 
    605 !  niterp  number of iteration done by the solver 
    606 &namspr 
    607    nmaxp   =   1000 
    608    epsp    =  1.e-3 
    609    niterp  =    400 
     313!       namsbc_rnf   runoffs namelist surface boundary condition 
     314!----------------------------------------------------------------------- 
     315!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
     316!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
     317!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
     318!cn_dir      Root directory for location of ssr files 
     319!sn_rnf      informations about the runoff file to be read 
     320!sn_cnf      informations about the runoff mouth file to be read 
     321!  THE ORDER OF THE FILES MATTER: 
     322!               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
     323!               !                    !   (if <0  months) !   name     !     (logical)      ! 
     324&namsbc_rnf 
     325   nn_runoff  = 0 
     326   rn_hrnf    = 0.e0 
     327   rn_avt_rnf = 1.e-3 
     328   cn_dir     = './' 
     329   sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0      
     330   sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0               
     331/ 
     332!----------------------------------------------------------------------- 
     333!       namsbc_ssr   surface boundary condition : sea surface restoring 
     334!----------------------------------------------------------------------- 
     335!  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
     336!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
     337!             (=1, salt flux, concentration/dillution emps only) 
     338!             (=2, volume flux, both emp and emps are updated)  
     339!  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
     340!  deds       magnitude of the damping on salinity [mm/day] 
     341!  THE ORDER OF THE FILES MATTER: 
     342!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     343!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     344&namsbc_ssr 
     345   nn_sstr      = 0 
     346   nn_sssr      = 0 
     347   dqdt         = -40. 
     348   deds         =  27.7 
     349   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0 
     350   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0 
     351/       
     352!----------------------------------------------------------------------- 
     353!       namalb   albedo parameters 
     354!----------------------------------------------------------------------- 
     355!  cgren    correction of the snow or ice albedo to take into account 
     356!  albice   albedo of melting ice in the arctic and antarctic 
     357!  alphd    coefficients for linear interpolation used to compute albedo 
     358!           between two extremes values (Pyane, 1972) 
     359!  alphc     "                                         " 
     360!  alphdi    "                                         " 
     361&namalb 
     362   cgren    =      0.06 
     363   albice   =      0.5 
     364   alphd    =      0.80 
     365   alphc    =      0.65 
     366   alphdi   =      0.72 
     367/ 
     368 
     369!!====================================================================== 
     370!!               ***  Lateral boundary condition  *** 
     371!!====================================================================== 
     372!!   namlbc        lateral momentum boundary condition 
     373!!   namcla        cross land advection 
     374!!   namobc        open boundaries parameters                           ("key_obc") 
     375!!====================================================================== 
     376 
     377!----------------------------------------------------------------------- 
     378!       namlbc   lateral momentum boundary condition 
     379!----------------------------------------------------------------------- 
     380!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip 
     381!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip 
     382&namlbc 
     383   shlat  =      0. 
     384/ 
     385!----------------------------------------------------------------------- 
     386!       namcla   cross land advection 
     387!----------------------------------------------------------------------- 
     388!  n_cla   advection between 2 ocean pts separates by land  
     389&namcla 
     390   n_cla   = 0 
    610391/ 
    611392!----------------------------------------------------------------------- 
     
    639420    ln_vol_cst  = .false. 
    640421/ 
    641 !----------------------------------------------------------------------- 
    642 !       namflo    float parameters (#ifdef key_float) 
     422 
     423!!====================================================================== 
     424!!                 ***  Bottom boundary condition  *** 
     425!!====================================================================== 
     426!!   nambfr        bottom friction 
     427!!   nambbc        bottom temperature boundary condition                ("key_trabbc") 
     428!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv") 
     429!!====================================================================== 
     430!----------------------------------------------------------------------- 
     431!       nambfr   bottom friction 
     432!----------------------------------------------------------------------- 
     433!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction 
     434!                                     = 3 : free slip  ;  = 2 : non linear friction 
     435!  bfri1    bottom drag coefficient (linear case) 
     436!  bfri2    bottom drag coefficient (non linear case) 
     437!  bfeb2    bottom turbulent kinetic energy  (m2/s2) 
     438&nambfr 
     439   nbotfr =       2 
     440   bfri1  =   4.e-4 
     441   bfri2  =   1.e-3 
     442   bfeb2  =  2.5e-3 
     443/ 
     444!----------------------------------------------------------------------- 
     445!       nambbc   bottom temperature boundary condition 
     446!----------------------------------------------------------------------- 
     447!  ngeo_flux  = 0 no geothermal heat flux 
     448!             = 1 constant geothermal heat flux 
     449!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 
     450!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
     451&nambbc 
     452   ngeo_flux =  0 
     453   ngeo_flux_const = 86.4e-3 
     454/ 
     455!----------------------------------------------------------------------- 
     456!       nambbl   bottom boundary layer scheme 
     457!----------------------------------------------------------------------- 
     458!                                !  diffusive bbl                       ("key_trabbl") 
     459!                                !  advective bbl                       ("key_trabbl_adv") 
     460!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s) 
     461&nambbl 
     462   atrbbl = 10000. 
     463/ 
     464 
     465!!====================================================================== 
     466!!                        Tracer (T & S ) namelists 
     467!!====================================================================== 
     468!!   nameos        equation of state 
     469!!   nam_traadv    advection scheme 
     470!!   nam_traldf    lateral diffusion scheme 
     471!!   namtdp        tracer newtonian damping                             ("key_tradmp") 
     472!!====================================================================== 
     473 
     474!----------------------------------------------------------------------- 
     475!       nameos   ocean physical parameters 
     476!----------------------------------------------------------------------- 
     477!  neos    type of equation of state and Brunt-Vaisala frequency 
     478!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
     479!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
     480!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
     481!  ralpha  thermal expension coefficient (neos= 1 or 2) 
     482!  rbeta   saline  expension coefficient (neos= 2) 
     483&nameos 
     484   neos   =      2 
     485   ralpha =  2.e-4 
     486   rbeta  =  7.7e-4 
     487/ 
     488!----------------------------------------------------------------------- 
     489!       nam_traadv   advection scheme for tracer 
     490!----------------------------------------------------------------------- 
     491!  ln_traadv_cen2     2nd order centered scheme          (default T) 
     492!  ln_traadv_tvd      TVD scheme                         (default F) 
     493!  ln_traadv_muscl    MUSCL scheme                       (default F) 
     494!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F) 
     495!  ln_traadv_ubs      UBS scheme                         (default F) 
     496!  ln_traadv_qck      QUICKEST scheme                    (default F) 
     497&nam_traadv 
     498   ln_traadv_cen2   =  .false. 
     499   ln_traadv_tvd    =  .true. 
     500   ln_traadv_muscl  =  .false. 
     501   ln_traadv_muscl2 =  .false. 
     502   ln_traadv_ubs    =  .false. 
     503   ln_traadv_qck    =  .false. 
     504/ 
     505!----------------------------------------------------------------------- 
     506!       nam_traldf   lateral diffusion scheme for tracer 
     507!----------------------------------------------------------------------- 
     508!  Type of the operator :  
     509!     ln_traldf_lap    laplacian operator          (default T) 
     510!     ln_traldf_bilap  bilaplacian operator        (default F) 
     511!  Direction of action  :  
     512!     ln_traldf_level  iso-level                   (default F) 
     513!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T) 
     514!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
     515!  Coefficient 
     516!     aht0    horizontal eddy diffusivity for tracers (m2/s) 
     517!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
     518!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv") 
     519&nam_traldf 
     520   ln_traldf_lap    =  .true. 
     521   ln_traldf_bilap  =  .false. 
     522   ln_traldf_level  =  .false. 
     523   ln_traldf_hor    =  .false. 
     524   ln_traldf_iso    =  .true. 
     525   aht0    =  1000. 
     526   ahtb0   =     0. 
     527   aeiv0   =     0. 
     528/ 
     529!----------------------------------------------------------------------- 
     530!       namtdp   tracer newtonian damping ('key_tradmp') 
     531!----------------------------------------------------------------------- 
     532!  ndmp    type of damping in temperature and salinity  
     533!          !   ='latitude', damping poleward of 'ndmp' degrees and function  
     534!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1 
     535!          !   =-1          damping only in Med and Red Seas 
     536!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
     537!  nmldmp  type of damping in the mixed layer 
     538!          !   =0 damping throughout the water column 
     539!          !   =1 no damping in the mixed layer defined by avt >5cm2/s  
     540!          !   =2 no damping in the mixed layer defined rho<rho(surf)+.01  
     541!  sdmp    surface time scale for internal damping (days) 
     542!  bdmp    bottom time scale for internal damping (days) 
     543!  hdmp    depth of transition between sdmp and bdmp (meters) 
     544&namtdp 
     545   ndmp   =   -1 
     546   ndmpf  =    1 
     547   nmldmp =    1 
     548   sdmp   =  50. 
     549   bdmp   = 360. 
     550   hdmp   = 800. 
     551/ 
     552 
     553!!====================================================================== 
     554!!                      ***  Dynamics namelists  *** 
     555!!====================================================================== 
     556!!   nam_dynadv    formulation of the momentum advection 
     557!!   nam_dynvor    advection scheme 
     558!!   nam_dynhpg    hydrostatic pressure gradient 
     559!!   namflg        hydrostatic pressure gradient time stepping 
     560!!   nam_dynspg    surface pressure gradient                            (CPP key only) 
     561!!   nam_dynldf    lateral diffusion scheme 
     562!!====================================================================== 
     563 
     564!----------------------------------------------------------------------- 
     565!       nam_dynadv   formulation of the momentum advection 
     566!----------------------------------------------------------------------- 
     567!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T) 
     568!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
     569!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
     570&nam_dynadv    
     571   ln_dynadv_vec   = .true. 
     572   ln_dynadv_cen2  = .false. 
     573   ln_dynadv_ubs   = .false. 
     574 
     575!----------------------------------------------------------------------- 
     576!       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
     577!----------------------------------------------------------------------- 
     578!  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
     579!  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
     580!  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
     581!  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
     582&nam_dynvor 
     583   ln_dynvor_ene = .true. 
     584   ln_dynvor_ens = .false. 
     585   ln_dynvor_mix = .false. 
     586   ln_dynvor_een = .false. 
     587/ 
     588!----------------------------------------------------------------------- 
     589!       nam_dynhpg   Hydrostatic pressure gradient option 
     590!----------------------------------------------------------------------- 
     591!  type of pressure gradient scheme (choose one only!) 
     592!     ln_hpg_zco    z-coordinate - full steps                   (default T) 
     593!     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
     594!     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
     595!     ln_hpg_hel    s-coordinate (helsinki modification) 
     596!     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
     597!     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
     598!     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
     599!  parameters 
     600!     gamm          weighting coefficient (wdj scheme) 
     601&nam_dynhpg 
     602   ln_hpg_zco = .true. 
     603   ln_hpg_zps = .false. 
     604   ln_hpg_sco = .false. 
     605   ln_hpg_hel = .false. 
     606   ln_hpg_wdj = .false. 
     607   ln_hpg_djc = .false. 
     608   ln_hpg_rot = .false. 
     609   gamm       = 0.e0 
     610/ 
     611!----------------------------------------------------------------------- 
     612!       namflg   algorithm flags (algorithm not control by CPP keys) 
     613!----------------------------------------------------------------------- 
     614!  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
     615!                                                  centered     time scheme  (F) 
     616!   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
     617&namflg 
     618   ln_dynhpg_imp   =  .false. 
     619   nn_dynhpg_rst   =  0 
     620/ 
     621!----------------------------------------------------------------------- 
     622!       nam_dynspg   surface pressure gradient   (CPP key only) 
     623!----------------------------------------------------------------------- 
     624!                                !  explicit free surface               ("key_dynspg_exp") 
     625!                                !  filtered free surface               ("key_dynspg_flt") 
     626!                                !  split-explicit free surface         ("key_dynspg_ts") 
     627!                                !  rigid-lid                           ("key_dynspg_rl") 
     628 
     629!----------------------------------------------------------------------- 
     630!       nam_dynldf   lateral diffusion on momentum 
     631!----------------------------------------------------------------------- 
     632!  Type of the operator :  
     633!     ln_dynldf_lap    laplacian operator          (default T) 
     634!     ln_dynldf_bilap  bilaplacian operator        (default F) 
     635!  Direction of action  :  
     636!     ln_dynldf_level  iso-level                   (default F) 
     637!     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T) 
     638!     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
     639!  Coefficient 
     640!     ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
     641!     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
     642&nam_dynldf 
     643   ln_dynldf_lap    =  .true. 
     644   ln_dynldf_bilap  =  .false. 
     645   ln_dynldf_level  =  .false. 
     646   ln_dynldf_hor    =  .true. 
     647   ln_dynldf_iso    =  .false. 
     648   ahm0    = 100000. 
     649   ahmb0   =      0. 
     650/ 
     651 
     652!!====================================================================== 
     653!!             Tracers & Dynamics vertical physics namelists 
     654!!====================================================================== 
     655!!       namzdf        vertical physics 
     656!!       namnpc        non penetrative convection                        
     657!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      ) 
     658!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      ) 
     659!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      ) 
     660!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      ) 
     661!!====================================================================== 
     662 
     663!----------------------------------------------------------------------- 
     664!       namzdf   vertical physics 
     665!----------------------------------------------------------------------- 
     666!  ln_zdfevd  enhanced vertical diffusion         (default T) 
     667!  ln_zdfnpc  Non-Penetrative Convection          (default T) 
     668!  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst") 
     669!  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst") 
     670!  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
     671!  n_evdm     = 0  apply enhanced mixing on tracer only 
     672!             = 1  apply enhanced mixing on both tracer and momentum 
     673!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
     674!                               (=F)  euler backward (F) 
     675!  n_zdfexp   number of sub-timestep for time splitting scheme 
     676&namzdf 
     677   ln_zdfevd = .true. 
     678   ln_zdfnpc = .false. 
     679   avm0      = 1.2e-4 
     680   avt0      = 1.2e-5 
     681   avevd     =   100. 
     682   n_evdm    =     1 
     683   ln_zdfexp =  .false. 
     684   n_zdfexp  =     3 
     685/ 
     686!----------------------------------------------------------------------- 
     687!       namnpc   non penetrative convection 
     688!----------------------------------------------------------------------- 
     689!  nnpc1   non penetrative convective scheme frequency    
     690!  nnpc2   non penetrative convective scheme print frequency 
     691&namnpc 
     692   nnpc1  =      1 
     693   nnpc2  =    365 
     694/ 
     695!----------------------------------------------------------------------- 
     696!       namric   richardson number dependent vertical diffusion         ("key_zdfric" ) 
     697!----------------------------------------------------------------------- 
     698!  avmri   maximum value of the vertical viscosity 
     699!  alp     coefficient of the parameterization 
     700!  nric    coefficient of the parameterization 
     701&namric 
     702   avmri = 100.e-4 
     703   alp   =      5. 
     704   nric  =       2 
     705/ 
     706!----------------------------------------------------------------------- 
     707!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
     708!----------------------------------------------------------------------- 
     709!  ln_rstke flag to restart with tke from a run without tke (default F) 
     710!  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
     711!  ediss    coef. of the Kolmogoroff dissipation   
     712!  ebb      coef. of the surface input of tke 
     713!  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
     714!  emin     minimum value of tke (m^2/s^2) 
     715!  emin0    surface minimum value of tke (m^2/s^2) 
     716!  nitke    number of restart iterative loops 
     717!  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
     718!  nmxl     length used   = 0 bounded by the distance to surface and bottom 
     719!           !             = 1 bounded by the local vertical scale factor 
     720!           !             = 2 first vertical derivative of mixing length bounded by 1 
     721!  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
     722!           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
     723!           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
     724!  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
     725!  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
     726&namtke 
     727   ln_rstke = .false. 
     728   ediff =       0.1 
     729   ediss =       0.7 
     730   ebb   =      3.75 
     731   efave =        1. 
     732   emin  =     1.e-5 
     733   emin0 =     1.e-4 
     734   nitke =        50 
     735   nmxl  =         2 
     736   npdl  =         1 
     737   navb  =         0 
     738/ 
     739!----------------------------------------------------------------------- 
     740!       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" ) 
     741!----------------------------------------------------------------------- 
     742!  ln_kpprimix   shear instability mixing  (default T) 
     743!  difmiw        constant internal wave viscosity (m2/s) 
     744!  difsiw        constant internal wave diffusivity (m2/s) 
     745!  Riinfty       local Richardson Number limit for shear instability 
     746!  difri         maximum shear mixing at Rig = 0    (m2/s) 
     747!  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
     748!  difcon        maximum mixing in interior convection (m2/s) 
     749!  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
     750!  navb          = 0/1 flag for constant or profile background avt 
     751&namkpp 
     752   ln_kpprimix  = .true. 
     753   difmiw       =  1.e-04 
     754   difsiw       =  0.1e-04 
     755   Riinfty      =  0.8 
     756   difri        =  0.0050 
     757   bvsqcon      = -0.01e-07 
     758   difcon       =  1. 
     759   navb         =  0 
     760   nave         =  1 
     761/ 
     762!----------------------------------------------------------------------- 
     763!       namddm   double diffusive mixing parameterization               ("key_zdfddm") 
     764!----------------------------------------------------------------------- 
     765!   avts    maximum avs for dd mixing  
     766!   hsbfr   heat/salt buoyancy flux ratio 
     767&namddm 
     768      avts  = 1.e-4 
     769      hsbfr = 1.6 
     770/ 
     771 
     772!!====================================================================== 
     773!!                  ***  Miscelaneous namelists  *** 
     774!!====================================================================== 
     775!!   namsol      elliptic solver / island / free surface  
     776!!====================================================================== 
     777 
     778!----------------------------------------------------------------------- 
     779!       namsol   elliptic solver / island / free surface  
     780!----------------------------------------------------------------------- 
     781!  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
     782!                            (=2 successive-over-relaxation: sor) 
     783!                            (=3 FETI currently it does not work!       ("key_feti") 
     784!  nsol_arp  absolute/relative (0/1) precision convergence test 
     785!  nmin      minimum of iterations for the SOR solver 
     786!  nmax      maximum of iterations for the SOR solver 
     787!  nmod      frequency of test for the SOR solver 
     788!  eps       absolute precision of the solver 
     789!  resmax    absolute precision for the SOR solver 
     790!  sor       optimal coefficient for SOR solver 
     791!  epsisl    absolute precision on stream function solver 
     792!  nmisl     maximum pcg iterations for island 
     793!  rnu       strength of the additional force used in free surface b.c. 
     794&namsol 
     795   nsolv     =      2 
     796   nsol_arp  =      0 
     797   nmin      =    100 
     798   nmax      =    800 
     799   nmod      =     10 
     800   eps       =  1.E-6 
     801   resmax    = 1.E-10 
     802   sor       =   1.96 
     803   epsisl    = 1.e-10 
     804   nmisl     =   4000 
     805   rnu       =     1. 
     806/ 
     807!!====================================================================== 
     808!!                       ***  Diagnostics namelists  *** 
     809!!====================================================================== 
     810!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld") 
     811!!   namgap       level mean model-data gap                             ("key_diagap") 
     812!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr") 
     813!!   namflo       float parameters                                      ("key_float") 
     814!!   namptr       Poleward Transport Diagnostics 
     815!!====================================================================== 
     816 
     817!----------------------------------------------------------------------- 
     818!       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra") 
     819!                          or mixed-layer trends                        ('key_trdmld') 
     820!                          or barotropic vorticity                      ("key_trdvor") 
     821!----------------------------------------------------------------------- 
     822!  ntrd                time step frequency dynamics and tracers trends 
     823!  nctls               control surface type in mixed-layer trends (0,1 or n<jpk) 
     824!  ln_trdmld_restart   restart for ML diagnostics 
     825!  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day) 
     826!  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S 
     827&namtrd 
     828   ntrd              = 365 
     829   nctls             =   0 
     830   ln_trdmld_restart = .false. 
     831   ucf               =  1. 
     832   ln_trdmld_instant = .false. 
     833/ 
     834!----------------------------------------------------------------------- 
     835!       namgap    level mean model-data gap ('key_diagap') 
     836!----------------------------------------------------------------------- 
     837!  ngap    time-step frequency of model-data gap computation 
     838!  nprg    time-step frequency of gap print in model output 
     839&namgap 
     840   ngap =  15 
     841   nprg =  10 
     842/ 
     843!----------------------------------------------------------------------- 
     844!       namspr  surface pressure diagnostic 
     845!----------------------------------------------------------------------- 
     846!  nmaxp   maximum of iterations for the solver 
     847!  epsp    absolute precision of the solver 
     848!  niterp  number of iteration done by the solver 
     849&namspr 
     850   nmaxp   =   1000 
     851   epsp    =  1.e-3 
     852   niterp  =    400 
     853/ 
     854!----------------------------------------------------------------------- 
     855!       namflo    float parameters                                      ("key_float") 
    643856!----------------------------------------------------------------------- 
    644857!  ln_rstflo   boolean term for float restart (true or false) 
     
    655868    ln_flork4 = .false. 
    656869/ 
    657 !!====================================================================== 
    658 !!              Surface Boundary Condition namelists 
    659 !! 
    660 !!       namsbc    surface boundary condition 
    661 !!       namsbc_ana    analytical         formulation of sbc 
    662 !!       namsbc_flx    flux               formulation of sbc 
    663 !!       namsbc_core   CORE bulk formulea formulation of sbc 
    664 !!       namsbc_clio   CLIO bulk formulea formulation of sbc 
    665 !!       namsbc_cpl    CouPLed            formulation of sbc 
    666 !!       namsbc_ssr    sea surface restoring term (for T and/or S) 
    667 !!====================================================================== 
    668 !----------------------------------------------------------------------- 
    669 !       namsbc       Surface Boundary Condition (surface module) 
    670 !----------------------------------------------------------------------- 
    671 !  nn_fsbc          frequency of surface boundary condition computation  
    672 !                   (also correspond to the frequency of sea-ice model call) 
    673 ! 
    674 !  ln_ana        analytical formulation (fill namsbc_ana )  
    675 !  ln_flx        flux formulation       (fill namsbc_flx ) (overwritten by key_sbc_flux     ) 
    676 !  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core) (overwritten by key_sbc_bulk_clio) 
    677 !  ln_blk_core   CORE bulk formulation  (fill namsbc_clio) (overwritten by key_sbc_bulk_core) 
    678 !  ln_cpl        Coupled formulation    (fill namsbc_cpl ) (overwritten by key_sbc_coupled  ) 
    679 ! 
    680 !  nn_ice        =0 no ice boundary condition   , 
    681 !                =1 observed ice-cover          , 
    682 !                =2 ice-model ice-cover (overwritten by key_ice_lim) 
    683 ! 
    684 !  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
    685 ! 
    686 !  ln_ssr       Sea Surface restoring on SST and/or SSS (fill namsbc_ssr) 
    687 !               (overwritten by key_dta_sst and/or key_dtasss) 
    688 ! 
    689 !  nn_fwb       FreshWater Budget: =0 no control                              ,  
    690 !                                  =1 annual global mean of emp set to zero   , 
    691 !                                  =2 global mean of emp set to zero at each nn_fsbc time step 
    692 &namsbc 
    693    nn_fsbc     = 0 
    694    ln_ana      = .true. 
    695    ln_flx      = .false. 
    696    ln_blk_clio = .false. 
    697    ln_blk_core = .false. 
    698    ln_cpl      = .false. 
    699    nn_ice      = 0 
    700    ln_dm2dc    = .false. 
    701    ln_rnf      = .false. 
    702    ln_ssr      = .false. 
    703    nn_fwb      = 0 
    704 / 
    705 !----------------------------------------------------------------------- 
    706 !       namsbc_rnf   runoffs namelist surface boundary condition 
    707 !----------------------------------------------------------------------- 
    708 !nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
    709 !rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
    710 !rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
    711 !cn_dir      Root directory for location of ssr files 
    712 !sn_rnf      informations about the runoff file to be read 
    713 !sn_cnf      informations about the runoff mouth file to be read 
    714 !  THE ORDER OF THE FILES MATTER: 
    715 !               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
    716 !               !                    !   (if <0  months) !   name     !     (logical)      ! 
    717 &namsbc_rnf 
    718    nn_runoff  = 0 
    719    rn_hrnf    = 0.e0 
    720    rn_avt_rnf = 1.e-3 
    721    cn_dir     = './' 
    722    sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .TRUE.   ,  1  ,  0      
    723    sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .FALSE.  ,  1  ,  0               
    724 / 
    725 !----------------------------------------------------------------------- 
    726 !       namsbc_ana   analytical surface boundary condition 
    727 !----------------------------------------------------------------------- 
    728 !  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
    729 !  rn_utau0    uniform value used as default wind stress -u 
    730 !  rn_vtau0    uniform value used as default wind stress -v 
    731 !  rn_q0       uniform value used as default surface heat flux 
    732 !  rn_qsr0     uniform value used as default solar radiation flux 
    733 !  rn_emp0     uniform value used as default surface freswater budget (E-P) 
    734 &namsbc_ana 
    735    nn_tau000  =   100 
    736    rn_utau0   =   0.1e0 
    737    rn_vtau0   =   0.e0 
    738    rn_q0      =   0.e0 
    739    rn_qsr0    =   0.e0 
    740    rn_emp0    =   0.e0 
    741 / 
    742 !----------------------------------------------------------------------- 
    743 !       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
    744 !----------------------------------------------------------------------- 
    745 !  cn_dir   directory in which the model is executed 
    746 ! 
    747 !  THE ORDER OF THE FILES MATTER: 
    748 !               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    749 !               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    750 &namsbc_flx 
    751    cn_dir       = './' 
    752    sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .FALSE.       ,    0   ,     0 
    753    sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .FALSE.       ,    0   ,     0 
    754    sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .FALSE.       ,    0   ,     0 
    755    sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .FALSE.       ,    0   ,     0 
    756    sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .FALSE.       ,    0   ,     0 
    757  
    758 /       
    759 !----------------------------------------------------------------------- 
    760 !       namsbc_clio  CLIO bulk formulea 
    761 !----------------------------------------------------------------------- 
    762 !  cn_dir   directory in which the model is executed 
    763 ! 
    764 !  THE ORDER OF THE FILES MATTER: 
    765 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    766 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    767 &namsbc_clio 
    768    cn_dir       = './' 
    769    sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .FALSE.        ,    1   ,    0     
    770    sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .FALSE.        ,    1   ,    0     
    771    sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .FALSE.        ,    1   ,    0       
    772    sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .FALSE.        ,    1   ,    0      
    773    sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .FALSE.        ,    1   ,    0     
    774    sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .FALSE.        ,    1   ,    0       
    775    sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .FALSE.        ,    1   ,    0        
    776 / 
    777 !----------------------------------------------------------------------- 
    778 !       namsbc_core  CORE bulk formulea 
    779 !----------------------------------------------------------------------- 
    780 !  cn_dir         directory in which the model is executed 
    781 !  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
    782 !  alpha_precip   multiplication factor for precipitation (total & snow) 
    783 ! 
    784 !  THE ORDER OF THE FILES MATTER: 
    785 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    786 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    787 &namsbc_core 
    788    cn_dir       = './' 
    789    ln_2m        = .FALSE. 
    790    alpha_precip = 1. 
    791    sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .FALSE.        ,    1   ,    0     
    792    sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .FALSE.        ,    1   ,    0     
    793    sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .FALSE.        ,    1   ,    0       
    794    sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .FALSE.        ,    1   ,    0       
    795    sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .FALSE.        ,    1   ,    0      
    796    sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .FALSE.        ,    1   ,    0     
    797    sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .FALSE.        ,    1   ,    0       
    798    sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .FALSE.        ,    1   ,    0        
    799 / 
    800 !----------------------------------------------------------------------- 
    801 !       namsbc_cpl    coupled ocean/atmosphere model  (#ifdef "key_coupled") 
    802 !----------------------------------------------------------------------- 
    803 &namsbc_cpl 
    804 / 
    805 !----------------------------------------------------------------------- 
    806 !       namsbc_ssr   surface boundary condition : sea surface restoring 
    807 !----------------------------------------------------------------------- 
    808 !  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
    809 !  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
    810 !             (=1, salt flux, concentration/dillution emps only) 
    811 !             (=2, volume flux, both emp and emps are updated)  
    812 !  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
    813 !  deds       magnitude of the damping on salinity [mm/day] 
    814 !  THE ORDER OF THE FILES MATTER: 
    815 !               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    816 !               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    817 &namsbc_ssr 
    818    nn_sstr      = 0 
    819    nn_sssr      = 0 
    820    dqdt         = -40. 
    821    deds         =  27.7 
    822    sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .FALSE.       ,    0   ,     0 
    823    sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .TRUE.        ,    0   ,     0 
    824 /       
     870!----------------------------------------------------------------------- 
     871!       namptr   Poleward Transport Diagnostic 
     872!----------------------------------------------------------------------- 
     873!  ln_diaptr  logical flag for Poleward transport computation 
     874!  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
     875!             need input basins mask file named "subbasins.nc" 
     876!  nf_ptr     Frequency of computation 
     877&namptr 
     878   ln_diaptr = .false. 
     879   ln_subbas = .false. 
     880   nf_ptr    =  15 
     881/ 
  • trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist

    r889 r907  
    1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! OPA namelist :  model option and parameter input 
    3 ! ------------- 
    4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    5 ! 
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/OPA  :       1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 
     3!! namelists         2 - Domain           (nam_zgr, nam_zgr_sco, namdom) 
     4!!                   3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
     5!!                                         namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 
     6!!                   4 - lateral boundary (namlbc, namcla, namobc) 
     7!!                   5 - bottom  boundary (nambfr, nambbc, nambbl) 
     8!!                   6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp) 
     9!!                   7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) 
     10!!                   8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm) 
     11!!                   9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr) 
     12!!                   9 - miscellaneous    (namsol) 
     13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     14!  CAUTION: some scripts does not support CAPITALs for logical definition  
     15!  *******       use .true. or .false. and NOT .TRUE. or .FALSE. 
     16 
     17!!====================================================================== 
     18!!                   ***  Run management namelists  *** 
     19!!====================================================================== 
     20!!   namrun            parameters of the run 
     21!!   nam_mpp           Massively Parallel Processing 
     22!!   nam_mpp_dyndist    
     23!!   nam_ctl           Control prints & Benchmark 
     24!!====================================================================== 
     25 
    626!----------------------------------------------------------------------- 
    727!       namrun   parameters of the run 
     
    1030!  cexper     experience name for vairmer format 
    1131!  ln_rstart  boolean term for restart (true or false) 
    12 !  nrstdt     control of the restart timestep: 
    13 !                = 0 restart, do not control nit000 in the restart file. 
    14 !                = 1 restart, control nit000 in the restart file. Do not 
    15 !                    use the date in the restart file (use ndate0 in namelist) 
    16 !                = 2 restart, control nit000 in the restart file, use the date 
    17 !                    in the restart file. ndate0 in the namelist is ignored. 
     32!  nrstdt     restart control = 0 restart, do not control nit000 in the restart file. 
     33!             !               = 1 restart, control nit000 in the restart file. Do not 
     34!             !                   use the date in the restart file (use ndate0 in namelist) 
     35!             !               = 2 restart, control nit000 in the restart file, use the date 
     36!             !                   in the restart file. ndate0 in the namelist is ignored. 
    1837!  nit000     number of the first time step 
    1938!  nitend     number of the last time step 
     
    2443!  nwrite     frequency of OUTPUT file 
    2544!  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 
    26 ! 
    27 !  CAUTION: for usual run scripts, logical value of 
    28 !  *******  ln_rstart must be .true. or .false. 
    29 !                     and NOT .TRUE. or .FALSE. 
    3045&namrun 
    3146   no         =       0 
     
    3853   nleapy     =       0 
    3954   ninist     =       0 
    40    nstock     =   10950 
    41    nwrite     =   10950 
    42    ln_dimgnnn = .FALSE. 
     55   nstock     =   10950  
     56   nwrite     =   10950  
     57   ln_dimgnnn = .false. 
     58/ 
     59!----------------------------------------------------------------------- 
     60!       nam_mpp      Massively Parallel Processing 
     61!----------------------------------------------------------------------- 
     62!  c_mpi_send         mpi send/recieve type 
     63!                      = 'S'  : standard blocking send 
     64!                      = 'B'  : buffer blocking send 
     65!                      = 'I'  : immediate non-blocking send 
     66&nam_mpp 
     67   c_mpi_send =  'S' 
     68/ 
     69!----------------------------------------------------------------------- 
     70!       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist") 
     71!----------------------------------------------------------------------- 
     72!  jpni   number of processors following i 
     73!  jpnj   number of processors following j 
     74!  jpnij  number of local domains 
     75&nam_mpp_dyndist 
     76   jpni  = 1 
     77   jpnj  = 1 
     78   jpnij = 1 
    4379/ 
    4480!----------------------------------------------------------------------- 
     
    69105   nbit_cmp =     0 
    70106/ 
    71 !----------------------------------------------------------------------- 
    72 !       nam_mpp      Massively Parallel Processing 
    73 !----------------------------------------------------------------------- 
    74 !  c_mpi_send         mpi send/recieve type 
    75 !                      = 'S'  : standard blocking send 
    76 !                      = 'B'  : buffer blocking send 
    77 !                      = 'I'  : immediate non-blocking send 
    78 &nam_mpp 
    79    c_mpi_send =  'S' 
    80 / 
    81 !----------------------------------------------------------------------- 
    82 !       nam_mpp_dyndist      Massively Parallel Distribution 
    83 !              #ifdef ( key_agrif && key_mpp_dyndist ) 
    84 !----------------------------------------------------------------------- 
    85 !  jpni   number of processors following i 
    86 !  jpnj   number of processors following j 
    87 !  jpnij  number of local domains 
    88 &nam_mpp_dyndist 
    89    jpni  = 1 
    90    jpnj  = 1 
    91    jpnij = 1 
    92 / 
     107 
     108!!====================================================================== 
     109!!                      ***  Domain namelists  *** 
     110!!====================================================================== 
     111!!   nam_zgr       vertical coordinate 
     112!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate 
     113!!   namdom        space and time domain (bathymetry, mesh, timestep) 
     114!!====================================================================== 
     115 
    93116!----------------------------------------------------------------------- 
    94117!       nam_zgr       vertical coordinate 
    95118!----------------------------------------------------------------------- 
    96 !  ln_zco     z-coordinate - full steps      (T/F) 
     119!  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined) 
    97120!  ln_zps     z-coordinate - partial steps   (T/F) 
    98121!  ln_sco     s- or hybrid z-s-coordinate    (T/F) 
     
    118141/ 
    119142!----------------------------------------------------------------------- 
    120 !       nam_traadv   advection scheme for tracer (option not control by CPP keys) 
    121 !----------------------------------------------------------------------- 
    122 !  ln_traadv_cen2     2nd order centered scheme    (default T) 
    123 !  ln_traadv_tvd      TVD scheme                   (default F) 
    124 !  ln_traadv_muscl    MUSCL scheme                 (default F) 
    125 !  ln_traadv_muscl2   MUSCL2 scheme                (default F) 
    126 !  ln_traadv_ubs      UBS scheme                   (default F) 
    127 !  ln_traadv_qck      QUICKEST scheme              (default F) 
    128 &nam_traadv 
    129    ln_traadv_cen2   =  .false. 
    130    ln_traadv_tvd    =  .true. 
    131    ln_traadv_muscl  =  .false. 
    132    ln_traadv_muscl2 =  .false. 
    133    ln_traadv_ubs    =  .false. 
    134    ln_traadv_qck    =  .false. 
    135 / 
    136 !----------------------------------------------------------------------- 
    137 !       nam_traldf   lateral diffusion scheme for tracer (option not control by CPP keys) 
    138 !----------------------------------------------------------------------- 
    139 !  Type of the operator :  
    140 !     ln_traldf_lap    laplacian operator          (default T) 
    141 !     ln_traldf_bilap  bilaplacian operator        (default F) 
    142 !  Direction of action  :  
    143 !     ln_traldf_level  iso-level                   (default F) 
    144 !     ln_traldf_hor    horizontal (geopotential)   (default F)^** 
    145 !     ln_traldf_iso    iso-neutral                 (default T)^* 
    146 !  Coefficient 
    147 !     aht0    horizontal eddy diffusivity for tracers (m2/s) 
    148 !     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
    149 !     aeiv0   eddy induced velocity coefficient (m2/s) 
    150 ! ^* require key_ldfslp to compute the direction of the lateral diffusion 
    151 ! ^** require key_ldfslp in s-coordinate 
    152 &nam_traldf 
    153    ln_traldf_lap    =  .true. 
    154    ln_traldf_bilap  =  .false. 
    155    ln_traldf_level  =  .false. 
    156    ln_traldf_hor    =  .false. 
    157    ln_traldf_iso    =  .true. 
    158    aht0    =  1000. 
    159    ahtb0   =     0. 
    160    aeiv0   =     0. 
    161 / 
    162 !----------------------------------------------------------------------- 
    163 !       nam_dynldf   lateral diffusion on momentum 
    164 !----------------------------------------------------------------------- 
    165 !  Type of the operator :  
    166 !     ln_dynldf_lap    laplacian operator          (default T) 
    167 !     ln_dynldf_bilap  bilaplacian operator        (default F) 
    168 !  Direction of action  :  
    169 !     ln_dynldf_level  iso-level                   (default F) 
    170 !     ln_dynldf_hor    horizontal (geopotential)   (default F)^** 
    171 !     ln_dynldf_iso    iso-neutral                 (default T)^* 
    172 !  Coefficient 
    173 !  ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
    174 !  ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
    175 &nam_dynldf 
    176    ln_dynldf_lap    =  .false. 
    177    ln_dynldf_bilap  =  .true. 
    178    ln_dynldf_level  =  .false. 
    179    ln_dynldf_hor    =  .true. 
    180    ln_dynldf_iso    =  .false. 
    181    ahm0    = -8.5e+11 
    182    ahmb0   =     0. 
    183 / 
    184 !----------------------------------------------------------------------- 
    185 !       namflg   algorithm flags (algorithm not control by CPP keys) 
    186 !----------------------------------------------------------------------- 
    187 !  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
    188 !                                                  centered      time scheme  (F) 
    189 !   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
    190 &namflg 
    191    ln_dynhpg_imp   =  .false. 
    192    nn_dynhpg_rst   =  0 
    193 / 
    194 !----------------------------------------------------------------------- 
    195 !       nam_dynhpg   Hydrostatic pressure gradient option 
    196 !----------------------------------------------------------------------- 
    197 !  type of pressure gradient scheme (choose one only!) 
    198 !     ln_hpg_zco    z-coordinate - full steps                   (default T) 
    199 !     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
    200 !     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
    201 !     ln_hpg_hel    s-coordinate (helsinki modification) 
    202 !     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
    203 !     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
    204 !     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
    205 !  parameters 
    206 !    gamm          weighting coefficient (wdj scheme) 
    207 &nam_dynhpg 
    208    ln_hpg_zco = .false. 
    209    ln_hpg_zps = .true. 
    210    ln_hpg_sco = .false. 
    211    ln_hpg_hel = .false. 
    212    ln_hpg_wdj = .false. 
    213    ln_hpg_djc = .false. 
    214    ln_hpg_rot = .false. 
    215    gamm       = 0.e0 
    216 / 
    217 !----------------------------------------------------------------------- 
    218 !       nam_dynadv   option of physics/algorithm (not control by CPP keys) 
    219 !----------------------------------------------------------------------- 
    220 !  ln_dynadv_vec      vector form flag 
    221 !  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
    222 !  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
    223 &nam_dynadv    
    224    ln_dynadv_vec   = .TRUE. 
    225    ln_dynadv_cen2  = .FALSE. 
    226    ln_dynadv_ubs   = .FALSE. 
    227  
    228 !----------------------------------------------------------------------- 
    229 !       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
    230 !----------------------------------------------------------------------- 
    231 !  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
    232 !  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
    233 !  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
    234 !  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
    235 &nam_dynvor 
    236    ln_dynvor_ene = .FALSE. 
    237    ln_dynvor_ens = .FALSE. 
    238    ln_dynvor_mix = .FALSE. 
    239    ln_dynvor_een = .TRUE. 
    240 / 
    241 !----------------------------------------------------------------------- 
    242 !       namalb   albedo parameters 
    243 !----------------------------------------------------------------------- 
    244 !  cgren    correction of the snow or ice albedo to take into account 
    245 !  albice   albedo of melting ice in the arctic and antarctic 
    246 !  alphd    coefficients for linear interpolation used to compute albedo 
    247 !           between two extremes values (Pyane, 1972) 
    248 !  alphc     "                                         " 
    249 !  alphdi    "                                         " 
    250 &namalb 
    251    cgren    =      0.06 
    252    albice   =      0.5 
    253    alphd    =      0.80 
    254    alphc    =      0.65 
    255    alphdi   =      0.72 
    256 / 
    257 !----------------------------------------------------------------------- 
    258143!       namdom   space and time domain (bathymetry, mesh, timestep) 
    259144!----------------------------------------------------------------------- 
    260 !  ntopo      = 0/1 ,compute/read the bathymetry file 
    261 !               (mbathy, nb of T-ocean levels) 
     145!  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 
    262146!  e3zps_min  the thickness of the partial step is set larger than the 
    263 !  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t 
    264 !                (N.B. 0<e3zps_rat<1) 
    265 !                  (coordinates, scale factors) 
     147!  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 
    266148!  nmsh       =1 create a mesh file (coordinates, scale factors, masks) 
    267 !  nacc       the acceleration of convergence method 
    268 !             = 0, no acceleration, rdt = rdttra 
    269 !             = 1, acceleration used, rdt < rdttra(k) 
     149!  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra 
     150!                                                  = 1,     used, rdt < rdttra(k) 
    270151!  atfp       asselin time filter parameter 
    271152!  rdt        time step for the dynamics (and tracer if nacc=0) 
     
    274155!  rdth       depth variation of tracer time step 
    275156!  rdtbt      barotropic time step (for the time splitting algorithm) 
    276 !  nclosea    0 no closed sea  
    277 !             1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
     157!  nclosea    =0 remove the closed sea from the global domain (orca configuration) 
     158!             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
    278159&namdom 
    279160   ntopo     =     1 
     
    290171   nclosea   =     0 
    291172/ 
    292 !----------------------------------------------------------------------- 
    293 !       namfwb   freshwater budget correction 
    294 !----------------------------------------------------------------------- 
    295 !  ln_fwb     logical flag for freshwater budget correction (0 annual mean) 
    296 &namfwb 
    297    ln_fwb    = .false. 
    298 / 
    299 !----------------------------------------------------------------------- 
    300 !       namptr   Poleward Transport Diagnostic 
    301 !----------------------------------------------------------------------- 
    302 !  ln_diaptr  logical flag for Poleward transport computation 
    303 !  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
    304 !             need input basins mask file named "subbasins.nc" 
    305 !  nf_ptr     Frequency of computation 
    306 &namptr 
    307    ln_diaptr = .false. 
    308    ln_subbas = .false. 
    309    nf_ptr    =  15 
    310 / 
    311 !----------------------------------------------------------------------- 
    312 !       namcro   cross land advection 
    313 !----------------------------------------------------------------------- 
    314 !  n_cla   advection between 2 ocean pts separates by land  
    315 &namcla 
    316    n_cla   = 0 
    317 / 
    318 !----------------------------------------------------------------------- 
    319 !       namzdf   vertical physics 
    320 !----------------------------------------------------------------------- 
    321 !  ln_zdfevd  enhanced vertical diffusion         (default T) 
    322 !  ln_zdfnpc  Non-Penetrative Convection          (default T) 
    323 !  avm0       vertical eddy viscosity for the dynamic (m2/s) 
    324 !  avt0       vertical eddy diffusivity for tracers (m2/s) 
    325 !  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
    326 !  nevdm      = 0  apply enhanced mixing on tracer only 
    327 !             = 1  apply enhanced mixing on both tracer and momentum 
    328 !  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
    329 !                               (=F)  euler backward (F) 
    330 !  n_zdfexp   number of sub-timestep for time splitting scheme 
    331 &namzdf 
    332    ln_zdfevd = .true. 
    333    ln_zdfnpc = .false. 
    334    avm0      = 1.2e-4 
    335    avt0      = 1.2e-5 
    336    avevd     =   100. 
    337    n_evdm    =     1 
    338    ln_zdfexp =  .false. 
    339    n_zdfexp  =     3 
    340 / 
    341 !----------------------------------------------------------------------- 
    342 !       namnpc   vnon penetrative convection 
    343 !----------------------------------------------------------------------- 
    344 !  nnpc1   non penetrative convective scheme frequency 
    345 !  nnpc2   non penetrative convective scheme print frequency 
    346 &namnpc 
    347    nnpc1  =      1 
    348    nnpc2  =    365 
    349 / 
    350 !----------------------------------------------------------------------- 
    351 !       nambbl   bottom boundary layer scheme 
    352 !----------------------------------------------------------------------- 
    353 !  atrbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s) 
    354 &nambbl 
    355    atrbbl = 10000. 
    356 / 
    357 !----------------------------------------------------------------------- 
    358 !       namric   richardson number dependent vertical diffusion 
    359 !                ( #ifdef "key_zdfrichardson" ) 
    360 !----------------------------------------------------------------------- 
    361 !  avmri   maximum value of the vertical viscosity 
    362 !  alp     coefficient of the parameterization 
    363 !  nric    coefficient of the parameterization 
    364 &namric 
    365    avmri = 100.e-4 
    366    alp   =      5. 
    367    nric  =       2 
    368 / 
    369 !----------------------------------------------------------------------- 
    370 !       namtke   turbulent eddy kinetic dependent vertical diffusion 
    371 !                ( #ifdef "key_zdftke" ) 
    372 !----------------------------------------------------------------------- 
    373 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    374 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    375 !  ediss    coef. of the Kolmogoroff dissipation   
    376 !  ebb      coef. of the surface input of tke 
    377 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    378 !  emin     minimum value of tke (m^2/s^2) 
    379 !  emin0    surface minimum value of tke (m^2/s^2) 
    380 !  nitke    number of restart iterative loops 
    381 !  ri_c     critic richardson number 
    382 !  nmxl     flag on mixing length used 
    383 !           = 0 bounded by the distance to surface and bottom 
    384 !           = 1 bounded by the local vertical scale factor 
    385 !           = 2 first vertical derivative of mixing length bounded by 1 
    386 !  npdl     flag on prandtl number 
    387 !           = 0 no vertical prandtl number (avt=avm) 
    388 !           = 1 prandtl number function of richarson number (avt=pdl*avm) 
    389 !           = 2 same as = 1 but a shapiro filter is applied on pdl 
    390 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    391 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    392 &namtke 
    393    ln_rstke = .false. 
    394    ediff =       0.1 
    395    ediss =       0.7 
    396    ebb   =       60.  
    397    efave =        1. 
    398    emin  =     1.e-6 
    399    emin0 =     1.e-4 
    400    nitke =        50 
    401    nmxl  =         2 
    402    npdl  =         1 
    403    navb  =         0 
    404 / 
    405 !----------------------------------------------------------------------- 
    406 !       namkpp   K-Profile Parameterization dependent vertical diffusion 
    407 !                ( #ifdef "key_zdfkpp" ) 
    408 !----------------------------------------------------------------------- 
    409 !  ln_kpprimix   shear instability mixing  (default T) 
    410 !  difmiw        constant internal wave viscosity (m2/s) 
    411 !  difsiw        constant internal wave diffusivity (m2/s) 
    412 !  Riinfty       local Richardson Number limit for shear instability 
    413 !  difri         maximum shear mixing at Rig = 0    (m2/s) 
    414 !  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
    415 !  difcon        maximum mixing in interior convection (m2/s) 
    416 !  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
    417 !  navb          = 0/1 flag for constant or profile background avt 
    418 &namkpp 
    419    ln_kpprimix  = .true. 
    420    difmiw       =  1.e-04 
    421    difsiw       =  0.1e-04 
    422    Riinfty      =  0.8 
    423    difri        =  0.0050 
    424    bvsqcon      = -0.01e-07 
    425    difcon       =  1. 
    426    navb         =  0 
    427    nave         =  1 
    428 / 
    429 !----------------------------------------------------------------------- 
    430 !       namddm   double diffusive mixing parameterization 
    431 !----------------------------------------------------------------------- 
    432 !   avts    maximum avs for dd mixing  
    433 !   hsbfr   heat/salt buoyancy flux ratio 
    434 &namddm 
    435       avts  = 1.e-4 
    436       hsbfr = 1.6 
    437 / 
    438 !----------------------------------------------------------------------- 
    439 !       namlbc   lateral momentum boundary condition 
    440 !----------------------------------------------------------------------- 
    441 !  shlat   lateral boundary condition on velocity 
    442 !                   shlat = 0 , free slip 
    443 !               0 < shlat < 2 , partial slip 
    444 !                   shlat = 2 , no slip 
    445 !               2 < shlat     , strong slip 
    446 &namlbc 
    447    shlat  =      2. 
    448 / 
    449 !----------------------------------------------------------------------- 
    450 !       nambfr   bottom friction 
    451 !----------------------------------------------------------------------- 
    452 !  nbotfr  type of bottom friction  
    453 !                  nbotfr = 0 , no slip  
    454 !                  nbotfr = 1 , linear friction 
    455 !                  nbotfr = 2 , nonlinear friction 
    456 !                  nbotfr = 3 , free slip 
    457 !  bfri1   bottom drag coefficient (linear case) 
    458 !  bfri2   bottom drag coefficient (non linear case) 
    459 !  bfeb2   bottom turbulent kinetic energy  (m^2/s^2) 
    460 &nambfr 
    461    nbotfr =       1 
    462    bfri1  =   4.e-4 
    463    bfri2  =   1.e-3 
    464    bfeb2  =  2.5e-3 
    465 / 
    466 !----------------------------------------------------------------------- 
    467 !       nambbc   bottom temperature boundary condition 
    468 !----------------------------------------------------------------------- 
    469 !  ngeo_flux  = 0 no geothermal heat flux 
    470 !             = 1 constant geothermal heat flux 
    471 !             = 2 variable geothermal heat flux (read in geothermal_heating.nc)  
    472 !                 ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 
    473 !  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
    474 &nambbc 
    475    ngeo_flux =  2 
    476    ngeo_flux_const = 86.4e-3 
     173 
     174!!====================================================================== 
     175!!            ***  Surface Boundary Condition namelists  *** 
     176!!====================================================================== 
     177!!   namsbc        surface boundary condition 
     178!!   namsbc_ana    analytical         formulation 
     179!!   namsbc_flx    flux               formulation 
     180!!   namsbc_clio   CLIO bulk formulea formulation 
     181!!   namsbc_core   CORE bulk formulea formulation 
     182!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled") 
     183!!   namqsr        penetrative solar radiation 
     184!!   namsbc_rnf    river runoffs 
     185!!   namsbc_ssr    sea surface restoring term (for T and/or S) 
     186!!   namalb        albedo parameters 
     187!!====================================================================== 
     188 
     189!----------------------------------------------------------------------- 
     190!       namsbc       Surface Boundary Condition (surface module) 
     191!----------------------------------------------------------------------- 
     192!  nn_fsbc       frequency of surface boundary condition computation  
     193!                (= the frequency of sea-ice model call) 
     194!  ln_ana        analytical formulation (fill namsbc_ana )  
     195!  ln_flx        flux formulation       (fill namsbc_flx )  
     196!  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core)  
     197!  ln_blk_core   CORE bulk formulation  (fill namsbc_clio)  
     198!  ln_cpl        Coupled formulation    (fill namsbc_cpl )  
     199!  nn_ice        =0 no ice boundary condition   , 
     200!                =1 observed ice-cover          , 
     201!                =2 ice-model used                                      ("key_lim3" or "key_lim2) 
     202!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
     203!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr) 
     204!  nn_fwb       FreshWater Budget: =0 no control                              ,  
     205!                                  =1 annual global mean of e-p-r set to zero   , 
     206!                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step 
     207&namsbc 
     208   nn_fsbc     = 5 
     209   ln_ana      = .false. 
     210   ln_flx      = .false. 
     211   ln_blk_clio = .true. 
     212   ln_blk_core = .false. 
     213   ln_cpl      = .false. 
     214   nn_ice      = 2 
     215   ln_dm2dc    = .false. 
     216   ln_rnf      = .false. 
     217   ln_ssr      = .false. 
     218   nn_fwb      = 0 
     219/ 
     220!----------------------------------------------------------------------- 
     221!       namsbc_ana   analytical surface boundary condition 
     222!----------------------------------------------------------------------- 
     223!  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
     224!  rn_utau0    uniform value for the i-stress 
     225!  rn_vtau0    uniform value for the j-stress 
     226!  rn_q0       uniform value for the total heat flux 
     227!  rn_qsr0     uniform value for the solar radiation 
     228!  rn_emp0     uniform value for the freswater budget (E-P) 
     229&namsbc_ana 
     230   nn_tau000  =   0 
     231   rn_utau0   =   0.5 
     232   rn_vtau0   =   0.e0 
     233   rn_q0      =   0.e0 
     234   rn_qsr0    =   0.e0 
     235   rn_emp0    =   0.e0 
     236/ 
     237!----------------------------------------------------------------------- 
     238!       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
     239!----------------------------------------------------------------------- 
     240!  cn_dir   directory in which the model is executed 
     241! 
     242!  THE ORDER OF THE FILES MATTER: 
     243!               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     244!               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     245&namsbc_flx 
     246   cn_dir       = './' 
     247   sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0 
     248   sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0 
     249   sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0 
     250   sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0 
     251   sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0 
     252/       
     253!----------------------------------------------------------------------- 
     254!       namsbc_clio  CLIO bulk formulea 
     255!----------------------------------------------------------------------- 
     256!  cn_dir   directory in which the model is executed 
     257! 
     258!  THE ORDER OF THE FILES MATTER: 
     259!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     260!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     261&namsbc_clio 
     262   cn_dir       = './' 
     263   sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0     
     264   sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0     
     265   sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0       
     266   sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0      
     267   sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0     
     268   sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0       
     269   sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0        
     270/ 
     271!----------------------------------------------------------------------- 
     272!       namsbc_core  CORE bulk formulea 
     273!----------------------------------------------------------------------- 
     274!  cn_dir         directory in which the model is executed 
     275!  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
     276!  alpha_precip   multiplication factor for precipitation (total & snow) 
     277! 
     278!  THE ORDER OF THE FILES MATTER: 
     279!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     280!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     281&namsbc_core 
     282   cn_dir       = './' 
     283   ln_2m        = .false. 
     284   alpha_precip = 1. 
     285   sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0     
     286   sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0     
     287   sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0       
     288   sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0       
     289   sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0      
     290   sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0     
     291   sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0       
     292   sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0        
     293/ 
     294!----------------------------------------------------------------------- 
     295!       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled") 
     296!----------------------------------------------------------------------- 
     297&namsbc_cpl 
    477298/ 
    478299!----------------------------------------------------------------------- 
     
    490311/ 
    491312!----------------------------------------------------------------------- 
    492 !       namtdp   tracer newtonian damping ('key_tradmp') 
    493 !----------------------------------------------------------------------- 
    494 !  ndmp    type of damping in temperature and salinity  
    495 !          (='latitude', damping poleward of 'ndmp' degrees and function  
    496 !             of the distance-to-coast. Red and Med Seas as ndmp=-1) 
    497 !          (=-1 damping only in Med and Red Seas) 
    498 !  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
    499 !  nmldmp  type of damping in the mixed layer 
    500 !          (=0 damping throughout the water column) 
    501 !     (=1 no damping in the mixed layer defined by avt >5cm2/s ) 
    502 !     (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
    503 !  sdmp    surface time scale for internal damping (days) 
    504 !  bdmp    bottom time scale for internal damping (days) 
    505 !  hdmp    depth of transition between sdmp and bdmp (meters) 
    506 &namtdp 
    507    ndmp   =   90  
    508    ndmpf  =    1 
    509    nmldmp =    1 
    510    sdmp   =  50. 
    511    bdmp   = 360. 
    512    hdmp   = 800. 
    513 / 
    514 !----------------------------------------------------------------------- 
    515 !       nameos   ocean physical parameters 
    516 !----------------------------------------------------------------------- 
    517 !  neos    type of equation of state and Brunt-Vaisala frequency 
    518 !          = 0, UNESCO (formulation of Jackett and McDougall (1994) 
    519 !                                         and of McDougall (1987) ) 
    520 !          = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    521 !          = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    522 !                               with rau0=1020 set in parcst routine 
    523 !  ralpha  thermal expension coefficient (linear equation of state) 
    524 !  rbeta   saline  expension coefficient (linear equation of state) 
    525 &nameos 
    526    neos   =      0 
    527    ralpha =  2.e-4 
    528    rbeta  =  0.001 
    529 / 
    530 !----------------------------------------------------------------------- 
    531 !       namsol   elliptic solver / island / free surface  
    532 !----------------------------------------------------------------------- 
    533 !  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
    534 !                            (=2 successive-over-relaxation: sor) 
    535 !                            (=3 FETI: fet, all require "key_feti" defined) 
    536 !  nsol_arp  absolute/relative (0/1) precision convergence test 
    537 !  nmin      minimum of iterations for the SOR solver 
    538 !  nmax      maximum of iterations for the SOR solver 
    539 !  nmod      frequency of test for the SOR solver 
    540 !  eps       absolute precision of the solver 
    541 !  resmax    absolute precision for the SOR solver 
    542 !  sor       optimal coefficient for SOR solver 
    543 !  epsisl    absolute precision on stream function solver 
    544 !  nmisl     maximum pcg iterations for island 
    545 !  rnu       strength of the additional force used in free surface b.c. 
    546 &namsol 
    547    nsolv     =      1 
    548    nsol_arp  =      0 
    549    nmin      =    300 
    550    nmax      =    800 
    551    nmod      =     10 
    552    eps       =  1.E-6 
    553    resmax    = 1.E-10 
    554    sor       =   1.92 
    555    epsisl    = 1.e-10 
    556    nmisl     =   4000 
    557    rnu       =     1. 
    558 / 
    559 !======================================================================= 
    560 !   Diagnostics namelists 
    561 !       namtrd    dynamics and/or tracer trends 
    562 !       namgap    level mean model-data gap 
    563 !       namznl    zonal mean heat & freshwater fluxes computation 
    564 !       namspr    surface pressure in rigid-lid 
    565 !======================================================================= 
    566 !----------------------------------------------------------------------- 
    567 !       namtrd    diagnostics on dynamics and/or tracer trends 
    568 !                         ('key_trdyn' and/or 'key_trdtra') 
    569 !                 or mixed-layer trends ('key_trdmld') 
    570 !----------------------------------------------------------------------- 
    571 !  ntrd              time step frequency dynamics and tracers trends 
    572 !  nctls             control surface type in mixed-layer trends (0,1 or n<jpk) 
    573 !  ln_trdmld_restart restart for ML diagnostics 
    574 !  ucf               unit conversion factor (=1 -> /seconds | =86400. -> /day) 
    575 !  ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 
    576 &namtrd 
    577    ntrd  = 365 
    578    nctls =   0 
    579    ln_trdmld_restart = .false. 
    580    ucf   =  1. 
    581    ln_trdmld_instant = .false. 
    582 / 
    583 !----------------------------------------------------------------------- 
    584 !       namgap    level mean model-data gap ('key_diagap') 
    585 !----------------------------------------------------------------------- 
    586 !  ngap    time-step frequency of model-data gap computation 
    587 !  nprg    time-step frequency of gap print in model output 
    588 &namgap 
    589    ngap =  15 
    590    nprg =  10 
    591 / 
    592 !----------------------------------------------------------------------- 
    593 !       namznl    zonal mean heat & freshwater fluxes computation 
    594 !                 (#ifdef "key_diaznl") 
    595 !----------------------------------------------------------------------- 
    596 !  nfznl   time-step frequency of zonal mean fluxes computation 
    597 &namznl 
    598    nfznl =  15 
    599 / 
    600 !----------------------------------------------------------------------- 
    601 !       namspr  surface pressure diagnostic 
    602 !----------------------------------------------------------------------- 
    603 !  nmaxp   maximum of iterations for the solver 
    604 !  epsp    absolute precision of the solver 
    605 !  niterp  number of iteration done by the solver 
    606 &namspr 
    607    nmaxp   =   1000 
    608    epsp    =  1.e-3 
    609    niterp  =    400 
     313!       namsbc_rnf   runoffs namelist surface boundary condition 
     314!----------------------------------------------------------------------- 
     315!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
     316!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
     317!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
     318!cn_dir      Root directory for location of ssr files 
     319!sn_rnf      informations about the runoff file to be read 
     320!sn_cnf      informations about the runoff mouth file to be read 
     321!  THE ORDER OF THE FILES MATTER: 
     322!               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
     323!               !                    !   (if <0  months) !   name     !     (logical)      ! 
     324&namsbc_rnf 
     325   nn_runoff  = 2 
     326   rn_hrnf    = 0.e0 
     327   rn_avt_rnf = 1.e-3 
     328   cn_dir     = './' 
     329   sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0      
     330   sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0               
     331/ 
     332!----------------------------------------------------------------------- 
     333!       namsbc_ssr   surface boundary condition : sea surface restoring 
     334!----------------------------------------------------------------------- 
     335!  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
     336!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
     337!             (=1, salt flux, concentration/dillution emps only) 
     338!             (=2, volume flux, both emp and emps are updated)  
     339!  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
     340!  deds       magnitude of the damping on salinity [mm/day] 
     341!  THE ORDER OF THE FILES MATTER: 
     342!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     343!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     344&namsbc_ssr 
     345   nn_sstr      = 0 
     346   nn_sssr      = 0 
     347   dqdt         = -40. 
     348   deds         =  27.7 
     349   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0 
     350   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0 
     351/       
     352!----------------------------------------------------------------------- 
     353!       namalb   albedo parameters 
     354!----------------------------------------------------------------------- 
     355!  cgren    correction of the snow or ice albedo to take into account 
     356!  albice   albedo of melting ice in the arctic and antarctic 
     357!  alphd    coefficients for linear interpolation used to compute albedo 
     358!           between two extremes values (Pyane, 1972) 
     359!  alphc     "                                         " 
     360!  alphdi    "                                         " 
     361&namalb 
     362   cgren    =      0.06 
     363   albice   =      0.5 
     364   alphd    =      0.80 
     365   alphc    =      0.65 
     366   alphdi   =      0.72 
     367/ 
     368 
     369!!====================================================================== 
     370!!               ***  Lateral boundary condition  *** 
     371!!====================================================================== 
     372!!   namlbc        lateral momentum boundary condition 
     373!!   namcla        cross land advection 
     374!!   namobc        open boundaries parameters                           ("key_obc") 
     375!!====================================================================== 
     376 
     377!----------------------------------------------------------------------- 
     378!       namlbc   lateral momentum boundary condition 
     379!----------------------------------------------------------------------- 
     380!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip 
     381!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip 
     382&namlbc 
     383   shlat  =      2. 
     384/ 
     385!----------------------------------------------------------------------- 
     386!       namcla   cross land advection 
     387!----------------------------------------------------------------------- 
     388!  n_cla   advection between 2 ocean pts separates by land  
     389&namcla 
     390   n_cla   = 0 
    610391/ 
    611392!----------------------------------------------------------------------- 
     
    639420    ln_vol_cst  = .false. 
    640421/ 
    641 !----------------------------------------------------------------------- 
    642 !       namflo    float parameters (#ifdef key_float) 
     422 
     423!!====================================================================== 
     424!!                 ***  Bottom boundary condition  *** 
     425!!====================================================================== 
     426!!   nambfr        bottom friction 
     427!!   nambbc        bottom temperature boundary condition                ("key_trabbc") 
     428!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv") 
     429!!====================================================================== 
     430!----------------------------------------------------------------------- 
     431!       nambfr   bottom friction 
     432!----------------------------------------------------------------------- 
     433!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction 
     434!                                     = 3 : free slip  ;  = 2 : non linear friction 
     435!  bfri1    bottom drag coefficient (linear case) 
     436!  bfri2    bottom drag coefficient (non linear case) 
     437!  bfeb2    bottom turbulent kinetic energy  (m2/s2) 
     438&nambfr 
     439   nbotfr =       1 
     440   bfri1  =   4.e-4 
     441   bfri2  =   1.e-3 
     442   bfeb2  =  2.5e-3 
     443/ 
     444!----------------------------------------------------------------------- 
     445!       nambbc   bottom temperature boundary condition 
     446!----------------------------------------------------------------------- 
     447!  ngeo_flux  = 0 no geothermal heat flux 
     448!             = 1 constant geothermal heat flux 
     449!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 
     450!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
     451&nambbc 
     452   ngeo_flux =  2 
     453   ngeo_flux_const = 86.4e-3 
     454/ 
     455!----------------------------------------------------------------------- 
     456!       nambbl   bottom boundary layer scheme 
     457!----------------------------------------------------------------------- 
     458!                                !  diffusive bbl                       ("key_trabbl") 
     459!                                !  advective bbl                       ("key_trabbl_adv") 
     460!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s) 
     461&nambbl 
     462   atrbbl = 10000. 
     463/ 
     464 
     465!!====================================================================== 
     466!!                        Tracer (T & S ) namelists 
     467!!====================================================================== 
     468!!   nameos        equation of state 
     469!!   nam_traadv    advection scheme 
     470!!   nam_traldf    lateral diffusion scheme 
     471!!   namtdp        tracer newtonian damping                             ("key_tradmp") 
     472!!====================================================================== 
     473 
     474!----------------------------------------------------------------------- 
     475!       nameos   ocean physical parameters 
     476!----------------------------------------------------------------------- 
     477!  neos    type of equation of state and Brunt-Vaisala frequency 
     478!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
     479!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
     480!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
     481!  ralpha  thermal expension coefficient (neos= 1 or 2) 
     482!  rbeta   saline  expension coefficient (neos= 2) 
     483&nameos 
     484   neos   =      0 
     485   ralpha =  2.e-4 
     486   rbeta  =  0.001 
     487/ 
     488!----------------------------------------------------------------------- 
     489!       nam_traadv   advection scheme for tracer 
     490!----------------------------------------------------------------------- 
     491!  ln_traadv_cen2     2nd order centered scheme          (default T) 
     492!  ln_traadv_tvd      TVD scheme                         (default F) 
     493!  ln_traadv_muscl    MUSCL scheme                       (default F) 
     494!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F) 
     495!  ln_traadv_ubs      UBS scheme                         (default F) 
     496!  ln_traadv_qck      QUICKEST scheme                    (default F) 
     497&nam_traadv 
     498   ln_traadv_cen2   =  .false. 
     499   ln_traadv_tvd    =  .true. 
     500   ln_traadv_muscl  =  .false. 
     501   ln_traadv_muscl2 =  .false. 
     502   ln_traadv_ubs    =  .false. 
     503   ln_traadv_qck    =  .false. 
     504/ 
     505!----------------------------------------------------------------------- 
     506!       nam_traldf   lateral diffusion scheme for tracer 
     507!----------------------------------------------------------------------- 
     508!  Type of the operator :  
     509!     ln_traldf_lap    laplacian operator          (default T) 
     510!     ln_traldf_bilap  bilaplacian operator        (default F) 
     511!  Direction of action  :  
     512!     ln_traldf_level  iso-level                   (default F) 
     513!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T) 
     514!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
     515!  Coefficient 
     516!     aht0    horizontal eddy diffusivity for tracers (m2/s) 
     517!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
     518!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv") 
     519&nam_traldf 
     520   ln_traldf_lap    =  .true. 
     521   ln_traldf_bilap  =  .false. 
     522   ln_traldf_level  =  .false. 
     523   ln_traldf_hor    =  .false. 
     524   ln_traldf_iso    =  .true. 
     525   aht0    =  1000. 
     526   ahtb0   =     0. 
     527   aeiv0   =     0. 
     528/ 
     529!----------------------------------------------------------------------- 
     530!       namtdp   tracer newtonian damping ('key_tradmp') 
     531!----------------------------------------------------------------------- 
     532!  ndmp    type of damping in temperature and salinity  
     533!          !   ='latitude', damping poleward of 'ndmp' degrees and function  
     534!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1 
     535!          !   =-1          damping only in Med and Red Seas 
     536!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
     537!  nmldmp  type of damping in the mixed layer 
     538!          !   =0 damping throughout the water column 
     539!          !   =1 no damping in the mixed layer defined by avt >5cm2/s  
     540!          !   =2 no damping in the mixed layer defined rho<rho(surf)+.01  
     541!  sdmp    surface time scale for internal damping (days) 
     542!  bdmp    bottom time scale for internal damping (days) 
     543!  hdmp    depth of transition between sdmp and bdmp (meters) 
     544&namtdp 
     545   ndmp   =   90 
     546   ndmpf  =    1 
     547   nmldmp =    1 
     548   sdmp   =  50. 
     549   bdmp   = 360. 
     550   hdmp   = 800. 
     551/ 
     552 
     553!!====================================================================== 
     554!!                      ***  Dynamics namelists  *** 
     555!!====================================================================== 
     556!!   nam_dynadv    formulation of the momentum advection 
     557!!   nam_dynvor    advection scheme 
     558!!   nam_dynhpg    hydrostatic pressure gradient 
     559!!   namflg        hydrostatic pressure gradient time stepping 
     560!!   nam_dynspg    surface pressure gradient                            (CPP key only) 
     561!!   nam_dynldf    lateral diffusion scheme 
     562!!====================================================================== 
     563 
     564!----------------------------------------------------------------------- 
     565!       nam_dynadv   formulation of the momentum advection 
     566!----------------------------------------------------------------------- 
     567!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T) 
     568!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
     569!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
     570&nam_dynadv    
     571   ln_dynadv_vec   = .true. 
     572   ln_dynadv_cen2  = .false. 
     573   ln_dynadv_ubs   = .false. 
     574 
     575!----------------------------------------------------------------------- 
     576!       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
     577!----------------------------------------------------------------------- 
     578!  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
     579!  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
     580!  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
     581!  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
     582&nam_dynvor 
     583   ln_dynvor_ene = .false. 
     584   ln_dynvor_ens = .false. 
     585   ln_dynvor_mix = .false. 
     586   ln_dynvor_een = .true. 
     587/ 
     588!----------------------------------------------------------------------- 
     589!       nam_dynhpg   Hydrostatic pressure gradient option 
     590!----------------------------------------------------------------------- 
     591!  type of pressure gradient scheme (choose one only!) 
     592!     ln_hpg_zco    z-coordinate - full steps                   (default T) 
     593!     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
     594!     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
     595!     ln_hpg_hel    s-coordinate (helsinki modification) 
     596!     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
     597!     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
     598!     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
     599!  parameters 
     600!     gamm          weighting coefficient (wdj scheme) 
     601&nam_dynhpg 
     602   ln_hpg_zco = .false. 
     603   ln_hpg_zps = .true. 
     604   ln_hpg_sco = .false. 
     605   ln_hpg_hel = .false. 
     606   ln_hpg_wdj = .false. 
     607   ln_hpg_djc = .false. 
     608   ln_hpg_rot = .false. 
     609   gamm       = 0.e0 
     610/ 
     611!----------------------------------------------------------------------- 
     612!       namflg   algorithm flags (algorithm not control by CPP keys) 
     613!----------------------------------------------------------------------- 
     614!  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
     615!                                                  centered     time scheme  (F) 
     616!   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
     617&namflg 
     618   ln_dynhpg_imp   =  .false. 
     619   nn_dynhpg_rst   =  0 
     620/ 
     621!----------------------------------------------------------------------- 
     622!       nam_dynspg   surface pressure gradient   (CPP key only) 
     623!----------------------------------------------------------------------- 
     624!                                !  explicit free surface               ("key_dynspg_exp") 
     625!                                !  filtered free surface               ("key_dynspg_flt") 
     626!                                !  split-explicit free surface         ("key_dynspg_ts") 
     627!                                !  rigid-lid                           ("key_dynspg_rl") 
     628 
     629!----------------------------------------------------------------------- 
     630!       nam_dynldf   lateral diffusion on momentum 
     631!----------------------------------------------------------------------- 
     632!  Type of the operator :  
     633!     ln_dynldf_lap    laplacian operator          (default T) 
     634!     ln_dynldf_bilap  bilaplacian operator        (default F) 
     635!  Direction of action  :  
     636!     ln_dynldf_level  iso-level                   (default F) 
     637!     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T) 
     638!     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
     639!  Coefficient 
     640!     ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
     641!     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
     642&nam_dynldf 
     643   ln_dynldf_lap    =  .false. 
     644   ln_dynldf_bilap  =  .true. 
     645   ln_dynldf_level  =  .false. 
     646   ln_dynldf_hor    =  .true. 
     647   ln_dynldf_iso    =  .false. 
     648   ahm0    = -8.5e+11 
     649   ahmb0   =     0. 
     650/ 
     651 
     652!!====================================================================== 
     653!!             Tracers & Dynamics vertical physics namelists 
     654!!====================================================================== 
     655!!       namzdf        vertical physics 
     656!!       namnpc        non penetrative convection                        
     657!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      ) 
     658!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      ) 
     659!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      ) 
     660!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      ) 
     661!!====================================================================== 
     662 
     663!----------------------------------------------------------------------- 
     664!       namzdf   vertical physics 
     665!----------------------------------------------------------------------- 
     666!  ln_zdfevd  enhanced vertical diffusion         (default T) 
     667!  ln_zdfnpc  Non-Penetrative Convection          (default T) 
     668!  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst") 
     669!  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst") 
     670!  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
     671!  n_evdm     = 0  apply enhanced mixing on tracer only 
     672!             = 1  apply enhanced mixing on both tracer and momentum 
     673!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
     674!                               (=F)  euler backward (F) 
     675!  n_zdfexp   number of sub-timestep for time splitting scheme 
     676&namzdf 
     677   ln_zdfevd = .true. 
     678   ln_zdfnpc = .false. 
     679   avm0      = 1.2e-4 
     680   avt0      = 1.2e-5 
     681   avevd     =   100. 
     682   n_evdm    =     1 
     683   ln_zdfexp =  .false. 
     684   n_zdfexp  =     3 
     685/ 
     686!----------------------------------------------------------------------- 
     687!       namnpc   non penetrative convection 
     688!----------------------------------------------------------------------- 
     689!  nnpc1   non penetrative convective scheme frequency    
     690!  nnpc2   non penetrative convective scheme print frequency 
     691&namnpc 
     692   nnpc1  =      1 
     693   nnpc2  =    365 
     694/ 
     695!----------------------------------------------------------------------- 
     696!       namric   richardson number dependent vertical diffusion         ("key_zdfric" ) 
     697!----------------------------------------------------------------------- 
     698!  avmri   maximum value of the vertical viscosity 
     699!  alp     coefficient of the parameterization 
     700!  nric    coefficient of the parameterization 
     701&namric 
     702   avmri = 100.e-4 
     703   alp   =      5. 
     704   nric  =       2 
     705/ 
     706!----------------------------------------------------------------------- 
     707!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
     708!----------------------------------------------------------------------- 
     709!  ln_rstke flag to restart with tke from a run without tke (default F) 
     710!  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
     711!  ediss    coef. of the Kolmogoroff dissipation   
     712!  ebb      coef. of the surface input of tke 
     713!  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
     714!  emin     minimum value of tke (m^2/s^2) 
     715!  emin0    surface minimum value of tke (m^2/s^2) 
     716!  nitke    number of restart iterative loops 
     717!  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
     718!  nmxl     length used   = 0 bounded by the distance to surface and bottom 
     719!           !             = 1 bounded by the local vertical scale factor 
     720!           !             = 2 first vertical derivative of mixing length bounded by 1 
     721!  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
     722!           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
     723!           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
     724!  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
     725!  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
     726&namtke 
     727   ln_rstke = .false. 
     728   ediff =       0.1 
     729   ediss =       0.7 
     730   ebb   =       60. 
     731   efave =        1. 
     732   emin  =     1.e-6 
     733   emin0 =     1.e-4 
     734   nitke =        50 
     735   nmxl  =         2 
     736   npdl  =         1 
     737   navb  =         0 
     738/ 
     739!----------------------------------------------------------------------- 
     740!       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" ) 
     741!----------------------------------------------------------------------- 
     742!  ln_kpprimix   shear instability mixing  (default T) 
     743!  difmiw        constant internal wave viscosity (m2/s) 
     744!  difsiw        constant internal wave diffusivity (m2/s) 
     745!  Riinfty       local Richardson Number limit for shear instability 
     746!  difri         maximum shear mixing at Rig = 0    (m2/s) 
     747!  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
     748!  difcon        maximum mixing in interior convection (m2/s) 
     749!  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
     750!  navb          = 0/1 flag for constant or profile background avt 
     751&namkpp 
     752   ln_kpprimix  = .true. 
     753   difmiw       =  1.e-04 
     754   difsiw       =  0.1e-04 
     755   Riinfty      =  0.8 
     756   difri        =  0.0050 
     757   bvsqcon      = -0.01e-07 
     758   difcon       =  1. 
     759   navb         =  0 
     760   nave         =  1 
     761/ 
     762!----------------------------------------------------------------------- 
     763!       namddm   double diffusive mixing parameterization               ("key_zdfddm") 
     764!----------------------------------------------------------------------- 
     765!   avts    maximum avs for dd mixing  
     766!   hsbfr   heat/salt buoyancy flux ratio 
     767&namddm 
     768      avts  = 1.e-4 
     769      hsbfr = 1.6 
     770/ 
     771 
     772!!====================================================================== 
     773!!                  ***  Miscelaneous namelists  *** 
     774!!====================================================================== 
     775!!   namsol      elliptic solver / island / free surface  
     776!!====================================================================== 
     777 
     778!----------------------------------------------------------------------- 
     779!       namsol   elliptic solver / island / free surface  
     780!----------------------------------------------------------------------- 
     781!  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
     782!                            (=2 successive-over-relaxation: sor) 
     783!                            (=3 FETI currently it does not work!       ("key_feti") 
     784!  nsol_arp  absolute/relative (0/1) precision convergence test 
     785!  nmin      minimum of iterations for the SOR solver 
     786!  nmax      maximum of iterations for the SOR solver 
     787!  nmod      frequency of test for the SOR solver 
     788!  eps       absolute precision of the solver 
     789!  resmax    absolute precision for the SOR solver 
     790!  sor       optimal coefficient for SOR solver 
     791!  epsisl    absolute precision on stream function solver 
     792!  nmisl     maximum pcg iterations for island 
     793!  rnu       strength of the additional force used in free surface b.c. 
     794&namsol 
     795   nsolv     =      1 
     796   nsol_arp  =      0 
     797   nmin      =    300 
     798   nmax      =    800 
     799   nmod      =     10 
     800   eps       =  1.E-6 
     801   resmax    = 1.E-10 
     802   sor       =   1.92 
     803   epsisl    = 1.e-10 
     804   nmisl     =   4000 
     805   rnu       =     1. 
     806/ 
     807!!====================================================================== 
     808!!                       ***  Diagnostics namelists  *** 
     809!!====================================================================== 
     810!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld") 
     811!!   namgap       level mean model-data gap                             ("key_diagap") 
     812!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr") 
     813!!   namflo       float parameters                                      ("key_float") 
     814!!   namptr       Poleward Transport Diagnostics 
     815!!====================================================================== 
     816 
     817!----------------------------------------------------------------------- 
     818!       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra") 
     819!                          or mixed-layer trends                        ('key_trdmld') 
     820!                          or barotropic vorticity                      ("key_trdvor") 
     821!----------------------------------------------------------------------- 
     822!  ntrd                time step frequency dynamics and tracers trends 
     823!  nctls               control surface type in mixed-layer trends (0,1 or n<jpk) 
     824!  ln_trdmld_restart   restart for ML diagnostics 
     825!  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day) 
     826!  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S 
     827&namtrd 
     828   ntrd              = 365 
     829   nctls             =   0 
     830   ln_trdmld_restart = .false. 
     831   ucf               =  1. 
     832   ln_trdmld_instant = .false. 
     833/ 
     834!----------------------------------------------------------------------- 
     835!       namgap    level mean model-data gap ('key_diagap') 
     836!----------------------------------------------------------------------- 
     837!  ngap    time-step frequency of model-data gap computation 
     838!  nprg    time-step frequency of gap print in model output 
     839&namgap 
     840   ngap =  15 
     841   nprg =  10 
     842/ 
     843!----------------------------------------------------------------------- 
     844!       namspr  surface pressure diagnostic 
     845!----------------------------------------------------------------------- 
     846!  nmaxp   maximum of iterations for the solver 
     847!  epsp    absolute precision of the solver 
     848!  niterp  number of iteration done by the solver 
     849&namspr 
     850   nmaxp   =   1000 
     851   epsp    =  1.e-3 
     852   niterp  =    400 
     853/ 
     854!----------------------------------------------------------------------- 
     855!       namflo    float parameters                                      ("key_float") 
    643856!----------------------------------------------------------------------- 
    644857!  ln_rstflo   boolean term for float restart (true or false) 
     
    655868    ln_flork4 = .false. 
    656869/ 
    657 !!====================================================================== 
    658 !!              Surface Boundary Condition namelists 
    659 !! 
    660 !!       namsbc    surface boundary condition 
    661 !!       namsbc_ana    analytical         formulation of sbc 
    662 !!       namsbc_flx    flux               formulation of sbc 
    663 !!       namsbc_core   CORE bulk formulea formulation of sbc 
    664 !!       namsbc_clio   CLIO bulk formulea formulation of sbc 
    665 !!       namsbc_cpl    CouPLed            formulation of sbc 
    666 !!       namsbc_ssr    sea surface restoring term (for T and/or S) 
    667 !!====================================================================== 
    668 !----------------------------------------------------------------------- 
    669 !       namsbc       Surface Boundary Condition (surface module) 
    670 !----------------------------------------------------------------------- 
    671 !  nn_fsbc          frequency of surface boundary condition computation  
    672 !                   (also correspond to the frequency of sea-ice model call) 
    673 ! 
    674 !  ln_ana        analytical formulation (fill namsbc_ana )  
    675 !  ln_flx        flux formulation       (fill namsbc_flx ) (overwritten by key_sbc_flux     ) 
    676 !  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core) (overwritten by key_sbc_bulk_clio) 
    677 !  ln_blk_core   CORE bulk formulation  (fill namsbc_clio) (overwritten by key_sbc_bulk_core) 
    678 !  ln_cpl        Coupled formulation    (fill namsbc_cpl ) (overwritten by key_sbc_coupled  ) 
    679 ! 
    680 !  nn_ice        =0 no ice boundary condition   , 
    681 !                =1 observed ice-cover          , 
    682 !                =2 ice-model ice-cover (overwritten by key_ice_lim) 
    683 ! 
    684 !  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
    685 ! 
    686 !  ln_ssr       Sea Surface restoring on SST and/or SSS (fill namsbc_ssr) 
    687 !               (overwritten by key_dta_sst and/or key_dtasss) 
    688 ! 
    689 !  nn_fwb       FreshWater Budget: =0 no control                              ,  
    690 !                                  =1 annual global mean of emp set to zero   , 
    691 !                                  =2 global mean of emp set to zero at each nn_fsbc time step 
    692 &namsbc 
    693    nn_fsbc     = 5 
    694    ln_ana      = .false. 
    695    ln_flx      = .false. 
    696    ln_blk_clio = .false. 
    697    ln_blk_core = .true. 
    698    ln_cpl      = .false. 
    699    nn_ice      = 2 
    700    ln_dm2dc    = .false. 
    701    ln_rnf      = .false. 
    702    ln_ssr      = .false. 
    703    nn_fwb      = 0 
    704 / 
    705 !----------------------------------------------------------------------- 
    706 !       namsbc_rnf   runoffs namelist surface boundary condition 
    707 !----------------------------------------------------------------------- 
    708 !nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
    709 !rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
    710 !rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
    711 !cn_dir      Root directory for location of ssr files 
    712 !sn_rnf      informations about the runoff file to be read 
    713 !sn_cnf      informations about the runoff mouth file to be read 
    714 !  THE ORDER OF THE FILES MATTER: 
    715 !               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
    716 !               !                    !   (if <0  months) !   name     !     (logical)      ! 
    717 &namsbc_rnf 
    718    nn_runoff  = 2 
    719    rn_hrnf    = 0.e0 
    720    rn_avt_rnf = 1.e-3 
    721    cn_dir     = './' 
    722    sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .TRUE.   ,  1  ,  0      
    723    sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .FALSE.  ,  1  ,  0               
    724 / 
    725 !----------------------------------------------------------------------- 
    726 !       namsbc_ana   analytical surface boundary condition 
    727 !----------------------------------------------------------------------- 
    728 !  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
    729 !  rn_utau0    uniform value used as default surface heat flux 
    730 !  rn_vtau0    uniform value used as default solar radiation flux 
    731 !  rn_q0       uniform value used as default surface heat flux 
    732 !  rn_qsr0     uniform value used as default solar radiation flux 
    733 !  rn_emp0     uniform value used as default surface freswater budget (E-P) 
    734 &namsbc_ana 
    735    nn_tau000  =   0 
    736    rn_utau0   =   0.5 
    737    rn_vtau0   =   0.e0 
    738    rn_q0      =   0.e0 
    739    rn_qsr0    =   0.e0 
    740    rn_emp0    =   0.e0 
    741 / 
    742 !----------------------------------------------------------------------- 
    743 !       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
    744 !----------------------------------------------------------------------- 
    745 !  cn_dir   directory in which the model is executed 
    746 ! 
    747 !  THE ORDER OF THE FILES MATTER: 
    748 !               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    749 !               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    750 &namsbc_flx 
    751    cn_dir       = './' 
    752    sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .FALSE.       ,    0   ,     0 
    753    sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .FALSE.       ,    0   ,     0 
    754    sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .FALSE.       ,    0   ,     0 
    755    sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .FALSE.       ,    0   ,     0 
    756    sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .FALSE.       ,    0   ,     0 
    757  
    758 /       
    759 !----------------------------------------------------------------------- 
    760 !       namsbc_clio  CLIO bulk formulea 
    761 !----------------------------------------------------------------------- 
    762 !  cn_dir   directory in which the model is executed 
    763 ! 
    764 !  THE ORDER OF THE FILES MATTER: 
    765 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    766 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    767 &namsbc_clio 
    768    cn_dir       = './' 
    769    sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .FALSE.        ,    1   ,    0     
    770    sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .FALSE.        ,    1   ,    0     
    771    sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .FALSE.        ,    1   ,    0       
    772    sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .FALSE.        ,    1   ,    0      
    773    sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .FALSE.        ,    1   ,    0     
    774    sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .FALSE.        ,    1   ,    0       
    775    sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .FALSE.        ,    1   ,    0        
    776 / 
    777 !----------------------------------------------------------------------- 
    778 !       namsbc_core  CORE bulk formulea 
    779 !----------------------------------------------------------------------- 
    780 !  cn_dir         directory in which the model is executed 
    781 !  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
    782 !  alpha_precip   multiplication factor for precipitation (total & snow) 
    783 ! 
    784 !  THE ORDER OF THE FILES MATTER: 
    785 !               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
    786 !               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
    787 &namsbc_core 
    788    cn_dir       = './' 
    789    ln_2m        = .FALSE. 
    790    alpha_precip = 1. 
    791    sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .FALSE.        ,    1   ,    0     
    792    sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .FALSE.        ,    1   ,    0     
    793    sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .FALSE.        ,    1   ,    0       
    794    sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .FALSE.        ,    1   ,    0       
    795    sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .FALSE.        ,    1   ,    0      
    796    sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .FALSE.        ,    1   ,    0     
    797    sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .FALSE.        ,    1   ,    0       
    798    sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .FALSE.        ,    1   ,    0        
    799 / 
    800 !----------------------------------------------------------------------- 
    801 !       namsbc_cpl    coupled ocean/atmosphere model  (#ifdef "key_coupled") 
    802 !----------------------------------------------------------------------- 
    803 &namsbc_cpl 
    804 / 
    805 !----------------------------------------------------------------------- 
    806 !       namsbc_ssr   surface boundary condition : sea surface restoring 
    807 !----------------------------------------------------------------------- 
    808 !  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
    809 !  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
    810 !             (=1, salt flux, concentration/dillution emps only) 
    811 !             (=2, volume flux, both emp and emps are updated)  
    812 !  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
    813 !  deds       magnitude of the damping on salinity [mm/day] 
    814 !  THE ORDER OF THE FILES MATTER: 
    815 !               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
    816 !               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
    817 &namsbc_ssr 
    818    nn_sstr      = 0 
    819    nn_sssr      = 0 
    820    dqdt         = -40. 
    821    deds         =  27.7 
    822    sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .FALSE.       ,    0   ,     0 
    823    sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .TRUE.        ,    0   ,     0 
    824 /       
    825 !----------------------------------------------------------------------- 
    826 !       namagrif    agrif parameters (#ifdef key_agrif) 
    827 !----------------------------------------------------------------------- 
    828 ! nbclineupdate baroclinic update frequency  
    829 ! ln_spc_dyn use 0 as special value for dynamics 
    830 ! visc_tra viscosity coeeficient for tracers sponge layer 
    831 ! visc_dyn viscosity coeeficient for dynamics sponge layer 
    832 &namagrif 
    833     nbclineupdate = 3 
    834     ln_spc_dyn = .TRUE. 
    835     visc_tra = 2880. 
    836     visc_dyn = 2880. 
    837 / 
     870!----------------------------------------------------------------------- 
     871!       namptr   Poleward Transport Diagnostic 
     872!----------------------------------------------------------------------- 
     873!  ln_diaptr  logical flag for Poleward transport computation 
     874!  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
     875!             need input basins mask file named "subbasins.nc" 
     876!  nf_ptr     Frequency of computation 
     877&namptr 
     878   ln_diaptr = .false. 
     879   ln_subbas = .false. 
     880   nf_ptr    =  15 
     881/ 
  • trunk/CONFIG/ORCA2_LIM/EXP00/namelist

    r889 r907  
    1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! OPA namelist :  model option and parameter input 
    3 ! ------------- 
    4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    5 ! 
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/OPA  :  1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl) 
     3!! namelists    2 - Domain           (nam_zgr, nam_zgr_sco, namdom) 
     4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
     5!!                                    namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) 
     6!!              4 - lateral boundary (namlbc, namcla, namobc) 
     7!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
     8!!              6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp) 
     9!!              7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) 
     10!!              8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm) 
     11!!              9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr) 
     12!!              9 - miscellaneous    (namsol) 
     13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     14!  CAUTION: some scripts does not support CAPITALs for logical definition  
     15!  *******       use .true. or .false. and NOT .TRUE. or .FALSE. 
     16 
     17!!====================================================================== 
     18!!                   ***  Run management namelists  *** 
     19!!====================================================================== 
     20!!   namrun            parameters of the run 
     21!!   nam_mpp           Massively Parallel Processing 
     22!!   nam_mpp_dyndist    
     23!!   nam_ctl           Control prints & Benchmark 
     24!!====================================================================== 
     25 
    626!----------------------------------------------------------------------- 
    727!       namrun   parameters of the run 
     
    1030!  cexper     experience name for vairmer format 
    1131!  ln_rstart  boolean term for restart (true or false) 
    12 !  nrstdt     control of the restart timestep: 
    13 !                = 0 restart, do not control nit000 in the restart file. 
    14 !                = 1 restart, control nit000 in the restart file. Do not 
    15 !                    use the date in the restart file (use ndate0 in namelist) 
    16 !                = 2 restart, control nit000 in the restart file, use the date 
    17 !                    in the restart file. ndate0 in the namelist is ignored. 
     32!  nrstdt     restart control = 0 restart, do not control nit000 in the restart file. 
     33!             !               = 1 restart, control nit000 in the restart file. Do not 
     34!             !                   use the date in the restart file (use ndate0 in namelist) 
     35!             !               = 2 restart, control nit000 in the restart file, use the date 
     36!             !                   in the restart file. ndate0 in the namelist is ignored. 
    1837!  nit000     number of the first time step 
    1938!  nitend     number of the last time step 
     
    2443!  nwrite     frequency of OUTPUT file 
    2544!  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc) 
    26 ! 
    27 !  CAUTION: for usual run scripts, logical value of 
    28 !  *******  ln_rstart must be .true. or .false. 
    29 !                     and NOT .TRUE. or .FALSE. 
    3045&namrun 
    3146   no         =       0 
     
    4055   nstock     =    5475  
    4156   nwrite     =    5475  
    42    ln_dimgnnn = .FALSE. 
     57   ln_dimgnnn = .false. 
     58/ 
     59!----------------------------------------------------------------------- 
     60!       nam_mpp      Massively Parallel Processing 
     61!----------------------------------------------------------------------- 
     62!  c_mpi_send         mpi send/recieve type 
     63!                      = 'S'  : standard blocking send 
     64!                      = 'B'  : buffer blocking send 
     65!                      = 'I'  : immediate non-blocking send 
     66&nam_mpp 
     67   c_mpi_send =  'S' 
     68/ 
     69!----------------------------------------------------------------------- 
     70!       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist") 
     71!----------------------------------------------------------------------- 
     72!  jpni   number of processors following i 
     73!  jpnj   number of processors following j 
     74!  jpnij  number of local domains 
     75&nam_mpp_dyndist 
     76   jpni  = 1 
     77   jpnj  = 1 
     78   jpnij = 1 
    4379/ 
    4480!----------------------------------------------------------------------- 
     
    69105   nbit_cmp =     0 
    70106/ 
    71 !----------------------------------------------------------------------- 
    72 !       nam_mpp      Massively Parallel Processing 
    73 !----------------------------------------------------------------------- 
    74 !  c_mpi_send         mpi send/recieve type 
    75 !                      = 'S'  : standard blocking send 
    76 !                      = 'B'  : buffer blocking send 
    77 !                      = 'I'  : immediate non-blocking send 
    78 &nam_mpp 
    79    c_mpi_send =  'S' 
    80 / 
    81 !----------------------------------------------------------------------- 
    82 !       nam_mpp_dyndist      Massively Parallel Distribution 
    83 !              #ifdef ( key_agrif && key_mpp_dyndist ) 
    84 !----------------------------------------------------------------------- 
    85 !  jpni   number of processors following i 
    86 !  jpnj   number of processors following j 
    87 !  jpnij  number of local domains 
    88 &nam_mpp_dyndist 
    89    jpni  = 1 
    90    jpnj  = 1 
    91    jpnij = 1 
    92 / 
     107 
     108!!====================================================================== 
     109!!                      ***  Domain namelists  *** 
     110!!====================================================================== 
     111!!   nam_zgr       vertical coordinate 
     112!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate 
     113!!   namdom        space and time domain (bathymetry, mesh, timestep) 
     114!!====================================================================== 
     115 
    93116!----------------------------------------------------------------------- 
    94117!       nam_zgr       vertical coordinate 
    95118!----------------------------------------------------------------------- 
    96 !  ln_zco     z-coordinate - full steps      (T/F) 
     119!  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined) 
    97120!  ln_zps     z-coordinate - partial steps   (T/F) 
    98121!  ln_sco     s- or hybrid z-s-coordinate    (T/F) 
     
    118141/ 
    119142!----------------------------------------------------------------------- 
    120 !       nam_traadv   advection scheme for tracer (option not control by CPP keys) 
    121 !----------------------------------------------------------------------- 
    122 !  ln_traadv_cen2     2nd order centered scheme    (default T) 
    123 !  ln_traadv_tvd      TVD scheme                   (default F) 
    124 !  ln_traadv_muscl    MUSCL scheme                 (default F) 
    125 !  ln_traadv_muscl2   MUSCL2 scheme                (default F) 
    126 !  ln_traadv_ubs      UBS scheme                   (default F) 
    127 !  ln_traadv_qck      QUICKEST scheme              (default F) 
    128 &nam_traadv 
    129    ln_traadv_cen2   =  .true. 
    130    ln_traadv_tvd    =  .false. 
    131    ln_traadv_muscl  =  .false. 
    132    ln_traadv_muscl2 =  .false. 
    133    ln_traadv_ubs    =  .false. 
    134    ln_traadv_qck    =  .false. 
    135 / 
    136 !----------------------------------------------------------------------- 
    137 !       nam_traldf   lateral diffusion scheme for tracer (option not control by CPP keys) 
    138 !----------------------------------------------------------------------- 
    139 !  Type of the operator :  
    140 !     ln_traldf_lap    laplacian operator          (default T) 
    141 !     ln_traldf_bilap  bilaplacian operator        (default F) 
    142 !  Direction of action  :  
    143 !     ln_traldf_level  iso-level                   (default F) 
    144 !     ln_traldf_hor    horizontal (geopotential)   (default F)^** 
    145 !     ln_traldf_iso    iso-neutral                 (default T)^* 
    146 !  Coefficient 
    147 !     aht0    horizontal eddy diffusivity for tracers (m2/s) 
    148 !     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
    149 !     aeiv0   eddy induced velocity coefficient (m2/s) 
    150 ! ^* require key_ldfslp to compute the direction of the lateral diffusion 
    151 ! ^** require key_ldfslp in s-coordinate 
    152 &nam_traldf 
    153    ln_traldf_lap    =  .true. 
    154    ln_traldf_bilap  =  .false. 
    155    ln_traldf_level  =  .false. 
    156    ln_traldf_hor    =  .false. 
    157    ln_traldf_iso    =  .true. 
    158    aht0    =  2000. 
    159    ahtb0   =     0. 
    160    aeiv0   =  2000. 
    161 / 
    162 !----------------------------------------------------------------------- 
    163 !       nam_dynldf   lateral diffusion on momentum 
    164 !----------------------------------------------------------------------- 
    165 !  Type of the operator :  
    166 !     ln_dynldf_lap    laplacian operator          (default T) 
    167 !     ln_dynldf_bilap  bilaplacian operator        (default F) 
    168 !  Direction of action  :  
    169 !     ln_dynldf_level  iso-level                   (default F) 
    170 !     ln_dynldf_hor    horizontal (geopotential)   (default F)^** 
    171 !     ln_dynldf_iso    iso-neutral                 (default T)^* 
    172 !  Coefficient 
    173 !  ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
    174 !  ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
    175 &nam_dynldf 
    176    ln_dynldf_lap    =  .true. 
    177    ln_dynldf_bilap  =  .false. 
    178    ln_dynldf_level  =  .false. 
    179    ln_dynldf_hor    =  .true. 
    180    ln_dynldf_iso    =  .false. 
    181    ahm0    = 40000. 
    182    ahmb0   =     0. 
    183 / 
    184 !----------------------------------------------------------------------- 
    185 !       namflg   algorithm flags (algorithm not control by CPP keys) 
    186 !----------------------------------------------------------------------- 
    187 !  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
    188 !                                                  centered      time scheme  (F) 
    189 !   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0) 
    190 &namflg 
    191    ln_dynhpg_imp   =  .false. 
    192    nn_dynhpg_rst   =  0 
    193 / 
    194 !----------------------------------------------------------------------- 
    195 !       nam_dynhpg   Hydrostatic pressure gradient option 
    196 !----------------------------------------------------------------------- 
    197 !  type of pressure gradient scheme (choose one only!) 
    198 !     ln_hpg_zco    z-coordinate - full steps                   (default T) 
    199 !     ln_hpg_zps    z-coordinate - partial steps (interpolation) 
    200 !     ln_hpg_sco    s-coordinate (standard jacobian formulation) 
    201 !     ln_hpg_hel    s-coordinate (helsinki modification) 
    202 !     ln_hpg_wdj    s-coordinate (weighted density jacobian) 
    203 !     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial) 
    204 !     ln_hpg_rot    s-coordinate (ROTated axes scheme) 
    205 !  parameters 
    206 !    gamm          weighting coefficient (wdj scheme) 
    207 &nam_dynhpg 
    208    ln_hpg_zco = .false. 
    209    ln_hpg_zps = .true. 
    210    ln_hpg_sco = .false. 
    211    ln_hpg_hel = .false. 
    212    ln_hpg_wdj = .false. 
    213    ln_hpg_djc = .false. 
    214    ln_hpg_rot = .false. 
    215    gamm       = 0.e0 
    216 / 
    217 !----------------------------------------------------------------------- 
    218 !       nam_dynadv   option of physics/algorithm (not control by CPP keys) 
    219 !----------------------------------------------------------------------- 
    220 !  ln_dynadv_vec      vector form flag 
    221 !  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
    222 !  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F) 
    223 &nam_dynadv    
    224    ln_dynadv_vec   = .TRUE. 
    225    ln_dynadv_cen2  = .FALSE. 
    226    ln_dynadv_ubs   = .FALSE. 
    227  
    228 !----------------------------------------------------------------------- 
    229 !       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
    230 !----------------------------------------------------------------------- 
    231 !  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
    232 !  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
    233 !  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
    234 !  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
    235 &nam_dynvor 
    236    ln_dynvor_ene = .FALSE. 
    237    ln_dynvor_ens = .TRUE. 
    238    ln_dynvor_mix = .FALSE. 
    239    ln_dynvor_een = .FALSE. 
    240 / 
    241 !----------------------------------------------------------------------- 
    242 !       namalb   albedo parameters 
    243 !----------------------------------------------------------------------- 
    244 !  cgren    correction of the snow or ice albedo to take into account 
    245 !  albice   albedo of melting ice in the arctic and antarctic 
    246 !  alphd    coefficients for linear interpolation used to compute albedo 
    247 !           between two extremes values (Pyane, 1972) 
    248 !  alphc     "                                         " 
    249 !  alphdi    "                                         " 
    250 &namalb 
    251    cgren    =      0.06 
    252    albice   =      0.5 
    253    alphd    =      0.80 
    254    alphc    =      0.65 
    255    alphdi   =      0.72 
    256 / 
    257 !----------------------------------------------------------------------- 
    258143!       namdom   space and time domain (bathymetry, mesh, timestep) 
    259144!----------------------------------------------------------------------- 
    260 !  ntopo      = 0/1 ,compute/read the bathymetry file 
    261 !               (mbathy, nb of T-ocean levels) 
     145!  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels) 
    262146!  e3zps_min  the thickness of the partial step is set larger than the 
    263 !  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t 
    264 !                (N.B. 0<e3zps_rat<1) 
     147!  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1) 
    265148!  nmsh       =1 create a mesh file (coordinates, scale factors, masks) 
    266 !  nacc       the acceleration of convergence method 
    267 !             = 0, no acceleration, rdt = rdttra 
    268 !             = 1, acceleration used, rdt < rdttra(k) 
     149!  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra 
     150!                                                  = 1,     used, rdt < rdttra(k) 
    269151!  atfp       asselin time filter parameter 
    270152!  rdt        time step for the dynamics (and tracer if nacc=0) 
     
    273155!  rdth       depth variation of tracer time step 
    274156!  rdtbt      barotropic time step (for the time splitting algorithm) 
    275 !  nclosea    0 no closed sea  
    276 !             1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
     157!  nclosea    =0 remove the closed sea from the global domain (orca configuration) 
     158!             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)  
    277159&namdom 
    278160   ntopo     =     1 
     
    289171   nclosea   =     0 
    290172/ 
    291 !----------------------------------------------------------------------- 
    292 !       namfwb   freshwater budget correction 
    293 !----------------------------------------------------------------------- 
    294 !  ln_fwb     logical flag for freshwater budget correction (0 annual mean) 
    295 &namfwb 
    296    ln_fwb    = .true. 
    297 / 
    298 !----------------------------------------------------------------------- 
    299 !       namptr   Poleward Transport Diagnostic 
    300 !----------------------------------------------------------------------- 
    301 !  ln_diaptr  logical flag for Poleward transport computation 
    302 !  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation 
    303 !             need input basins mask file named "subbasins.nc" 
    304 !  nf_ptr     Frequency of computation 
    305 &namptr 
    306    ln_diaptr = .false. 
    307    ln_subbas = .false. 
    308    nf_ptr    =  15 
    309 / 
    310 !----------------------------------------------------------------------- 
    311 !       namcro   cross land advection 
    312 !----------------------------------------------------------------------- 
    313 !  n_cla   advection between 2 ocean pts separates by land  
    314 &namcla 
    315    n_cla   = 1 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 !       namzdf   vertical physics 
    319 !----------------------------------------------------------------------- 
    320 !  ln_zdfevd  enhanced vertical diffusion         (default T) 
    321 !  ln_zdfnpc  Non-Penetrative Convection          (default T) 
    322 !  avm0       vertical eddy viscosity for the dynamic (m2/s) 
    323 !  avt0       vertical eddy diffusivity for tracers (m2/s) 
    324 !  avevd      vertical coefficient for enhanced diffusion scheme (m2/s) 
    325 !  nevdm      = 0  apply enhanced mixing on tracer only 
    326 !             = 1  apply enhanced mixing on both tracer and momentum 
    327 !  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
    328 !                               (=F)  euler backward (F) 
    329 !  n_zdfexp   number of sub-timestep for time splitting scheme 
    330 &namzdf 
    331    ln_zdfevd = .true. 
    332    ln_zdfnpc = .false. 
    333    avm0      = 1.2e-4 
    334    avt0      = 1.2e-5 
    335    avevd     =   100. 
    336    n_evdm    =     1 
    337    ln_zdfexp =  .false. 
    338    n_zdfexp  =     3 
    339 / 
    340 !----------------------------------------------------------------------- 
    341 !       namnpc   vnon penetrative convection 
    342 !----------------------------------------------------------------------- 
    343 !  nnpc1   non penetrative convective scheme frequency 
    344 !  nnpc2   non penetrative convective scheme print frequency 
    345 &namnpc 
    346    nnpc1  =      1 
    347    nnpc2  =    365 
    348 / 
    349 !----------------------------------------------------------------------- 
    350 !       nambbl   bottom boundary layer scheme 
    351 !----------------------------------------------------------------------- 
    352 !  atrbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s) 
    353 &nambbl 
    354    atrbbl = 10000. 
    355 / 
    356 !----------------------------------------------------------------------- 
    357 !       namric   richardson number dependent vertical diffusion 
    358 !                ( #ifdef "key_zdfrichardson" ) 
    359 !----------------------------------------------------------------------- 
    360 !  avmri   maximum value of the vertical viscosity 
    361 !  alp     coefficient of the parameterization 
    362 !  nric    coefficient of the parameterization 
    363 &namric 
    364    avmri = 100.e-4 
    365    alp   =      5. 
    366    nric  =       2 
    367 / 
    368 !----------------------------------------------------------------------- 
    369 !       namtke   turbulent eddy kinetic dependent vertical diffusion 
    370 !                ( #ifdef "key_zdftke" ) 
    371 !----------------------------------------------------------------------- 
    372 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    373 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    374 !  ediss    coef. of the Kolmogoroff dissipation   
    375 !  ebb      coef. of the surface input of tke 
    376 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    377 !  emin     minimum value of tke (m^2/s^2) 
    378 !  emin0    surface minimum value of tke (m^2/s^2) 
    379 !  nitke    number of restart iterative loops 
    380 !  ri_c     critic richardson number 
    381 !  nmxl     flag on mixing length used 
    382 !           = 0 bounded by the distance to surface and bottom 
    383 !           = 1 bounded by the local vertical scale factor 
    384 !           = 2 first vertical derivative of mixing length bounded by 1 
    385 !  npdl     flag on prandtl number 
    386 !           = 0 no vertical prandtl number (avt=avm) 
    387 !           = 1 prandtl number function of richarson number (avt=pdl*avm) 
    388 !           = 2 same as = 1 but a shapiro filter is applied on pdl 
    389 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    390 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    391 &namtke 
    392    ln_rstke = .false. 
    393    ediff =       0.1 
    394    ediss =       0.7 
    395    ebb   =      3.75 
    396    efave =        1. 
    397    emin  =     1.e-6 
    398    emin0 =     1.e-4 
    399    nitke =        50 
    400    nmxl  =         2 
    401    npdl  =         1 
    402    navb  =         0 
    403 / 
    404 !----------------------------------------------------------------------- 
    405 !       namkpp   K-Profile Parameterization dependent vertical diffusion 
    406 !                ( #ifdef "key_zdfkpp" ) 
    407 !----------------------------------------------------------------------- 
    408 !  ln_kpprimix   shear instability mixing  (default T) 
    409 !  difmiw        constant internal wave viscosity (m2/s) 
    410 !  difsiw        constant internal wave diffusivity (m2/s) 
    411 !  Riinfty       local Richardson Number limit for shear instability 
    412 !  difri         maximum shear mixing at Rig = 0    (m2/s) 
    413 !  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection 
    414 !  difcon        maximum mixing in interior convection (m2/s) 
    415 !  nave          = 0/1 flag for horizontal average on avt, avmu, avmv 
    416 !  navb          = 0/1 flag for constant or profile background avt 
    417 &namkpp 
    418    ln_kpprimix  = .true. 
    419    difmiw       =  1.e-04 
    420    difsiw       =  0.1e-04 
    421    Riinfty      =  0.8 
    422    difri        =  0.0050 
    423    bvsqcon      = -0.01e-07 
    424    difcon       =  1. 
    425    navb         =  0 
    426    nave         =  1 
    427 / 
    428 !----------------------------------------------------------------------- 
    429 !       namddm   double diffusive mixing parameterization 
    430 !----------------------------------------------------------------------- 
    431 !   avts    maximum avs for dd mixing  
    432 !   hsbfr   heat/salt buoyancy flux ratio 
    433 &namddm 
    434       avts  = 1.e-4 
    435       hsbfr = 1.6 
    436 / 
    437 !----------------------------------------------------------------------- 
    438 !       namlbc   lateral momentum boundary condition 
    439 !----------------------------------------------------------------------- 
    440 !  shlat   lateral boundary condition on velocity 
    441 !                   shlat = 0 , free slip 
    442 !               0 < shlat < 2 , partial slip 
    443 !                   shlat = 2 , no slip 
    444 !               2 < shlat     , strong slip 
    445 &namlbc 
    446    shlat  =      2. 
    447 / 
    448 !----------------------------------------------------------------------- 
    449 !       nambfr   bottom friction 
    450 !----------------------------------------------------------------------- 
    451 !  nbotfr  type of bottom friction  
    452 !                  nbotfr = 0 , no slip  
    453 !                  nbotfr = 1 , linear friction 
    454 !                  nbotfr = 2 , nonlinear friction 
    455 !                  nbotfr = 3 , free slip 
    456 !  bfri1   bottom drag coefficient (linear case) 
    457 !  bfri2   bottom drag coefficient (non linear case) 
    458 !  bfeb2   bottom turbulent kinetic energy  (m^2/s^2) 
    459 &nambfr 
    460    nbotfr =       1 
    461    bfri1  =   4.e-4 
    462    bfri2  =   1.e-3 
    463    bfeb2  =  2.5e-3 
    464 / 
    465 !----------------------------------------------------------------------- 
    466 !       nambbc   bottom temperature boundary condition 
    467 !----------------------------------------------------------------------- 
    468 !  ngeo_flux  = 0 no geothermal heat flux 
    469 !             = 1 constant geothermal heat flux 
    470 !             = 2 variable geothermal heat flux (read in geothermal_heating.nc)  
    471 !                 ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 
    472 !  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
    473 &nambbc 
    474    ngeo_flux =  2 
    475    ngeo_flux_const = 86.4e-3 
     173 
     174!!====================================================================== 
     175!!            ***  Surface Boundary Condition namelists  *** 
     176!!====================================================================== 
     177!!   namsbc        surface boundary condition 
     178!!   namsbc_ana    analytical         formulation 
     179!!   namsbc_flx    flux               formulation 
     180!!   namsbc_clio   CLIO bulk formulea formulation 
     181!!   namsbc_core   CORE bulk formulea formulation 
     182!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled") 
     183!!   namqsr        penetrative solar radiation 
     184!!   namsbc_rnf    river runoffs 
     185!!   namsbc_ssr    sea surface restoring term (for T and/or S) 
     186!!   namalb        albedo parameters 
     187!!====================================================================== 
     188 
     189!----------------------------------------------------------------------- 
     190!       namsbc       Surface Boundary Condition (surface module) 
     191!----------------------------------------------------------------------- 
     192!  nn_fsbc       frequency of surface boundary condition computation  
     193!                (= the frequency of sea-ice model call) 
     194!  ln_ana        analytical formulation (fill namsbc_ana )  
     195!  ln_flx        flux formulation       (fill namsbc_flx )  
     196!  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core)  
     197!  ln_blk_core   CORE bulk formulation  (fill namsbc_clio)  
     198!  ln_cpl        Coupled formulation    (fill namsbc_cpl )  
     199!  nn_ice        =0 no ice boundary condition   , 
     200!                =1 observed ice-cover          , 
     201!                =2 ice-model used                                      ("key_lim3" or "key_lim2) 
     202!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
     203!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr) 
     204!  nn_fwb       FreshWater Budget: =0 no control                              ,  
     205!                                  =1 annual global mean of e-p-r set to zero   , 
     206!                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step 
     207&namsbc 
     208   nn_fsbc     = 5 
     209   ln_ana      = .false. 
     210   ln_flx      = .false. 
     211   ln_blk_clio = .true. 
     212   ln_blk_core = .false. 
     213   ln_cpl      = .false. 
     214   nn_ice      = 2 
     215   ln_dm2dc    = .false. 
     216   ln_rnf      = .false. 
     217   ln_ssr      = .false. 
     218   nn_fwb      = 0 
     219/ 
     220!----------------------------------------------------------------------- 
     221!       namsbc_ana   analytical surface boundary condition 
     222!----------------------------------------------------------------------- 
     223!  nn_tau000   gently increase the stress over the first ntau_rst time-steps 
     224!  rn_utau0    uniform value for the i-stress 
     225!  rn_vtau0    uniform value for the j-stress 
     226!  rn_q0       uniform value for the total heat flux 
     227!  rn_qsr0     uniform value for the solar radiation 
     228!  rn_emp0     uniform value for the freswater budget (E-P) 
     229&namsbc_ana 
     230   nn_tau000  =   0 
     231   rn_utau0   =   0.5 
     232   rn_vtau0   =   0.e0 
     233   rn_q0      =   0.e0 
     234   rn_qsr0    =   0.e0 
     235   rn_emp0    =   0.e0 
     236/ 
     237!----------------------------------------------------------------------- 
     238!       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux") 
     239!----------------------------------------------------------------------- 
     240!  cn_dir   directory in which the model is executed 
     241! 
     242!  THE ORDER OF THE FILES MATTER: 
     243!               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     244!               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     245&namsbc_flx 
     246   cn_dir       = './' 
     247   sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0 
     248   sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0 
     249   sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0 
     250   sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0 
     251   sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0 
     252/       
     253!----------------------------------------------------------------------- 
     254!       namsbc_clio  CLIO bulk formulea 
     255!----------------------------------------------------------------------- 
     256!  cn_dir   directory in which the model is executed 
     257! 
     258!  THE ORDER OF THE FILES MATTER: 
     259!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     260!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     261&namsbc_clio 
     262   cn_dir       = './' 
     263   sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0     
     264   sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0     
     265   sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0       
     266   sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0      
     267   sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0     
     268   sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0       
     269   sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0        
     270/ 
     271!----------------------------------------------------------------------- 
     272!       namsbc_core  CORE bulk formulea 
     273!----------------------------------------------------------------------- 
     274!  cn_dir         directory in which the model is executed 
     275!  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m 
     276!  alpha_precip   multiplication factor for precipitation (total & snow) 
     277! 
     278!  THE ORDER OF THE FILES MATTER: 
     279!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting ! 
     280!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  ! 
     281&namsbc_core 
     282   cn_dir       = './' 
     283   ln_2m        = .false. 
     284   alpha_precip = 1. 
     285   sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0     
     286   sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0     
     287   sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0       
     288   sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0       
     289   sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0      
     290   sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0     
     291   sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0       
     292   sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0        
     293/ 
     294!----------------------------------------------------------------------- 
     295!       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled") 
     296!----------------------------------------------------------------------- 
     297&namsbc_cpl 
    476298/ 
    477299!----------------------------------------------------------------------- 
     
    489311/ 
    490312!----------------------------------------------------------------------- 
    491 !       namtdp   tracer newtonian damping ('key_tradmp') 
    492 !----------------------------------------------------------------------- 
    493 !  ndmp    type of damping in temperature and salinity  
    494 !          (='latitude', damping poleward of 'ndmp' degrees and function  
    495 !             of the distance-to-coast. Red and Med Seas as ndmp=-1) 
    496 !          (=-1 damping only in Med and Red Seas) 
    497 !  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
    498 !  nmldmp  type of damping in the mixed layer 
    499 !          (=0 damping throughout the water column) 
    500 !     (=1 no damping in the mixed layer defined by avt >5cm2/s ) 
    501 !     (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
    502 !  sdmp    surface time scale for internal damping (days) 
    503 !  bdmp    bottom time scale for internal damping (days) 
    504 !  hdmp    depth of transition between sdmp and bdmp (meters) 
    505 &namtdp 
    506    ndmp   =   -1 
    507    ndmpf  =    1 
    508    nmldmp =    1 
    509    sdmp   =  50. 
    510    bdmp   = 360. 
    511    hdmp   = 800. 
    512 / 
    513 !----------------------------------------------------------------------- 
    514 !       nameos   ocean physical parameters 
    515 !----------------------------------------------------------------------- 
    516 !  neos    type of equation of state and Brunt-Vaisala frequency 
    517 !          = 0, UNESCO (formulation of Jackett and McDougall (1994) 
    518 !                                         and of McDougall (1987) ) 
    519 !          = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    520 !          = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    521 !                               with rau0=1020 set in parcst routine 
    522 !  ralpha  thermal expension coefficient (linear equation of state) 
    523 !  rbeta   saline  expension coefficient (linear equation of state) 
    524 &nameos 
    525    neos   =      0 
    526    ralpha =  2.e-4 
    527    rbeta  =  0.001 
    528 / 
    529 !----------------------------------------------------------------------- 
    530 !       namsol   elliptic solver / island / free surface  
    531 !----------------------------------------------------------------------- 
    532 !  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg) 
    533 !                            (=2 successive-over-relaxation: sor) 
    534 !                            (=3 FETI: fet, all require "key_feti" defined) 
    535 !  nsol_arp  absolute/relative (0/1) precision convergence test 
    536 !  nmin      minimum of iterations for the SOR solver 
    537 !  nmax      maximum of iterations for the SOR solver 
    538 !  nmod      frequency of test for the SOR solver 
    539 !  eps       absolute precision of the solver 
    540 !  resmax    absolute precision for the SOR solver 
    541 !  sor       optimal coefficient for SOR solver 
    542 !  epsisl    absolute precision on stream function solver 
    543 !  nmisl     maximum pcg iterations for island 
    544 !  rnu       strength of the additional force used in free surface b.c. 
    545 &namsol 
    546    nsolv     =      1 
    547    nsol_arp  =      0 
    548    nmin      =    300 
    549    nmax      =    800 
    550    nmod      =     10 
    551    eps       =  1.E-6 
    552    resmax    = 1.E-10 
    553    sor       =   1.92 
    554    epsisl    = 1.e-10 
    555    nmisl     =   4000 
    556    rnu       =     1. 
    557 / 
    558 !======================================================================= 
    559 !   Diagnostics namelists 
    560 !       namtrd    dynamics and/or tracer trends 
    561 !       namgap    level mean model-data gap 
    562 !       namznl    zonal mean heat & freshwater fluxes computation 
    563 !       namspr    surface pressure in rigid-lid 
    564 !======================================================================= 
    565 !----------------------------------------------------------------------- 
    566 !       namtrd    diagnostics on dynamics and/or tracer trends 
    567 !                         ('key_trdyn' and/or 'key_trdtra') 
    568 !                 or mixed-layer trends ('key_trdmld') 
    569 !----------------------------------------------------------------------- 
    570 !  ntrd              time step frequency dynamics and tracers trends 
    571 !  nctls             control surface type in mixed-layer trends (0,1 or n<jpk) 
    572 !  ln_trdmld_restart restart for ML diagnostics 
    573 !  ucf               unit conversion factor (=1 -> /seconds | =86400. -> /day) 
    574 !  ln_trdmld_instant flag to diagnose trends of instantantaneous or mean ML T/S 
    575 &namtrd 
    576    ntrd  = 365 
    577    nctls =   0 
    578    ln_trdmld_restart = .false. 
    579    ucf   =  1. 
    580    ln_trdmld_instant = .false. 
    581 / 
    582 !----------------------------------------------------------------------- 
    583 !       namgap    level mean model-data gap ('key_diagap') 
    584 !----------------------------------------------------------------------- 
    585 !  ngap    time-step frequency of model-data gap computation 
    586 !  nprg    time-step frequency of gap print in model output 
    587 &namgap 
    588    ngap =  15 
    589    nprg =  10 
    590 / 
    591 !----------------------------------------------------------------------- 
    592 !       namznl    zonal mean heat & freshwater fluxes computation 
    593 !                 (#ifdef "key_diaznl") 
    594 !----------------------------------------------------------------------- 
    595 !  nfznl   time-step frequency of zonal mean fluxes computation 
    596 &namznl 
    597    nfznl =  15 
    598 / 
    599 !----------------------------------------------------------------------- 
    600 !       namspr  surface pressure diagnostic 
    601 !----------------------------------------------------------------------- 
    602 !  nmaxp   maximum of iterations for the solver 
    603 !  epsp    absolute precision of the solver 
    604 !  niterp  number of iteration done by the solver 
    605 &namspr 
    606    nmaxp   =   1000 
    607    epsp    =  1.e-3 
    608    niterp  =    400 
     313!       namsbc_rnf   runoffs namelist surface boundary condition 
     314!----------------------------------------------------------------------- 
     315!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv 
     316!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used 
     317!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s] 
     318!cn_dir      Root directory for location of ssr files 
     319!sn_rnf      informations about the runoff file to be read 
     320!sn_cnf      informations about the runoff mouth file to be read 
     321!  THE ORDER OF THE FILES MATTER: 
     322!               !   file name        ! frequency (hours) ! variable   ! time interpolation ! 
     323!               !                    !   (if <0  months) !   name     !     (logical)      ! 
     324&namsbc_rnf 
     325   nn_runoff  = 2 
     326   rn_hrnf    = 0.e0 
     327   rn_avt_rnf = 1.e-3 
     328   cn_dir     = './' 
     329   sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0      
     330   sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0               
     331/ 
     332!----------------------------------------------------------------------- 
     333!       namsbc_ssr   surface boundary condition : sea surface restoring 
     334!----------------------------------------------------------------------- 
     335!  nn_sstr   =0/1   add a retroaction term in the surface heat flux 
     336!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux 
     337!             (=1, salt flux, concentration/dillution emps only) 
     338!             (=2, volume flux, both emp and emps are updated)  
     339!  dqdt       magnitude of the retroaction on temperature [W/m2/K] 
     340!  deds       magnitude of the damping on salinity [mm/day] 
     341!  THE ORDER OF THE FILES MATTER: 
     342!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting ! 
     343!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  ! 
     344&namsbc_ssr 
     345   nn_sstr      = 0 
     346   nn_sssr      = 0 
     347   dqdt         = -40. 
     348   deds         =  27.7 
     349   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0 
     350   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0 
     351/       
     352!----------------------------------------------------------------------- 
     353!       namalb   albedo parameters 
     354!----------------------------------------------------------------------- 
     355!  cgren    correction of the snow or ice albedo to take into account 
     356!  albice   albedo of melting ice in the arctic and antarctic 
     357!  alphd    coefficients for linear interpolation used to compute albedo 
     358!           between two extremes values (Pyane, 1972) 
     359!  alphc     "                                         " 
     360!  alphdi    "                                         " 
     361&namalb 
     362   cgren    =      0.06 
     363   albice   =      0.5 
     364   alphd    =      0.80 
     365   alphc    =      0.65 
     366   alphdi   =      0.72 
     367/ 
     368 
     369!!====================================================================== 
     370!!               ***  Lateral boundary condition  *** 
     371!!====================================================================== 
     372!!   namlbc        lateral momentum boundary condition 
     373!!   namcla        cross land advection 
     374!!   namobc        open boundaries parameters                           ("key_obc") 
     375!!====================================================================== 
     376 
     377!----------------------------------------------------------------------- 
     378!       namlbc   lateral momentum boundary condition 
     379!----------------------------------------------------------------------- 
     380!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip 
     381!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip 
     382&namlbc 
     383   shlat  =      2. 
     384/ 
     385!----------------------------------------------------------------------- 
     386!       namcla   cross land advection 
     387!----------------------------------------------------------------------- 
     388!  n_cla   advection between 2 ocean pts separates by land  
     389&namcla 
     390   n_cla   = 1 
    609391/ 
    610392!----------------------------------------------------------------------- 
     
    638420    ln_vol_cst  = .false. 
    639421/ 
    640 !----------------------------------------------------------------------- 
    641 !       namflo    float parameters (#ifdef key_float) 
     422 
     423!!====================================================================== 
     424!!                 ***  Bottom boundary condition  *** 
     425!!====================================================================== 
     426!!   nambfr        bottom friction 
     427!!   nambbc        bottom temperature boundary condition                ("key_trabbc") 
     428!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv") 
     429!!====================================================================== 
     430!----------------------------------------------------------------------- 
     431!       nambfr   bottom friction 
     432!----------------------------------------------------------------------- 
     433!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction 
     434!                                     = 3 : free slip  ;  = 2 : non linear friction 
     435!  bfri1    bottom drag coefficient (linear case) 
     436!  bfri2    bottom drag coefficient (non linear case) 
     437!  bfeb2    bottom turbulent kinetic energy  (m2/s2) 
     438&nambfr 
     439   nbotfr =       1 
     440   bfri1  =   4.e-4 
     441   bfri2  =   1.e-3 
     442   bfeb2  =  2.5e-3 
     443/ 
     444!----------------------------------------------------------------------- 
     445!       nambbc   bottom temperature boundary condition 
     446!----------------------------------------------------------------------- 
     447!  ngeo_flux  = 0 no geothermal heat flux 
     448!             = 1 constant geothermal heat flux 
     449!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 
     450!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
     451&nambbc 
     452   ngeo_flux =  2 
     453   ngeo_flux_const = 86.4e-3 
     454/ 
     455!----------------------------------------------------------------------- 
     456!       nambbl   bottom boundary layer scheme 
     457!----------------------------------------------------------------------- 
     458!                                !  diffusive bbl                       ("key_trabbl") 
     459!                                !  advective bbl                       ("key_trabbl_adv") 
     460!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s) 
     461&nambbl 
     462   atrbbl = 10000. 
     463/ 
     464 
     465!!====================================================================== 
     466!!                        Tracer (T & S ) namelists 
     467!!====================================================================== 
     468!!   nameos        equation of state 
     469!!   nam_traadv    advection scheme 
     470!!   nam_traldf    lateral diffusion scheme 
     471!!   namtdp        tracer newtonian damping                             ("key_tradmp") 
     472!!====================================================================== 
     473 
     474!----------------------------------------------------------------------- 
     475!       nameos   ocean physical parameters 
     476!----------------------------------------------------------------------- 
     477!  neos    type of equation of state and Brunt-Vaisala frequency 
     478!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
     479!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
     480!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
     481!  ralpha  thermal expension coefficient (neos= 1 or 2) 
     482!  rbeta   saline  expension coefficient (neos= 2) 
     483&nameos 
     484   neos   =      0 
     485   ralpha =  2.e-4 
     486   rbeta  =  0.001 
     487/ 
     488!----------------------------------------------------------------------- 
     489!       nam_traadv   advection scheme for tracer 
     490!----------------------------------------------------------------------- 
     491!  ln_traadv_cen2     2nd order centered scheme          (default T) 
     492!  ln_traadv_tvd      TVD scheme                         (default F) 
     493!  ln_traadv_muscl    MUSCL scheme                       (default F) 
     494!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F) 
     495!  ln_traadv_ubs      UBS scheme                         (default F) 
     496!  ln_traadv_qck      QUICKEST scheme                    (default F) 
     497&nam_traadv 
     498   ln_traadv_cen2   =  .true. 
     499   ln_traadv_tvd    =  .false. 
     500   ln_traadv_muscl  =  .false. 
     501   ln_traadv_muscl2 =  .false. 
     502   ln_traadv_ubs    =  .false. 
     503   ln_traadv_qck    =  .false. 
     504/ 
     505!----------------------------------------------------------------------- 
     506!       nam_traldf   lateral diffusion scheme for tracer 
     507!----------------------------------------------------------------------- 
     508!  Type of the operator :  
     509!     ln_traldf_lap    laplacian operator          (default T) 
     510!     ln_traldf_bilap  bilaplacian operator        (default F) 
     511!  Direction of action  :  
     512!     ln_traldf_level  iso-level                   (default F) 
     513!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T) 
     514!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp") 
     515!  Coefficient 
     516!     aht0    horizontal eddy diffusivity for tracers (m2/s) 
     517!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
     518!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv") 
     519&nam_traldf 
     520   ln_traldf_lap    =  .true. 
     521   ln_traldf_bilap  =  .false. 
     522   ln_traldf_level  =  .false. 
     523   ln_traldf_hor    =  .false. 
     524   ln_traldf_iso    =  .true. 
     525   aht0    =  2000. 
     526   ahtb0   =     0. 
     527   aeiv0   =  2000. 
     528/ 
     529!----------------------------------------------------------------------- 
     530!       namtdp   tracer newtonian damping ('key_tradmp') 
     531!----------------------------------------------------------------------- 
     532!  ndmp    type of damping in temperature and salinity  
     533!          !   ='latitude', damping poleward of 'ndmp' degrees and function  
     534!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1 
     535!          !   =-1          damping only in Med and Red Seas 
     536!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array) 
     537!  nmldmp  type of damping in the mixed layer 
     538!          !   =0 damping throughout the water column 
     539!          !   =1 no damping in the mixed layer defined by avt >5cm2/s  
     540!          !   =2 no damping in the mixed layer defined rho<rho(surf)+.01  
     541!  sdmp    surface time scale for internal damping (days) 
     542!  bdmp    bottom time scale for internal damping (days) 
     543!  hdmp    depth of transition between sdmp and bdmp (meters) 
     544&namtdp 
     545   ndmp   =   -1 
     546   ndmpf  =    1 
     547   nmldmp =    1 
     548   sdmp   =  50. 
     549   bdmp   = 360. 
     550   hdmp   = 800. 
     551/ 
     552 
     553!!====================================================================== 
     554!!                      ***  Dynamics namelists  *** 
     555!!====================================================================== 
     556!!   nam_dynadv    formulation of the momentum advection 
     557!!   nam_dynvor    advection scheme 
     558!!   nam_dynhpg    hydrostatic pressure gradient 
     559!!   namflg        hydrostatic pressure gradient time stepping 
     560!!   nam_dynspg    surface pressure gradient                            (CPP key only) 
     561!!   nam_dynldf    lateral diffusion scheme 
     562!!====================================================================== 
     563 
     564!----------------------------------------------------------------------- 
     565!       nam_dynadv   formulation of the momentum advection 
     566!----------------------------------------------------------------------- 
     567!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T) 
     568!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T) 
     569!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F)