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 7646 for trunk/NEMOGCM/CONFIG/SHARED/namelist_ref – NEMO

Ignore:
Timestamp:
2017-02-06T10:25:03+01:00 (7 years ago)
Author:
timgraham
Message:

Merge of dev_merge_2016 into trunk. UPDATE TO ARCHFILES NEEDED for XIOS2.
LIM_SRC_s/limrhg.F90 to follow in next commit due to change of kind (I'm unable to do it in this commit).
Merged using the following steps:

1) svn merge --reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk .
2) Resolve minor conflicts in sette.sh and namelist_cfg for ORCA2LIM3 (due to a change in trunk after branch was created)
3) svn commit
4) svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
5) svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2016/dev_merge_2016 .
6) At this stage I checked out a clean copy of the branch to compare against what is about to be committed to the trunk.
6) svn commit #Commit code to the trunk

In this commit I have also reverted a change to Fcheck_archfile.sh which was causing problems on the Paris machine.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/SHARED/namelist_ref

    r6497 r7646  
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!! NEMO/OPA  :  1 - run manager      (namrun) 
    5 !! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
    6 !!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
     5!! namelists    2 - Domain           (namcfg, namzgr, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
     6!!              3 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_sas) 
    77!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    88!!                                    namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave) 
     
    1111!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
    1212!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    13 !!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 
     13!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx, namzdf_tmx_new) 
    1414!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    1515!!             10 - miscellaneous    (nammpp, namctl) 
     
    6060!!   namcfg       parameters of the configuration 
    6161!!   namzgr       vertical coordinate                                   (default: NO selection) 
    62 !!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    6362!!   namdom       space and time domain (bathymetry, mesh, timestep) 
    6463!!   namwad       Wetting and drying                                    (default F) 
     
    7372&namcfg        !   parameters of the configuration 
    7473!----------------------------------------------------------------------- 
    75    cp_cfg      = "default" !  name of the configuration 
    76    cp_cfz      = "no zoom" !  name of the zoom of configuration 
    77    jp_cfg      =      0    !  resolution of the configuration 
    78    jpidta      =     10    !  1st lateral dimension ( >= jpi ) 
    79    jpjdta      =     12    !  2nd    "         "    ( >= jpj ) 
    80    jpkdta      =     31    !  number of levels      ( >= jpk ) 
    81    jpiglo      =     10    !  1st dimension of global domain --> i =jpidta 
    82    jpjglo      =     12    !  2nd    -                  -    --> j =jpjdta 
    83    jpizoom     =      1    !  left bottom (i,j) indices of the zoom 
    84    jpjzoom     =      1    !  in data domain indices 
    85    jperio      =      0    !  lateral cond. type (between 0 and 6) 
    86                                  !  = 0 closed                 ;   = 1 cyclic East-West 
    87                                  !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot 
    88                                  !  = 4 cyclic East-West AND North fold T-point pivot 
    89                                  !  = 5 North fold F-point pivot 
    90                                  !  = 6 cyclic East-West AND North fold F-point pivot 
     74   ln_read_cfg = .false.   !  (=T) read the domain configuration file 
     75      !                    !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules 
     76      cn_domcfg = "domain_cfg"         ! domain configuration filename 
     77      ! 
     78   ln_write_cfg= .false.   !  (=T) create the domain configuration file 
     79      cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename 
     80      ! 
    9181   ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present 
    92                            !  in netcdf input files, as the start j-row for reading 
    93 / 
    94 !----------------------------------------------------------------------- 
    95 &namzgr        !   vertical coordinate                                  (default: NO selection) 
    96 !----------------------------------------------------------------------- 
    97    ln_zco      = .false.   !  z-coordinate - full    steps 
    98    ln_zps      = .false.   !  z-coordinate - partial steps 
    99    ln_sco      = .false.   !  s- or hybrid z-s-coordinate 
    100    ln_isfcav   = .false.   !  ice shelf cavity 
    101    ln_linssh   = .false.   !  linear free surface 
    102 / 
    103 !----------------------------------------------------------------------- 
    104 &namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default F) 
    105 !----------------------------------------------------------------------- 
    106    ln_s_sh94   = .false.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
    107    ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
    108    ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
    109                            !  stretching coefficients for all functions 
    110    rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m) 
    111    rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
    112    rn_hc       =  150.0    !  critical depth for transition to stretched coordinates 
    113                         !!!!!!!  Envelop bathymetry 
    114    rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1) 
    115                         !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
    116    rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
    117    rn_bb       =    0.8    !  stretching with SH94 s-sigma 
    118                         !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
    119    rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
    120    rn_efold    =    0.0    !  efold length scale for transition to stretched coord 
    121    rn_zs       =    1.0    !  depth of surface grid box 
    122                            !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 
    123    rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb 
    124    rn_zb_b     =   -0.2    !  offset for calculating Zb 
    125                         !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
    126    rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1) 
     82   !                       !  in netcdf input files, as the start j-row for reading 
    12783/ 
    12884!----------------------------------------------------------------------- 
    12985&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    13086!----------------------------------------------------------------------- 
    131    nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file 
    132    rn_bathy    =    0.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
     87   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    13388   nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    134    nn_msh      =    1      !  create (=1) a mesh file or not (=0) 
    135    rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
     89   ! 
     90   nn_msh      =    0      !  create (>0) a mesh file or not (=0) 
    13691   rn_isfhmin  =    1.00   !  treshold (m) to discriminate grounding ice to floating ice 
    137    rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
    138    rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
    139                            ! 
     92   ! 
    14093   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
    14194   rn_atfp     =    0.1    !  asselin time filter parameter 
    142    ln_crs      = .false.      !  Logical switch for coarsening module 
    143    jphgr_msh   =       0               !  type of horizontal mesh 
    144                                        !  = 0 curvilinear coordinate on the sphere read in coordinate.nc 
    145                                        !  = 1 geographical mesh on the sphere with regular grid-spacing 
    146                                        !  = 2 f-plane with regular grid-spacing 
    147                                        !  = 3 beta-plane with regular grid-spacing 
    148                                        !  = 4 Mercator grid with T/U point at the equator 
    149    ppglam0     =       0.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
    150    ppgphi0     =     -35.0             ! latitude  of first raw and column T-point (jphgr_msh = 1) 
    151    ppe1_deg    =       1.0             !  zonal      grid-spacing (degrees) 
    152    ppe2_deg    =       0.5             !  meridional grid-spacing (degrees) 
    153    ppe1_m      =    5000.0             !  zonal      grid-spacing (degrees) 
    154    ppe2_m      =    5000.0             !  meridional grid-spacing (degrees) 
    155    ppsur       =    -4762.96143546300  !  ORCA r4, r2 and r05 coefficients 
    156    ppa0        =      255.58049070440  ! (default coefficients) 
    157    ppa1        =      245.58132232490  ! 
    158    ppkth       =       21.43336197938  ! 
    159    ppacr       =        3.0            ! 
    160    ppdzmin     =       10.             !  Minimum vertical spacing 
    161    pphmax      =     5000.             !  Maximum depth 
    162    ldbletanh   =    .TRUE.             !  Use/do not use double tanf function for vertical coordinates 
    163    ppa2        =      100.760928500000 !  Double tanh function parameters 
    164    ppkth2      =       48.029893720000 ! 
    165    ppacr2      =       13.000000000000 ! 
    166 / 
    167 !----------------------------------------------------------------------- 
    168 &namwad        !   Wetting and drying                                   (default F) 
    169 !----------------------------------------------------------------------- 
    170    ln_wd       = .false.   !  T/F activation of wetting and drying 
    171    rn_wdmin1   =  0.1      !  Minimum wet depth on dried cells 
    172    rn_wdmin2   =  0.01     !  Tolerance of min wet depth on dried cells 
    173    rn_wdld     =  20.0     !  Land elevation below which wetting/drying is allowed 
    174    nn_wdit     =  10       !  Max iterations for W/D limiter 
     95   ! 
     96   ln_crs      = .false.   !  Logical switch for coarsening module 
    17597/ 
    17698!----------------------------------------------------------------------- 
     
    185107   ln_tsd_init = .true.    !  Initialisation of ocean T & S with T & S input data (T) or not (F) 
    186108   ln_tsd_tradmp = .true.  !  damping of ocean T & S toward T & S input data (T) or not (F) 
     109/ 
     110!----------------------------------------------------------------------- 
     111&namwad        !   Wetting and drying                                   (default F) 
     112!----------------------------------------------------------------------- 
     113   ln_wd       = .false.   !  T/F activation of wetting and drying 
     114   rn_wdmin1   =  0.1      !  Minimum wet depth on dried cells 
     115   rn_wdmin2   =  0.01     !  Tolerance of min wet depth on dried cells 
     116   rn_wdld     =  20.0     !  Land elevation below which wetting/drying is allowed 
     117   nn_wdit     =  10       !  Max iterations for W/D limiter 
    187118/ 
    188119!----------------------------------------------------------------------- 
     
    230161!!====================================================================== 
    231162!!   namsbc          surface boundary condition 
    232 !!   namsbc_ana      analytical         formulation                     (ln_ana     =T) 
    233163!!   namsbc_flx      flux               formulation                     (ln_flx     =T) 
    234 !!   namsbc_clio     CLIO bulk formulae formulation                     (ln_blk_clio=T) 
    235 !!   namsbc_core     CORE bulk formulae formulation                     (ln_blk_core=T) 
    236 !!   namsbc_mfs      MFS  bulk formulae formulation                     (ln_blk_mfs =T) 
     164!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T) 
    237165!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" ) 
    238 !!   namsbc_sas      StAndalone Surface module 
     166!!   namsbc_sas      Stand-Alone Surface module 
    239167!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T) 
    240168!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T) 
     
    254182                           !     (also = the frequency of sea-ice & iceberg model call) 
    255183                     ! Type of air-sea fluxes  
    256    ln_ana      = .false.   !  analytical formulation                    (T => fill namsbc_ana ) 
     184   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc) 
    257185   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    258    ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio) 
    259    ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    260    ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
     186   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    261187                     ! Type of coupling (Ocean/Ice/Atmosphere) : 
    262188   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
     
    274200   nn_ice      = 2         !  =0 no ice boundary condition   , 
    275201                           !  =1 use observed ice-cover      , 
    276                            !  =2 ice-model used                         ("key_lim3", "key_lim2", "key_cice") 
     202                           !  =2 to 4 :  ice-model used (LIM2, LIM3 or CICE)                         ("key_lim3", "key_lim2", or "key_cice") 
    277203   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    278204                           !  =1 levitating ice with mass and salt exchange but no presure effect 
    279205                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    280206                     ! Misc. options of sbc :  
    281    ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr ) 
     207   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr) 
    282208   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    283209   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
     
    288214   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    289215   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf) 
    290    ln_wave     = .false.   !  coupling with surface wave                (T => fill namsbc_wave) 
     216   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
     217   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
     218   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)  
     219   ln_tauoc    = .false.   !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave) 
     220   ln_stcor    = .false.   !  Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) 
    291221   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
    292222                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
    293 / 
    294 !----------------------------------------------------------------------- 
    295 &namsbc_ana    !   analytical surface boundary condition 
    296 !----------------------------------------------------------------------- 
    297    nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps 
    298    rn_utau0    =   0.5     !  uniform value for the i-stress 
    299    rn_vtau0    =   0.e0    !  uniform value for the j-stress 
    300    rn_qns0     =   0.e0    !  uniform value for the total heat flux 
    301    rn_qsr0     =   0.e0    !  uniform value for the solar radiation 
    302    rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P) 
    303223/ 
    304224!----------------------------------------------------------------------- 
     
    316236/ 
    317237!----------------------------------------------------------------------- 
    318 &namsbc_clio   !   namsbc_clio  CLIO bulk formulae 
    319 !----------------------------------------------------------------------- 
    320 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    321 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    322    sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    323    sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    324    sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    325    sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    326    sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    327    sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    328    sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    329  
    330    cn_dir      = './'      !  root directory for the location of the bulk files are 
    331 / 
    332 !----------------------------------------------------------------------- 
    333 &namsbc_core   !   namsbc_core  CORE bulk formulae 
    334 !----------------------------------------------------------------------- 
    335 !              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
    336 !              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
     238&namsbc_blk   !   namsbc_blk  generic Bulk formula                      (ln_blk = T) 
     239!----------------------------------------------------------------------- 
     240!              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                              ! rotation ! land/sea mask ! 
     241!              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                             ! pairing  ! filename      ! 
    337242   sn_wndi     = 'u_10.15JUNE2009_fill'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , '' 
    338243   sn_wndj     = 'v_10.15JUNE2009_fill'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , '' 
     
    343248   sn_prec     = 'ncar_precip.15JUNE2009_fill',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    344249   sn_snow     = 'ncar_precip.15JUNE2009_fill',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     250   sn_slp      = 'slp.15JUNE2009_fill'        ,         6         , 'SLP'     ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    345251   sn_tdif     = 'taudif_core'                ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    346  
     252   !                    !  bulk algorithm : 
     253   ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008) 
     254   ln_COARE_3p0= .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003) 
     255   ln_COARE_3p5= .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013) 
     256   ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31) 
     257   ! 
    347258   cn_dir      = './'      !  root directory for the location of the bulk files 
    348259   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
     
    353264   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity 
    354265                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    355 / 
    356 !----------------------------------------------------------------------- 
    357 &namsbc_mfs   !   namsbc_mfs  MFS bulk formulae 
    358 !----------------------------------------------------------------------- 
    359 !              !  file name  ! frequency (hours) ! variable ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask ! 
    360 !              !             !  (if <0  months)  !   name   !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      ! 
    361    sn_wndi     =   'ecmwf'   ,        6          , 'u10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    362    sn_wndj     =   'ecmwf'   ,        6          , 'v10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    363    sn_clc      =   'ecmwf'   ,        6          , 'clc'    ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,   '' 
    364    sn_msl      =   'ecmwf'   ,        6          , 'msl'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    365    sn_tair     =   'ecmwf'   ,        6          , 't2'     ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    366    sn_rhm      =   'ecmwf'   ,        6          , 'rh'     ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,   '' 
    367    sn_prec     =   'ecmwf'   ,        6          , 'precip' ,    .true.    , .true. , 'daily'   ,'bicubic.nc' ,   ''     ,   '' 
    368  
    369    cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files 
     266   ln_Cd_L12   = .false.   !  Modify the drag ice-atm and oce-atm depending on ice concentration 
     267                           !  This parameterization is from Lupkes et al. (JGR 2012) 
    370268/ 
    371269!----------------------------------------------------------------------- 
     
    380278   sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    381279   sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     280   sn_snd_crtw   =       'none'                 ,    'no'    ,     ''      ,         ''           , 'U,V' 
     281   sn_snd_ifrac  =       'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     282   sn_snd_wlev   =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    382283! receive 
    383284   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     
    391292   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    392293   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     294   sn_rcv_hsig   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     295   sn_rcv_iceflx =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     296   sn_rcv_mslp   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     297   sn_rcv_phioc  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     298   sn_rcv_sdrfx  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     299   sn_rcv_sdrfy  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     300   sn_rcv_wper   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     301   sn_rcv_wnum   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     302   sn_rcv_wstrf  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     303   sn_rcv_wdrag  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    393304! 
    394305   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
     
    397308/ 
    398309!----------------------------------------------------------------------- 
    399 &namsbc_sas    !   analytical surface boundary condition 
     310&namsbc_sas    !   Stand Alone Surface boundary condition 
    400311!----------------------------------------------------------------------- 
    401312!              !  file name  ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    402313!              !             !  (if <0  months)  !   name    !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     314   l_sasread   = .TRUE.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
    403315   sn_usp      = 'sas_grid_U',     120           , 'vozocrtx',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
    404316   sn_vsp      = 'sas_grid_V',     120           , 'vomecrty',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    '' 
     
    495407/ 
    496408!----------------------------------------------------------------------- 
    497 &namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk 
     409&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T) 
    498410!----------------------------------------------------------------------- 
    499411!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     
    537449&namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
    538450!----------------------------------------------------------------------- 
    539 !              ! file name ! frequency (hours) ! variable    ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    540 !              !           !  (if <0  months)  !   name      !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    541    sn_cdg      = 'cdg_wave',        1          , 'drag_coeff',   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    542    sn_usd      = 'sdw_wave',        1          , 'u_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    543    sn_vsd      = 'sdw_wave',        1          , 'v_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    544    sn_wn       = 'sdw_wave',        1          , 'wave_num'  ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    545 ! 
    546    cn_dir_cdg  = './'      !  root directory for the location of drag coefficient files 
    547    ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model 
    548    ln_sdw      = .false.   !  Computation of 3D stokes drift                
     451!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     452!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     453   sn_cdg      =  'sdw_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     454   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     455   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     456   sn_hsw      =  'sdw_wave' ,        1          , 'hs'         ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     457   sn_wmp      =  'sdw_wave' ,        1          , 'wmp'        ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     458   sn_wnum     =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     459   sn_tauoc    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     460! 
     461   cn_dir  = './'  !  root directory for the location of drag coefficient files 
    549462/ 
    550463!----------------------------------------------------------------------- 
    551464&namberg       !   iceberg parameters                                   (default: No iceberg) 
    552465!----------------------------------------------------------------------- 
    553       ln_icebergs              = .false.              ! iceberg floats or not 
    554       ln_bergdia               = .true.               ! Calculate budgets 
    555       nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
    556       nn_verbose_write         = 15                   ! Timesteps between verbose messages 
    557       nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
    558                                                       ! Initial mass required for an iceberg of each class 
    559       rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
    560                                                       ! Proportion of calving mass to apportion to each class 
    561       rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    562                                                       ! Ratio between effective and real iceberg mass (non-dim) 
    563                                                       ! i.e. number of icebergs represented at a point 
    564       rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    565                                                       ! thickness of newly calved bergs (m) 
    566       rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
    567       rn_rho_bergs             = 850.                 ! Density of icebergs 
    568       rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
    569       ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
    570       rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
    571       rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
    572       ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
    573       nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
    574                                                       ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
    575       rn_test_box              = 108.0,  116.0, -66.0, -58.0 
    576       rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
    577  
    578 !            ! file name ! frequency (hours) !   variable   ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    579 !            !           !  (if <0  months)  !     name     !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
    580       sn_icb =  'calving',       -1          , 'calvingmask',   .true.    , .true.  , 'yearly'  ,    ''    ,    ''    ,    '' 
    581  
    582       cn_dir = './' 
     466   ln_icebergs              = .false.              ! iceberg floats or not 
     467   ln_bergdia               = .true.               ! Calculate budgets 
     468   nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
     469   nn_verbose_write         = 15                   ! Timesteps between verbose messages 
     470   nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     471                                                   ! Initial mass required for an iceberg of each class 
     472   rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
     473                                                   ! Proportion of calving mass to apportion to each class 
     474   rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
     475                                                   ! Ratio between effective and real iceberg mass (non-dim) 
     476                                                   ! i.e. number of icebergs represented at a point 
     477   rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
     478                                                   ! thickness of newly calved bergs (m) 
     479   rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
     480   rn_rho_bergs             = 850.                 ! Density of icebergs 
     481   rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
     482   ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
     483   rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
     484   rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
     485   ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
     486   nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
     487                                                   ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
     488   rn_test_box              = 108.0,  116.0, -66.0, -58.0 
     489   rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
     490 
     491!         ! file name ! frequency (hours) !   variable   ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     492!         !           !  (if <0  months)  !     name     !   (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     493   sn_icb =  'calving',       -1          , 'calvingmask',  .true.      , .true.  , 'yearly'  , ''       , ''       , '' 
     494 
     495   cn_dir = './' 
    583496/ 
    584497 
     
    589502!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
    590503!!   nam_tide      Tidal forcing  
    591 !!   nambdy        Unstructured open boundaries                         ("key_bdy") 
    592 !!   nambdy_dta    Unstructured open boundaries - external data         ("key_bdy") 
    593 !!   nambdy_tide   tidal forcing at open boundaries                     ("key_bdy_tides") 
     504!!   nambdy        Unstructured open boundaries                          
     505!!   nambdy_dta    Unstructured open boundaries - external data          
     506!!   nambdy_tide   tidal forcing at open boundaries                      
    594507!!====================================================================== 
    595508! 
     
    611524/ 
    612525!----------------------------------------------------------------------- 
    613 &nam_tide      !   tide parameters                                      ("key_tide") 
    614 !----------------------------------------------------------------------- 
     526&nam_tide      !   tide parameters 
     527!----------------------------------------------------------------------- 
     528   ln_tide     = .false. 
    615529   ln_tide_pot = .true.    !  use tidal potential forcing 
    616530   ln_tide_ramp= .false.   ! 
     
    619533/ 
    620534!----------------------------------------------------------------------- 
    621 &nambdy        !  unstructured open boundaries                          ("key_bdy") 
    622 !----------------------------------------------------------------------- 
     535&nambdy        !  unstructured open boundaries                           
     536!----------------------------------------------------------------------- 
     537    ln_bdy         = .false.              !  Use unstructured open boundaries 
    623538    nb_bdy         = 0                    !  number of open boundary sets 
    624539    ln_coords_file = .true.               !  =T : read bdy coordinates from file 
     
    651566    ln_vol        = .false.               !  total volume correction (see nn_volctl parameter) 
    652567    nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero 
    653 / 
    654 !----------------------------------------------------------------------- 
    655 &nambdy_dta    !  open boundaries - external data                       ("key_bdy") 
     568    nb_jpk_bdy    = -1                    ! number of levels in the bdy data (set < 0 if consistent with planned run) 
     569/ 
     570!----------------------------------------------------------------------- 
     571&nambdy_dta    !  open boundaries - external data                        
    656572!----------------------------------------------------------------------- 
    657573!              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     
    758674   !                     ! S-EOS coefficients (ln_seos=T): 
    759675   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    760    rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
    761    rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     676   rn_a0       =  1.6550e-1      !  thermal expension coefficient 
     677   rn_b0       =  7.6554e-1      !  saline  expension coefficient 
    762678   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
    763679   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
     
    900816!----------------------------------------------------------------------- 
    901817   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    902    ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
     818   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation) 
    903819   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
    904820   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf 
     
    940856   !                                !  = 30  F(i,j,k)=c2d*c1d 
    941857   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     858   !                                !  = 32  F(i,j,k)=F(local gridscale and deformation rate) 
     859   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    942860   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
    943861   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
    944862   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    945    ! 
    946    ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     863   !                       !  Smagorinsky settings (nn_ahm_ijk_t  = 32) : 
     864   rn_csmc       = 3.5         !  Smagorinsky constant of proportionality 
     865   rn_minfac     = 1.0         !  multiplier of theorectical lower limit 
     866   rn_maxfac     = 1.0         !  multiplier of theorectical upper limit 
    947867/ 
    948868 
     
    973893   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping 
    974894      nn_zdfexp   =    3        ! number of sub-timestep for ln_zdfexp=T 
     895   ln_zdfqiao  = .false.   !  Enhanced wave vertical mixing Qiao (2010) (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) 
    975896/ 
    976897!----------------------------------------------------------------------- 
     
    1027948   rn_hsro       =  0.02   !  Minimum surface roughness 
    1028949   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2) 
    1029    nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2) 
     950   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3) 
     951   !                             ! =3 requires ln_wave=T 
    1030952   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum) 
    1031953   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum) 
     
    1056978   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    1057979/ 
    1058  
    1059  
    1060980!!====================================================================== 
    1061981!!                  ***  Miscellaneous namelists  *** 
    1062982!!====================================================================== 
    1063983!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    1064 !!   namctl            Control prints & Benchmark 
     984!!   namctl            Control prints  
    1065985!!   namsto            Stochastic parametrization of EOS 
    1066986!!====================================================================== 
     
    1078998/ 
    1079999!----------------------------------------------------------------------- 
    1080 &namctl        !   Control prints & Benchmark 
     1000&namctl        !   Control prints  
    10811001!----------------------------------------------------------------------- 
    10821002   ln_ctl      = .false.   !  trends control print (expensive!) 
     
    10881008   nn_isplt    =    1      !  number of processors in i-direction 
    10891009   nn_jsplt    =    1      !  number of processors in j-direction 
    1090    nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench 
    1091                            !     (no physical validity of the results) 
    10921010   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    10931011   nn_diacfl   =    0      !  Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 
     
    11171035!!   namhsb       Heat and salt budgets                                 (default F) 
    11181036!!   namdiu       Cool skin and warm layer models                       (default F) 
     1037!!   namdiu       Cool skin and warm layer models                       (default F) 
    11191038!!   namflo       float parameters                                      ("key_float") 
    11201039!!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
     
    11591078&namdiu        !   Cool skin and warm layer models                       (default F) 
    11601079!----------------------------------------------------------------------- 
    1161    ln_diurnal      = .false.   !  
     1080   ln_diurnal      = .false.   ! 
    11621081   ln_diurnal_only = .false.   ! 
    11631082/ 
     
    11961115!----------------------------------------------------------------------- 
    11971116&nam_diatmb    !  Top Middle Bottom Output                               (default F) 
     1117!----------------------------------------------------------------------- 
     1118   ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
     1119/ 
     1120!----------------------------------------------------------------------- 
     1121&nam_dia25h    !  25h Mean Output                                        (default F) 
     1122!----------------------------------------------------------------------- 
     1123   ln_dia25h   = .false.   ! Choose 25h mean output or not 
     1124/ 
     1125!----------------------------------------------------------------------- 
     1126&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4") 
    11981127!----------------------------------------------------------------------- 
    11991128   ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not 
Note: See TracChangeset for help on using the changeset viewer.