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 12377 for NEMO/trunk/cfgs/SHARED/namelist_ref – NEMO

Ignore:
Timestamp:
2020-02-12T15:39:06+01:00 (4 years ago)
Author:
acc
Message:

The big one. Merging all 2019 developments from the option 1 branch back onto the trunk.

This changeset reproduces 2019/dev_r11943_MERGE_2019 on the trunk using a 2-URL merge
onto a working copy of the trunk. I.e.:

svn merge --ignore-ancestry \

svn+ssh://acc@forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/NEMO/trunk \
svn+ssh://acc@forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/NEMO/branches/2019/dev_r11943_MERGE_2019 ./

The --ignore-ancestry flag avoids problems that may otherwise arise from the fact that
the merge history been trunk and branch may have been applied in a different order but
care has been taken before this step to ensure that all applicable fixes and updates
are present in the merge branch.

The trunk state just before this step has been branched to releases/release-4.0-HEAD
and that branch has been immediately tagged as releases/release-4.0.2. Any fixes
or additions in response to tickets on 4.0, 4.0.1 or 4.0.2 should be done on
releases/release-4.0-HEAD. From now on future 'point' releases (e.g. 4.0.2) will
remain unchanged with periodic releases as needs demand. Note release-4.0-HEAD is a
transitional naming convention. Future full releases, say 4.2, will have a release-4.2
branch which fulfills this role and the first point release (e.g. 4.2.0) will be made
immediately following the release branch creation.

2020 developments can be started from any trunk revision later than this one.

Location:
NEMO/trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
  • NEMO/trunk/cfgs/SHARED/namelist_ref

    r12288 r12377  
    55!! namelists    2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 
    66!!                                    namsbc_sas, namtra_qsr, namsbc_rnf, 
    7 !!                                    namsbc_isf, namsbc_iscpl, namsbc_apr,  
     7!!                                    namisf, namsbc_apr,  
    88!!                                    namsbc_ssr, namsbc_wave, namberg) 
    99!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
     
    5151      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output) 
    5252      cn_ocerst_outdir = "."        !  directory in which to write output ocean restarts 
    53    ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model 
    5453   nn_istate   =       0   !  output the initial state (1) or not (0) 
    5554   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) 
     
    6564   ln_clobber  = .true.    !  clobber (overwrite) an existing file 
    6665   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    67    ln_xios_read = .FALSE.  !  use XIOS to read restart file (only for a single file restart) 
     66   ln_xios_read = .false.  !  use XIOS to read restart file (only for a single file restart) 
    6867   nn_wxios = 0      !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 
    6968/ 
     
    7271!----------------------------------------------------------------------- 
    7372   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    74    rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice 
    7573   ! 
    7674   rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     
    7977   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    8078   ! 
    81    ln_meshmask = .false.   !  =T create a mesh file 
     79   ln_meshmask = .true.   !  =T create a mesh file 
    8280/ 
    8381!----------------------------------------------------------------------- 
    8482&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg) 
    8583!----------------------------------------------------------------------- 
    86    ln_read_cfg = .false.   !  (=T) read the domain configuration file 
    87       !                    !  (=F) user defined configuration           (F => create/check namusr_def) 
     84   ln_read_cfg = .false.     !  (=T) read the domain configuration file 
     85      !                      !  (=F) user defined configuration           (F => create/check namusr_def) 
    8886      cn_domcfg = "domain_cfg"  ! domain configuration filename 
    8987      ! 
    90       ln_closea    = .false.    !  T => keep closed seas (defined by closea_mask field) in the   
    91       !                         !       domain and apply special treatment of freshwater fluxes. 
    92       !                         !  F => suppress closed seas (defined by closea_mask field)  
    93       !                         !       from the bathymetry at runtime. 
    94       !                         !  If closea_mask field doesn't exist in the domain_cfg file 
    95       !                         !       then this logical does nothing. 
    96    ln_write_cfg = .false.   !  (=T) create the domain configuration file 
     88      ln_closea    = .false. !  (=T => fill namclo)  
     89      !                      !  (=F) no control of net precip/evap over closed sea 
     90      ! 
     91   ln_write_cfg = .false.    !  (=T) create the domain configuration file 
    9792      cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename 
    9893      ! 
    99    ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present 
    100    !                       !  in netcdf input files, as the start j-row for reading 
     94   ln_use_jattr = .false.    !  use (T) the file attribute: open_ocean_jstart, if present 
     95   !                         !  in netcdf input files, as the start j-row for reading 
     96/ 
     97!----------------------------------------------------------------------- 
     98&namclo        !   parameters of the closed sea (cs) behavior                (default: OFF) 
     99!----------------------------------------------------------------------- 
     100   ln_maskcs = .false.        ! (=T) cs are masked ; So, in this case ln_mask_csundef and ln_clo_rnf have no effect. 
     101      !                       ! (=F => set ln_mask_csundef and ln_clo_rnf)  
     102      !                       ! cs masks are read and net evap/precip over closed sea spread out depending on domain_cfg.nc masks. 
     103      !                       ! See ln_mask_csundef and ln_clo_rnf for specific option related to this case 
     104      ! 
     105      ln_mask_csundef = .true.   ! (=T) undefined closed seas are masked ;  
     106      !                          ! (=F) undefined closed seas are kept and no specific treatment is done for these closed seas 
     107      ! 
     108      ln_clo_rnf = .true.        ! (=T) river mouth specified in domain_cfg.nc masks (rnf and emp case) are added to the runoff mask. 
     109      !                          !      allow the treatment of closed sea outflow grid-points to be the same as river mouth grid-points 
    101110/ 
    102111!----------------------------------------------------------------------- 
     
    106115   ln_tsd_init = .false.         !  ocean initialisation 
    107116   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp) 
    108     
     117 
    109118   cn_dir      = './'      !  root directory for the T-S data location 
    110119   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     
    184193!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
    185194!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T) 
    186 !!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 
    187 !!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T) 
    188195!!   namsbc_wave     external fields from wave model                    (ln_wave    =T) 
    189196!!   namberg         iceberg floats                                     (ln_icebergs=T) 
     
    195202   nn_fsbc     = 2         !  frequency of SBC module call 
    196203      !                    !  (control sea-ice & iceberg model call) 
    197                      ! Type of air-sea fluxes  
     204                     ! Type of air-sea fluxes 
    198205   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc) 
    199206   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    200207   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk ) 
     208   ln_abl      = .false.   !  ABL  formulation                          (T => fill namsbc_abl ) 
    201209      !              ! Type of coupling (Ocean/Ice/Atmosphere) : 
    202210   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
     
    205213      !                    !  =0 no opa-sas OASIS coupling: default single executable config. 
    206214      !                    !  =1 opa-sas OASIS coupling: multi executable config., OPA component 
    207       !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component  
     215      !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component 
    208216                     ! Sea-ice : 
    209    nn_ice      = 0         !  =0 no ice boundary condition     
     217   nn_ice      = 0         !  =0 no ice boundary condition 
    210218      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif ) 
    211219      !                    !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice") 
     
    213221   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges) 
    214222      !                    !  =F levitating ice (no pressure, mass and salt exchanges) 
    215                      ! Misc. options of sbc :  
     223                     ! Misc. options of sbc : 
    216224   ln_traqsr   = .false.   !  Light penetration in the ocean            (T => fill namtra_qsr) 
    217225   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
     
    222230   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    223231   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    224    ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf & namsbc_iscpl) 
    225232   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
    226233   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
    227    ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)  
     234   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 
    228235   nn_sdrift   =  0        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 
    229236      !                    !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] 
     
    250257/ 
    251258!----------------------------------------------------------------------- 
    252 &namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T) 
     259&namsbc_blk    !   namsbc_blk  generic Bulk formula          (ln_blk =T) 
    253260!----------------------------------------------------------------------- 
    254261   !                    !  bulk algorithm : 
    255    ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
     262   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008) 
    256263   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
    257    ln_COARE_3p5 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
    258    ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     264   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013) 
     265   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1) 
    259266      ! 
    260       rn_zqt      = 10.       !  Air temperature & humidity reference height (m) 
    261       rn_zu       = 10.       !  Wind vector reference height (m) 
    262       ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012) 
    263       ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015) 
    264       ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    265       rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    266       rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
    267       rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to 
    268       !                       !  calculate the wind stress (0.=absolute or 1.=relative winds) 
    269  
     267      rn_zqt     = 10.      !  Air temperature & humidity reference height (m) 
     268      rn_zu      = 10.      !  Wind vector reference height (m) 
     269      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012) 
     270      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015) 
     271      !                     !  - module of the mean stress" data 
     272      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow) 
     273      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.) 
     274      rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity  
     275      !                     !  used to calculate the wind stress 
     276      !                     ! (0. => absolute or 1. => relative winds) 
     277      ln_skin_cs = .false.  !  use the cool-skin parameterization 
     278      ln_skin_wl = .false.  !  use the warm-layer parameterization 
     279      !                     !   ==> only available in ECMWF and COARE algorithms 
     280      ln_humi_sph = .true.  !  humidity "sn_humi" is specific humidity  [kg/kg] 
     281      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K] 
     282      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%] 
     283   ! 
    270284   cn_dir      = './'      !  root directory for the bulk data location 
    271285   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 
     
    278292   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    279293   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     294   sn_hpgi     = 'NONE'                       ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , '' 
     295   sn_hpgj     = 'NONE'                       ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , '' 
    280296   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    281297   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    282298   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
    283    sn_tdif     = 'taudif_core'                ,   24         , 'taudif'  ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , '' 
     299/ 
     300!----------------------------------------------------------------------- 
     301&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T) 
     302!----------------------------------------------------------------------- 
     303   cn_dir           = './'      !  root directory for the location of the ABL grid file 
     304   cn_dom           = 'dom_cfg_abl.nc' 
     305 
     306   cn_ablrst_in     = "restart_abl"   !  suffix of abl restart name (input) 
     307   cn_ablrst_out    = "restart_abl"   !  suffix of abl restart name (output) 
     308   cn_ablrst_indir  = "."             !  directory to read   input abl restarts 
     309   cn_ablrst_outdir = "."             !  directory to write output abl restarts 
     310 
     311   ln_hpgls_frc   = .false. 
     312   ln_geos_winds  = .false. 
     313   nn_dyn_restore = 2         ! restoring option for dynamical ABL variables: = 0 no restoring 
     314                              !                                               = 1 equatorial restoring 
     315                              !                                               = 2 global restoring 
     316   rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour] 
     317   rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour] 
     318   rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour] 
     319   rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour] 
     320   nn_amxl       =  0         ! mixing length: = 0 Deardorff 80 length-scale 
     321                              !                = 1 length-scale based on the distance to the PBL height 
     322                              !                = 2 Bougeault & Lacarrere 89 length-scale 
     323   rn_Cm         = 0.0667     ! 0.126 in MesoNH 
     324   rn_Ct         = 0.1667     ! 0.143 in MesoNH 
     325   rn_Ce         = 0.4        ! 0.4   in MesoNH 
     326   rn_Ceps       = 0.7        ! 0.85  in MesoNH 
     327   rn_Rod        = 0.15       ! c0 in RMCA17 mixing length formulation (not yet implemented) 
     328   rn_Ric        = 0.139      !  Critical Richardson number (to compute PBL height and diffusivities) 
    284329/ 
    285330!----------------------------------------------------------------------- 
     
    375420   nn_chldta   =      0       !  RGB : Chl data (=1) or cst value (=0) 
    376421   rn_si1      =   23.0       !  2BD : longest depth of extinction 
    377     
     422 
    378423   cn_dir      = './'      !  root directory for the chlorophyl data location 
    379424   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     
    394439      nn_sssr_ice =   1       ! control of sea surface restoring under sea-ice 
    395440                              ! 0 = no restoration under ice : * (1-icefrac) 
    396                               ! 1 = restoration everywhere 
     441                              ! 1 = restoration everywhere  
    397442                              ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac 
     443 
    398444   cn_dir      = './'      !  root directory for the SST/SSS data location 
    399445   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
     
    442488/ 
    443489!----------------------------------------------------------------------- 
    444 &namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)  
    445 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    446    !                 ! type of top boundary layer  
    447    nn_isf      = 1         !  ice shelf melting/freezing 
    448                            !  1 = presence of ISF   ;  2 = bg03 parametrisation  
    449                            !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux 
    450                            !  options 1 and 4 need ln_isfcav = .true. (domzgr) 
    451       !              !  nn_isf = 1 or 2 cases: 
    452       rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula 
    453       rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula 
    454       !              !  nn_isf = 1 or 4 cases: 
    455       rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
    456       !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
    457       !              ! nn_isf = 1 case 
    458       nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006) 
    459       !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 
    460       nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s) 
    461       !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
    462       !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999) 
    463  
    464    !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
    465    !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    466    !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
    467 !* nn_isf = 4 case 
    468    sn_fwfisf   = 'rnfisf'    ,         -12.      ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    469 !* nn_isf = 3 case 
    470    sn_rnfisf   = 'rnfisf'    ,         -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    471 !* nn_isf = 2 and 3 cases  
    472    sn_depmax_isf ='rnfisf'   ,         -12.      ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    473    sn_depmin_isf ='rnfisf'   ,         -12.      ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    474 !* nn_isf = 2 case 
    475    sn_Leff_isf = 'rnfisf'    ,         -12.      ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
    476 / 
    477 !----------------------------------------------------------------------- 
    478 &namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T)  
    479 !-----------------------------------------------------------------------             or set or usr_def_zgr ) 
    480    nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
    481    ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl) 
    482    nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 
     490&namisf       !  Top boundary layer (ISF)                               (default: OFF) 
     491!----------------------------------------------------------------------- 
     492   ! 
     493   ! ---------------- ice shelf load ------------------------------- 
     494   ! 
     495   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 
     496      rn_isfload_T = -1.9 
     497      rn_isfload_S = 34.4 
     498   ! 
     499   ! ---------------- ice shelf melt formulation ------------------------------- 
     500   ! 
     501   ln_isf = .false.           ! activate ice shelf module 
     502      ln_isfdebug = .false.      ! add debug print in ISF code (global min/max/sum of specific variable) 
     503      cn_isfdir   = './'         ! directory for all ice shelf input file 
     504      ! 
     505      ! ---------------- cavities opened ------------------------------- 
     506      ! 
     507      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 
     508         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis) 
     509         !                       ! spe = fwfisf is read from a forcing field 
     510         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description) 
     511         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 
     512         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 
     513         !              !  cn_isfcav_mlt = 2eq or 3eq cases: 
     514         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99) 
     515         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0) 
     516         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 
     517         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 
     518         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method 
     519         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method 
     520         ! 
     521         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008) 
     522         !                       ! 0 => thickness of the tbl = thickness of the first wet cell 
     523         ! 
     524         !* 'spe' and 'oasis' case 
     525         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     526         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     527         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     528         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     529      ! 
     530      ! ---------------- cavities parametrised ------------------------------- 
     531      ! 
     532      ln_isfpar_mlt = .false.   ! ice shelf melting parametrised 
     533         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis) 
     534         !                      ! spe   = fwfisf is read from a forcing field 
     535         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation 
     536         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 
     537         ! 
     538         !* all cases 
     539         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 
     540         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     541         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     542         sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     543         sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    '' 
     544         !* 'spe' and 'oasis' case 
     545         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     546         !* 'bg03' case 
     547         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    '' 
     548      ! 
     549      ! ---------------- ice sheet coupling ------------------------------- 
     550      ! 
     551      ln_isfcpl = .false. 
     552         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells) 
     553         ln_isfcpl_cons = .false. 
    483554/ 
    484555!----------------------------------------------------------------------- 
     
    562633&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    563634!----------------------------------------------------------------------- 
     635   ln_agrif_2way = .true.  !  activate two way nesting 
    564636   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics 
    565637   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    566638   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
     639   rn_trelax_tra = 0.01    !  inverse of relaxation time (in steps) for tracers [] 
     640   rn_trelax_dyn = 0.01    !  inverse of relaxation time (in steps) for dynamics [] 
    567641   ln_chk_bathy  = .false. !  =T  check the parent bathymetry 
    568642/ 
     
    571645!----------------------------------------------------------------------- 
    572646   ln_tide     = .false.      ! Activate tides 
    573       ln_tide_pot   = .true.                !  use tidal potential forcing 
     647      nn_tide_var   = 1          !  Variant of tidal parameter set and tide-potential computation 
     648      !                          !     (1: default; 0: compatibility with previous versions) 
     649      ln_tide_dia   = .false.    !  Enable tidal diagnostic output 
     650      ln_tide_pot   = .false.               !  use tidal potential forcing 
     651         rn_tide_gamma = 0.7                   ! Tidal tilt factor 
    574652         ln_scal_load  = .false.               ! Use scalar approximation for 
    575653            rn_scal_load = 0.094               !     load potential 
    576654         ln_read_load  = .false.               ! Or read load potential from file 
    577655            cn_tide_load = 'tide_LOAD_grid_T.nc'  ! filename for load potential 
    578             !       
     656            ! 
    579657      ln_tide_ramp  = .false.               !  Use linear ramp for tides at startup 
    580          rdttideramp   =    0.                 !  ramp duration in days 
    581       clname(1)    = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg 
     658         rn_tide_ramp_dt = 0.               !  ramp duration in days 
     659      sn_tide_cnames(1) = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg 
    582660/ 
    583661!----------------------------------------------------------------------- 
     
    655733   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
    656734   ln_bdytide_2ddta = .false.                   ! 
    657    ln_bdytide_conj  = .false.                   !  
    658735/ 
    659736 
     
    682759!----------------------------------------------------------------------- 
    683760   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    684    rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     761   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0) 
    685762   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
    686763   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
     
    693770!----------------------------------------------------------------------- 
    694771   rn_Cd0      =  1.e-3    !  drag coefficient [-] 
    695    rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)  
     772   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0) 
    696773   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag) 
    697774   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases) 
     
    760837      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
    761838   ln_traadv_fct = .false. !  FCT scheme 
    762       nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order  
    763       nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order  
     839      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order 
     840      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order 
    764841   ln_traadv_mus = .false. !  MUSCL scheme 
    765842      ln_mus_ups = .false.       !  use upstream scheme near river mouths 
     
    782859   ln_traldf_triad = .false.   !  iso-neutral (triad    operator) 
    783860   ! 
    784    !                       !  iso-neutral options:         
     861   !                       !  iso-neutral options: 
    785862   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators) 
    786863   rn_slpmax       =  0.01     !  slope limit                           (both operators) 
     
    792869   nn_aht_ijk_t    = 0         !  space/time variation of eddy coefficient: 
    793870      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    794       !                             !   =  0           constant  
    795       !                             !   = 10 F(k)      =ldf_c1d  
    796       !                             !   = 20 F(i,j)    =ldf_c2d  
     871      !                             !   =  0           constant 
     872      !                             !   = 10 F(k)      =ldf_c1d 
     873      !                             !   = 20 F(i,j)    =ldf_c2d 
    797874      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    798875      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    799876      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    800       !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)  
     877      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case) 
    801878      !                             !                           or   = 1/12 Ud*Ld^3 (blp case) 
    802879      rn_Ud        = 0.01           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
     
    824901      nn_aei_ijk_t    = 0           !  space/time variation of eddy coefficient: 
    825902      !                             !   =-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  
     903      !                             !   =  0           constant 
     904      !                             !   = 10 F(k)      =ldf_c1d 
     905      !                             !   = 20 F(i,j)    =ldf_c2d 
    829906      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    830907      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
    831       !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le  
     908      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le 
    832909      rn_Ue        = 0.02           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 
    833910      rn_Le        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10) 
     
    869946   rn_lf_cutoff  =  5.0             !  cutoff frequency for low-pass filter  [days] 
    870947   rn_zdef_max   =  0.9             !  maximum fractional e3t deformation 
    871    ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
     948   ln_vvl_dbg    = .false.          !  debug prints    (T/F) 
    872949/ 
    873950!----------------------------------------------------------------------- 
     
    889966   ln_dynvor_eeT = .false. !  energy conserving scheme (een using e3t) 
    890967   ln_dynvor_een = .false. !  energy & enstrophy scheme 
    891       nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4  
     968      nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4 
    892969      !                       ! =1  e3f = mi(mj(e3t))/mi(mj( tmask)) 
    893970   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T)        ==>>> PLEASE DO NOT ACTIVATE 
     
    9341011      !                             !  =-30  read in eddy_viscosity_3D.nc file 
    9351012      !                             !  =-20  read in eddy_viscosity_2D.nc file 
    936       !                             !  =  0  constant  
     1013      !                             !  =  0  constant 
    9371014      !                             !  = 10  F(k)=c1d 
    9381015      !                             !  = 20  F(i,j)=F(grid spacing)=c2d 
     
    9401017      !                             !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
    9411018      !                             !  = 32  F(i,j,k)=F(local gridscale and deformation rate) 
    942       !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case)  
     1019      !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case) 
    9431020      !                             !                            or  = 1/12 Uv*Lv^3 (blp case) 
    9441021      rn_Uv      = 0.1              !  lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) 
     
    10641141                              !        = 0  constant 10 m length scale 
    10651142                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    1066       rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4    
     1143      rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 
    10671144/ 
    10681145!----------------------------------------------------------------------- 
     
    11211198!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    11221199!!   namflo       float parameters                                      (default: OFF) 
    1123 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    11241200!!   nam_diadct   transports through some sections                      (default: OFF) 
    11251201!!   nam_dia25h   25h Mean Output                                       (default: OFF) 
     
    11471223!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    11481224!!gm 
    1149 !----------------------------------------------------------------------- 
    1150 &namptr        !   Poleward Transport Diagnostic                        (default: OFF) 
    1151 !----------------------------------------------------------------------- 
    1152    ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F) 
    1153    ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not 
    11541225/ 
    11551226!----------------------------------------------------------------------- 
     
    11781249      ln_ariane   = .true.    !    Input with Ariane tool convention(T) 
    11791250      ln_flo_ascii= .true.    !    Output with Ariane tool netcdf convention(F) or ascii file (T) 
    1180 / 
    1181 !----------------------------------------------------------------------- 
    1182 &nam_diaharm   !   Harmonic analysis of tidal constituents              (default: OFF) 
    1183 !----------------------------------------------------------------------- 
    1184     ln_diaharm = .false.   ! Choose tidal harmonic output or not 
    1185        nit000_han = 1      !    First time step used for harmonic analysis 
    1186        nitend_han = 75     !    Last time step used for harmonic analysis 
    1187        nstep_han  = 15     !    Time step frequency for harmonic analysis 
    1188        tname(1)   = 'M2'   !    Name of tidal constituents 
    1189        tname(2)   = 'K1'   !              --- 
    11901251/ 
    11911252!----------------------------------------------------------------------- 
     
    13191380&namctl        !   Control prints                                       (default: OFF) 
    13201381!----------------------------------------------------------------------- 
    1321    ln_ctl = .FALSE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 
     1382   sn_cfctl%l_glochk = .FALSE.    ! Range sanity checks are local (F) or global (T). Set T for debugging only 
     1383   sn_cfctl%l_allon  = .FALSE.    ! IF T activate all options. If F deactivate all unless l_config is T 
    13221384     sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following 
    1323        sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 
     1385       sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 
    13241386       sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 
    13251387       sn_cfctl%l_oceout  = .FALSE. ! that  all areas report. 
    13261388       sn_cfctl%l_layout  = .FALSE. ! 
    1327        sn_cfctl%l_mppout  = .FALSE. ! 
    1328        sn_cfctl%l_mpptop  = .FALSE. ! 
     1389       sn_cfctl%l_prtctl  = .FALSE. ! 
     1390       sn_cfctl%l_prttrc  = .FALSE. ! 
     1391       sn_cfctl%l_oasout  = .FALSE. ! 
    13291392       sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0] 
    13301393       sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000] 
Note: See TracChangeset for help on using the changeset viewer.