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

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

Ignore:
Timestamp:
2016-03-29T11:24:48+02:00 (8 years ago)
Author:
timgraham
Message:

First attempt at upgrading branch to the head of the trunk. This should include all of the simplification branch from the merge in Dec 2015.

Location:
branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG
Files:
39 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/AMM12/EXP00/iodef.xml

    r6401 r6404  
    2525      <file_group id="1ts" output_freq="1ts"  output_level="10" enabled=".TRUE."/> <!-- 1 time step files --> 
    2626 
    27       <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."/> <!-- 1h files --> 
     27<!--old      <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."/> old --> <!-- 1h files -->  
     28   <!-- TMB files --> 
     29      <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."  >                       
     30 
     31   <file id="file9" name_suffix="_shelftmb_grid_T" description="ocean T grid variables" enabled=".TRUE." > 
     32     <field field_ref="top_temp"           name="votemper_top"   operation="instant" enabled=".TRUE." /> 
     33     <field field_ref="mid_temp"           name="votemper_mid"   operation="instant" enabled=".TRUE." /> 
     34     <field field_ref="bot_temp"           name="votemper_bot"   operation="instant" enabled=".TRUE." /> 
     35     <field field_ref="sshnmasked"         name="sossheig"       operation="instant" enabled=".TRUE." />  
     36     <field field_ref="top_sal"            name="vosaline_top"   operation="instant" enabled=".TRUE." /> 
     37     <field field_ref="mid_sal"            name="vosaline_mid"   operation="instant" enabled=".TRUE." /> 
     38     <field field_ref="bot_sal"            name="vosaline_bot"   operation="instant" enabled=".TRUE." /> 
     39   </file> 
     40 
     41   <file id="file10" name_suffix="_shelftmb_grid_U" description="TMB ocean U grid variables" enabled=".TRUE." > 
     42     <field field_ref="top_u"         name="vozocrtx_top"   operation="instant" /> 
     43     <field field_ref="mid_u"         name="vozocrtx_mid"   operation="instant" /> 
     44     <field field_ref="bot_u"         name="vozocrtx_bot"   operation="instant" /> 
     45     <field field_ref="baro_u"        name="vobtcrtx"       operation="instant" /> 
     46   </file> 
     47 
     48   <file id="file11" name_suffix="_shelftmb_grid_V" description="TMB ocean V grid variables" enabled=".TRUE." > 
     49     <field field_ref="top_v"         name="vomecrty_top"   operation="instant" /> 
     50     <field field_ref="mid_v"         name="vomecrty_mid"   operation="instant" /> 
     51     <field field_ref="bot_v"         name="vomecrty_bot"   operation="instant" /> 
     52     <field field_ref="baro_v"        name="vobtcrty"       operation="instant" /> 
     53   </file> 
     54 
     55      </file_group> 
     56 
    2857      <file_group id="2h" output_freq="2h"  output_level="10" enabled=".TRUE."/> <!-- 2h files --> 
    2958      <file_group id="3h" output_freq="3h"  output_level="10" enabled=".TRUE."/> <!-- 3h files -->      
    3059      <file_group id="4h" output_freq="4h"  output_level="10" enabled=".TRUE."/> <!-- 4h files --> 
    3160      <file_group id="6h" output_freq="6h"  output_level="10" enabled=".TRUE."/> <!-- 6h files --> 
     61 
     62      <file_group id="25h_mean" output_freq="1d" output_level="10" enabled=".TRUE."> 
     63 
     64 
     65     <file id="file12" name_suffix="_25hourm_grid_T" description="ocean T grid variables, 25h meaned" enabled=".TRUE." > 
     66     <field field_ref="temper25h"   name="votemper"  operation="instant" enabled=".TRUE."/> 
     67     <field field_ref="tempis25h"   name="votempis"  operation="instant" enabled=".TRUE."/> 
     68     <field field_ref="salin25h"    name="vosaline"  operation="instant" enabled=".TRUE."/> 
     69     <field field_ref="ssh25h"      name="sossheig"  operation="instant" enabled=".TRUE."/> 
     70    </file> 
     71 
     72     <file id="file13" name_suffix="_25hourm_grid_U" description="ocean U grid variables, 25h meaned" enabled=".TRUE." > 
     73     <field field_ref="vozocrtx25h" name="vozocrtx"  operation="instant" enabled=".TRUE."/> 
     74    </file> 
     75 
     76     <file id="file14" name_suffix="_25hourm_grid_V" description="ocean V grid variables, 25h meaned" enabled=".TRUE." > 
     77     <field field_ref="vomecrty25h" name="vomecrty"  operation="instant" enabled=".TRUE."/> 
     78    </file> 
     79 
     80     <file id="file15" name_suffix="_25hourm_grid_W" description="ocean W grid variables, 25h meaned" enabled=".TRUE." > 
     81     <field field_ref="vomecrtz25h"          name="vomerctz"  operation="instant" enabled=".TRUE."/> 
     82     <field field_ref="woce"                 name="vovecrtz"  operation="instant" enabled=".TRUE."/> 
     83     <field field_ref="avt25h"               name="votkeavt"  operation="instant" enabled=".TRUE."/> 
     84     <field field_ref="avm25h"               name="votkeavm"  operation="instant" enabled=".TRUE."/> 
     85     <field field_ref="tke25h"               name="votke"  operation="instant" enabled=".TRUE."/> 
     86     <field field_ref="mxln25h"              name="mxln"   operation="instant" enabled=".TRUE."/> 
     87    </file> 
     88      </file_group> 
    3289      
    3390      <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."> <!-- 1d files --> 
    … …  
    98155      <file_group id="5y"  output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files --> 
    99156      <file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files --> 
     157 
    100158 
    101159   </file_definition> 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r6401 r6404  
    3434&namzgr        !   vertical coordinate 
    3535!----------------------------------------------------------------------- 
    36    ln_zps      = .false.   !  z-coordinate - partial steps   (T/F) 
    37    ln_sco      = .true.    !  s- or hybrid z-s-coordinate    (T/F) 
     36   ln_sco      = .true.    !  s- or hybrid z-s-coordinate 
    3837/ 
    3938!----------------------------------------------------------------------- 
    … …  
    5049!----------------------------------------------------------------------- 
    5150   rn_rdt      =   600.    !  time step for the dynamics (and tracer if nn_acc=0) 
    52    rn_rdtmin   =   600.          !  minimum time step on tracers (used if nn_acc=1) 
    53    rn_rdtmax   =   600.          !  maximum time step on tracers (used if nn_acc=1) 
    54    rn_rdth     =   600.          !  depth variation of tracer time step  (used if nn_acc=1) 
    5551   ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
    5652   ppgphi0     =  999999.0             ! latitude  of first raw and column T-point (jphgr_msh = 1) 
    … …  
    7268/ 
    7369!----------------------------------------------------------------------- 
    74 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    75 !----------------------------------------------------------------------- 
    76    ln_bt_nn_auto =    .FALSE.          !  Set nn_baro automatically to be just below 
    77                                        !  a user defined maximum courant number (rn_bt_cmax) 
    78    nn_baro       =    30               !  Number of iterations of barotropic mode 
    79 / 
    80 !----------------------------------------------------------------------- 
    8170&namcrs        !   Grid coarsening for dynamics output and/or 
    8271               !   passive tracer coarsened online simulations 
    … …  
    10089   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
    10190   ln_apr_dyn  = .false.    !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
     91   ln_traqsr   = .false.   !  Light penetration (T) or not (F) 
    10292 
    10393/ 
    … …  
    137127&namtra_qsr    !   penetrative solar radiation 
    138128!----------------------------------------------------------------------- 
    139    ln_traqsr   = .false.   !  Light penetration (T) or not (F) 
    140129   nn_chldta   =      0    !  RGB : Chl data (=1) or cst value (=0) 
    141130/ 
    … …  
    194183   rn_shlat    =     0     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    195184                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    196 / 
    197 !----------------------------------------------------------------------- 
    198 &namcla        !   cross land advection 
    199 !----------------------------------------------------------------------- 
    200 / 
    201 !----------------------------------------------------------------------- 
    202 &namobc        !   open boundaries parameters                           ("key_obc") 
    203 !----------------------------------------------------------------------- 
    204185/ 
    205186!----------------------------------------------------------------------- 
    … …  
    265246/ 
    266247!----------------------------------------------------------------------- 
    267 &nambbc        !   bottom temperature boundary condition 
    268 !----------------------------------------------------------------------- 
    269    ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
     248&nambbc        !   bottom temperature boundary condition                (default: NO) 
     249!----------------------------------------------------------------------- 
    270250/ 
    271251!----------------------------------------------------------------------- 
    … …  
    281261&namtra_adv    !   advection scheme for tracer 
    282262!----------------------------------------------------------------------- 
     263   ln_traadv_fct =  .true.   !  FCT scheme 
     264      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     265      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     266      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     267      !                             !        (number of sub-timestep = nn_fct_zts) 
    283268/ 
    284269!----------------------------------------------------------------------- 
    … …  
    289274&namtra_ldf    !   lateral diffusion scheme for tracers 
    290275!---------------------------------------------------------------------------------- 
    291    ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    292    ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    293    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    294    rn_aht_0         =    50.    !  horizontal eddy diffusivity for tracers [m2/s] 
     276   !                       !  Operator type: 
     277   ln_traldf_lap   =  .true.   !    laplacian operator 
     278   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     279   !                       !  Direction of action: 
     280   ln_traldf_lev   =  .false.  !  iso-level 
     281   ln_traldf_hor   =  .true.   !  horizontal (geopotential) 
     282   ln_traldf_iso   =  .false.  !  iso-neutral 
     283   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     284   ! 
     285   !                       !  iso-neutral options:         
     286   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     287   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     288   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     289   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     290   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     291   ! 
     292   !                       !  Coefficients: 
     293   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     294   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     295   !                                !   =  0           constant  
     296   !                                !   = 10 F(k)      =ldf_c1d  
     297   !                                !   = 20 F(i,j)    =ldf_c2d  
     298   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     299   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     300   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     301   rn_aht_0        =   50.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     302   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     303/ 
     304!---------------------------------------------------------------------------------- 
     305&namtra_ldfeiv !   eddy induced velocity param. 
     306!---------------------------------------------------------------------------------- 
     307   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    295308/ 
    296309!----------------------------------------------------------------------- 
    … …  
    306319&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    307320!----------------------------------------------------------------------- 
     321   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     322   ln_dynvor_ens = .false. !  energy conserving scheme 
     323   ln_dynvor_mix = .false. !  mixed scheme 
     324   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     325      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    308326/ 
    309327!----------------------------------------------------------------------- 
    … …  
    314332/ 
    315333!----------------------------------------------------------------------- 
    316 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    317 !----------------------------------------------------------------------- 
    318 !                          !  explicit free surface                     ("key_dynspg_exp") 
    319 !                          !  filtered free surface                     ("key_dynspg_flt") 
    320 !                          !  split-explicit free surface               ("key_dynspg_ts") 
    321  
     334&namdyn_spg    !   surface pressure gradient  
     335!----------------------------------------------------------------------- 
     336   ln_dynspg_ts = .true.    ! split-explicit free surface 
     337      ln_bt_auto    = .false.  ! Number of sub-step defined from: 
     338         nn_baro       = 30       ! =F : the number of sub-step in rn_rdt seconds 
     339/ 
    322340!----------------------------------------------------------------------- 
    323341&namdyn_ldf    !   lateral diffusion on momentum 
    324342!----------------------------------------------------------------------- 
    325343   !                       !  Type of the operator : 
    326    ln_dynldf_bilap  =  .true.   !  bilaplacian operator 
    327    ln_dynldf_lap    =  .false.  !  bilaplacian operator 
     344   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     345   ln_dynldf_lap =  .false.    !    laplacian operator 
     346   ln_dynldf_blp =  .true.     !  bilaplacian operator 
    328347   !                       !  Direction of action  : 
    329    ln_dynldf_level  =  .true.   !  iso-level 
    330    ln_dynldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    331                            !  Coefficient 
    332    rn_ahm_0_lap     = 60.0      !  horizontal laplacian eddy viscosity   [m2/s] 
    333    rn_ahm_0_blp     = -1.0e+10  !  horizontal bilaplacian eddy viscosity [m4/s] 
     348   ln_dynldf_lev =  .true.     !  iso-level 
     349   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     350   ln_dynldf_iso =  .false.    !  iso-neutral 
     351   !                       !  Coefficient 
     352   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     353   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     354   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     355   !                                !  =  0  constant  
     356   !                                !  = 10  F(k)=c1d 
     357   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     358   !                                !  = 30  F(i,j,k)=c2d*c1d 
     359   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     360   rn_ahm_0      =     60.     !  horizontal laplacian eddy viscosity   [m2/s] 
     361   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     362   rn_bhm_0      = 1.0e+10     !  horizontal bilaplacian eddy viscosity [m4/s] 
     363   ! 
     364   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    334365/ 
    335366!----------------------------------------------------------------------- 
    … …  
    349380!----------------------------------------------------------------------- 
    350381/ 
    351 !------------------------------------------------------------------------ 
    352 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    353 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    354 / 
    355382!----------------------------------------------------------------------- 
    356383&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    … …  
    369396/ 
    370397!----------------------------------------------------------------------- 
    371 &namsol        !   elliptic solver / island / free surface 
    372 !----------------------------------------------------------------------- 
    373 / 
    374 !----------------------------------------------------------------------- 
    375398&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    376399!----------------------------------------------------------------------- 
    … …  
    421444/ 
    422445!----------------------------------------------------------------------- 
    423 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    424 !----------------------------------------------------------------------- 
    425 / 
     446&nam_diatmb  !  Top Middle Bottom Output 
     447!----------------------------------------------------------------------- 
     448   ln_diatmb  = .true.    !  Choose Top Middle and Bottom output or not 
     449/ 
     450!----------------------------------------------------------------------- 
     451&nam_dia25h  !  25h Mean Output 
     452!----------------------------------------------------------------------- 
     453   ln_dia25h  = .true.    ! Choose 25h mean output or not 
     454/ 
     455 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/AMM12/cpp_AMM12.fcm

    r6401 r6404  
    1  bld::tool::fppkeys  key_bdy key_tide key_dynspg_ts key_ldfslp  key_zdfgls  key_vvl key_diainstant key_mpp_mpi key_iomput 
     1 bld::tool::fppkeys  key_bdy key_tide key_zdfgls key_diainstant key_mpp_mpi key_iomput 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r6401 r6404  
    3030&namzgr        !   vertical coordinate 
    3131!----------------------------------------------------------------------- 
     32   ln_zps      = .true.    !  z-coordinate - partial steps 
    3233/ 
    3334!----------------------------------------------------------------------- 
    … …  
    4142   rn_bathy    = 4198.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
    4243   nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
    43    rn_rdt      =  360.     !  time step for the dynamics (and tracer if nn_acc=0) 
    44    rn_rdtmin   =  360.     !  minimum time step on tracers (used if nn_acc=1) 
    45    rn_rdtmax   =  360.     !  maximum time step on tracers (used if nn_acc=1) 
     44   rn_rdt      =  360.     !  time step for the dynamics  
    4645   jphgr_msh   =       1                 !  type of horizontal mesh 
    4746   ppglam0     =    -150.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
    … …  
    6463/ 
    6564!----------------------------------------------------------------------- 
    66 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    67 !----------------------------------------------------------------------- 
    68 / 
    69 !----------------------------------------------------------------------- 
    7065&namcrs        !   Grid coarsening for dynamics output and/or 
    7166               !   passive tracer coarsened online simulations 
    … …  
    141136&namtra_qsr    !   penetrative solar radiation 
    142137!----------------------------------------------------------------------- 
     138!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     139!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     140   sn_chl      ='chlorophyll_PAPASTATION', -1    , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    143141/ 
    144142!----------------------------------------------------------------------- 
    … …  
    172170/ 
    173171!----------------------------------------------------------------------- 
    174 &namcla        !   cross land advection 
    175 !----------------------------------------------------------------------- 
    176 / 
    177 !----------------------------------------------------------------------- 
    178 &namobc        !   open boundaries parameters                           ("key_obc") 
    179 !----------------------------------------------------------------------- 
    180 / 
    181 !----------------------------------------------------------------------- 
    182172&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    183173!----------------------------------------------------------------------- 
    … …  
    205195/ 
    206196!----------------------------------------------------------------------- 
    207 &nambbc        !   bottom temperature boundary condition 
    208 !----------------------------------------------------------------------- 
    209    ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    210    nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux 
     197&nambbc        !   bottom temperature boundary condition                (default: NO) 
     198!----------------------------------------------------------------------- 
    211199/ 
    212200!----------------------------------------------------------------------- 
    … …  
    226214&namtra_adv    !   advection scheme for tracer 
    227215!----------------------------------------------------------------------- 
     216! C1D : no advection scheme  
    228217/ 
    229218!----------------------------------------------------------------------- 
    … …  
    231220!----------------------------------------------------------------------- 
    232221/ 
    233 !---------------------------------------------------------------------------------- 
     222!----------------------------------------------------------------------- 
    234223&namtra_ldf    !   lateral diffusion scheme for tracers 
    235 !---------------------------------------------------------------------------------- 
    236 !---------------------------------------------------------------------------------- 
    237    ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    238    ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    239    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    240    rn_aht_0         =     0.    !  horizontal eddy diffusivity for tracers [m2/s] 
     224!----------------------------------------------------------------------- 
     225! C1D : no lateral diffusion   
     226/ 
     227!----------------------------------------------------------------------- 
     228&namtra_ldfeiv !   eddy induced velocity param. 
     229!----------------------------------------------------------------------- 
     230! C1D : no eiv   
    241231/ 
    242232!----------------------------------------------------------------------- 
    … …  
    248238&namdyn_adv    !   formulation of the momentum advection 
    249239!----------------------------------------------------------------------- 
     240! C1D : no advection scheme  
    250241/ 
    251242!----------------------------------------------------------------------- 
    … …  
    266257/ 
    267258!----------------------------------------------------------------------- 
    268 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    269 !----------------------------------------------------------------------- 
     259&namdyn_spg    !   surface pressure gradient                            (default: NO spg) 
     260!----------------------------------------------------------------------- 
     261/ 
    270262!----------------------------------------------------------------------- 
    271263&namdyn_ldf    !   lateral diffusion on momentum 
    … …  
    286278!----------------------------------------------------------------------- 
    287279/ 
    288 !------------------------------------------------------------------------ 
    289 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    291 / 
    292280!----------------------------------------------------------------------- 
    293281&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    … …  
    302290!----------------------------------------------------------------------- 
    303291   ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    304 / 
    305 !----------------------------------------------------------------------- 
    306 &namsol        !   elliptic solver / island / free surface 
    307 !----------------------------------------------------------------------- 
    308    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    309    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    310    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    311292/ 
    312293!----------------------------------------------------------------------- 
    … …  
    358339!----------------------------------------------------------------------- 
    359340/ 
    360 !----------------------------------------------------------------------- 
    361 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    362 !----------------------------------------------------------------------- 
    363    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    364 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r6401 r6404  
    3232&namzgr        !   vertical coordinate 
    3333!----------------------------------------------------------------------- 
    34    ln_zco      = .true.    !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
    35    ln_zps      = .false.   !  z-coordinate - partial steps   (T/F) 
     34   ln_zco      = .true.    !  z-coordinate - full    steps 
     35   ln_linssh   = .true.    !  linear free surface 
    3636/ 
    3737!----------------------------------------------------------------------- 
    … …  
    4343!----------------------------------------------------------------------- 
    4444   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    45    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    46    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    47    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
     45   rn_rdt      = 7200.     !  time step for the dynamics  
    4846   jphgr_msh   =       5                 !  type of horizontal mesh 
    4947   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
    … …  
    6664/ 
    6765!----------------------------------------------------------------------- 
    68 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    69 !----------------------------------------------------------------------- 
    70 / 
    71 !----------------------------------------------------------------------- 
    7266&namcrs        !   Grid coarsening for dynamics output and/or 
    7367               !   passive tracer coarsened online simulations 
    … …  
    154148!----------------------------------------------------------------------- 
    155149   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    156 / 
    157 !----------------------------------------------------------------------- 
    158 &namcla        !   cross land advection 
    159 !----------------------------------------------------------------------- 
    160 / 
    161 !----------------------------------------------------------------------- 
    162 &namobc        !   open boundaries parameters                           ("key_obc") 
    163 !----------------------------------------------------------------------- 
    164150/ 
    165151!----------------------------------------------------------------------- 
    … …  
    223209&namtra_adv    !   advection scheme for tracer 
    224210!----------------------------------------------------------------------- 
    225    ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
    226    ln_traadv_tvd    =  .true.    !  TVD scheme 
    227    ln_traadv_muscl  =  .false.   !  MUSCL scheme 
    228    ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
    229    ln_traadv_ubs    =  .false.   !  UBS scheme 
    230    ln_traadv_qck    =  .false.   !  QUICKEST scheme 
    231    ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
     211   ln_traadv_fct =  .true.   !  FCT scheme 
     212      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     213      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     214      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     215      !                             !        (number of sub-timestep = nn_fct_zts) 
    232216/ 
    233217!----------------------------------------------------------------------- 
    … …  
    238222&namtra_ldf    !   lateral diffusion scheme for tracers 
    239223!---------------------------------------------------------------------------------- 
    240    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    241    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     224   !                       !  Operator type: 
     225   ln_traldf_lap   =  .true.   !    laplacian operator 
     226   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     227   !                       !  Direction of action: 
     228   ln_traldf_lev   =  .false.  !  iso-level 
     229   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     230   ln_traldf_iso   =  .true.   !  iso-neutral 
     231   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     232   ! 
     233   !                       !  iso-neutral options:         
     234   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
     235   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     236   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     237   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     238   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     239   ! 
     240   !                       !  Coefficients: 
     241   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     242   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     243   !                                !   =  0           constant  
     244   !                                !   = 10 F(k)      =ldf_c1d  
     245   !                                !   = 20 F(i,j)    =ldf_c2d  
     246   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     247   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     248   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     249   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     250   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     251/ 
     252!---------------------------------------------------------------------------------- 
     253&namtra_ldfeiv !   eddy induced velocity param. 
     254!---------------------------------------------------------------------------------- 
     255   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    242256/ 
    243257!----------------------------------------------------------------------- 
    … …  
    253267&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    254268!----------------------------------------------------------------------- 
    255    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    256    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
     269   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
     270   ln_dynvor_ens = .false. !  energy conserving scheme 
     271   ln_dynvor_mix = .false. !  mixed scheme 
    257272   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     273      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    258274/ 
    259275!----------------------------------------------------------------------- 
    … …  
    264280/ 
    265281!----------------------------------------------------------------------- 
    266 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    267 !----------------------------------------------------------------------- 
    268  
     282&namdyn_spg    !   surface pressure gradient  
     283!----------------------------------------------------------------------- 
     284   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     285/ 
    269286!----------------------------------------------------------------------- 
    270287&namdyn_ldf    !   lateral diffusion on momentum 
    271288!----------------------------------------------------------------------- 
    272    rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
     289   !                       !  Type of the operator : 
     290   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     291   ln_dynldf_lap =  .true.     !    laplacian operator 
     292   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     293   !                       !  Direction of action  : 
     294   ln_dynldf_lev =  .true.     !  iso-level 
     295   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     296   ln_dynldf_iso =  .false.    !  iso-neutral 
     297   !                       !  Coefficient 
     298   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     299   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     300   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     301   !                                !  =  0  constant  
     302   !                                !  = 10  F(k)=c1d 
     303   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     304   !                                !  = 30  F(i,j,k)=c2d*c1d 
     305   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     306   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     307   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     308   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
    273309/ 
    274310!----------------------------------------------------------------------- 
    … …  
    286322   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    287323/ 
    288 !------------------------------------------------------------------------ 
    289 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    291 / 
    292324!----------------------------------------------------------------------- 
    293325&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    … …  
    302334!----------------------------------------------------------------------- 
    303335   ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    304 / 
    305 !----------------------------------------------------------------------- 
    306 &namsol        !   elliptic solver / island / free surface 
    307 !----------------------------------------------------------------------- 
    308    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    309    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    310    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    311336/ 
    312337!----------------------------------------------------------------------- 
    … …  
    374399!----------------------------------------------------------------------- 
    375400/ 
    376 !----------------------------------------------------------------------- 
    377 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    378 !----------------------------------------------------------------------- 
    379    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    380 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm

    r6401 r6404  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi key_nosignedzero 
     1 bld::tool::fppkeys key_zdftke key_iomput key_mpp_mpi 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r6401 r6404  
    3535&namzgr        !   vertical coordinate 
    3636!----------------------------------------------------------------------- 
    37    ln_zco      = .true.    !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
    38    ln_zps      = .false.   !  z-coordinate - partial steps   (T/F) 
     37   ln_zco      = .true.    !  z-coordinate - full    steps 
     38   ln_linssh   = .true.   !  linear free surface 
    3939/ 
    4040!----------------------------------------------------------------------- 
    … …  
    4646!----------------------------------------------------------------------- 
    4747   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    48    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    49    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    50    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
     48   rn_rdt      = 7200.     !  time step for the dynamics  
    5149   jphgr_msh   =       5                 !  type of horizontal mesh 
    5250   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
    … …  
    6967/ 
    7068!----------------------------------------------------------------------- 
    71 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    72 !----------------------------------------------------------------------- 
    73 / 
    74 !----------------------------------------------------------------------- 
    7569&namcrs        !   Grid coarsening for dynamics output and/or 
    7670               !   passive tracer coarsened online simulations 
    … …  
    161155/ 
    162156!----------------------------------------------------------------------- 
    163 &namcla        !   cross land advection 
    164 !----------------------------------------------------------------------- 
    165 / 
    166 !----------------------------------------------------------------------- 
    167 &namobc        !   open boundaries parameters                           ("key_obc") 
    168 !----------------------------------------------------------------------- 
    169 / 
    170 !----------------------------------------------------------------------- 
    171157&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    172158!----------------------------------------------------------------------- 
    … …  
    194180/ 
    195181!----------------------------------------------------------------------- 
    196 &nambbc        !   bottom temperature boundary condition 
    197 !----------------------------------------------------------------------- 
    198    ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    199    nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux 
     182&nambbc        !   bottom temperature boundary condition                (default: NO) 
     183!----------------------------------------------------------------------- 
    200184/ 
    201185!----------------------------------------------------------------------- 
    … …  
    228212&namtra_adv    !   advection scheme for tracer 
    229213!----------------------------------------------------------------------- 
    230    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     214   ln_traadv_fct =  .true.   !  FCT scheme 
     215      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     216      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     217      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     218      !                             !        (number of sub-timestep = nn_fct_zts) 
    231219/ 
    232220!----------------------------------------------------------------------- 
    … …  
    237225&namtra_ldf    !   lateral diffusion scheme for tracers 
    238226!---------------------------------------------------------------------------------- 
    239    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    240    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     227   !                       !  Operator type: 
     228   ln_traldf_lap   =  .true.   !    laplacian operator 
     229   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     230   !                       !  Direction of action: 
     231   ln_traldf_lev   =  .false.  !  iso-level 
     232   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     233   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     234   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     235   ! 
     236   !                       !  iso-neutral options:         
     237   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     238   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     239   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     240   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     241   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     242   ! 
     243   !                       !  Coefficients: 
     244   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     245   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     246   !                                !   =  0           constant  
     247   !                                !   = 10 F(k)      =ldf_c1d  
     248   !                                !   = 20 F(i,j)    =ldf_c2d  
     249   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     250   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     251   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     252   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     253   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     254/ 
     255!---------------------------------------------------------------------------------- 
     256&namtra_ldfeiv !   eddy induced velocity param. 
     257!---------------------------------------------------------------------------------- 
     258   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    241259/ 
    242260!----------------------------------------------------------------------- 
    … …  
    252270&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    253271!----------------------------------------------------------------------- 
    254    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    255    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
    256    ln_dynvor_een = .false. !  energy & enstrophy scheme 
     272   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
    257273/ 
    258274!----------------------------------------------------------------------- 
    … …  
    263279/ 
    264280!----------------------------------------------------------------------- 
    265 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    266 !----------------------------------------------------------------------- 
     281&namdyn_spg    !   surface pressure gradient 
     282!----------------------------------------------------------------------- 
     283   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     284/ 
    267285!----------------------------------------------------------------------- 
    268286&namdyn_ldf    !   lateral diffusion on momentum 
    269287!----------------------------------------------------------------------- 
     288   !                       !  Type of the operator : 
     289   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     290   ln_dynldf_lap =  .true.     !    laplacian operator 
     291   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     292   !                       !  Direction of action  : 
     293   ln_dynldf_lev =  .true.     !  iso-level 
     294   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     295   ln_dynldf_iso =  .false.    !  iso-neutral 
     296   !                       !  Coefficient 
     297   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     298   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     299   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     300   !                                !  =  0  constant  
     301   !                                !  = 10  F(k)=c1d 
     302   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     303   !                                !  = 30  F(i,j,k)=c2d*c1d 
     304   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     305   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     306   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     307   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
     308   ! 
     309   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     310/ 
    270311   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    271312/ 
    … …  
    284325   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    285326/ 
    286 !------------------------------------------------------------------------ 
    287 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    288 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    289 / 
    290327!----------------------------------------------------------------------- 
    291328&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    … …  
    300337!----------------------------------------------------------------------- 
    301338   ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    302 / 
    303 !----------------------------------------------------------------------- 
    304 &namsol        !   elliptic solver / island / free surface 
    305 !----------------------------------------------------------------------- 
    306    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    307    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    308    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    309339/ 
    310340!----------------------------------------------------------------------- 
    … …  
    356386!----------------------------------------------------------------------- 
    357387/ 
    358 !----------------------------------------------------------------------- 
    359 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    360 !----------------------------------------------------------------------- 
    361    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    362 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_top_cfg

    r6401 r6404  
    2323!----------------------------------------------------------------------- 
    2424&namtrc_adv    !   advection scheme for passive tracer  
    25 !-----------------------------------------------------------------------   
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     25!----------------------------------------------------------------------- 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm

    r6401 r6404  
    1 bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_vectopt_loop key_top key_my_trc key_mpp_mpi key_iomput 
     1bld::tool::fppkeys key_zdftke key_top key_my_trc key_mpp_mpi key_iomput 
    22inc $BFMDIR/src/nemo/bfm.fcm 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r6401 r6404  
    2929&namzgr        !   vertical coordinate 
    3030!----------------------------------------------------------------------- 
    31    ln_zco      = .true.    !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
    32    ln_zps      = .false.   !  z-coordinate - partial steps   (T/F) 
     31   ln_zco      = .true.    !  z-coordinate - full    steps 
     32   ln_linssh   = .true.    !  linear free surface 
    3333/ 
    3434!----------------------------------------------------------------------- 
    … …  
    3636!----------------------------------------------------------------------- 
    3737   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    38    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    39    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    40    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
     38   rn_rdt      = 7200.     !  time step for the dynamics  
    4139   jphgr_msh   =       5                 !  type of horizontal mesh 
    4240   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
    … …  
    5957/ 
    6058!----------------------------------------------------------------------- 
    61 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    62 !----------------------------------------------------------------------- 
    63 / 
    64 !----------------------------------------------------------------------- 
    6559&namcrs        !   Grid coarsening for dynamics output and/or 
    6660               !   passive tracer coarsened online simulations 
    … …  
    10498/ 
    10599!----------------------------------------------------------------------- 
    106 &namcla        !   cross land advection 
    107 !----------------------------------------------------------------------- 
    108 / 
    109 !----------------------------------------------------------------------- 
    110100&nambfr        !   bottom friction 
    111101!----------------------------------------------------------------------- 
    … …  
    113103/ 
    114104!----------------------------------------------------------------------- 
    115 &nambbc        !   bottom temperature boundary condition 
    116 !----------------------------------------------------------------------- 
    117    ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    118    nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux 
     105&nambbc        !   bottom temperature boundary condition                (default: NO) 
     106!----------------------------------------------------------------------- 
    119107/ 
    120108!----------------------------------------------------------------------- 
    … …  
    143131&namtra_adv    !   advection scheme for tracer 
    144132!----------------------------------------------------------------------- 
    145    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     133   ln_traadv_fct =  .true.   !  FCT scheme 
     134      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     135      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     136      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     137      !                             !        (number of sub-timestep = nn_fct_zts) 
    146138/ 
    147139!---------------------------------------------------------------------------------- 
    148140&namtra_ldf    !   lateral diffusion scheme for tracers 
    149141!---------------------------------------------------------------------------------- 
    150    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    151    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     142   !                       !  Operator type: 
     143   ln_traldf_lap   =  .true.   !    laplacian operator 
     144   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     145   !                       !  Direction of action: 
     146   ln_traldf_lev   =  .false.  !  iso-level 
     147   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     148   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     149   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     150   ! 
     151   !                       !  iso-neutral options:         
     152   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     153   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     154   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     155   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     156   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     157   ! 
     158   !                       !  Coefficients: 
     159   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     160   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     161   !                                !   =  0           constant  
     162   !                                !   = 10 F(k)      =ldf_c1d  
     163   !                                !   = 20 F(i,j)    =ldf_c2d  
     164   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     165   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     166   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     167   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     168   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     169/ 
     170!---------------------------------------------------------------------------------- 
     171&namtra_ldfeiv !   eddy induced velocity param. 
     172!---------------------------------------------------------------------------------- 
     173   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    152174/ 
    153175!----------------------------------------------------------------------- 
    … …  
    163185&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    164186!----------------------------------------------------------------------- 
    165    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    166    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
     187   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
     188   ln_dynvor_ens = .false. !  energy conserving scheme 
     189   ln_dynvor_mix = .false. !  mixed scheme 
    167190   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     191      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    168192/ 
    169193!----------------------------------------------------------------------- 
    … …  
    174198/ 
    175199!----------------------------------------------------------------------- 
     200&namdyn_spg    !   surface pressure gradient 
     201!----------------------------------------------------------------------- 
     202   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     203/ 
     204!----------------------------------------------------------------------- 
    176205&namdyn_ldf    !   lateral diffusion on momentum 
    177206!----------------------------------------------------------------------- 
     207   !                       !  Type of the operator : 
     208   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     209   ln_dynldf_lap =  .true.    !    laplacian operator 
     210   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     211   !                       !  Direction of action  : 
     212   ln_dynldf_lev =  .true.    !  iso-level 
     213   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     214   ln_dynldf_iso =  .false.    !  iso-neutral 
     215   !                       !  Coefficient 
     216   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     217   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     218   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     219   !                                !  =  0  constant  
     220   !                                !  = 10  F(k)=c1d 
     221   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     222   !                                !  = 30  F(i,j,k)=c2d*c1d 
     223   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     224   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     225   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     226   rn_bhm_0      =      0.      !  horizontal bilaplacian eddy viscosity [m4/s] 
     227   ! 
     228   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     229/ 
    178230   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    179231/ 
    … …  
    189241/ 
    190242!----------------------------------------------------------------------- 
    191 &namsol        !   elliptic solver / island / free surface 
    192 !----------------------------------------------------------------------- 
    193    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    194    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    195    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    196 / 
    197 !----------------------------------------------------------------------- 
    198243&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    199244!----------------------------------------------------------------------- 
    … …  
    211256!----------------------------------------------------------------------- 
    212257/ 
    213 !----------------------------------------------------------------------- 
    214 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    215 !----------------------------------------------------------------------- 
    216    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    217 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_top_cfg

    r6401 r6404  
    2424&namtrc_adv    !   advection scheme for passive tracer  
    2525!-----------------------------------------------------------------------    
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

    r6401 r6404  
    1 bld::tool::fppkeys  key_dynspg_flt key_ldfslp key_zdftke key_top key_pisces_reduced key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_zdftke key_top key_pisces_reduced key_mpp_mpi 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_XIOS/EXP00/namelist_cfg

    r6401 r6404  
    2929&namzgr        !   vertical coordinate 
    3030!----------------------------------------------------------------------- 
    31    ln_zco      = .true.    !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
    32    ln_zps      = .false.   !  z-coordinate - partial steps   (T/F) 
     31   ln_zco      = .true.    !  z-coordinate - full    steps 
     32   ln_linssh   = .true.    !  linear free surface 
    3333/ 
    3434!----------------------------------------------------------------------- 
    … …  
    4040!----------------------------------------------------------------------- 
    4141   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    42    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
     42   rn_rdt      = 7200.     !  time step for the dynamics  
    4343!   nn_baro     =   60      !  number of barotropic time step            ("key_dynspg_ts") 
    44    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    45    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
    4644   jphgr_msh   =       5                 !  type of horizontal mesh 
    4745   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
    … …  
    150148/ 
    151149!----------------------------------------------------------------------- 
    152 &namcla        !   cross land advection 
    153 !----------------------------------------------------------------------- 
    154 / 
    155 !----------------------------------------------------------------------- 
    156 &namobc        !   open boundaries parameters                           ("key_obc") 
    157 !----------------------------------------------------------------------- 
    158 / 
    159 !----------------------------------------------------------------------- 
    160150&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    161151!----------------------------------------------------------------------- 
    … …  
    183173/ 
    184174!----------------------------------------------------------------------- 
    185 &nambbc        !   bottom temperature boundary condition 
    186 !----------------------------------------------------------------------- 
    187    ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    188    nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux 
     175&nambbc        !   bottom temperature boundary condition                (default: NO) 
     176!----------------------------------------------------------------------- 
    189177/ 
    190178!----------------------------------------------------------------------- 
    … …  
    200188&namtra_adv    !   advection scheme for tracer 
    201189!----------------------------------------------------------------------- 
    202    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     190   ln_traadv_fct =  .true.    !  FCT scheme 
     191      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     192      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     193      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     194      !                             !        (number of sub-timestep = nn_fct_zts) 
    203195/ 
    204196!----------------------------------------------------------------------- 
    … …  
    209201&namtra_ldf    !   lateral diffusion scheme for tracers 
    210202!---------------------------------------------------------------------------------- 
    211    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    212    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     203   !                       !  Operator type: 
     204   ln_traldf_lap   =  .true.   !    laplacian operator 
     205   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     206   !                       !  Direction of action: 
     207   ln_traldf_lev   =  .false.  !  iso-level 
     208   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     209   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     210   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     211   ! 
     212   !                       !  iso-neutral options:         
     213   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     214   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     215   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     216   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     217   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     218   ! 
     219   !                       !  Coefficients: 
     220   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     221   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     222   !                                !   =  0           constant  
     223   !                                !   = 10 F(k)      =ldf_c1d  
     224   !                                !   = 20 F(i,j)    =ldf_c2d  
     225   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     226   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     227   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     228   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     229   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     230/ 
     231!---------------------------------------------------------------------------------- 
     232&namtra_ldfeiv !   eddy induced velocity param. 
     233!---------------------------------------------------------------------------------- 
     234   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    213235/ 
    214236!----------------------------------------------------------------------- 
    … …  
    235257/ 
    236258!----------------------------------------------------------------------- 
    237 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    238 !----------------------------------------------------------------------- 
     259&namdyn_spg    !   surface pressure gradient 
     260!----------------------------------------------------------------------- 
     261   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     262/ 
    239263!----------------------------------------------------------------------- 
    240264&namdyn_ldf    !   lateral diffusion on momentum 
    … …  
    274298/ 
    275299!----------------------------------------------------------------------- 
    276 &namsol        !   elliptic solver / island / free surface 
    277 !----------------------------------------------------------------------- 
    278    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    279    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    280    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    281 / 
    282 !----------------------------------------------------------------------- 
    283300&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    284301!----------------------------------------------------------------------- 
    … …  
    328345!----------------------------------------------------------------------- 
    329346/ 
    330 !----------------------------------------------------------------------- 
    331 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    332 !----------------------------------------------------------------------- 
    333    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    334 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/GYRE_XIOS/cpp_GYRE_XIOS.fcm

    r6401 r6404  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_zdftke key_iomput key_mpp_mpi 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r6401 r6404  
    44!----------------------------------------------------------------------- 
    55&namrun        !   parameters of the run 
     6 nn_it000=1 
    67!----------------------------------------------------------------------- 
    78   cn_exp      = "Agulhas" !  experience name  
    8    nn_itend    =   10950   !  last  time step 
     9   nn_itend    =   10950 
    910   nn_stock    =   10950   !  frequency of creation of a restart file (modulo referenced to 1) 
    1011   nn_write    =   10950   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    11    ln_clobber  = .true.    !  clobber (overwrite) an existing file 
     12   ln_clobber  =  .true. 
    1213/ 
    1314!----------------------------------------------------------------------- 
    … …  
    2829&namzgr        !   vertical coordinate 
    2930!----------------------------------------------------------------------- 
     31   ln_zps      = .true.    !  z-coordinate - partial steps 
     32   ln_linssh   = .true.    !  linear free surface 
    3033/ 
    3134!----------------------------------------------------------------------- 
    … …  
    5154   ppacr2      =  999999.              ! 
    5255   rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
    53    rn_rdtmin   = 14400.          !  minimum time step on tracers (used if nn_acc=1) 
    54    rn_rdtmax   = 14400.          !  maximum time step on tracers (used if nn_acc=1) 
    55 / 
    56 !----------------------------------------------------------------------- 
    57 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    58 !----------------------------------------------------------------------- 
    5956/ 
    6057!----------------------------------------------------------------------- 
    … …  
    116113/ 
    117114!----------------------------------------------------------------------- 
    118 &namcla        !   cross land advection 
    119 !----------------------------------------------------------------------- 
    120 / 
    121 !----------------------------------------------------------------------- 
    122115&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    123116!----------------------------------------------------------------------- 
    … …  
    128121/ 
    129122!----------------------------------------------------------------------- 
    130 &nambbc        !   bottom temperature boundary condition 
    131 !----------------------------------------------------------------------- 
     123&nambbc        !   bottom temperature boundary condition                (default: NO) 
     124!----------------------------------------------------------------------- 
     125   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    132126/ 
    133127!----------------------------------------------------------------------- 
    … …  
    145139/ 
    146140!----------------------------------------------------------------------- 
    147 &namtra_adv    !   advection scheme for tracer  
    148 !----------------------------------------------------------------------- 
    149 / 
    150 !----------------------------------------------------------------------- 
    151 &namtra_ldf    !   lateral diffusion scheme for tracers  
    152 !----------------------------------------------------------------------- 
    153  
    154    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    155    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     141&namtra_adv    !   advection scheme for tracer 
     142!----------------------------------------------------------------------- 
     143   ln_traadv_fct =  .true.    !  FCT scheme 
     144      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order 
     145      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order 
     146      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     147      !                             !        (number of sub-timestep = nn_fct_zts) 
     148/ 
     149!----------------------------------------------------------------------- 
     150&namtra_ldf    !   lateral diffusion scheme for tracers 
     151!----------------------------------------------------------------------- 
     152   !                       !  Operator type: 
     153   ln_traldf_lap   =  .true.   !    laplacian operator 
     154   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     155   !                       !  Direction of action: 
     156   ln_traldf_lev   =  .false.  !  iso-level 
     157   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     158   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     159   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     160   ! 
     161   !                       !  iso-neutral options:         
     162   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     163   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     164   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     165   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     166   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     167   ! 
     168   !                       !  Coefficients: 
     169   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     170   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     171   !                                !   =  0           constant 
     172   !                                !   = 10 F(k)      =ldf_c1d 
     173   !                                !   = 20 F(i,j)    =ldf_c2d 
     174   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     175   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     176   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     177   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     178   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     179/ 
     180!---------------------------------------------------------------------------------- 
     181&namtra_ldfeiv !   eddy induced velocity param. 
     182!---------------------------------------------------------------------------------- 
     183   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    156184/ 
    157185!----------------------------------------------------------------------- 
    … …  
    162190&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    163191!----------------------------------------------------------------------- 
     192   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     193   ln_dynvor_ens = .false. !  energy conserving scheme 
     194   ln_dynvor_mix = .false. !  mixed scheme 
     195   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     196      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    164197/ 
    165198!----------------------------------------------------------------------- 
    … …  
    168201/ 
    169202!----------------------------------------------------------------------- 
     203&namdyn_spg    !   Surface pressure gradient 
     204!----------------------------------------------------------------------- 
     205   ln_dynspg_ts = .true.   !  Split-explicit free surface 
     206/ 
     207!----------------------------------------------------------------------- 
    170208&namdyn_ldf    !   lateral diffusion on momentum 
    171209!----------------------------------------------------------------------- 
    172    !                       !  Type of the operator :  
    173    ln_dynldf_lap    =  .false.  !  laplacian operator          
    174    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    175    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     210   !                       !  Type of the operator : 
     211   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F 
     212   ln_dynldf_lap =  .false.    !    laplacian operator 
     213   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     214   !                       !  Direction of action  : 
     215   ln_dynldf_lev =  .true.     !  iso-level 
     216   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     217   ln_dynldf_iso =  .false.    !  iso-neutral 
     218   !                       !  Coefficient 
     219   nn_ahm_ijk_t  =   0         !  space/time variation of eddy coef 
     220   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     221   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     222   !                                !  =  0  constant 
     223   !                                !  = 10  F(k)=c1d 
     224   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     225   !                                !  = 30  F(i,j,k)=c2d*c1d 
     226   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     227   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     228   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     229   rn_bhm_0      = 8.5e+11      !  horizontal bilaplacian eddy viscosity [m4/s] 
    176230/ 
    177231!----------------------------------------------------------------------- 
    … …  
    193247/ 
    194248!----------------------------------------------------------------------- 
    195 &namsol        !   elliptic solver / island / free surface  
    196 !----------------------------------------------------------------------- 
    197 / 
    198 !----------------------------------------------------------------------- 
    199249&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    200250!----------------------------------------------------------------------- 
    … …  
    217267/ 
    218268!----------------------------------------------------------------------- 
    219 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    220 !----------------------------------------------------------------------- 
    221 / 
    222 !----------------------------------------------------------------------- 
    223269&namobs       !  observation usage                                      ('key_diaobs') 
    224270!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r6401 r6404  
    66&namrun        !   parameters of the run 
    77!----------------------------------------------------------------------- 
     8   nn_no       =       0   !  job number (no more used...) 
     9   cn_exp      =  "ORCA2"  !  experience name 
     10   nn_it000    =       1   !  first time step 
     11   nn_itend    =     300   !  last  time step (std 5475) 
    812/ 
    913!----------------------------------------------------------------------- 
    … …  
    2428&namzgr        !   vertical coordinate 
    2529!----------------------------------------------------------------------- 
     30   ln_zps      = .true.    !  z-coordinate - partial steps 
     31   ln_linssh   = .true.    !  linear free surface 
    2632/ 
    2733!----------------------------------------------------------------------- 
    … …  
    4854/ 
    4955!----------------------------------------------------------------------- 
    50 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    51 !----------------------------------------------------------------------- 
    52 / 
    53 !----------------------------------------------------------------------- 
    5456&namcrs        !   Grid coarsening for dynamics output and/or 
    5557               !   passive tracer coarsened online simulations 
    … …  
    9395/ 
    9496!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9997&nambfr        !   bottom friction 
    10098!----------------------------------------------------------------------- 
    10199/ 
    102100!----------------------------------------------------------------------- 
    103 &nambbc        !   bottom temperature boundary condition 
    104 !----------------------------------------------------------------------- 
     101&nambbc        !   bottom temperature boundary condition                (default: NO) 
     102!----------------------------------------------------------------------- 
     103   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    105104/ 
    106105!----------------------------------------------------------------------- 
    … …  
    115114&namtra_adv    !   advection scheme for tracer 
    116115!----------------------------------------------------------------------- 
     116   ln_traadv_fct =  .true.    !  FCT scheme 
     117      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     118      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     119      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     120      !                             !        (number of sub-timestep = nn_fct_zts) 
    117121/ 
    118122!----------------------------------------------------------------------- 
    … …  
    120124!----------------------------------------------------------------------- 
    121125/ 
    122 !----------------------------------------------------------------------- 
    123 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    124 !----------------------------------------------------------------------- 
    125 / 
    126126!---------------------------------------------------------------------------------- 
    127127&namtra_ldf    !   lateral diffusion scheme for tracers 
    128128!---------------------------------------------------------------------------------- 
     129   !                       !  Operator type: 
     130   ln_traldf_lap   =  .true.   !    laplacian operator 
     131   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     132   !                       !  Direction of action: 
     133   ln_traldf_lev   =  .false.  !  iso-level 
     134   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     135   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     136   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     137   ! 
     138   !                       !  iso-neutral options:         
     139   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     140   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     141   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     142   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     143   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     144   ! 
     145   !                       !  Coefficients: 
     146   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     147   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     148   !                                !   =  0           constant 
     149   !                                !   = 10 F(k)      =ldf_c1d 
     150   !                                !   = 20 F(i,j)    =ldf_c2d 
     151   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     152   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     153   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     154   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     155   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     156/ 
     157!---------------------------------------------------------------------------------- 
     158&namtra_ldfeiv !   eddy induced velocity param. 
     159!---------------------------------------------------------------------------------- 
     160   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     161   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     162   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     163   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     164   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     165   !                                !   =  0           constant 
     166   !                                !   = 10 F(k)      =ldf_c1d 
     167   !                                !   = 20 F(i,j)    =ldf_c2d 
     168   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     169   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    129170/ 
    130171!----------------------------------------------------------------------- 
    … …  
    139180&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    140181!----------------------------------------------------------------------- 
     182   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     183   ln_dynvor_ens = .false. !  energy conserving scheme 
     184   ln_dynvor_mix = .false. !  mixed scheme 
     185   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     186      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    141187/ 
    142188!----------------------------------------------------------------------- 
    … …  
    145191/ 
    146192!----------------------------------------------------------------------- 
     193&namdyn_spg    !   Surface pressure gradient 
     194!----------------------------------------------------------------------- 
     195   ln_dynspg_ts = .true.   !  Split-explicit free surface 
     196/ 
     197!----------------------------------------------------------------------- 
    147198&namdyn_ldf    !   lateral diffusion on momentum 
    148199!----------------------------------------------------------------------- 
     200   !                       !  Type of the operator : 
     201   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     202   ln_dynldf_lap =  .true.     !    laplacian operator 
     203   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     204   !                       !  Direction of action  : 
     205   ln_dynldf_lev =  .true.     !  iso-level 
     206   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     207   ln_dynldf_iso =  .false.    !  iso-neutral 
     208   !                       !  Coefficient 
     209   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     210   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     211   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     212   !                                !  =  0  constant  
     213   !                                !  = 10  F(k)=c1d 
     214   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     215   !                                !  = 30  F(i,j,k)=c2d*c1d 
     216   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     217   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     218   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     219   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    149220/ 
    150221!----------------------------------------------------------------------- 
    … …  
    165236/ 
    166237!----------------------------------------------------------------------- 
    167 &namsol        !   elliptic solver / island / free surface 
    168 !----------------------------------------------------------------------- 
    169 / 
    170 !----------------------------------------------------------------------- 
    171238&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    172239!----------------------------------------------------------------------- 
    … …  
    185252/ 
    186253!----------------------------------------------------------------------- 
    187 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    188 !----------------------------------------------------------------------- 
    189 / 
    190 !----------------------------------------------------------------------- 
    191 &namobs       !  observation usage                                      ('key_diaobs') 
     254&namobs       !  observation usage 
    192255!----------------------------------------------------------------------- 
    193256/ 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

    r6401 r6404  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc  
     1bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_iomput key_nosignedzero 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg

    r6401 r6404  
    2828&namzgr        !   vertical coordinate 
    2929!----------------------------------------------------------------------- 
     30   ln_zps      = .true.    !  z-coordinate - partial steps 
    3031/ 
    3132!----------------------------------------------------------------------- 
    … …  
    5152   ppacr2      =  999999.              ! 
    5253   rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
    53    rn_rdtmin   = 14400.          !  minimum time step on tracers (used if nn_acc=1) 
    54    rn_rdtmax   = 14400.          !  maximum time step on tracers (used if nn_acc=1) 
    55 / 
    56 !----------------------------------------------------------------------- 
    57 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    58 !----------------------------------------------------------------------- 
    5954/ 
    6055!----------------------------------------------------------------------- 
    … …  
    116111/ 
    117112!----------------------------------------------------------------------- 
    118 &namcla        !   cross land advection 
    119 !----------------------------------------------------------------------- 
    120 / 
    121 !----------------------------------------------------------------------- 
    122113&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    123114!----------------------------------------------------------------------- 
    … …  
    130121&nambbc        !   bottom temperature boundary condition 
    131122!----------------------------------------------------------------------- 
     123   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    132124/ 
    133125!----------------------------------------------------------------------- 
    … …  
    145137/ 
    146138!----------------------------------------------------------------------- 
    147 &namtra_adv    !   advection scheme for tracer  
    148 !----------------------------------------------------------------------- 
    149 / 
    150 !----------------------------------------------------------------------- 
    151 &namtra_ldf    !   lateral diffusion scheme for tracers  
    152 !----------------------------------------------------------------------- 
    153  
    154    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    155    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     139&namtra_adv    !   advection scheme for tracer 
     140!----------------------------------------------------------------------- 
     141   ln_traadv_fct =  .true.    !  FCT scheme 
     142      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     143      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     144      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     145      !                             !        (number of sub-timestep = nn_fct_zts) 
     146/ 
     147!----------------------------------------------------------------------- 
     148&namtra_ldf    !   lateral diffusion scheme for tracers 
     149!---------------------------------------------------------------------------------- 
     150   !                       !  Operator type: 
     151   ln_traldf_lap   =  .true.   !    laplacian operator 
     152   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     153   !                       !  Direction of action: 
     154   ln_traldf_lev   =  .false.  !  iso-level 
     155   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     156   ln_traldf_iso   =  .true.   !  iso-neutral (Standard operator) 
     157   ln_traldf_triad =  .false.  !  iso-neutral (Triads   operator) 
     158   ! 
     159   !                       !  iso-neutral options:         
     160   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     161   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     162   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     163   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     164   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     165   ! 
     166   !                       !  Coefficients: 
     167   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     168   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     169   !                                !   =  0           constant  
     170   !                                !   = 10 F(k)      =ldf_c1d  
     171   !                                !   = 20 F(i,j)    =ldf_c2d  
     172   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     173   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     174   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     175   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     176   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     177/ 
     178!---------------------------------------------------------------------------------- 
     179&namtra_ldfeiv !   eddy induced velocity param. 
     180!---------------------------------------------------------------------------------- 
     181   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    156182/ 
    157183!----------------------------------------------------------------------- 
    … …  
    162188&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    163189!----------------------------------------------------------------------- 
     190   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     191   ln_dynvor_ens = .false. !  energy conserving scheme 
     192   ln_dynvor_mix = .false. !  mixed scheme 
     193   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     194      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     195   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    164196/ 
    165197!----------------------------------------------------------------------- 
    … …  
    168200/ 
    169201!----------------------------------------------------------------------- 
     202&namdyn_spg    !   surface pressure gradient 
     203!----------------------------------------------------------------------- 
     204   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     205/ 
     206!----------------------------------------------------------------------- 
    170207&namdyn_ldf    !   lateral diffusion on momentum 
    171208!----------------------------------------------------------------------- 
    172    !                       !  Type of the operator :  
    173    ln_dynldf_lap    =  .false.  !  laplacian operator          
    174    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    175    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     209   !                       !  Type of the operator : 
     210   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     211   ln_dynldf_lap =  .false.    !    laplacian operator 
     212   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     213   !                       !  Direction of action  : 
     214   ln_dynldf_lev =  .true.     !  iso-level 
     215   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     216   ln_dynldf_iso =  .false.    !  iso-neutral 
     217   !                       !  Coefficient 
     218   nn_ahm_ijk_t  = 20          !  space/time variation of eddy coef 
     219   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     220   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     221   !                                !  =  0  constant  
     222   !                                !  = 10  F(k)=c1d 
     223   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     224   !                                !  = 30  F(i,j,k)=c2d*c1d 
     225   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     226   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     227   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     228   rn_bhm_0      = 8.5e+11     !  horizontal bilaplacian eddy viscosity [m4/s] 
    176229/ 
    177230!----------------------------------------------------------------------- 
    … …  
    193246/ 
    194247!----------------------------------------------------------------------- 
    195 &namsol        !   elliptic solver / island / free surface  
    196 !----------------------------------------------------------------------- 
    197 / 
    198 !----------------------------------------------------------------------- 
    199248&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    200249!----------------------------------------------------------------------- 
    … …  
    217266/ 
    218267!----------------------------------------------------------------------- 
    219 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    220 !----------------------------------------------------------------------- 
    221 / 
    222 !----------------------------------------------------------------------- 
    223268&namobs       !  observation usage                                      ('key_diaobs') 
    224269!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r6401 r6404  
    66&namrun        !   parameters of the run 
    77!----------------------------------------------------------------------- 
     8   nn_no       =       0   !  job number (no more used...) 
     9   cn_exp      =  "ORCA2"  !  experience name 
     10   nn_it000    =       1   !  first time step 
     11   nn_itend    =     300   !  last  time step (std 5475) 
    812/ 
    913!----------------------------------------------------------------------- 
    … …  
    2428&namzgr        !   vertical coordinate 
    2529!----------------------------------------------------------------------- 
     30   ln_zps      = .true.    !  z-coordinate - partial steps 
    2631/ 
    2732!----------------------------------------------------------------------- 
    … …  
    4853/ 
    4954!----------------------------------------------------------------------- 
    50 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    51 !----------------------------------------------------------------------- 
    52 / 
    53 !----------------------------------------------------------------------- 
    5455&namcrs        !   Grid coarsening for dynamics output and/or 
    5556               !   passive tracer coarsened online simulations 
    … …  
    9394/ 
    9495!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9996&nambfr        !   bottom friction 
    10097!----------------------------------------------------------------------- 
    10198/ 
    10299!----------------------------------------------------------------------- 
    103 &nambbc        !   bottom temperature boundary condition 
    104 !----------------------------------------------------------------------- 
     100&nambbc        !   bottom temperature boundary condition                (default: NO) 
     101!----------------------------------------------------------------------- 
     102   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    105103/ 
    106104!----------------------------------------------------------------------- 
    … …  
    115113&namtra_adv    !   advection scheme for tracer 
    116114!----------------------------------------------------------------------- 
     115   ln_traadv_fct =  .true.    !  FCT scheme 
     116      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     117      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     118      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     119      !                             !        (number of sub-timestep = nn_fct_zts) 
    117120/ 
    118121!----------------------------------------------------------------------- 
    … …  
    123126&namtra_ldf    !   lateral diffusion scheme for tracers 
    124127!---------------------------------------------------------------------------------- 
     128   !                       !  Operator type: 
     129   ln_traldf_lap   =  .true.   !    laplacian operator 
     130   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     131   !                       !  Direction of action: 
     132   ln_traldf_lev   =  .false.  !  iso-level 
     133   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     134   ln_traldf_iso   =  .true.   !  iso-neutral (Standard operator) 
     135   ln_traldf_triad =  .false.  !  iso-neutral (Triads   operator) 
     136   ! 
     137   !                       !  iso-neutral options:         
     138   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     139   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     140   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     141   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     142   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     143   ! 
     144   !                       !  Coefficients: 
     145   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     146   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     147   !                                !   =  0           constant  
     148   !                                !   = 10 F(k)      =ldf_c1d  
     149   !                                !   = 20 F(i,j)    =ldf_c2d  
     150   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     151   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     152   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     153   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     154   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     155/ 
     156!---------------------------------------------------------------------------------- 
     157&namtra_ldfeiv !   eddy induced velocity param. 
     158!---------------------------------------------------------------------------------- 
     159   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     160   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     161   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     162   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     163   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     164   !                                !   =  0           constant  
     165   !                                !   = 10 F(k)      =ldf_c1d  
     166   !                                !   = 20 F(i,j)    =ldf_c2d  
     167   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     168   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    125169/ 
    126170!----------------------------------------------------------------------- 
    … …  
    135179&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    136180!----------------------------------------------------------------------- 
     181   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     182   ln_dynvor_ens = .false. !  energy conserving scheme 
     183   ln_dynvor_mix = .false. !  mixed scheme 
     184   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     185      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    137186/ 
    138187!----------------------------------------------------------------------- 
    139188&namdyn_hpg    !   Hydrostatic pressure gradient option 
    140189!----------------------------------------------------------------------- 
    141    ln_hpg_zco  = .false.   !  z-coordinate - full steps 
     190   ln_hpg_sco  = .true.    !  s-coordinate (standard jacobian formulation) 
    142191   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation) 
    143    ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
    144    ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    145    ln_hpg_prj  = .true.    !  s-coordinate (Pressure Jacobian scheme) 
    146    ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    147                                  !           centered      time scheme  (F) 
     192/ 
     193!----------------------------------------------------------------------- 
     194&namdyn_spg    !   surface pressure gradient 
     195!----------------------------------------------------------------------- 
     196   ln_dynspg_ts  = .true.  !  split-explicit free surface 
    148197/ 
    149198!----------------------------------------------------------------------- 
    150199&namdyn_ldf    !   lateral diffusion on momentum 
    151200!----------------------------------------------------------------------- 
     201   !                       !  Type of the operator : 
     202   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     203   ln_dynldf_lap =  .true.     !    laplacian operator 
     204   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     205   !                       !  Direction of action  : 
     206   ln_dynldf_lev =  .true.     !  iso-level 
     207   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     208   ln_dynldf_iso =  .false.    !  iso-neutral 
     209   !                       !  Coefficient 
     210   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     211   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     212   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     213   !                                !  =  0  constant  
     214   !                                !  = 10  F(k)=c1d 
     215   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     216   !                                !  = 30  F(i,j,k)=c2d*c1d 
     217   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     218   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     219   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     220   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     221   ! 
     222   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    152223/ 
    153224!----------------------------------------------------------------------- 
    … …  
    168239/ 
    169240!----------------------------------------------------------------------- 
    170 &namsol        !   elliptic solver / island / free surface 
    171 !----------------------------------------------------------------------- 
    172 / 
    173 !----------------------------------------------------------------------- 
    174241&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    175242!----------------------------------------------------------------------- 
    … …  
    184251/ 
    185252!----------------------------------------------------------------------- 
    186 &namhsb       !  Heat and salt budgets 
    187 !----------------------------------------------------------------------- 
    188 / 
    189 !----------------------------------------------------------------------- 
    190 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    191 !----------------------------------------------------------------------- 
    192 / 
     253&namhsb       !  Heat and salt budgets                                  (default F) 
     254!----------------------------------------------------------------------- 
     255   ln_diahsb  = .true.    !  check the heat and salt budgets (T) or not (F) 
     256/ 
     257!----------------------------------------------------------------------- 
     258&namobs       !  observation usage                                      ('key_diaobs') 
     259!----------------------------------------------------------------------- 
     260/ 
     261!----------------------------------------------------------------------- 
     262&nam_asminc   !   assimilation increments                               ('key_asminc') 
     263!----------------------------------------------------------------------- 
     264/ 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM3/cpp_ORCA2_LIM3.fcm

    r6401 r6404  
    1  bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc  
     1bld::tool::fppkeys key_trabbl key_lim3 key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_asminc key_nosignedzero 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_cfg

    r6401 r6404  
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,  
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     7!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    88!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    99!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
    … …  
    5959&namzgr        !   vertical coordinate 
    6060!----------------------------------------------------------------------- 
    61    ln_zco      = .false.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
    6261   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    63    ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
     62   ln_linssh   = .true.    !  linear free surface 
    6463/ 
    6564!----------------------------------------------------------------------- 
    … …  
    8685                           ! 
    8786   rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
    88    nn_baro     =   64      !  number of barotropic time step            ("key_dynspg_ts") 
    8987   rn_atfp     =    0.1    !  asselin time filter parameter 
    90    nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
    91                                  !                          =0, not used, rdt = rdttra 
    92    rn_rdtmin   = 14400.          !  minimum time step on tracers (used if nn_acc=1) 
    93    rn_rdtmax   = 14400.          !  maximum time step on tracers (used if nn_acc=1) 
    94    rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
    9588/ 
    9689!----------------------------------------------------------------------- 
    … …  
    297290   rn_alphdi   =    0.72   !  (Pyane, 1972) 
    298291/ 
     292!----------------------------------------------------------------------- 
     293&namsbc_wave   ! External fields from wave model 
     294!----------------------------------------------------------------------- 
     295!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     296!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     297   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' , .true.   , .false. , 'daily'  ,''         , '' 
     298! 
     299   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
     300/ 
    299301 
    300302!!====================================================================== 
    … …  
    302304!!====================================================================== 
    303305!!   namlbc        lateral momentum boundary condition 
    304 !!   namcla        cross land advection 
    305 !!   namobc        open boundaries parameters                           ("key_obc") 
    306306!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")  
    307307!!   nambdy        Unstructured open boundaries                         ("key_bdy") 
    … …  
    314314   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    315315                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 &namcla        !   cross land advection 
    319 !----------------------------------------------------------------------- 
    320    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    321 / 
    322 !----------------------------------------------------------------------- 
    323 &namobc        !   open boundaries parameters                           ("key_obc") 
    324 !----------------------------------------------------------------------- 
    325    ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    326    ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F) 
    327    ln_obc_fla  = .false.   !  Flather open boundary condition  
    328    nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    329                            !  = 1 the obc data are read in 'obc.dta' files 
    330    cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data 
    331                            !  set to monthly if obc datafile hold 1 month of data 
    332    rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary 
    333    rn_dpwin    =    1.     !     -           -         -       west    -      - 
    334    rn_dpnin    =    1.     !     -           -         -       north   -      - 
    335    rn_dpsin    =    1.     !     -           -         -       south   -      - 
    336    rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary 
    337    rn_dpwob    =   15.     !     -           -         -     west    -      - 
    338    rn_dpnob    = 3000.     !     -           -         -     north   -      - 
    339    rn_dpsob    =   15.     !     -           -         -     south   -      - 
    340    rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R) 
    341                            !  = 1 the total volume remains constant 
    342316/ 
    343317!----------------------------------------------------------------------- 
    … …  
    394368&nambdy_tide     ! tidal forcing at open boundaries               
    395369!----------------------------------------------------------------------- 
    396    filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
     370    filtide       = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
    397371    tide_cpt(1)   ='Q1'  !  names of tidal components used 
    398372    tide_cpt(2)   ='O1'  !  names of tidal components used 
    … …  
    485459/ 
    486460!----------------------------------------------------------------------- 
    487 &namtra_adv    !   advection scheme for tracer  
    488 !----------------------------------------------------------------------- 
    489    ln_traadv_cen2   =  .false.  !  2nd order centered scheme    
    490    ln_traadv_tvd    =  .true.   !  TVD scheme                 
    491    ln_traadv_muscl  =  .false.  !  MUSCL scheme              
    492    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries   
    493    ln_traadv_ubs    =  .false.  !  UBS scheme                  
    494    ln_traadv_qck    =  .false.  !  QUCIKEST scheme                  
     461&namtra_adv    !   advection scheme for tracer 
     462!----------------------------------------------------------------------- 
     463   ln_traadv_fct =  .true.    !  FCT scheme 
     464      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     465      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     466      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     467      !                             !        (number of sub-timestep = nn_fct_zts) 
    495468/ 
    496469!----------------------------------------------------------------------- 
    497470&namtra_ldf    !   lateral diffusion scheme for tracer  
    498471!----------------------------------------------------------------------- 
    499    !                       !  Type of the operator :  
    500    ln_traldf_lap    =  .true.   !  laplacian operator        
    501    ln_traldf_bilap  =  .false.  !  bilaplacian operator      
    502    !                       !  Direction of action  : 
    503    ln_traldf_level  =  .false.  !  iso-level                 
    504    ln_traldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" when ln_sco=T) 
    505    ln_traldf_iso    =  .true.   !  iso-neutral                          (require "key_ldfslp") 
    506    ln_traldf_grif   =  .false.  !  griffies skew flux formulation       (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    507    ln_traldf_gdia   =  .false.  !  griffies operator strfn diagnostics  (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    508    !                       !  Coefficient 
    509    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    510    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    511    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     472   !                       !  Operator type: 
     473   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     474   ln_traldf_lap   =  .true.   !    laplacian operator 
     475   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     476   !                       !  Direction of action: 
     477   ln_traldf_lev   =  .false.  !  iso-level 
     478   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     479   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     480   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     481   ! 
     482   !                       !  iso-neutral options:         
     483   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     484   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     485   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     486   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     487   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     488   ! 
     489   !                       !  Coefficients: 
     490   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     491   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     492   !                                !   =  0           constant  
     493   !                                !   = 10 F(k)      =ldf_c1d  
     494   !                                !   = 20 F(i,j)    =ldf_c2d  
     495   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     496   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     497   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     498   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     499   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     500/ 
     501!---------------------------------------------------------------------------------- 
     502&namtra_ldfeiv !   eddy induced velocity param. 
     503!---------------------------------------------------------------------------------- 
     504   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    512505/ 
    513506!----------------------------------------------------------------------- 
    … …  
    546539&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    547540!----------------------------------------------------------------------- 
    548    ln_dynvor_ene = .false. !  enstrophy conserving scheme   
    549    ln_dynvor_ens = .false. !  energy conserving scheme     
    550    ln_dynvor_mix = .false. !  mixed scheme                
    551    ln_dynvor_een = .true.  !  energy & enstrophy scheme   
     541   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     542   ln_dynvor_ens = .false. !  energy conserving scheme 
     543   ln_dynvor_mix = .false. !  mixed scheme 
     544   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     545      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     546   ln_dynvor_msk = .true.  !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    552547/ 
    553548!----------------------------------------------------------------------- 
    … …  
    559554   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    560555   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    561    ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    562                                  !           centered      time scheme  (F) 
    563 / 
    564 !----------------------------------------------------------------------- 
    565 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    566 !----------------------------------------------------------------------- 
    567 !                          !  explicit free surface                     ("key_dynspg_exp") 
    568 !                          !  filtered free surface                     ("key_dynspg_flt") 
    569 !                          !  split-explicit free surface               ("key_dynspg_ts") 
    570  
     556/ 
     557!----------------------------------------------------------------------- 
     558&namdyn_spg    !   surface pressure gradient 
     559!----------------------------------------------------------------------- 
     560   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     561/ 
    571562!----------------------------------------------------------------------- 
    572563&namdyn_ldf    !   lateral diffusion on momentum 
    573564!----------------------------------------------------------------------- 
    574    !                       !  Type of the operator :  
    575    ln_dynldf_lap    =  .false.  !  laplacian operator          
    576    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    577    !                       !  Direction of action  :  
    578    ln_dynldf_level  =  .false.  !  iso-level                
    579    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    580    ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     565   !                       !  Type of the operator : 
     566   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     567   ln_dynldf_lap =  .false.    !    laplacian operator 
     568   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     569   !                       !  Direction of action  : 
     570   ln_dynldf_lev =  .true.     !  iso-level 
     571   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     572   ln_dynldf_iso =  .false.    !  iso-neutral 
    581573   !                       !  Coefficient 
    582    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    583    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    584    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     574   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     575   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     576   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     577   !                                !  =  0  constant  
     578   !                                !  = 10  F(k)=c1d 
     579   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     580   !                                !  = 30  F(i,j,k)=c2d*c1d 
     581   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     582   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     583   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     584   rn_bhm_0      = 8.5e+11     !  horizontal bilaplacian eddy viscosity [m4/s] 
    585585/ 
    586586 
    … …  
    700700!!   nammpp_dyndist    Massively Parallel domain decomposition          ("key_agrif" && "key_mpp_dyndist") 
    701701!!   namctl            Control prints & Benchmark 
    702 !!   namsol            elliptic solver / island / free surface  
    703 !!====================================================================== 
    704 ! 
    705 !----------------------------------------------------------------------- 
    706 &namsol        !   elliptic solver / island / free surface  
    707 !----------------------------------------------------------------------- 
    708    nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    709                            !                   =2 successive-over-relaxation (sor) 
    710    nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test 
    711    rn_eps      =  1.e-6    !  absolute precision of the solver 
    712    nn_nmin     =    300    !  minimum of iterations for the SOR solver 
    713    nn_nmax     =    800    !  maximum of iterations for the SOR solver 
    714    nn_nmod     =     10    !  frequency of test for the SOR solver 
    715    rn_resmax   =  1.e-10   !  absolute precision for the SOR solver 
    716    rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    717 / 
     702!!====================================================================== 
     703! 
    718704!----------------------------------------------------------------------- 
    719705&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    … …  
    877863   ln_ignmis  = .true.     ! Logical switch for ignoring missing files    
    878864                           !     endailyavtypes   ENACT daily average types                     
    879    ln_grid_global = .true. 
    880865   ln_grid_search_lookup = .false. 
    881866/  
    … …  
    897882    salfixmin = -9999      !  Minimum salinity after applying the increments 
    898883/ 
    899 !----------------------------------------------------------------------- 
    900 &namsbc_wave   ! External fields from wave model 
    901 !----------------------------------------------------------------------- 
    902 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    903 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    904    sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' , .true.   , .false. , 'daily'  ,''         , '' 
    905 ! 
    906    cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
    907 / 
    908 !----------------------------------------------------------------------- 
    909 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    910 !----------------------------------------------------------------------- 
    911    ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 
    912    ln_neptsimp       = .false.  ! yes/no use simplified neptune 
    913  
    914    ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep 
    915    rn_tslse          =  1.2e4   ! value of lengthscale L at the equator 
    916    rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole 
    917    ! Specify whether to ramp down the Neptune velocity in shallow 
    918    ! water, and if so the depth range controlling such ramping down 
    919    ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water 
    920    rn_htrmin         =  100.0   ! min. depth of transition range 
    921    rn_htrmax         =  200.0   ! max. depth of transition range 
    922 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_top_cfg

    r6401 r6404  
    3939&namtrc_adv    !   advection scheme for passive tracer  
    4040!----------------------------------------------------------------------- 
    41    ln_trcadv_cen2   =  .false.  !  2nd order centered scheme    
    42    ln_trcadv_tvd    =  .true.  !  TVD scheme 
    43    ln_trcadv_muscl  =  .false.   !  MUSCL scheme 
    44    ln_trcadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    45    ln_trcadv_ubs    =  .false.  !  UBS scheme 
    46    ln_trcadv_qck    =  .false.  !  QUICKEST scheme 
     41   ln_trcadv_fct =  .true.   !  FCT scheme 
     42      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     43      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     44      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     45      !                             !        (number of sub-timestep = nn_fct_zts) 
    4746/ 
    4847!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_cfg

    r6401 r6404  
    3131&namzgr        !   vertical coordinate 
    3232!----------------------------------------------------------------------- 
     33   ln_zps      = .true.    !  z-coordinate - partial steps 
     34   ln_linssh   = .true.    !  linear free surface 
    3335/ 
    3436!----------------------------------------------------------------------- 
    … …  
    132134/ 
    133135!----------------------------------------------------------------------- 
    134 &namcla        !   cross land advection 
    135 !----------------------------------------------------------------------- 
    136 / 
    137 !----------------------------------------------------------------------- 
    138 &namobc        !   open boundaries parameters                           ("key_obc") 
    139 !----------------------------------------------------------------------- 
    140 / 
    141 !----------------------------------------------------------------------- 
    142136&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    143137!----------------------------------------------------------------------- 
    … …  
    164158/ 
    165159!----------------------------------------------------------------------- 
    166 &nambbc        !   bottom temperature boundary condition 
     160&nambbc        !   bottom temperature boundary condition                (default: NO) 
    167161!-----------------------------------------------------------------------  
     162   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    168163/ 
    169164!----------------------------------------------------------------------- 
    … …  
    176171/ 
    177172!----------------------------------------------------------------------- 
    178 &namtra_adv    !   advection scheme for tracer  
    179 !-----------------------------------------------------------------------      
     173&namtra_adv    !   advection scheme for tracer 
     174!----------------------------------------------------------------------- 
     175   ln_traadv_fct =  .true.    !  FCT scheme 
     176      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     177      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     178      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     179      !                             !        (number of sub-timestep = nn_fct_zts) 
    180180/ 
    181181!---------------------------------------------------------------------------------- 
    182182&namtra_ldf    !   lateral diffusion scheme for tracers 
    183183!---------------------------------------------------------------------------------- 
     184   !                       !  Operator type: 
     185   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     186   ln_traldf_lap   =  .true.   !    laplacian operator 
     187   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     188   !                       !  Direction of action: 
     189   ln_traldf_lev   =  .false.  !  iso-level 
     190   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     191   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     192   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     193   ! 
     194   !                       !  iso-neutral options:         
     195   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     196   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     197   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     198   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     199   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     200   ! 
     201   !                       !  Coefficients: 
     202   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     203   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     204   !                                !   =  0           constant  
     205   !                                !   = 10 F(k)      =ldf_c1d  
     206   !                                !   = 20 F(i,j)    =ldf_c2d  
     207   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     208   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     209   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     210   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     211   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     212/ 
     213!---------------------------------------------------------------------------------- 
     214&namtra_ldfeiv !   eddy induced velocity param. 
     215!---------------------------------------------------------------------------------- 
     216   ln_ldfeiv     =.true.    ! use eddy induced velocity parameterization 
     217   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     218   rn_aeiv_0     = 2000.    ! eddy induced velocity coefficient   [m2/s] 
     219   nn_aei_ijk_t  = 21       ! space/time variation of the eiv coeficient 
     220   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     221   !                                !   =  0           constant  
     222   !                                !   = 10 F(k)      =ldf_c1d  
     223   !                                !   = 20 F(i,j)    =ldf_c2d  
     224   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     225   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    184226/ 
    185227!----------------------------------------------------------------------- 
    … …  
    188230/ 
    189231!----------------------------------------------------------------------- 
     232&namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
     233!----------------------------------------------------------------------- 
     234/ 
     235!----------------------------------------------------------------------- 
    190236&namdyn_adv    !   formulation of the momentum advection 
    191237!----------------------------------------------------------------------- 
    192238/   
    193239!----------------------------------------------------------------------- 
    194 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    195 !----------------------------------------------------------------------- 
    196 / 
    197 !----------------------------------------------------------------------- 
    198240&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    199241!----------------------------------------------------------------------- 
     242   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     243   ln_dynvor_ens = .false. !  energy conserving scheme 
     244   ln_dynvor_mix = .false. !  mixed scheme 
     245   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     246      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    200247/ 
    201248!----------------------------------------------------------------------- 
    … …  
    209256&namdyn_ldf    !   lateral diffusion on momentum 
    210257!-----------------------------------------------------------------------       
     258   !                       !  Type of the operator : 
     259   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     260   ln_dynldf_lap =  .true.     !    laplacian operator 
     261   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     262   !                       !  Direction of action  : 
     263   ln_dynldf_lev =  .true.     !  iso-level 
     264   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     265   ln_dynldf_iso =  .false.    !  iso-neutral 
     266   !                       !  Coefficient 
     267   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     268   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     269   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     270   !                                !  =  0  constant  
     271   !                                !  = 10  F(k)=c1d 
     272   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     273   !                                !  = 30  F(i,j,k)=c2d*c1d 
     274   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     275   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     276   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     277   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    211278/ 
    212279!----------------------------------------------------------------------- 
    … …  
    221288&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    222289!----------------------------------------------------------------------- 
    223 / 
    224 !------------------------------------------------------------------------ 
    225 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    226 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    227290/ 
    228291!----------------------------------------------------------------------- 
    … …  
    291354!----------------------------------------------------------------------- 
    292355/ 
    293 !----------------------------------------------------------------------- 
    294 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    295 !----------------------------------------------------------------------- 
    296 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_top_cfg

    r6401 r6404  
    2424&namtrc_adv    !   advection scheme for passive tracer  
    2525!-----------------------------------------------------------------------    
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/cpp_ORCA2_LIM_CFC_C14b.fcm

    r6401 r6404  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_cfg

    r6401 r6404  
    2424&namzgr        !   vertical coordinate 
    2525!----------------------------------------------------------------------- 
     26   ln_zps      = .true.    !  z-coordinate - partial steps 
     27   ln_linssh   = .true.    !  linear free surface 
    2628/ 
    2729!----------------------------------------------------------------------- 
    … …  
    4850/ 
    4951!----------------------------------------------------------------------- 
    50 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    51 !----------------------------------------------------------------------- 
    52 / 
    53 !----------------------------------------------------------------------- 
    5452&namcrs        !   Grid coarsening for dynamics output and/or 
    5553               !   passive tracer coarsened online simulations 
    … …  
    9391/ 
    9492!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9993&nambfr        !   bottom friction 
    10094!----------------------------------------------------------------------- 
    10195/ 
    10296!----------------------------------------------------------------------- 
    103 &nambbc        !   bottom temperature boundary condition 
    104 !----------------------------------------------------------------------- 
     97&nambbc        !   bottom temperature boundary condition                (default: NO) 
     98!----------------------------------------------------------------------- 
     99   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    105100/ 
    106101!----------------------------------------------------------------------- 
    … …  
    115110&namtra_adv    !   advection scheme for tracer 
    116111!----------------------------------------------------------------------- 
     112   ln_traadv_fct =  .true.    !  FCT scheme 
     113      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     114      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     115      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     116      !                             !        (number of sub-timestep = nn_fct_zts) 
    117117/ 
    118118!----------------------------------------------------------------------- 
    … …  
    120120!----------------------------------------------------------------------- 
    121121/ 
    122 !----------------------------------------------------------------------- 
    123 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    124 !----------------------------------------------------------------------- 
    125 / 
    126122!---------------------------------------------------------------------------------- 
    127123&namtra_ldf    !   lateral diffusion scheme for tracers 
    128124!---------------------------------------------------------------------------------- 
     125   !                       !  Operator type: 
     126   ln_traldf_lap   =  .true.   !    laplacian operator 
     127   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     128   !                       !  Direction of action: 
     129   ln_traldf_lev   =  .false.  !  iso-level 
     130   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     131   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     132   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     133   ! 
     134   !                       !  iso-neutral options:         
     135   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     136   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     137   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     138   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     139   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     140   ! 
     141   !                       !  Coefficients: 
     142   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     143   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     144   !                                !   =  0           constant 
     145   !                                !   = 10 F(k)      =ldf_c1d 
     146   !                                !   = 20 F(i,j)    =ldf_c2d 
     147   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     148   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     149   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     150   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     151   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     152/ 
     153!---------------------------------------------------------------------------------- 
     154&namtra_ldfeiv !   eddy induced velocity param. 
     155!---------------------------------------------------------------------------------- 
     156   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     157   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     158   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     159   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     160   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     161   !                                !   =  0           constant 
     162   !                                !   = 10 F(k)      =ldf_c1d 
     163   !                                !   = 20 F(i,j)    =ldf_c2d 
     164   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     165   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    129166/ 
    130167!----------------------------------------------------------------------- 
    … …  
    139176&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    140177!----------------------------------------------------------------------- 
     178   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     179   ln_dynvor_ens = .false. !  energy conserving scheme 
     180   ln_dynvor_mix = .false. !  mixed scheme 
     181   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     182      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    141183/ 
    142184!----------------------------------------------------------------------- 
    … …  
    145187/ 
    146188!----------------------------------------------------------------------- 
     189&namdyn_spg    !   surface pressure gradient 
     190!----------------------------------------------------------------------- 
     191   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     192/ 
     193!----------------------------------------------------------------------- 
    147194&namdyn_ldf    !   lateral diffusion on momentum 
    148195!----------------------------------------------------------------------- 
     196   !                       !  Type of the operator : 
     197   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     198   ln_dynldf_lap =  .true.     !    laplacian operator 
     199   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     200   !                       !  Direction of action  : 
     201   ln_dynldf_lev =  .true.     !  iso-level 
     202   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     203   ln_dynldf_iso =  .false.    !  iso-neutral 
     204   !                       !  Coefficient 
     205   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     206   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     207   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     208   !                                !  =  0  constant  
     209   !                                !  = 10  F(k)=c1d 
     210   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     211   !                                !  = 30  F(i,j,k)=c2d*c1d 
     212   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     213   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     214   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     215   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    149216/ 
    150217!----------------------------------------------------------------------- 
    … …  
    165232/ 
    166233!----------------------------------------------------------------------- 
    167 &namsol        !   elliptic solver / island / free surface 
    168 !----------------------------------------------------------------------- 
    169 / 
    170 !----------------------------------------------------------------------- 
    171234&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    172235!----------------------------------------------------------------------- 
    … …  
    184247!----------------------------------------------------------------------- 
    185248/ 
    186 !----------------------------------------------------------------------- 
    187 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    188 !----------------------------------------------------------------------- 
    189 / 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_top_cfg

    r6401 r6404  
    6262&namtrc_adv    !   advection scheme for passive tracer  
    6363!----------------------------------------------------------------------- 
    64    ln_trcadv_tvd     =  .false.  !  TVD scheme 
    65    ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     64   ln_trcadv_mus =  .true.  !  MUSCL scheme 
     65      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
    6666/ 
    6767!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm

    r6401 r6404  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_mpp_mpi key_iomput 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r6401 r6404  
    2727&namzgr        !   vertical coordinate 
    2828!----------------------------------------------------------------------- 
     29   ln_zps      = .true.    !  z-coordinate - partial steps 
     30   ln_linssh   = .true.    !  linear free surface 
    2931/ 
    3032!----------------------------------------------------------------------- 
    … …  
    3234!----------------------------------------------------------------------- 
    3335   nn_msh      =    1      !  create (=1) a mesh file or not (=0) 
    34    rn_rdt      = 21600.    !  time step for the dynamics (and tracer if nn_acc=0) 
    35    rn_rdtmin   = 21600.          !  minimum time step on tracers (used if nn_acc=1) 
    36    rn_rdtmax   = 21600.          !  maximum time step on tracers (used if nn_acc=1) 
     36   rn_rdt      = 21600.    !  time step for the dynamics  
    3737   jphgr_msh   =       0               !  type of horizontal mesh 
    3838   ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
    … …  
    6868   nn_fsbc     = 1         !  frequency of surface boundary condition computation  
    6969   ln_rnf      = .false.   !  runoffs 
     70   ln_traqsr   = .false.   !  Light penetration (T) or not (F) 
    7071/   
    7172!----------------------------------------------------------------------- 
    7273&namtra_qsr    !   penetrative solar radiation 
    7374!----------------------------------------------------------------------- 
    74    ln_traqsr   = .false.   !  Light penetration (T) or not (F) 
    7575/ 
    7676!----------------------------------------------------------------------- 
    … …  
    7878!----------------------------------------------------------------------- 
    7979/       
    80  
    81 !----------------------------------------------------------------------- 
    82 &namcla        !   cross land advection 
    83 !----------------------------------------------------------------------- 
    84 / 
    8580!----------------------------------------------------------------------- 
    8681&nambbl        !   bottom boundary layer scheme 
    … …  
    9186!----------------------------------------------------------------------- 
    9287/ 
    93 !----------------------------------------------------------------------- 
    94 &namtra_ldf    !   lateral diffusion scheme for tracer  
    95 !----------------------------------------------------------------------- 
    96    ln_triad_iso     =  .true.  !  griffies operator calculates triads twice => pure lateral mixing in ML (require "key_ldfslp") 
     88!---------------------------------------------------------------------------------- 
     89&namtra_ldf    !   lateral diffusion scheme for tracers 
     90!---------------------------------------------------------------------------------- 
     91   !                       !  Operator type: 
     92   ln_traldf_lap   =  .true.   !    laplacian operator 
     93   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     94   !                       !  Direction of action: 
     95   ln_traldf_lev   =  .false.  !  iso-level 
     96   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     97   ln_traldf_iso   =  .true.   !  iso-neutral 
     98   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     99   ! 
     100   !                       !  iso-neutral options:         
     101   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     102   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     103   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     104   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     105   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     106   ! 
     107   !                       !  Coefficients: 
     108   nn_aht_ijk_t    = 0        !  space/time variation of eddy coef 
     109   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     110   !                                !   =  0           constant  
     111   !                                !   = 10 F(k)      =ldf_c1d  
     112   !                                !   = 20 F(i,j)    =ldf_c2d  
     113   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     114   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     115   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     116   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     117   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     118/ 
     119!---------------------------------------------------------------------------------- 
     120&namtra_ldfeiv !   eddy induced velocity param. 
     121!---------------------------------------------------------------------------------- 
     122   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     123   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     124   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     125   nn_aei_ijk_t  = 0      ! space/time variation of the eiv coeficient 
     126   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     127   !                                !   =  0           constant 
     128   !                                !   = 10 F(k)      =ldf_c1d 
     129   !                                !   = 20 F(i,j)    =ldf_c2d 
     130   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     131   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    97132/ 
    98133!----------------------------------------------------------------------- 
    … …  
    115150   sn_ubl  = 'dyna_grid_U' ,    120            , 'sobblcox' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    116151   sn_vbl  = 'dyna_grid_V' ,    120            , 'sobblcoy' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    117    sn_ahu  = 'dyna_grid_U' ,    120            , 'vozoahtu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    118    sn_ahv  = 'dyna_grid_V' ,    120            , 'vomeahtv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    119    sn_ahw  = 'dyna_grid_W' ,    120            , 'voveahtz' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    120    sn_eiu  = 'dyna_grid_U' ,    120            , 'vozoaeiu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    121    sn_eiv  = 'dyna_grid_V' ,    120            , 'vomeaeiv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    122    sn_eiw  = 'dyna_grid_W' ,    120            , 'soleaeiw' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    123152! 
    124153   cn_dir      = './'      !  root directory for the location of the dynamical files 
    125    ln_degrad   =  .false.  !  flag for degradation -                requires ("key_degrad") 
    126154   ln_dynwzv   =  .true.   !  computation of vertical velocity instead of using the one read in file 
    127155   ln_dynbbl   =  .true.   !  bbl coef are in files, so read them - requires ("key_trabbl") 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_top_cfg

    r6401 r6404  
    6565&namtrc_adv    !   advection scheme for passive tracer  
    6666!----------------------------------------------------------------------- 
    67    ln_trcadv_tvd     =  .false.  !  TVD scheme 
    68    ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     67   ln_trcadv_mus =  .true.  !  MUSCL scheme 
     68      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
    6969/ 
    7070!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

    r6401 r6404  
    1 bld::tool::fppkeys key_trabbl key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_top key_offline key_pisces key_iomput key_mpp_mpi 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/namelist_cfg

    r6401 r6404  
    2626&namzgr        !   vertical coordinate 
    2727!----------------------------------------------------------------------- 
     28   ln_zco      = .true.    !  z-coordinate - full    steps 
     29   ln_linssh   = .true.    !  linear free surface 
    2830/ 
    2931!----------------------------------------------------------------------- 
    … …  
    4850   ppkth2      =  999999.              ! 
    4951   ppacr2      =  999999.              ! 
    50 / 
    51 !----------------------------------------------------------------------- 
    52 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    53 !----------------------------------------------------------------------- 
    5452/ 
    5553!----------------------------------------------------------------------- 
    … …  
    8785/ 
    8886!----------------------------------------------------------------------- 
    89 &namcla        !   cross land advection 
    90 !----------------------------------------------------------------------- 
    91 / 
    92 !----------------------------------------------------------------------- 
    9387&nameos        !   ocean physical parameters 
    9488!----------------------------------------------------------------------- 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/ORCA2_SAS_LIM/cpp_ORCA2_SAS_LIM.fcm

    r6401 r6404  
    1  bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/SHARED/domain_def.xml

    r6401 r6404  
    163163 
    164164      <domain id="grid_V" long_name="grid V"/> 
     165       
    165166      <domain_group id="grid_W"> 
    166167         <domain id="grid_W" long_name="grid W"/> 
    … …  
    168169         <domain id="EqW" zoom_ibegin="1" zoom_jbegin="0000" zoom_ni="0000" zoom_nj="1" /> 
    169170      </domain_group> 
     171 
     172      <domain id="grid_F" long_name="grid F"/> 
    170173 
    171174      <domain_group id="scalarpoint">  

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/SHARED/field_def.xml

    r6401 r6404  
    2323      <field_group id="grid_T" grid_ref="grid_T_2D" > 
    2424         <field id="e3t"          long_name="T-cell thickness"   standard_name="cell_thickness"   unit="m"   grid_ref="grid_T_3D"/> 
     25         <field id="e3t_0"        long_name="Initial T-cell thickness"   standard_name="ref_cell_thickness"   unit="m"   grid_ref="grid_T_3D"/> 
    2526 
    2627         <field id="toce"         long_name="temperature"         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
    27          <field id="toce_e3t"     long_name="temperature * e3t"                                                     unit="degC*m"   grid_ref="grid_T_3D" > toce * e3t </field > 
     28         <field id="toce_e3t"     long_name="temperature (thickness weighted)"                                      unit="degC"     grid_ref="grid_T_3D" > toce * e3t </field > 
    2829         <field id="soce"         long_name="salinity"            standard_name="sea_water_practical_salinity"      unit="1e-3"     grid_ref="grid_T_3D"/> 
    29          <field id="soce_e3t"     long_name="salinity * e3t"                                                        unit="1e-3*m"   grid_ref="grid_T_3D" > soce * e3t </field > 
     30         <field id="soce_e3t"     long_name="salinity    (thickness weighted)"                                      unit="1e-3"     grid_ref="grid_T_3D" > soce * e3t </field > 
     31 
     32         <!-- t-eddy viscosity coefficients (ldfdyn) --> 
     33         <field id="ahmt_2d"      long_name=" surface t-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
     34         <field id="ahmt_3d"      long_name=" 3D      t-eddy viscosity coefficient"   unit="m2/s or m4/s"                           grid_ref="grid_T_3D"/> 
    3035 
    3136         <field id="sst"          long_name="sea surface temperature"             standard_name="sea_surface_temperature"             unit="degC"     /> 
    … …  
    3641         <field id="sstgrad2"     long_name="square of module of sst gradient"                                                        unit="degC2/m2" /> 
    3742         <field id="sbt"          long_name="sea bottom temperature"                                                                  unit="degC"     /> 
    38      
     43         <field id="sst_wl"       long_name="Delta SST of warm layer"                                                                 unit="degC"     /> 
     44         <field id="sst_cs"       long_name="Delta SST of cool skin"                                                                  unit="degC"     /> 
     45    <field id="temp_3m"      long_name="temperature at 3m"                                                                       unit="degC"     /> 
     46          
    3947         <field id="sss"          long_name="sea surface salinity"             standard_name="sea_surface_salinity"   unit="1e-3" /> 
    4048         <field id="sss2"         long_name="square of sea surface salinity"                                          unit="1e-6"  > sss * sss </field > 
    … …  
    174182      <field_group id="SBC" grid_ref="grid_T_2D" > <!-- time step automaticaly defined based on nn_fsbc --> 
    175183         <field id="empmr"        long_name="Net Upward Water Flux"                standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   /> 
     184         <field id="empbmr"       long_name="Net Upward Water Flux at pre. tstep"  standard_name="water_flux_out_of_sea_ice_and_sea_water"                              unit="kg/m2/s"   /> 
    176185         <field id="saltflx"      long_name="Downward salt flux"                                                                                                        unit="1e-3/m2/s" /> 
    177186         <field id="fmmflx"       long_name="Water flux due to freezing/melting"                                                                                        unit="kg/m2/s"   /> 
    … …  
    193202 
    194203         <!-- * variable related to ice shelf forcing * --> 
    195          <field id="fwfisf"       long_name="Ice shelf melting"                             unit="kg/m2/s"  /> 
    196          <field id="qisf"         long_name="Ice Shelf Heat Flux"                           unit="W/m2"     /> 
    197          <field id="isfgammat"    long_name="transfert coefficient for isf (temperature)"   unit="m/s"      /> 
    198          <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)"      unit="m/s"      /> 
    199          <field id="stbl"         long_name="salinity in the Losh tbl"                      unit="1e-3"     /> 
    200          <field id="ttbl"         long_name="temperature in the Losh tbl"                   unit="degC"     /> 
     204         <field id="fwfisf"       long_name="Ice shelf melting"                                            unit="Kg/m2/s"  /> 
     205         <field id="qisf"         long_name="Ice Shelf Heat Flux"                                          unit="W/m2"     /> 
     206         <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "                 unit="m/s"      /> 
     207         <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "                 unit="m/s"      /> 
     208    <field id="stbl"         long_name="salinity in the Losh tbl                    "                 unit="PSU"      /> 
     209    <field id="ttbl"         long_name="temperature in the Losh tbl                 "                 unit="C"        /> 
     210    <field id="utbl"         long_name="zonal current in the Losh tbl at T point    "                 unit="m/s"      /> 
     211    <field id="vtbl"         long_name="merid current in the Losh tbl at T point    "                 unit="m/s"      /> 
     212    <field id="thermald"     long_name="thermal driving of ice shelf melting        "                 unit="C"        /> 
     213    <field id="tfrz"         long_name="top freezing point (used to compute melt)   "                 unit="C"        /> 
     214    <field id="tinsitu"      long_name="top insitu temperature (used to cmpt melt)  "                 unit="C"        /> 
     215    <field id="ustar"        long_name="ustar at T point used in ice shelf melting  "                 unit="m/s"      /> 
    201216 
    202217         <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
    … …  
    365380      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    366381         <field id="e3u"          long_name="U-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_U_3D" /> 
     382         <field id="e3u_0"        long_name="Initial U-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_U_3D"/> 
    367383         <field id="utau"         long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            /> 
    368384         <field id="uoce"         long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" /> 
    369          <field id="uoce_e3u"     long_name="ocean current along i-axis * e3u"                                                                   unit="m2/s"       grid_ref="grid_U_3D"  > uoce * e3u </field> 
     385         <field id="uoce_e3u"     long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
    370386         <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             /> 
    371387         <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"                             /> 
    … …  
    375391         <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" /> 
    376392 
     393         <!-- u-eddy coefficients (ldftra) --> 
     394         <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" /> 
     395         <field id="ahtu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s or m4/s"      grid_ref="grid_U_3D"/> 
     396         <field id="aeiu_2d"      long_name=" surface u-EIV coefficient"                unit="m2/s" /> 
     397         <field id="aeiu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s"              grid_ref="grid_U_3D"/> 
     398 
    377399         <!-- variables available with MLE --> 
    378400         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"   unit="m3/s"   grid_ref="grid_U_3D" /> 
    379401 
    380          <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     402         <!-- uoce_eiv: available EIV --> 
    381403         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"   standard_name="bolus_sea_water_x_velocity"   unit="m/s"   grid_ref="grid_U_3D" /> 
    382404 
    … …  
    400422      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
    401423         <field id="e3v"          long_name="V-cell thickness"                                       standard_name="cell_thickness"              unit="m"          grid_ref="grid_V_3D" /> 
     424         <field id="e3v_0"        long_name="Initial V-cell thickness"                               standard_name="ref_cell_thickness"          unit="m"          grid_ref="grid_V_3D"/> 
    402425         <field id="vtau"         long_name="Wind Stress along j-axis"                               standard_name="surface_downward_y_stress"   unit="N/m2"                            /> 
    403426         <field id="voce"         long_name="ocean current along j-axis"                             standard_name="sea_water_y_velocity"        unit="m/s"        grid_ref="grid_V_3D" /> 
    404          <field id="voce_e3v"     long_name="ocean current along j-axis * e3v"                                                                   unit="m2/s"       grid_ref="grid_V_3D"  > voce * e3v </field> 
     427         <field id="voce_e3v"     long_name="ocean current along j-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_V_3D"  > voce * e3v </field> 
    405428         <field id="ssv"          long_name="ocean surface current along j-axis"                                                                 unit="m/s"                             /> 
    406429         <field id="sbv"          long_name="ocean bottom current along j-axis"                                                                  unit="m/s"                             /> 
    … …  
    410433         <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" /> 
    411434 
     435         <!-- v-eddy coefficients (ldftra, ldfdyn) --> 
     436         <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" /> 
     437         <field id="ahtv_3d"      long_name=" 3D v-eddy diffusivity coefficient"          unit="m2/s or (m4/s)^1/2"           grid_ref="grid_V_3D"/> 
     438         <field id="aeiv_2d"      long_name=" surface v-EIV coefficient"                  unit="m2/s" /> 
     439         <field id="aeiv_3d"      long_name=" 3D v-EIV coefficient"                       unit="m2/s"                         grid_ref="grid_V_3D" /> 
     440 
    412441         <!-- variables available with MLE --> 
    413442         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"   unit="m3/s"   grid_ref="grid_V_3D" /> 
    414443 
    415          <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     444         <!-- voce_eiv: available with EIV --> 
    416445         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"   standard_name="bolus_sea_water_y_velocity"   unit="m/s"   grid_ref="grid_V_3D" /> 
    417446 
    … …  
    438467        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                               unit="m3/s" /> 
    439468 
    440         <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     469        <!-- woce_eiv: available with EIV --> 
    441470        <field id="woce_eiv"     long_name="EIV ocean vertical velocity"   standard_name="bolus_upward_sea_water_velocity"   unit="m/s" /> 
    442471 
    443         <!-- woce_eiv: available with key_trabbl_adv --> 
     472 
    444473        <field id="avt"          long_name="vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
     474        <field id="logavt"       long_name="logarithm of vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    445475        <field id="avm"          long_name="vertical eddy viscosity"     standard_name="ocean_vertical_momentum_diffusivity"   unit="m2/s" /> 
    446476 
    447477        <!-- avs: available with key_zdfddm --> 
    448478        <field id="avs"          long_name="salt vertical eddy diffusivity"   standard_name="ocean_vertical_salt_diffusivity"   unit="m2/s" /> 
     479        <field id="logavs"       long_name="logarithm of salt vertical eddy diffusivity"   standard_name="ocean_vertical_heat_diffusivity"       unit="m2/s" /> 
    449480 
    450481        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    451         <field id="avt_evd"      long_name="convective enhancement to vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
    452         <field id="avm_evd"      long_name="convective enhancement to vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
     482        <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
     483        <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    453484 
    454485        <!-- avt_tide: available with key_zdftmx --> 
    … …  
    459490        <field id="w_masstr2"    long_name="square of vertical mass transport"   standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" /> 
    460491 
    461         <!-- aht2d and  aht2d_eiv: available with key_traldf_eiv and key_traldf_c2d --> 
     492        <!-- aht2d and  aht2d_eiv --> 
    462493        <field id="aht2d"        long_name="lateral eddy diffusivity"       standard_name="ocean_tracer_xy_laplacian_diffusivity"      unit="m2/s"   grid_ref="grid_W_2D" /> 
    463494        <field id="aht2d_eiv"    long_name="EIV lateral eddy diffusivity"   standard_name="ocean_tracer_bolus_laplacian_diffusivity"   unit="m2/s"   grid_ref="grid_W_2D" /> 
    464495      </field_group> 
    465            
     496         
     497      <!-- F grid --> 
     498         <!-- f-eddy viscosity coefficients (ldfdyn) --> 
     499         <field id="ahmf_2d"      long_name=" surface f-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
     500         <field id="ahmf_3d"      long_name=" 3D      f-eddy viscosity coefficient"   unit="m2/s or m4/s"                           grid_ref="grid_T_3D"/> 
     501 
    466502      <!-- scalar variables available with key_diaar5 --> 
    467503 
    … …  
    10181054      <field field_ref="woce"         name="wo"       long_name="ocean vertical velocity"  /> 
    10191055    </field_group> 
    1020      
     1056 
     1057       <!-- TMB diagnostic output --> 
     1058    <field_group  id="1h_grid_T_tmb" grid_ref="grid_T_2D" operation="instant"> 
     1059     <field id="top_temp"           name="votemper_top"  unit="degC"  /> 
     1060     <field id="mid_temp"           name="votemper_mid"  unit="degC"  /> 
     1061     <field id="bot_temp"           name="votemper_bot"  unit="degC"  /> 
     1062     <field id="top_sal"            name="vosaline_top"  unit="psu"   /> 
     1063     <field id="mid_sal"            name="vosaline_mid"  unit="psu"   /> 
     1064     <field id="bot_sal"            name="vosaline_bot"  unit="psu"   /> 
     1065     <field id="sshnmasked"         name="sossheig"      unit="m"     />  
     1066   </field_group> 
     1067    <field_group  id="1h_grid_U_tmb" grid_ref="grid_U_2D" operation="instant"> 
     1068     <field id="top_u"           name="vozocrtx_top"  unit="m/s"  /> 
     1069     <field id="mid_u"           name="vozocrtx_mid"  unit="m/s"  /> 
     1070     <field id="bot_u"           name="vozocrtx_bot"  unit="m/s"  /> 
     1071     <field id="baro_u"          name="vobtcrtx"      unit="m/s"  /> 
     1072   </field_group> 
     1073    <field_group  id="1h_grid_V_tmb" grid_ref="grid_V_2D" operation="instant"> 
     1074     <field id="top_v"           name="vomecrty_top"  unit="m/s"  /> 
     1075     <field id="mid_v"           name="vomecrty_mid"  unit="m/s"  /> 
     1076     <field id="bot_v"           name="vomecrty_bot"  unit="m/s"  /> 
     1077     <field id="baro_v"          name="vobtcrty"      unit="m/s"  /> 
     1078   </field_group> 
     1079       <!-- 25h diagnostic output --> 
     1080      <field_group id="25h_grid_T" grid_ref="grid_T_3D" operation="instant"> 
     1081   <field id="temper25h"         name="potential temperature 25h mean"    unit="degC" /> 
     1082   <field id="tempis25h"         name="insitu temperature 25h mean"    unit="degC" /> 
     1083   <field id="salin25h"          name="salinity 25h mean"                 unit="psu"  /> 
     1084   <field id="ssh25h"            name="sea surface height 25h mean"  grid_ref="grid_T_2D"      unit="m"    /> 
     1085      </field_group> 
     1086 
     1087      <field_group id="25h_grid_U" grid_ref="grid_U_3D" operation="instant" > 
     1088   <field id="vozocrtx25h"         name="i current 25h mean"    unit="m/s"   /> 
     1089      </field_group> 
     1090 
     1091      <field_group id="25h_grid_V" grid_ref="grid_V_3D" operation="instant"> 
     1092   <field id="vomecrty25h"         name="j current 25h mean"    unit="m/s"    /> 
     1093      </field_group> 
     1094 
     1095      <field_group id="25h_grid_W" grid_ref="grid_W_3D" operation="instant"> 
     1096   <field id="vomecrtz25h"         name="k current 25h mean"                 unit="m/s"      /> 
     1097   <field id="avt25h"             name="vertical diffusivity25h mean"       unit="m2/s" /> 
     1098   <field id="avm25h"              name="vertical viscosity 25h mean"        unit="m2/s" /> 
     1099   <field id="tke25h"              name="turbulent kinetic energy 25h mean" /> 
     1100   <field id="mxln25h"             name="mixing length 25h mean"             unit="m" /> 
     1101      </field_group> 
     1102 
     1103 
    10211104 
    10221105   </field_definition> 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/SHARED/namelist_ref

    r6401 r6404  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!!                            namelist_ref 
    13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    24!! NEMO/OPA  :  1 - run manager      (namrun) 
    … …  
    46!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    57!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    6 !!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     8!!                                    namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave) 
     9!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    810!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    9 !!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     11!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
    1012!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    11 !!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
     13!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 
    1214!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    13 !!             10 - miscellaneous    (namsol, nammpp, namctl) 
     15!!             10 - miscellaneous    (nammpp, namctl) 
    1416!!             11 - Obs & Assim      (namobs, nam_asminc) 
    1517!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    … …  
    2931   nn_itend    =    5475   !  last  time step (std 5475) 
    3032   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
     33   nn_time0    =       0   !  initial time of day in hhmm 
    3134   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    3235   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    33    nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=T 
    34    nn_rstctl   =       0   !  restart control ==> activated only if ln_rstart=T 
    35                            !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    36                            !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
    37                            !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
    38    cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input) 
    39    cn_ocerst_indir = "."       !  directory from which to read input ocean restarts 
    40    cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
    41    cn_ocerst_outdir = "."      !  directory in which to write output ocean restarts 
     36      nn_euler    =    1            !  = 0 : start with forward time step if ln_rstart=T 
     37      nn_rstctl   =    0            !  restart control ==> activated only if ln_rstart=T 
     38      !                             !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
     39      !                             !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     40      !                             !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
     41      cn_ocerst_in    = "restart"   !  suffix of ocean restart name (input) 
     42      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts 
     43      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output) 
     44      cn_ocerst_outdir= "."         !  directory in which to write output ocean restarts 
     45   ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model 
    4246   nn_istate   =       0   !  output the initial state (1) or not (0) 
    4347   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
    … …  
    4549   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written 
    4650   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    47    ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4851   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4952   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard 
    50    ln_clobber  = .true.   !  clobber (overwrite) an existing file 
     53   ln_clobber  = .true.    !  clobber (overwrite) an existing file 
    5154   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    5255/ 
    … …  
    5962!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    6063!!   namdom       space and time domain (bathymetry, mesh, timestep) 
     64!!   namcrs       coarsened grid (for outputs and/or TOP)               ("key_crs") 
     65!!   namc1d       1D configuration options                              ("key_c1d") 
     66!!   namc1d_uvd   1D data (currents)                                    ("key_c1d") 
     67!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
    6168!!   namtsd       data: temperature & salinity 
    6269!!====================================================================== 
    6370! 
    6471!----------------------------------------------------------------------- 
    65 &namcfg     !   parameters of the configuration 
    66 !----------------------------------------------------------------------- 
    67    cp_cfg      =  "default"            !  name of the configuration 
    68    cp_cfz      =  "no zoom"            !  name of the zoom of configuration 
    69    jp_cfg      =       0               !  resolution of the configuration 
    70    jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
    71    jpjdta      =      12               !  2nd    "         "    ( >= jpj ) 
    72    jpkdta      =      31               !  number of levels      ( >= jpk ) 
    73    jpiglo      =      10               !  1st dimension of global domain --> i =jpidta 
    74    jpjglo      =      12               !  2nd    -                  -    --> j =jpjdta 
    75    jpizoom     =       1               !  left bottom (i,j) indices of the zoom 
    76    jpjzoom     =       1               !  in data domain indices 
    77    jperio      =       0               !  lateral cond. type (between 0 and 6) 
    78                                        !  = 0 closed                 ;   = 1 cyclic East-West 
    79                                        !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
    80                                        !  = 4 cyclic East-West AND North fold T-point pivot 
    81                                        !  = 5 North fold F-point pivot 
    82                                        !  = 6 cyclic East-West AND North fold F-point pivot 
    83    ln_use_jattr = .false.              !  use (T) the file attribute: open_ocean_jstart, if present 
    84                                        !  in netcdf input files, as the start j-row for reading 
    85 / 
    86 !----------------------------------------------------------------------- 
    87 &namzgr        !   vertical coordinate 
    88 !----------------------------------------------------------------------- 
    89    ln_zco      = .false.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
    90    ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    91    ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
    92    ln_isfcav   = .false.   !  ice shelf cavity               (T/F) 
     72&namcfg        !   parameters of the configuration 
     73!----------------------------------------------------------------------- 
     74   cp_cfg      = "default" !  name of the configuration 
     75   cp_cfz      = "no zoom" !  name of the zoom of configuration 
     76   jp_cfg      =      0    !  resolution of the configuration 
     77   jpidta      =     10    !  1st lateral dimension ( >= jpi ) 
     78   jpjdta      =     12    !  2nd    "         "    ( >= jpj ) 
     79   jpkdta      =     31    !  number of levels      ( >= jpk ) 
     80   jpiglo      =     10    !  1st dimension of global domain --> i =jpidta 
     81   jpjglo      =     12    !  2nd    -                  -    --> j =jpjdta 
     82   jpizoom     =      1    !  left bottom (i,j) indices of the zoom 
     83   jpjzoom     =      1    !  in data domain indices 
     84   jperio      =      0    !  lateral cond. type (between 0 and 6) 
     85                                 !  = 0 closed                 ;   = 1 cyclic East-West 
     86                                 !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
     87                                 !  = 4 cyclic East-West AND North fold T-point pivot 
     88                                 !  = 5 North fold F-point pivot 
     89                                 !  = 6 cyclic East-West AND North fold F-point pivot 
     90   ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present 
     91                           !  in netcdf input files, as the start j-row for reading 
     92/ 
     93!----------------------------------------------------------------------- 
     94&namzgr        !   vertical coordinate                                  (default: NO selection) 
     95!----------------------------------------------------------------------- 
     96   ln_zco      = .false.   !  z-coordinate - full    steps 
     97   ln_zps      = .false.   !  z-coordinate - partial steps 
     98   ln_sco      = .false.   !  s- or hybrid z-s-coordinate 
     99   ln_isfcav   = .false.   !  ice shelf cavity 
     100   ln_linssh   = .false.   !  linear free surface 
    93101/ 
    94102!----------------------------------------------------------------------- 
    95103&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate 
    96104!----------------------------------------------------------------------- 
    97    ln_s_sh94   = .true.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
     105   ln_s_sh94   = .false.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
    98106   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
    99107   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
    … …  
    125133   nn_msh      =    1      !  create (=1) a mesh file or not (=0) 
    126134   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
     135   rn_isfhmin  =    1.00   !  treshold (m) to discriminate grounding ice to floating ice 
    127136   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
    128137   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
    … …  
    130139   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
    131140   rn_atfp     =    0.1    !  asselin time filter parameter 
    132    nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
    133                                  !                          =0, not used, rdt = rdttra 
    134    rn_rdtmin   = 28800.          !  minimum time step on tracers (used if nn_acc=1) 
    135    rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nn_acc=1) 
    136    rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
    137141   ln_crs      = .false.      !  Logical switch for coarsening module 
    138142   jphgr_msh   =       0               !  type of horizontal mesh 
    … …  
    161165/ 
    162166!----------------------------------------------------------------------- 
    163 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    164 !----------------------------------------------------------------------- 
    165    ln_bt_fw      =    .TRUE.           !  Forward integration of barotropic equations 
    166    ln_bt_av      =    .TRUE.           !  Time filtering of barotropic variables 
    167    ln_bt_nn_auto =    .TRUE.           !  Set nn_baro automatically to be just below 
    168                                        !  a user defined maximum courant number (rn_bt_cmax) 
    169    nn_baro       =    30               !  Number of iterations of barotropic mode 
    170                                        !  during rn_rdt seconds. Only used if ln_bt_nn_auto=F 
    171    rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T 
    172    nn_bt_flt     =    1                !  Time filter choice 
    173                                        !  = 0 None 
    174                                        !  = 1 Boxcar over   nn_baro barotropic steps 
    175                                        !  = 2 Boxcar over 2*nn_baro     "        " 
    176 / 
    177 !----------------------------------------------------------------------- 
    178 &namcrs        !   Grid coarsening for dynamics output and/or 
    179                !   passive tracer coarsened online simulations 
     167&namcrs        !   coarsened grid (for outputs and/or TOP)              ("key_crs") 
    180168!----------------------------------------------------------------------- 
    181169   nn_factx    = 3         !  Reduction factor of x-direction 
    … …  
    199187/ 
    200188!----------------------------------------------------------------------- 
     189&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d") 
     190!----------------------------------------------------------------------- 
     191   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
     192/ 
     193!----------------------------------------------------------------------- 
     194&namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
     195!----------------------------------------------------------------------- 
     196!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     197!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     198   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , '' 
     199   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , '' 
     200! 
     201   cn_dir        = './'    !  root directory for the location of the files 
     202   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F) 
     203   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F) 
     204/ 
     205!----------------------------------------------------------------------- 
    201206&namtsd    !   data : Temperature  & Salinity 
    202 !----------------------------------------------------------------------- 
    203207!----------------------------------------------------------------------- 
    204208!          !  file name                            ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    … …  
    211215   ln_tsd_tradmp = .true.   !  damping of ocean T & S toward T &S input data (T) or not (F) 
    212216/ 
     217 
    213218!!====================================================================== 
    214219!!            ***  Surface Boundary Condition namelists  *** 
    215220!!====================================================================== 
    216221!!   namsbc          surface boundary condition 
    217 !!   namsbc_ana      analytical         formulation 
    218 !!   namsbc_flx      flux               formulation 
    219 !!   namsbc_clio     CLIO bulk formulae formulation 
    220 !!   namsbc_core     CORE bulk formulae formulation 
    221 !!   namsbc_mfs      MFS  bulk formulae formulation 
    222 !!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3") 
     222!!   namsbc_ana      analytical         formulation                     (ln_ana     =T) 
     223!!   namsbc_flx      flux               formulation                     (ln_flx     =T) 
     224!!   namsbc_clio     CLIO bulk formulae formulation                     (ln_blk_clio=T) 
     225!!   namsbc_core     CORE bulk formulae formulation                     (ln_blk_core=T) 
     226!!   namsbc_mfs      MFS  bulk formulae formulation                     (ln_blk_mfs =T) 
     227!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" ) 
    223228!!   namsbc_sas      StAndalone Surface module 
    224 !!   namtra_qsr      penetrative solar radiation 
    225 !!   namsbc_rnf      river runoffs 
    226 !!   namsbc_isf      ice shelf melting/freezing 
    227 !!   namsbc_apr      Atmospheric Pressure 
    228 !!   namsbc_ssr      sea surface restoring term (for T and/or S) 
     229!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T) 
     230!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
     231!!   namsbc_isf      ice shelf melting/freezing                         (nn_isf     >0) 
     232!!   namsbc_iscpl    coupling option between land ice model and ocean 
     233!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
     234!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T) 
    229235!!   namsbc_alb      albedo parameters 
     236!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
     237!!   namberg         iceberg floats                                     (ln_icebergs=T) 
    230238!!====================================================================== 
    231239! 
    … …  
    234242!----------------------------------------------------------------------- 
    235243   nn_fsbc     = 5         !  frequency of surface boundary condition computation 
    236                            !     (also = the frequency of sea-ice model call) 
     244                           !     (also = the frequency of sea-ice & iceberg model call) 
     245                     ! Type of air-sea fluxes  
    237246   ln_ana      = .false.   !  analytical formulation                    (T => fill namsbc_ana ) 
    238247   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    … …  
    240249   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    241250   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
     251                     ! Type of coupling (Ocean/Ice/Atmosphere) : 
    242252   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
    243253   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 ) 
    … …  
    246256                           !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component 
    247257                           !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component  
    248    ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    249    nn_ice      = 2         !  =0 no ice boundary condition   , 
    250                            !  =1 use observed ice-cover      , 
    251                            !  =2 ice-model used                         ("key_lim3" or "key_lim2") 
    252    nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    253                            !  =1 levitating ice with mass and salt exchange but no presure effect 
    254                            !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    255    ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    256    ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf) 
    257    nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf) 
    258                            !  0 =no isf                  1 = presence of ISF 
    259                            !  2 = bg03 parametrisation   3 = rnf file for isf 
    260                            !  4 = ISF fwf specified 
    261                            !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    262    ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    263    nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
    264                            !     =1 global mean of e-p-r set to zero at each time step 
    265                            !     =2 annual global mean of e-p-r set to zero 
    266    ln_wave = .false.       !  Activate coupling with wave (either Stokes Drift or Drag coefficient, or both)  (T => fill namsbc_wave) 
    267    ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => fill namsbc_wave) 
    268    ln_sdw  = .false.       !  Computation of 3D stokes drift                (T => fill namsbc_wave) 
    269    nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
    270                            !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
    271258   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
    272259                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
    … …  
    274261                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
    275262                           !  = 2  Redistribute a single flux over categories (coupled mode only) 
     263                     ! Sea-ice : 
     264   nn_ice      = 2         !  =0 no ice boundary condition   , 
     265                           !  =1 use observed ice-cover      , 
     266                           !  =2 ice-model used                         ("key_lim3", "key_lim2", "key_cice") 
     267   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
     268                           !  =1 levitating ice with mass and salt exchange but no presure effect 
     269                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
     270                     ! Misc. options of sbc :  
     271   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr ) 
     272   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
     273   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
     274   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
     275   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
     276                           !     =1 global mean of e-p-r set to zero at each time step 
     277                           !     =2 annual global mean of e-p-r set to zero 
     278   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
     279   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf) 
     280   ln_wave = .false.       !  coupling with surface wave                (T => fill namsbc_wave) 
     281   nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
     282                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
    276283/ 
    277284!----------------------------------------------------------------------- 
    … …  
    316323&namsbc_core   !   namsbc_core  CORE bulk formulae 
    317324!----------------------------------------------------------------------- 
    318 !              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
    319 !              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
    320    sn_wndi     = 'u_10.15JUNE2009_fill'        ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , '' 
    321    sn_wndj     = 'v_10.15JUNE2009_fill'        ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , '' 
    322    sn_qsr      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    323    sn_qlw      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    324    sn_tair     = 't_10.15JUNE2009_fill'        ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    325    sn_humi     = 'q_10.15JUNE2009_fill'        ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    326    sn_prec     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    327    sn_snow     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    328    sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     325!              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
     326!              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
     327   sn_wndi     = 'u_10.15JUNE2009_fill'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , '' 
     328   sn_wndj     = 'v_10.15JUNE2009_fill'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , '' 
     329   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     330   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     331   sn_tair     = 't_10.15JUNE2009_fill'       ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     332   sn_humi     = 'q_10.15JUNE2009_fill'       ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     333   sn_prec     = 'ncar_precip.15JUNE2009_fill',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     334   sn_snow     = 'ncar_precip.15JUNE2009_fill',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     335   sn_tdif     = 'taudif_core'                ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    329336 
    330337   cn_dir      = './'      !  root directory for the location of the bulk files 
    331338   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    332    rn_zqt      = 10.        !  Air temperature and humidity reference height (m) 
    333    rn_zu       = 10.        !  Wind vector reference height (m) 
     339   rn_zqt      = 10.       !  Air temperature and humidity reference height (m) 
     340   rn_zu       = 10.       !  Wind vector reference height (m) 
    334341   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    335342   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    … …  
    340347&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae 
    341348!----------------------------------------------------------------------- 
    342 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask ! 
    343 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      ! 
    344    sn_wndi     =   'ecmwf'   ,        6          , 'u10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
    345    sn_wndj     =   'ecmwf'   ,        6          , 'v10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
    346    sn_clc      =   'ecmwf'   ,        6          , 'clc'     ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , '' 
    347    sn_msl      =   'ecmwf'   ,        6          , 'msl'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
    348    sn_tair     =   'ecmwf'   ,        6          , 't2'      ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
    349    sn_rhm      =   'ecmwf'   ,        6          , 'rh'      ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , '' 
    350    sn_prec     =   'ecmwf'   ,        6          , 'precip'  ,    .true.    , .true.  , 'daily'  ,'bicubic.nc' , ''       , '' 
     349!              !  file name  ! frequency (hours) ! variable ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask ! 
     350!              !             !  (if <0  months)  !   name   !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      ! 
     351   sn_wndi     =   'ecmwf'   ,        6          , 'u10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,  '' 
     352   sn_wndj     =   'ecmwf'   ,        6          , 'v10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,  '' 
     353   sn_clc      =   'ecmwf'   ,        6          , 'clc'    ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,  '' 
     354   sn_msl      =   'ecmwf'   ,        6          , 'msl'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,  '' 
     355   sn_tair     =   'ecmwf'   ,        6          , 't2'     ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,  '' 
     356   sn_rhm      =   'ecmwf'   ,        6          , 'rh'     ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,  '' 
     357   sn_prec     =   'ecmwf'   ,        6          , 'precip' ,    .true.    , .true. , 'daily'   ,'bicubic.nc' ,   ''     ,  '' 
    351358 
    352359   cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files 
    … …  
    355362&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3") 
    356363!----------------------------------------------------------------------- 
    357 !                    !     description       !  multiple  !    vector   !      vector          ! vector ! 
    358 !                    !                       ! categories !  reference  !    orientation       ! grids  ! 
     364!                    !     description      !  multiple  !    vector   !      vector          ! vector ! 
     365!                    !                      ! categories !  reference  !    orientation       ! grids  ! 
    359366! send 
    360    sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    361    sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    362    sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
    363    sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    364    sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     367   sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     368   sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     369   sn_snd_thick  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     370   sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
     371   sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    365372! receive 
    366    sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    367    sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    368    sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
    369    sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    370    sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    371    sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    372    sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    373    sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    374    sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    375    sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     373   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     374   sn_rcv_taumod =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     375   sn_rcv_tau    =   'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     376   sn_rcv_dqnsdt =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     377   sn_rcv_qsr    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     378   sn_rcv_qns    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     379   sn_rcv_emp    =   'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     380   sn_rcv_rnf    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     381   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     382   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    376383! 
    377384   nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
    … …  
    397404/ 
    398405!----------------------------------------------------------------------- 
    399 &namtra_qsr    !   penetrative solar radiation 
     406&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr=T) 
    400407!----------------------------------------------------------------------- 
    401408!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    … …  
    404411 
    405412   cn_dir      = './'      !  root directory for the location of the runoff files 
    406    ln_traqsr   = .true.    !  Light penetration (T) or not (F) 
    407413   ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration 
    408414   ln_qsr_2bd  = .false.   !  2 bands              light penetration 
    … …  
    415421/ 
    416422!----------------------------------------------------------------------- 
    417 &namsbc_rnf    !   runoffs namelist surface boundary condition 
     423&namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf=T) 
    418424!----------------------------------------------------------------------- 
    419425!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    … …  
    439445/ 
    440446!----------------------------------------------------------------------- 
    441 &namsbc_isf    !  Top boundary layer (ISF) 
    442 !----------------------------------------------------------------------- 
    443 !              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    444 !              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
     447&namsbc_isf    !  Top boundary layer (ISF)                              (nn_isf >0) 
     448!----------------------------------------------------------------------- 
     449!              ! file name ! frequency (hours) ! variable ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     450!              !           !  (if <0  months)  !   name   !  (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    445451! nn_isf == 4 
    446    sn_qisf      = 'rnfisf' ,         -12      ,'sohflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
    447    sn_fwfisf    = 'rnfisf' ,         -12      ,'sowflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     452   sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    448453! nn_isf == 3 
    449    sn_rnfisf    = 'runoffs' ,         -12      ,'sofwfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,  '' 
     454   sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    450455! nn_isf == 2 and 3 
    451    sn_depmax_isf = 'runoffs' ,       -12        ,'sozisfmax' ,   .false.  , .true.  , 'yearly'  ,  ''      ,  '' 
    452    sn_depmin_isf = 'runoffs' ,       -12        ,'sozisfmin' ,   .false.  , .true.  , 'yearly'  ,  ''      ,  '' 
     456   sn_depmax_isf='rnfisf'  ,         -12       ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
     457   sn_depmin_isf='rnfisf'  ,         -12       ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
    453458! nn_isf == 2 
    454    sn_Leff_isf = 'rnfisf' ,       0          ,'Leff'         ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     459   sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , '' 
     460! 
    455461! for all case 
    456    ln_divisf   = .true.  ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) 
     462   nn_isf      = 1         !  ice shelf melting/freezing 
     463                           !  1 = presence of ISF    2 = bg03 parametrisation  
     464                           !  3 = rnf file for isf   4 = ISF fwf specified 
     465                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    457466! only for nn_isf = 1 or 2 
    458    rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula 
    459    rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula 
     467   rn_gammat0  = 1.e-4    ! gammat coefficient used in blk formula 
     468   rn_gammas0  = 1.e-4    ! gammas coefficient used in blk formula 
     469! only for nn_isf = 1 or 4 
     470   rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     471                          ! 0 => thickness of the tbl = thickness of the first wet cell 
    460472! only for nn_isf = 1 
    461    nn_isfblk   =  1       ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 
    462    rn_hisf_tbl =  30.      ! thickness of the top boundary layer           (Losh et al. 2008) 
    463                           ! 0 => thickness of the tbl = thickness of the first wet cell 
    464    ln_conserve = .true.   ! conservative case (take into account meltwater advection) 
     473   nn_isfblk   = 1        ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
     474                          ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    465475   nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s) 
    466476                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    467                           !     if you want to keep the cd as in global config, adjust rn_gammat0 to compensate 
    468                           ! 2 = velocity and stability dependent Gamma    Holland et al. 1999 
     477                          ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
     478/ 
     479!----------------------------------------------------------------------- 
     480&namsbc_iscpl  !   land ice / ocean coupling option                      
     481!----------------------------------------------------------------------- 
     482   nn_drown  = 10       ! number of iteration of the extrapolation loop (fill the new wet cells) 
     483   ln_hsb    = .false.  ! activate conservation module (conservation exact after a time of rn_fiscpl) 
     484   nn_fiscpl = 43800    ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
    469485/ 
    470486!----------------------------------------------------------------------- 
    471487&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk 
    472488!----------------------------------------------------------------------- 
    473 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    474 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    475    sn_apr      = 'patm'      ,         -1        ,'somslpre',    .true.     , .true. , 'yearly'  ,  ''      ,   ''     , '' 
     489!              !  file name ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     490!              !            !  (if <0  months)  !   name    !  (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     491   sn_apr      = 'patm'     ,         -1        ,'somslpre',    .true.     , .true. , 'yearly'  ,  ''      ,   ''     , '' 
    476492 
    477493   cn_dir      = './'       !  root directory for the location of the bulk files 
    … …  
    481497/ 
    482498!----------------------------------------------------------------------- 
    483 &namsbc_ssr    !   surface boundary condition : sea surface restoring 
     499&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr=T) 
    484500!----------------------------------------------------------------------- 
    485501!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    … …  
    507523/ 
    508524!----------------------------------------------------------------------- 
    509 &namberg       !   iceberg parameters 
    510 !----------------------------------------------------------------------- 
    511       ln_icebergs              = .false. 
     525&namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
     526!----------------------------------------------------------------------- 
     527!              !  file name  ! frequency (hours) ! variable    ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     528!              !             !  (if <0  months)  !   name      !  (logical)   !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     529   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff',     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
     530   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'    ,     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
     531   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'    ,     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
     532   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'  ,     .true.   , .false., 'daily'   ,  ''      , ''       , '' 
     533! 
     534   cn_dir_cdg  = './'      !  root directory for the location of drag coefficient files 
     535   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model 
     536   ln_sdw  = .false.       !  Computation of 3D stokes drift                
     537/ 
     538!----------------------------------------------------------------------- 
     539&namberg       !   iceberg parameters                                   (default: No iceberg) 
     540!----------------------------------------------------------------------- 
     541      ln_icebergs              = .false.              ! iceberg floats or not 
    512542      ln_bergdia               = .true.               ! Calculate budgets 
    513543      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
    … …  
    534564      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
    535565 
    536 !              ! file name ! frequency (hours) !   variable   ! time interp.  !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    537 !              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    538       sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , '' 
     566!            ! file name ! frequency (hours) !   variable   ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     567!            !           !  (if <0  months)  !     name     !   (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     568      sn_icb =  'calving',       -1          , 'calvingmask',  .true.      , .true.  , 'yearly'  , ''       , ''       , '' 
    539569 
    540570      cn_dir = './' 
    … …  
    545575!!====================================================================== 
    546576!!   namlbc        lateral momentum boundary condition 
    547 !!   namcla        cross land advection 
    548577!!   namobc        open boundaries parameters                           ("key_obc") 
    549578!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
    … …  
    555584&namlbc        !   lateral momentum boundary condition 
    556585!----------------------------------------------------------------------- 
     586   !                       !  free slip  !   partial slip  !   no slip   ! strong slip 
    557587   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    558                            !  free slip  !   partial slip  !   no slip   ! strong slip 
    559    ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs. 
    560 / 
    561 !----------------------------------------------------------------------- 
    562 &namcla        !   cross land advection 
    563 !----------------------------------------------------------------------- 
    564    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    565 / 
    566 !----------------------------------------------------------------------- 
    567 &namobc        !   open boundaries parameters                           ("key_obc") 
    568 !----------------------------------------------------------------------- 
    569    ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    570    ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F) 
    571    ln_obc_fla  = .false.   !  Flather open boundary condition 
    572    nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    573                            !  = 1 the obc data are read in 'obc.dta' files 
    574    cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data 
    575                            !  set to monthly if obc datafile hold 1 month of data 
    576    rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary 
    577    rn_dpwin    =    1.     !     -           -         -       west    -      - 
    578    rn_dpnin    =    1.     !     -           -         -       north   -      - 
    579    rn_dpsin    =    1.     !     -           -         -       south   -      - 
    580    rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary 
    581    rn_dpwob    =   15.     !     -           -         -     west    -      - 
    582    rn_dpnob    = 3000.     !     -           -         -     north   -      - 
    583    rn_dpsob    =   15.     !     -           -         -     south   -      - 
    584    rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R) 
    585                            !  = 1 the total volume remains constant 
     588   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs. 
    586589/ 
    587590!----------------------------------------------------------------------- 
    … …  
    595598/ 
    596599!----------------------------------------------------------------------- 
    597 &nam_tide      !   tide parameters (#ifdef key_tide) 
     600&nam_tide      !   tide parameters                                      ("key_tide") 
    598601!----------------------------------------------------------------------- 
    599602   ln_tide_pot   = .true.   !  use tidal potential forcing 
    … …  
    637640/ 
    638641!----------------------------------------------------------------------- 
    639 &nambdy_dta      !  open boundaries - external data           ("key_bdy") 
    640 !----------------------------------------------------------------------- 
    641 !              !  file name      ! frequency (hours) ! variable   ! time interp.  !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    642 !              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    643    bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    644    bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    645    bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    646    bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    647    bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    648    bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    649    bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     642&nambdy_dta    !  open boundaries - external data                       ("key_bdy") 
     643!----------------------------------------------------------------------- 
     644!              !  file name      ! frequency (hours) ! variable  ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     645!              !                 !  (if <0  months)  !   name    !  (logical)   !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     646   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     647   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     648   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     649   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     650   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     651   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     652   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    650653! for lim2 
    651 !   bn_frld  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    652 !   bn_hicif =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    653 !   bn_hsnif =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     654!   bn_frld  =   'amm12_bdyT_ice' ,         24        , 'ileadfra',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     655!   bn_hicif =   'amm12_bdyT_ice' ,         24        , 'iicethic',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     656!   bn_hsnif =   'amm12_bdyT_ice' ,         24        , 'isnowthi',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    654657! for lim3 
    655 !   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    656 !   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    657 !   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    658    cn_dir  =    'bdydta/' 
    659    ln_full_vel = .false. 
     658!   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     659!   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     660!   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi',     .true.   , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     661 
     662   cn_dir      =    'bdydta/'   !  root directory for the location of the bulk files 
     663   ln_full_vel = .false.        !   
    660664/ 
    661665!----------------------------------------------------------------------- 
    662666&nambdy_tide     ! tidal forcing at open boundaries 
    663667!----------------------------------------------------------------------- 
    664    filtide          = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
    665    ln_bdytide_2ddta = .false. 
    666    ln_bdytide_conj  = .false. 
     668   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
     669   ln_bdytide_2ddta = .false.                   ! 
     670   ln_bdytide_conj  = .false.                   !  
    667671/ 
    668672!!====================================================================== 
    … …  
    675679! 
    676680!----------------------------------------------------------------------- 
    677 &nambfr        !   bottom friction 
     681&nambfr        !   bottom friction                                      (default: linear) 
    678682!----------------------------------------------------------------------- 
    679683   nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    … …  
    698702/ 
    699703!----------------------------------------------------------------------- 
    700 &nambbc        !   bottom temperature boundary condition 
    701 !----------------------------------------------------------------------- 
    702 !              !                              !  (if <0  months)  !   
    703 !              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    704 !              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    705    sn_qgh      ='geothermal_heating.nc',  -12.  , 'heatflow'      ,   .false.      , .true.  , 'yearly'  , ''       , ''       , '' 
     704&nambbc        !   bottom temperature boundary condition                (default: NO) 
     705!----------------------------------------------------------------------- 
     706!              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     707!              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     708   sn_qgh      ='geothermal_heating.nc',  -12.       , 'heatflow',   .false.   , .true.  , 'yearly'  ,   ''     ,   ''     ,   '' 
    706709   ! 
    707    cn_dir      = './'      !  root directory for the location of the runoff files 
    708    ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
     710   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    709711   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux 
    710712                           !     = 1 constant flux 
    711713                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    712714   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    713  
    714 / 
    715 !----------------------------------------------------------------------- 
    716 &nambbl        !   bottom boundary layer scheme 
     715   cn_dir      = './'      !  root directory for the location of the runoff files  
     716/ 
     717!----------------------------------------------------------------------- 
     718&nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    717719!----------------------------------------------------------------------- 
    718720   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0) 
    … …  
    725727!!                        Tracer (T & S ) namelists 
    726728!!====================================================================== 
    727 !!   nameos        equation of state 
    728 !!   namtra_adv    advection scheme 
     729!!   nameos           equation of state 
     730!!   namtra_adv       advection scheme 
    729731!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    730 !!   namtra_ldf    lateral diffusion scheme 
    731 !!   namtra_dmp    T & S newtonian damping 
     732!!   namtra_ldf       lateral diffusion scheme 
     733!!   namtra_ldfeiv    eddy induced velocity param. 
     734!!   namtra_dmp       T & S newtonian damping 
    732735!!====================================================================== 
    733736! 
    … …  
    740743                                 !  = 1, S-EOS   (simplified eos) 
    741744   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
    742    !                             ! 
     745                                 ! 
    743746   !                     ! S-EOS coefficients : 
    744    !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     747                                 !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    745748   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
    746749   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
    … …  
    752755/ 
    753756!----------------------------------------------------------------------- 
    754 &namtra_adv    !   advection scheme for tracer 
    755 !----------------------------------------------------------------------- 
    756    ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
    757    ln_traadv_tvd    =  .true.    !  TVD scheme 
    758    ln_traadv_muscl  =  .false.   !  MUSCL scheme 
    759    ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
    760    ln_traadv_ubs    =  .false.   !  UBS scheme 
    761    ln_traadv_qck    =  .false.   !  QUICKEST scheme 
    762    ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
    763    ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection 
    764 / 
    765 !----------------------------------------------------------------------- 
    766 &namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
    767 !----------------------------------------------------------------------- 
    768    ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     757&namtra_adv    !   advection scheme for tracer                          (default: NO advection) 
     758!----------------------------------------------------------------------- 
     759   ln_traadv_cen =  .false.  !  2nd order centered scheme 
     760      nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     761      nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
     762   ln_traadv_fct =  .false.  !  FCT scheme 
     763      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     764      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     765      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     766      !                             !        (number of sub-timestep = nn_fct_zts) 
     767   ln_traadv_mus =  .false.  !  MUSCL scheme 
     768      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
     769   ln_traadv_ubs =  .false.  !  UBS scheme 
     770      nn_ubs_v   =  2               !  =2  , vertical 2nd order FCT / COMPACT 4th order 
     771   ln_traadv_qck =  .false.  !  QUICKEST scheme 
     772/ 
     773!----------------------------------------------------------------------- 
     774&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) 
     775!----------------------------------------------------------------------- 
     776   ln_mle    = .false.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
    769777   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
    770778   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
    … …  
    776784   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK 
    777785/ 
    778 !---------------------------------------------------------------------------------- 
    779 &namtra_ldf    !   lateral diffusion scheme for tracers 
    780 !---------------------------------------------------------------------------------- 
     786!----------------------------------------------------------------------- 
     787&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO diffusion) 
     788!----------------------------------------------------------------------- 
    781789   !                       !  Operator type: 
    782    ln_traldf_lap    =  .true.   !  laplacian operator 
    783    ln_traldf_bilap  =  .false.  !  bilaplacian operator 
     790   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     791   ln_traldf_lap   =  .false.  !    laplacian operator 
     792   ln_traldf_blp   =  .false.  !  bilaplacian operator 
    784793   !                       !  Direction of action: 
    785    ln_traldf_level  =  .false.  !  iso-level 
    786    ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    787    ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp") 
    788    !                 !  Griffies parameters              (all need "key_ldfslp") 
    789    ln_traldf_grif   =  .false.  !  use griffies triads 
    790    ln_traldf_gdia   =  .false.  !  output griffies eddy velocities 
    791    ln_triad_iso     =  .false.  !  pure lateral mixing in ML 
    792    ln_botmix_grif   =  .false.  !  lateral mixing on bottom 
    793    !                       !  Coefficients 
    794    ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 
    795    ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 
    796    !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 
    797    rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s] 
    798    rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    799    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    800    !                                           (normally=0; not used with Griffies) 
    801    rn_slpmax        =     0.01  !  slope limit 
    802    rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity 
    803    rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
    804    rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
    805 / 
    806 !----------------------------------------------------------------------- 
    807 &namtra_dmp    !   tracer: T & S newtonian damping 
     794   ln_traldf_lev   =  .false.  !  iso-level 
     795   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     796   ln_traldf_iso   =  .false.  !  iso-neutral (standard operator) 
     797   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     798   ! 
     799   !                       !  iso-neutral options:         
     800   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
     801   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     802   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     803   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     804   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     805   ! 
     806   !                       !  Coefficients: 
     807   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     808   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     809   !                                !   =  0           constant  
     810   !                                !   = 10 F(k)      =ldf_c1d  
     811   !                                !   = 20 F(i,j)    =ldf_c2d  
     812   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     813   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     814   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     815   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     816   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     817/ 
     818!----------------------------------------------------------------------- 
     819&namtra_ldfeiv !   eddy induced velocity param.                         (default: NO) 
     820!----------------------------------------------------------------------- 
     821   ln_ldfeiv     =.false.  ! use eddy induced velocity parameterization 
     822   ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities 
     823   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     824   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     825   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     826   !                                !   =  0           constant  
     827   !                                !   = 10 F(k)      =ldf_c1d  
     828   !                                !   = 20 F(i,j)    =ldf_c2d  
     829   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     830   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     831/ 
     832!----------------------------------------------------------------------- 
     833&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO) 
    808834!----------------------------------------------------------------------- 
    809835   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F) 
    … …  
    811837                           !                   =1 no damping in the mixing layer (kz  criteria) 
    812838                           !                   =2 no damping in the mixed  layer (rho crieria) 
    813    cn_resto    = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this) 
     839   cn_resto    ='resto.nc' !  Name of file containing restoration coeff. field (use dmp_tools to create this) 
    814840/ 
    815841 
    … …  
    820846!!   namdyn_vor    advection scheme 
    821847!!   namdyn_hpg    hydrostatic pressure gradient 
    822 !!   namdyn_spg    surface pressure gradient                            (CPP key only) 
     848!!   namdyn_spg    surface pressure gradient 
    823849!!   namdyn_ldf    lateral diffusion scheme 
    824850!!====================================================================== 
    825851! 
    826852!----------------------------------------------------------------------- 
    827 &namdyn_adv    !   formulation of the momentum advection 
     853&namdyn_adv    !   formulation of the momentum advection                (default: vector form) 
    828854!----------------------------------------------------------------------- 
    829855   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
    … …  
    834860/ 
    835861!----------------------------------------------------------------------- 
    836 &nam_vvl    !   vertical coordinate options 
     862&nam_vvl    !   vertical coordinate options                             (default: zstar) 
    837863!----------------------------------------------------------------------- 
    838864   ln_vvl_zstar  = .true.           !  zstar vertical coordinate 
    … …  
    848874/ 
    849875!----------------------------------------------------------------------- 
    850 &namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
     876&namdyn_vor    !   option of physics/algorithm                          (default: NO) 
    851877!----------------------------------------------------------------------- 
    852878   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
    853879   ln_dynvor_ens = .false. !  energy conserving scheme 
    854880   ln_dynvor_mix = .false. !  mixed scheme 
    855    ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    856    ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation 
    857 / 
    858 !----------------------------------------------------------------------- 
    859 &namdyn_hpg    !   Hydrostatic pressure gradient option 
     881   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     882      nn_een_e3f = 1          ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     883   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes)  ! PLEASE DO NOT ACTIVATE 
     884/ 
     885!----------------------------------------------------------------------- 
     886&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: zps) 
    860887!----------------------------------------------------------------------- 
    861888   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    … …  
    865892   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    866893   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    867    ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    868                                  !           centered      time scheme  (F) 
    869 / 
    870 !----------------------------------------------------------------------- 
    871 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    872 !----------------------------------------------------------------------- 
    873 !                          !  explicit free surface                     ("key_dynspg_exp") 
    874 !                          !  filtered free surface                     ("key_dynspg_flt") 
    875 !                          !  split-explicit free surface               ("key_dynspg_ts") 
    876  
    877 !----------------------------------------------------------------------- 
    878 &namdyn_ldf    !   lateral diffusion on momentum 
     894/ 
     895!----------------------------------------------------------------------- 
     896&namdyn_spg    !   surface pressure gradient                            (default: NO) 
     897!----------------------------------------------------------------------- 
     898   ln_dynspg_exp  = .false.   ! explicit free surface 
     899   ln_dynspg_ts   = .false.   ! split-explicit free surface 
     900      ln_bt_fw      = .true.     ! Forward integration of barotropic Eqs. 
     901      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
     902         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
     903         !                          !                     = 1 Boxcar over   nn_baro sub-steps 
     904         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     905      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
     906         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed 
     907         nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds 
     908/ 
     909!----------------------------------------------------------------------- 
     910&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO) 
    879911!----------------------------------------------------------------------- 
    880912   !                       !  Type of the operator : 
    881    ln_dynldf_lap    =  .true.   !  laplacian operator 
    882    ln_dynldf_bilap  =  .false.  !  bilaplacian operator 
     913   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     914   ln_dynldf_lap =  .false.    !    laplacian operator 
     915   ln_dynldf_blp =  .false.    !  bilaplacian operator 
    883916   !                       !  Direction of action  : 
    884    ln_dynldf_level  =  .false.  !  iso-level 
    885    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    886    ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     917   ln_dynldf_lev =  .false.    !  iso-level 
     918   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     919   ln_dynldf_iso =  .false.    !  iso-neutral 
    887920   !                       !  Coefficient 
    888    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    889    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    890    rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
    891    rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
    892    rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
    893    rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    894    rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    895    rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
     921   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     922   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     923   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     924   !                                !  =  0  constant  
     925   !                                !  = 10  F(k)=c1d 
     926   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     927   !                                !  = 30  F(i,j,k)=c2d*c1d 
     928   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     929   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     930   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     931   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     932   ! 
     933   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    896934/ 
    897935 
    … …  
    902940!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric") 
    903941!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke") 
    904 !!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp") 
     942!!    namzdf_gls    GLS vertical mixing                                 ("key_zdfgls") 
    905943!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm") 
    906944!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx") 
    … …  
    934972   rn_wtmix    =  10.0     !  vertical eddy viscosity coeff [m2/s] in the mixed-layer 
    935973   rn_wvmix    =  10.0     !  vertical eddy diffusion coeff [m2/s] in the mixed-layer 
    936    ln_mldw     = .true.    !  Flag to use or not the mized layer depth param. 
     974   ln_mldw     = .true.    !  Flag to use or not the mixed layer depth param. 
    937975/ 
    938976!----------------------------------------------------------------------- 
    … …  
    9641002                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    9651003/ 
    966 !------------------------------------------------------------------------ 
    967 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    968 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    969    ln_kpprimix = .true.    !  shear instability mixing 
    970    rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s] 
    971    rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s] 
    972    rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability 
    973    rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s] 
    974    rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2] 
    975    rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s] 
    976    nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv 
    977    nn_ave      =  1        !  constant (=0) or profile (=1) background on avt 
    978 / 
    9791004!----------------------------------------------------------------------- 
    9801005&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    … …  
    10151040!!                  ***  Miscellaneous namelists  *** 
    10161041!!====================================================================== 
    1017 !!   namsol            elliptic solver / island / free surface 
    10181042!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    10191043!!   namctl            Control prints & Benchmark 
    1020 !!   namc1d            1D configuration options                         ("key_c1d") 
    1021 !!   namc1d_uvd        data: U & V currents                             ("key_c1d") 
    1022 !!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d") 
    10231044!!   namsto            Stochastic parametrization of EOS 
    10241045!!====================================================================== 
    10251046! 
    1026 !----------------------------------------------------------------------- 
    1027 &namsol        !   elliptic solver / island / free surface 
    1028 !----------------------------------------------------------------------- 
    1029    nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    1030                            !                   =2 successive-over-relaxation (sor) 
    1031    nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test 
    1032    rn_eps      =  1.e-6    !  absolute precision of the solver 
    1033    nn_nmin     =    300    !  minimum of iterations for the SOR solver 
    1034    nn_nmax     =    800    !  maximum of iterations for the SOR solver 
    1035    nn_nmod     =     10    !  frequency of test for the SOR solver 
    1036    rn_resmax   =  1.e-10   !  absolute precision for the SOR solver 
    1037    rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    1038 / 
    10391047!----------------------------------------------------------------------- 
    10401048&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    … …  
    10621070                           !     (no physical validity of the results) 
    10631071   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    1064 / 
    1065 !----------------------------------------------------------------------- 
    1066 &namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
    1067 !----------------------------------------------------------------------- 
    1068 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    1069 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    1070    sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , '' 
    1071    sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , '' 
    1072 ! 
    1073    cn_dir        = './'    !  root directory for the location of the files 
    1074    ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F) 
    1075    ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F) 
    1076 / 
    1077 !----------------------------------------------------------------------- 
    1078 &namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d") 
    1079 !----------------------------------------------------------------------- 
    1080    ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    1081 / 
    1082 !----------------------------------------------------------------------- 
    1083 &namsto       ! Stochastic parametrization of EOS 
    1084 !----------------------------------------------------------------------- 
    1085    ln_rststo = .false.           ! start from mean parameter (F) or from restart file (T) 
     1072   nn_diacfl   =    0      !  Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 
     1073/ 
     1074!----------------------------------------------------------------------- 
     1075&namsto       ! Stochastic parametrization of EOS                       (default: NO) 
     1076!----------------------------------------------------------------------- 
     1077   ln_sto_eos   = .false.  ! stochastic equation of state 
     1078   nn_sto_eos   = 1        ! number of independent random walks 
     1079   rn_eos_stdxy = 1.4      ! random walk horz. standard deviation (in grid points) 
     1080   rn_eos_stdz  = 0.7      ! random walk vert. standard deviation (in grid points) 
     1081   rn_eos_tcor  = 1440.    ! random walk time correlation (in timesteps) 
     1082   nn_eos_ord   = 1        ! order of autoregressive processes 
     1083   nn_eos_flt   = 0        ! passes of Laplacian filter 
     1084   rn_eos_lim   = 2.0      ! limitation factor (default = 3.0) 
     1085   ln_rststo    = .false.  ! start from mean parameter (F) or from restart file (T) 
    10861086   ln_rstseed = .true.           ! read seed of RNG from restart file 
    10871087   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input) 
    10881088   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output) 
    1089  
    1090    ln_sto_eos = .false.          ! stochastic equation of state 
    1091    nn_sto_eos = 1                ! number of independent random walks 
    1092    rn_eos_stdxy = 1.4            ! random walk horz. standard deviation (in grid points) 
    1093    rn_eos_stdz  = 0.7            ! random walk vert. standard deviation (in grid points) 
    1094    rn_eos_tcor  = 1440.0         ! random walk time correlation (in timesteps) 
    1095    nn_eos_ord  = 1               ! order of autoregressive processes 
    1096    nn_eos_flt  = 0               ! passes of Laplacian filter 
    1097    rn_eos_lim  = 2.0             ! limitation factor (default = 3.0) 
    10981089/ 
    10991090 
    … …  
    11011092!!                  ***  Diagnostics namelists  *** 
    11021093!!====================================================================== 
    1103 !!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    11041094!!   namtrd       dynamics and/or tracer trends 
    11051095!!   namptr       Poleward Transport Diagnostics 
     1096!!   namhsb       Heat and salt budgets 
    11061097!!   namflo       float parameters                                      ("key_float") 
    1107 !!   namhsb       Heat and salt budgets 
    1108 !!====================================================================== 
    1109 ! 
    1110 !----------------------------------------------------------------------- 
    1111 &namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
    1112 !----------------------------------------------------------------------- 
    1113    nn_nchunks_i=   4       !  number of chunks in i-dimension 
    1114    nn_nchunks_j=   4       !  number of chunks in j-dimension 
    1115    nn_nchunks_k=   31      !  number of chunks in k-dimension 
    1116                            !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
    1117                            !  is optimal for postprocessing which works exclusively with horizontal slabs 
    1118    ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
    1119                            !  (F) ignore chunking information and produce netcdf3-compatible files 
    1120 / 
    1121 !----------------------------------------------------------------------- 
    1122 &namtrd        !   diagnostics on dynamics and/or tracer trends 
    1123 !              !       and/or mixed-layer trends and/or barotropic vorticity 
     1098!!   nam_diaharm  Harmonic analysis of tidal constituents               ('key_diaharm') 
     1099!!   namdct       transports through some sections 
     1100!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
     1101!!====================================================================== 
     1102! 
     1103!----------------------------------------------------------------------- 
     1104&namtrd        !   diagnostics on dynamics and/or tracer trends         (default F) 
     1105!              !   and/or mixed-layer trends and/or barotropic vorticity 
    11241106!----------------------------------------------------------------------- 
    11251107   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
    11261108   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
    1127    ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
    1128    ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet) 
     1109   ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
     1110   ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet) 
    11291111   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
    11301112   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
    1131    ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output 
     1113   ln_tra_trd  = .false.   ! (T) 3D tracer trend output 
    11321114   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
    11331115   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
    … …  
    11401122!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    11411123!!gm 
     1124!----------------------------------------------------------------------- 
     1125&namptr       !   Poleward Transport Diagnostic                         (default F) 
     1126!----------------------------------------------------------------------- 
     1127   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
     1128   ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not 
     1129/ 
     1130!----------------------------------------------------------------------- 
     1131&namhsb       !  Heat and salt budgets                                  (default F) 
     1132!----------------------------------------------------------------------- 
     1133   ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F) 
     1134/ 
    11421135!----------------------------------------------------------------------- 
    11431136&namflo       !   float parameters                                      ("key_float") 
    … …  
    11551148/ 
    11561149!----------------------------------------------------------------------- 
    1157 &namptr       !   Poleward Transport Diagnostic 
    1158 !----------------------------------------------------------------------- 
    1159    ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
    1160    ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    1161 / 
    1162 !----------------------------------------------------------------------- 
    1163 &namhsb       !  Heat and salt budgets 
    1164 !----------------------------------------------------------------------- 
    1165    ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F) 
    1166 / 
    1167 !----------------------------------------------------------------------- 
    1168 &nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm') 
     1150&nam_diaharm   !   Harmonic analysis of tidal constituents              ('key_diaharm') 
    11691151!----------------------------------------------------------------------- 
    11701152    nit000_han = 1         ! First time step used for harmonic analysis 
    … …  
    11751157/ 
    11761158!----------------------------------------------------------------------- 
    1177 &namdct        ! transports through sections 
     1159&namdct        ! transports through some sections 
    11781160!----------------------------------------------------------------------- 
    11791161    nn_dct      = 15       !  time step frequency for transports computing 
    … …  
    11831165                           !  0 < n : debug section number n 
    11841166/ 
    1185  
    1186 !!====================================================================== 
    1187 !!            ***  Observation & Assimilation namelists *** 
    1188 !!====================================================================== 
    1189 !!   namobs       observation and model comparison                      ('key_diaobs') 
     1167!----------------------------------------------------------------------- 
     1168&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
     1169!----------------------------------------------------------------------- 
     1170   nn_nchunks_i=   4       !  number of chunks in i-dimension 
     1171   nn_nchunks_j=   4       !  number of chunks in j-dimension 
     1172   nn_nchunks_k=   31      !  number of chunks in k-dimension 
     1173                           !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 
     1174                           !  is optimal for postprocessing which works exclusively with horizontal slabs 
     1175   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression 
     1176                           !  (F) ignore chunking information and produce netcdf3-compatible files 
     1177/ 
     1178 
     1179!!====================================================================== 
     1180!!               ***  Observation & Assimilation  *** 
     1181!!====================================================================== 
     1182!!   namobs       observation and model comparison 
    11901183!!   nam_asminc   assimilation increments                               ('key_asminc') 
    11911184!!====================================================================== 
    11921185! 
    11931186!----------------------------------------------------------------------- 
    1194 &namobs       !  observation usage switch                               ('key_diaobs') 
    1195 !----------------------------------------------------------------------- 
    1196    ln_t3d     = .false.    ! Logical switch for T profile observations 
    1197    ln_s3d     = .false.    ! Logical switch for S profile observations 
    1198    ln_ena     = .false.    ! Logical switch for ENACT insitu data set 
    1199    ln_cor     = .false.    ! Logical switch for Coriolis insitu data set 
    1200    ln_profb   = .false.    ! Logical switch for feedback insitu data set 
    1201    ln_sla     = .false.    ! Logical switch for SLA observations 
    1202    ln_sladt   = .false.    ! Logical switch for AVISO SLA data 
    1203    ln_slafb   = .false.    ! Logical switch for feedback SLA data 
    1204    ln_ssh     = .false.    ! Logical switch for SSH observations 
    1205    ln_sst     = .false.    ! Logical switch for SST observations 
    1206    ln_reysst  = .false.    ! Logical switch for Reynolds observations 
    1207    ln_ghrsst  = .false.    ! Logical switch for GHRSST observations 
    1208    ln_sstfb   = .false.    ! Logical switch for feedback SST data 
    1209    ln_sss     = .false.    ! Logical switch for SSS observations 
    1210    ln_seaice  = .false.    ! Logical switch for Sea Ice observations 
    1211    ln_vel3d   = .false.    ! Logical switch for velocity observations 
    1212    ln_velavcur= .false     ! Logical switch for velocity daily av. cur. 
    1213    ln_velhrcur= .false     ! Logical switch for velocity high freq. cur. 
    1214    ln_velavadcp = .false.  ! Logical switch for velocity daily av. ADCP 
    1215    ln_velhradcp = .false.  ! Logical switch for velocity high freq. ADCP 
    1216    ln_velfb   = .false.    ! Logical switch for feedback velocity data 
    1217    ln_grid_global = .false. ! Global distribtion of observations 
    1218    ln_grid_search_lookup = .false. !  Logical switch for obs grid search w/lookup table 
    1219    grid_search_file = 'grid_search'  !  Grid search lookup file header 
     1187&namobs       !  observation usage switch 
     1188!----------------------------------------------------------------------- 
     1189   ln_diaobs  = .false.             ! Logical switch for the observation operator 
     1190   ln_t3d     = .false.             ! Logical switch for T profile observations 
     1191   ln_s3d     = .false.             ! Logical switch for S profile observations 
     1192   ln_sla     = .false.             ! Logical switch for SLA observations 
     1193   ln_sst     = .false.             ! Logical switch for SST observations 
     1194   ln_sic     = .false.             ! Logical switch for Sea Ice observations 
     1195   ln_vel3d   = .false.             ! Logical switch for velocity observations 
     1196   ln_altbias = .false.             ! Logical switch for altimeter bias correction 
     1197   ln_nea     = .false.             ! Logical switch for rejection of observations near land 
     1198   ln_grid_global = .true.          ! Logical switch for global distribution of observations 
     1199   ln_grid_search_lookup = .false.  ! Logical switch for obs grid search w/lookup table 
     1200   ln_ignmis  = .true.              ! Logical switch for ignoring missing files 
     1201   ln_s_at_t  = .false.             ! Logical switch for computing model S at T obs if not there 
     1202   ln_sstnight = .false.            ! Logical switch for calculating night-time average for SST obs 
    12201203! All of the *files* variables below are arrays. Use namelist_cfg to add more files 
    1221    enactfiles = 'enact.nc' !  ENACT input observation file names (specify full array in namelist_cfg) 
    1222    coriofiles = 'corio.nc' !  Coriolis input observation file name 
    1223    profbfiles = 'profiles_01.nc' ! Profile feedback input observation file name 
    1224    ln_profb_enatim = .false !        Enact feedback input time setting switch 
    1225    slafilesact = 'sla_act.nc' !  Active SLA input observation file names 
    1226    slafilespas = 'sla_pass.nc' ! Passive SLA input observation file names 
    1227    slafbfiles = 'sla_01.nc' ! slafbfiles: Feedback SLA input observation file names 
    1228    sstfiles = 'ghrsst.nc'   ! GHRSST input observation file names 
    1229    sstfbfiles = 'sst_01.nc' ! Feedback SST input observation file names 
    1230    seaicefiles = 'seaice_01.nc' ! Sea Ice input observation file names 
    1231    velavcurfiles = 'velavcurfile.nc'  ! Vel. cur. daily av. input file name 
    1232    velhrcurfiles = 'velhrcurfile.nc'  ! Vel. cur. high freq. input file name 
    1233    velavadcpfiles = 'velavadcpfile.nc' ! Vel. ADCP daily av. input file name 
    1234    velhradcpfiles = 'velhradcpfile.nc' ! Vel. ADCP high freq. input file name 
    1235    velfbfiles = 'velfbfile.nc' ! Vel. feedback input observation file name 
    1236    dobsini = 20000101.000000  !  Initial date in window YYYYMMDD.HHMMSS 
    1237    dobsend = 20010101.000000  !  Final date in window YYYYMMDD.HHMMSS 
    1238    n1dint = 0  !               Type of vertical interpolation method 
    1239    n2dint = 0  !               Type of horizontal interpolation method 
    1240    ln_nea = .false.   !        Rejection of observations near land switch 
    1241    nmsshc     = 0     !        MSSH correction scheme 
    1242    mdtcorr = 1.61     !        MDT  correction 
    1243    mdtcutoff = 65.0   !        MDT cutoff for computed correction 
    1244    ln_altbias = .false.    ! Logical switch for alt bias 
    1245    ln_ignmis  = .true.     ! Logical switch for ignoring missing files 
    1246    endailyavtypes = 820    ! ENACT daily average types - array (use namelist_cfg to set more values) 
    1247    ln_grid_global = .true. 
    1248    ln_grid_search_lookup = .false. 
     1204   cn_profbfiles = 'profiles_01.nc'    ! Profile feedback input observation file names 
     1205   cn_slafbfiles = 'sla_01.nc'         ! SLA feedback input observation file names 
     1206   cn_sstfbfiles = 'sst_01.nc'         ! SST feedback input observation file names 
     1207   cn_sicfbfiles = 'sic_01.nc'         ! SIC feedback input observation file names 
     1208   cn_velfbfiles = 'vel_01.nc'         ! Velocity feedback input observation file names 
     1209   cn_altbiasfile = 'altbias.nc'       ! Altimeter bias input file name 
     1210   cn_gridsearchfile = 'gridsearch.nc' ! Grid search file name 
     1211   rn_gridsearchres = 0.5              ! Grid search resolution 
     1212   rn_dobsini = 00010101.000000        ! Initial date in window YYYYMMDD.HHMMSS 
     1213   rn_dobsend = 00010102.000000        ! Final date in window YYYYMMDD.HHMMSS 
     1214   nn_1dint = 0                        ! Type of vertical interpolation method 
     1215   nn_2dint = 0                        ! Type of horizontal interpolation method 
     1216   nn_msshc = 0                        ! MSSH correction scheme 
     1217   rn_mdtcorr = 1.61                   ! MDT  correction 
     1218   rn_mdtcutoff = 65.0                 ! MDT cutoff for computed correction 
     1219   nn_profdavtypes = -1                ! Profile daily average types - array 
     1220   ln_sstbias = .false.  
     1221   cn_sstbias_files = 'sstbias.nc' 
    12491222/ 
    12501223!----------------------------------------------------------------------- 
    … …  
    12671240/ 
    12681241!----------------------------------------------------------------------- 
    1269 &namsbc_wave   ! External fields from wave model 
    1270 !----------------------------------------------------------------------- 
    1271 !              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    1272 !              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    1273    sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    1274    sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    1275    sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    1276    sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    1277 ! 
    1278    cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
    1279 / 
    1280 !----------------------------------------------------------------------- 
    1281 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    1282 !----------------------------------------------------------------------- 
    1283    ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 
    1284    ln_neptsimp       = .false.  ! yes/no use simplified neptune 
    1285  
    1286    ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep 
    1287    rn_tslse          =  1.2e4   ! value of lengthscale L at the equator 
    1288    rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole 
    1289    ! Specify whether to ramp down the Neptune velocity in shallow 
    1290    ! water, and if so the depth range controlling such ramping down 
    1291    ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water 
    1292    rn_htrmin         =  100.0   ! min. depth of transition range 
    1293    rn_htrmax         =  200.0   ! max. depth of transition range 
    1294 / 
     1242&namdiu !   Cool skin and warm layer models 
     1243!----------------------------------------------------------------------- 
     1244   ln_diurnal      = .false.   !  
     1245   ln_diurnal_only = .false.   ! 
     1246/ 
     1247!----------------------------------------------------------------------- 
     1248&nam_diatmb  !  Top Middle Bottom Output 
     1249!----------------------------------------------------------------------- 
     1250   ln_diatmb  = .false.    !  Choose Top Middle and Bottom output or not 
     1251/ 
     1252!----------------------------------------------------------------------- 
     1253&namwad  !   Wetting and drying 
     1254!----------------------------------------------------------------------- 
     1255   ln_wd             = .false.  ! T/F activation of wetting and drying 
     1256   rn_wdmin1         =  0.1     ! Minimum wet depth on dried cells 
     1257   rn_wdmin2         =  0.01    ! Tolerance of min wet depth on dried cells 
     1258   rn_wdld           =  20.0    ! Land elevation below which wetting/drying is allowed 
     1259   nn_wdit           =  10      ! Max iterations for W/D limiter 
     1260/ 
     1261!----------------------------------------------------------------------- 
     1262&nam_dia25h  !  25h Mean Output 
     1263!----------------------------------------------------------------------- 
     1264   ln_dia25h  = .false.    ! Choose 25h mean output or not 
     1265/ 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/SHARED/namelist_top_ref

    r6401 r6404  
    1111!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1212!----------------------------------------------------------------------- 
    13 &namtrc_run     !   run information 
     13&namtrc_run      !   run information 
    1414!----------------------------------------------------------------------- 
    1515   nn_dttrc      =  1        !  time step frequency for passive sn_tracers 
    1616   nn_writetrc   =  5475     !  time step frequency for sn_tracer outputs 
    17    ln_top_euler  = .false.    !  use Euler time-stepping for TOP 
     17   ln_top_euler  = .false.   !  use Euler time-stepping for TOP 
    1818   ln_rsttr      = .false.   !  start from a restart file (T) or not (F) 
    1919   nn_rsttr      =   0       !  restart control = 0 initial time step is not compared to the restart file value 
    20                            !                  = 1 do not use the value in the restart file 
    21                            !                  = 2 calendar parameters read in the restart file 
     20                             !                  = 1 do not use the value in the restart file 
     21                             !                  = 2 calendar parameters read in the restart file 
    2222   cn_trcrst_in  = "restart_trc"   !  suffix of pass. sn_tracer restart name (input) 
    2323   cn_trcrst_indir = "."           !  directory from which to read input passive tracer restarts 
    … …  
    2626/ 
    2727!----------------------------------------------------------------------- 
    28 &namtrc     !   tracers definition 
     28&namtrc          !   tracers definition 
    2929!----------------------------------------------------------------------- 
    3030   ln_trcdta     =   .true.  !  Initialisation from data input file (T) or not (F) 
    … …  
    3535&namtrc_dta      !    Initialisation from data input file 
    3636!----------------------------------------------------------------------- 
    37 ! 
    38    cn_dir        =  './'      !  root directory for the location of the data files 
     37   cn_dir        =  './'     !  root directory for the location of the data files 
    3938/ 
    4039!----------------------------------------------------------------------- 
    41 &namtrc_adv    !   advection scheme for passive tracer  
     40&namtrc_adv      !   advection scheme for passive tracer  
    4241!----------------------------------------------------------------------- 
    43    ln_trcadv_cen2    =  .false.  !  2nd order centered scheme    
    44    ln_trcadv_tvd     =  .true.  !  TVD scheme 
    45    ln_trcadv_muscl   =  .false.   !  MUSCL scheme 
    46    ln_trcadv_muscl2  =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    47    ln_trcadv_ubs     =  .false.  !  UBS scheme 
    48    ln_trcadv_qck     =  .false.  !  QUICKEST scheme 
    49    ln_trcadv_msc_ups =  .false.  !  use upstream scheme within muscl 
     42   ln_trcadv_cen =  .false.  !  2nd order centered scheme 
     43      nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     44      nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
     45   ln_trcadv_fct =  .false.  !  FCT scheme 
     46      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     47      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     48      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     49      !                             !        (number of sub-timestep = nn_fct_zts) 
     50   ln_trcadv_mus =  .false.  !  MUSCL scheme 
     51      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
     52   ln_trcadv_ubs =  .false.  !  UBS scheme 
     53      nn_ubs_v   =  2               !  =2  , vertical 2nd order FCT 
     54   ln_trcadv_qck =  .false.  !  QUICKEST scheme 
    5055/ 
    5156!----------------------------------------------------------------------- 
    52 &namtrc_ldf    !   lateral diffusion scheme for passive tracer  
     57&namtrc_ldf      !   lateral diffusion scheme for passive tracer  
    5358!----------------------------------------------------------------------- 
    54 !                               !  Type of the operator :  
    55    ln_trcldf_lap    =  .true.   !     laplacian operator        
    56    ln_trcldf_bilap  =  .false.  !     bilaplacian operator      
    57                                 !  Direction of action  : 
    58    ln_trcldf_level  =  .false.  !     iso-level                 
    59    ln_trcldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    60    ln_trcldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
    61 !                               !  Coefficient 
    62    rn_ahtrc_0       =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    63    rn_ahtrb_0       =     0.    !     background eddy diffusivity for ldf_iso [m2/s] 
     59!                            !  Type of the operator:   
     60   ln_trcldf_lap   =  .true.     !    laplacian operator 
     61   ln_trcldf_blp   =  .false.    !  bilaplacian operator 
     62   !                         !  Direction of action: 
     63   ln_trcldf_lev   =  .false.    !  iso-level 
     64   ln_trcldf_hor   =  .false.    !  horizontal (geopotential) 
     65   ln_trcldf_iso   =  .true.     !  iso-neutral (standard operator) 
     66   ln_trcldf_triad =  .false.    !  iso-neutral (triad    operator) 
     67   !                         !  Coefficient  
     68   rn_ahtrc_0      = 2000.       !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     69   rn_bhtrc_0      = 1.e+12      !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    6470/ 
    6571!----------------------------------------------------------------------- 
    66 &namtrc_zdf        !   vertical physics 
     72&namtrc_zdf      !   vertical physics 
    6773!----------------------------------------------------------------------- 
    68    ln_trczdf_exp   =  .false.  !  split explicit (T) or implicit (F) time stepping 
    69    nn_trczdf_exp   =   3       !  number of sub-timestep for ln_trczdfexp=T 
     74   ln_trczdf_exp =  .false.  !  split explicit (T) or implicit (F) time stepping 
     75   nn_trczdf_exp =   3       !  number of sub-timestep for ln_trczdfexp=T 
    7076/ 
    7177!----------------------------------------------------------------------- 
    72 &namtrc_rad        !  treatment of negative concentrations  
     78&namtrc_rad      !  treatment of negative concentrations  
    7379!----------------------------------------------------------------------- 
    74    ln_trcrad   =  .true.  !  artificially correct negative concentrations (T) or not (F) 
     80   ln_trcrad     =  .true.   !  artificially correct negative concentrations (T) or not (F) 
    7581/ 
    7682!----------------------------------------------------------------------- 
    77 &namtrc_dmp    !   passive tracer newtonian damping    
     83&namtrc_dmp      !   passive tracer newtonian damping    
    7884!----------------------------------------------------------------------- 
    79    nn_zdmp_tr  =    1      !  vertical   shape =0    damping throughout the water column 
    80                            !                   =1 no damping in the mixing layer (kz  criteria) 
    81                            !                   =2 no damping in the mixed  layer (rho crieria) 
    82    cn_resto_tr  = 'resto_tr.nc'    !  create a damping.coeff NetCDF file (=1) or not (=0) 
     85   nn_zdmp_tr    =    1      !  vertical   shape =0    damping throughout the water column 
     86                             !                   =1 no damping in the mixing layer (kz  criteria) 
     87                             !                   =2 no damping in the mixed  layer (rho crieria) 
     88   cn_resto_tr   = 'resto_tr.nc'    !  create a damping.coeff NetCDF file (=1) or not (=0) 
    8389/ 
    8490!----------------------------------------------------------------------- 
    85 &namtrc_ice       !    Representation of sea ice growth & melt effects 
     91&namtrc_ice      !    Representation of sea ice growth & melt effects 
    8692!----------------------------------------------------------------------- 
    87    nn_ice_tr   =  -1        !  tracer concentration in sea ice  
    88                            !    =-1 (no vvl: identical cc in ice and ocean / vvl: cc_ice = 0) 
    89                            !    = 0 (no vvl: cc_ice = zero / vvl: cc_ice = ) 
    90                            !    = 1 prescribed to a namelist value (implemented in pisces only) 
     93   nn_ice_tr     =  -1       !  tracer concentration in sea ice  
     94                             !    =-1 (no vvl: identical cc in ice and ocean / vvl: cc_ice = 0) 
     95                             !    = 0 (no vvl: cc_ice = zero / vvl: cc_ice = ) 
     96                             !    = 1 prescribed to a namelist value (implemented in pisces only) 
    9197/ 
    9298!----------------------------------------------------------------------- 
    93 &namtrc_trd       !   diagnostics on tracer trends        ('key_trdtrc') 
    94 !                          or mixed-layer trends          ('key_trdmld_trc') 
     99&namtrc_trd      !   diagnostics on tracer trends                       ('key_trdtrc') 
     100!                          or mixed-layer trends                        ('key_trdmld_trc') 
    95101!---------------------------------------------------------------------- 
    96    nn_trd_trc  =  5475      !  time step frequency and tracers trends 
    97    nn_ctls_trc =   0        !  control surface type in mixed-layer trends (0,1 or n<jpk) 
    98    rn_ucf_trc  =   1        !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
    99    ln_trdmld_trc_restart = .false.  !  restart for ML diagnostics 
     102   nn_trd_trc    =  5475     !  time step frequency and tracers trends 
     103   nn_ctls_trc   =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     104   rn_ucf_trc    =   1       !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     105   ln_trdmld_trc_restart = .false. !  restart for ML diagnostics 
    100106   ln_trdmld_trc_instant = .true.  !  flag to diagnose trends of instantantaneous or mean ML T/S 
    101    ln_trdtrc(1)  =  .true. 
    102    ln_trdtrc(2)  =  .true. 
    103    ln_trdtrc(23) =   .true. 
     107   ln_trdtrc( 1) = .true. 
     108   ln_trdtrc( 2) = .true. 
     109   ln_trdtrc(23) = .true. 
    104110/ 
    105111!----------------------------------------------------------------------- 
    106 &namtrc_dia       !   parameters for passive tracer additional diagnostics 
     112&namtrc_dia      !   parameters for passive tracer additional diagnostics 
    107113!---------------------------------------------------------------------- 
    108114   ln_diatrc     =  .true.   !  save additional diag. (T) or not (F) 
    109115   ln_diabio     =  .true.   !  output biological trends 
    110116   nn_writedia   =  5475     !  time step frequency for diagnostics 
    111    nn_writebio   =    10     !: frequency of biological outputs 
     117   nn_writebio   =    10     !  frequency of biological outputs 
    112118/ 
    113119!---------------------------------------------------------------------- 
    114 ! namtrc_bc       !   data for boundary conditions 
     120&namtrc_bc       !   data for boundary conditions 
    115121!----------------------------------------------------------------------- 
    116 &namtrc_bc 
    117 ! 
    118    cn_dir        =  './'      !  root directory for the location of the data files 
     122   cn_dir_sbc    =  './'     !  root directory for the location of SURFACE data files 
     123   cn_dir_cbc    =  './'     !  root directory for the location of COASTAL data files 
     124   cn_dir_obc    =  './'     !  root directory for the location of OPEN data files 
    119125/ 
     126!---------------------------------------------------------------------- 
     127&namtrc_bdy      !   Setup of tracer boundary conditions 
     128!----------------------------------------------------------------------- 
     129   cn_trc_dflt   = 'neumann' !  OBC applied by default to all tracers 
     130   cn_trc        = 'none'    !  Boundary conditions used for tracers with data files (selected in namtrc) 
     131 
     132   nn_trcdmp_bdy = 0         !  Use damping timescales defined in nambdy of namelist 
     133                             !  = 0 NO damping of tracers at open boudaries 
     134                             !  = 1 Only for tracers forced with external data 
     135                             !  = 2 Damping applied to all tracers 
     136/ 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/cfg.txt

    r6401 r6404  
    77AMM12 OPA_SRC 
    88ORCA2_LIM_PISCES OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC 
    9 GYRE OPA_SRC 
    109ORCA2_LIM3 OPA_SRC LIM_SRC_3 NST_SRC 
    1110ORCA2_LIM OPA_SRC LIM_SRC_2 NST_SRC 
    1211ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC 
     12GYRE OPA_SRC 

  • branches/2015/dev_r5803_UKMO_AGRIF_Vert_interp/NEMOGCM/CONFIG/uspcfg.txt

    r6401 r6404  
    11ORCA1_CICE # ORCA2_LIM # OPA_SRC TOP_SRC  # http://gws-access.ceda.ac.uk/public/nemo/uspconfigs/ORCA1_CICE/v3.6.0/ORCA1_CICE_ctl.txt 
    2 ISOMIP     # GYRE      # OPA_SRC          # http://gws-access.ceda.ac.uk/public/nemo/uspconfigs/ISOMIP/v3.6.0/ISOMIP_ctl.txt 
     2ISOMIP     # GYRE      # OPA_SRC          # http://gws-access.ceda.ac.uk/public/nemo/uspconfigs/ISOMIP/trunk/ISOMIP_ctl.txt 
     3IRISHSEA   # GYRE      # OPA_SRC          # http://gws-access.ceda.ac.uk/public/nemo/uspconfigs/IRISHSEA/v3.7.0/IRISHSEA_ctl.txt 
Note: See TracChangeset for help on using the changeset viewer.