source: branches/2017/dev_r8624_ENHANCE4_FREESURFACE/NEMOGCM/CONFIG/SHARED/namelist_ref @ 8888

Last change on this file since 8888 was 8888, checked in by jchanut, 3 years ago

Enhance4-freesurface. Minor corrections after testing, update namelists - #1959

  • Property svn:keywords set to Id
  • Property svn:mime-type set to text/x-fortran
File size: 87.2 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!!                            namelist_ref
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!! NEMO/OPA  :  1 - run manager      (namrun)
5!! namelists    2 - Domain           (namcfg, namzgr, namdom, namtsd, namcrs, namc1d, namc1d_uvd)
6!!              3 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_sas)
7!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,
8!!                                    namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave)
9!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide)
10!!              5 - bottom  boundary (nambfr, nambbc, nambbl)
11!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp)
12!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
13!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx, namzdf_tmx_new)
14!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto)
15!!             10 - miscellaneous    (nammpp, namctl)
16!!             11 - Obs & Assim      (namobs, nam_asminc)
17!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
18
19!!======================================================================
20!!                   ***  Run management namelists  ***
21!!======================================================================
22!!   namrun       parameters of the run
23!!======================================================================
24!
25!-----------------------------------------------------------------------
26&namrun        !   parameters of the run
27!-----------------------------------------------------------------------
28   nn_no       =       0   !  job number (no more used...)
29   cn_exp      =  "ORCA2"  !  experience name
30   nn_it000    =       1   !  first time step
31   nn_itend    =    5475   !  last  time step (std 5475)
32   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
34   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
35   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
36      nn_euler    =    1            !  = 0 : start with forward time step if ln_rstart=T
37      nn_rstctl   =    0            !  restart control ==> activated only if ln_rstart=T
38      !                             !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
39      !                             !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
40      !                             !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart
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
46   nn_istate   =       0   !  output the initial state (1) or not (0)
47   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
48   nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1)
49   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
50   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000)
51   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%)
52   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard
53   ln_clobber  = .true.    !  clobber (overwrite) an existing file
54   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
55/
56!
57!!======================================================================
58!!                      ***  Domain namelists  ***
59!!======================================================================
60!!   namcfg       parameters of the configuration
61!!   namdom       space and time domain (bathymetry, mesh, timestep)
62!!   namwad       Wetting and drying                                    (default F)
63!!   namtsd       data: temperature & salinity
64!!   namcrs       coarsened grid (for outputs and/or TOP)               ("key_crs")
65!!   namc1d       1D configuration options                              ("key_c1d")
66!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
67!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
68!!======================================================================
69!
70!-----------------------------------------------------------------------
71&namcfg        !   parameters of the configuration
72!-----------------------------------------------------------------------
73   ln_read_cfg = .false.   !  (=T) read the domain configuration file
74      !                    !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules
75      cn_domcfg = "domain_cfg"         ! domain configuration filename
76      !
77   ln_write_cfg= .false.   !  (=T) create the domain configuration file
78      cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename
79      !
80   ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present
81   !                       !  in netcdf input files, as the start j-row for reading
82/
83!-----------------------------------------------------------------------
84&namdom        !   space and time domain (bathymetry, mesh, timestep)
85!-----------------------------------------------------------------------
86   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time
87   nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA)
88   !
89   nn_msh      =    0      !  create (>0) a mesh file or not (=0)
90   rn_isfhmin  =    1.00   !  treshold (m) to discriminate grounding ice to floating ice
91   !
92   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0)
93   rn_atfp     =    0.1    !  asselin time filter parameter
94   !
95   ln_crs      = .false.   !  Logical switch for coarsening module
96/
97!-----------------------------------------------------------------------
98&namtsd        !   data : Temperature  & Salinity
99!-----------------------------------------------------------------------
100!              !  file name                 ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
101!              !                            !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
102   sn_tem = 'data_1m_potential_temperature_nomask',     -1      ,'votemper',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
103   sn_sal = 'data_1m_salinity_nomask'             ,     -1      ,'vosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
104   !
105   cn_dir      = './'      !  root directory for the location of the runoff files
106   ln_tsd_init = .true.    !  Initialisation of ocean T & S with T & S input data (T) or not (F)
107   ln_tsd_tradmp = .true.  !  damping of ocean T & S toward T & S input data (T) or not (F)
108/
109!-----------------------------------------------------------------------
110&namwad        !   Wetting and drying                                   (default F)
111!-----------------------------------------------------------------------
112   ln_wd       = .false.   !  T/F activation of wetting and drying
113   rn_wdmin1   =  0.1      !  Minimum wet depth on dried cells
114   rn_wdmin2   =  0.01     !  Tolerance of min wet depth on dried cells
115   rn_wdld     =  20.0     !  Land elevation below which wetting/drying is allowed
116   nn_wdit     =  10       !  Max iterations for W/D limiter
117/
118!-----------------------------------------------------------------------
119&namcrs        !   coarsened grid (for outputs and/or TOP)              ("key_crs")
120!-----------------------------------------------------------------------
121   nn_factx    = 3         !  Reduction factor of x-direction
122   nn_facty    = 3         !  Reduction factor of y-direction
123   nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT
124                           !  0, coarse grid is binned with preferential treatment of the north fold
125                           !  1, coarse grid is binned with centering at the equator
126                           !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.
127   nn_msh_crs  = 1         !  create (=1) a mesh file or not (=0)
128   nn_crs_kz   = 0         ! 0, MEAN of volume boxes
129                           ! 1, MAX of boxes
130                           ! 2, MIN of boxes
131   ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F )
132/
133!-----------------------------------------------------------------------
134&namc1d        !   1D configuration options                             ("key_c1d")
135!-----------------------------------------------------------------------
136   rn_lat1d    =      50   !  Column latitude (default at PAPA station)
137   rn_lon1d    =    -145   !  Column longitude (default at PAPA station)
138   ln_c1d_locpt=  .true.   ! Localization of 1D config in a grid (T) or independant point (F)
139/
140!-----------------------------------------------------------------------
141&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d")
142!-----------------------------------------------------------------------
143   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
144/
145!-----------------------------------------------------------------------
146&namc1d_uvd    !   data: U & V currents                                 ("key_c1d")
147!-----------------------------------------------------------------------
148!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
149!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
150   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , ''
151   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , ''
152!
153   cn_dir        = './'    !  root directory for the location of the files
154   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F)
155   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F)
156/
157
158!!======================================================================
159!!            ***  Surface Boundary Condition namelists  ***
160!!======================================================================
161!!   namsbc          surface boundary condition
162!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
163!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
164!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
165!!   namsbc_sas      Stand-Alone Surface module
166!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
167!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
168!!   namsbc_isf      ice shelf melting/freezing                         (nn_isf     >0)
169!!   namsbc_iscpl    coupling option between land ice model and ocean
170!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
171!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
172!!   namsbc_alb      albedo parameters
173!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
174!!   namberg         iceberg floats                                     (ln_icebergs=T)
175!!======================================================================
176!
177!-----------------------------------------------------------------------
178&namsbc        !   Surface Boundary Condition (surface module)
179!-----------------------------------------------------------------------
180   nn_fsbc     = 5         !  frequency of surface boundary condition computation
181                           !     (also = the frequency of sea-ice & iceberg model call)
182                     ! Type of air-sea fluxes
183   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc)
184   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx )
185   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
186                     ! Type of coupling (Ocean/Ice/Atmosphere) :
187   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 )
188   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 )
189   nn_components = 0       !  configuration of the opa-sas OASIS coupling
190                           !  =0 no opa-sas OASIS coupling: default single executable configuration
191                           !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component
192                           !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component
193   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used)
194                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled
195                           !  = 0  Average per-category fluxes (forced and coupled mode)
196                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled
197                           !  = 2  Redistribute a single flux over categories (coupled mode only)
198                     ! Sea-ice :
199   nn_ice      = 3         !  =0 no ice boundary condition   ,
200                           !  =1 use observed ice-cover      ,
201                           !  =2 to 4 :  ice-model used (LIM2, LIM3 or CICE)                         ("key_lim3", "key_lim2", or "key_cice")
202   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect)
203                           !  =1 levitating ice with mass and salt exchange but no presure effect
204                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure)
205                     ! Misc. options of sbc :
206   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr)
207   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave
208   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf)
209   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
210   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked
211                           !     =1 global mean of e-p-r set to zero at each time step
212                           !     =2 annual global mean of e-p-r set to zero
213   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr )
214   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf)
215   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave)
216   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)
217   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)
218   ln_tauoc    = .false.   !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)
219   ln_stcor    = .false.   !  Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
220   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
221                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
222/
223!-----------------------------------------------------------------------
224&namsbc_flx    !   surface boundary condition : flux formulation
225!-----------------------------------------------------------------------
226!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
227!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
228   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
229   sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
230   sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
231   sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
232   sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
233
234   cn_dir      = './'      !  root directory for the location of the flux files
235/
236!-----------------------------------------------------------------------
237&namsbc_blk   !   namsbc_blk  generic Bulk formula                      (ln_blk = T)
238!-----------------------------------------------------------------------
239!              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                              ! rotation ! land/sea mask !
240!              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                             ! pairing  ! filename      !
241   sn_wndi     = 'u_10.15JUNE2009_fill'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , ''
242   sn_wndj     = 'v_10.15JUNE2009_fill'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , ''
243   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
244   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
245   sn_tair     = 't_10.15JUNE2009_fill'       ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
246   sn_humi     = 'q_10.15JUNE2009_fill'       ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
247   sn_prec     = 'ncar_precip.15JUNE2009_fill',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
248   sn_snow     = 'ncar_precip.15JUNE2009_fill',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
249   sn_slp      = 'slp.15JUNE2009_fill'        ,         6         , 'SLP'     ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
250   sn_tdif     = 'taudif_core'                ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
251   !                    !  bulk algorithm :
252   ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008)
253   ln_COARE_3p0= .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003)
254   ln_COARE_3p5= .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013)
255   ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31)
256   !
257   cn_dir      = './'      !  root directory for the location of the bulk files
258   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data
259   rn_zqt      = 10.       !  Air temperature and humidity reference height (m)
260   rn_zu       = 10.       !  Wind vector reference height (m)
261   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow)
262   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.)
263   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity
264                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)
265   ln_Cd_L12   = .false.   !  Modify the drag ice-atm and oce-atm depending on ice concentration
266                           !  This parameterization is from Lupkes et al. (JGR 2012)
267/
268!-----------------------------------------------------------------------
269&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
270!-----------------------------------------------------------------------
271!                    !     description      !  multiple  !    vector   !      vector          ! vector !
272!                    !                      ! categories !  reference  !    orientation       ! grids  !
273! send
274   sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   ''
275   sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
276   sn_snd_thick  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
277   sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T'
278   sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
279   sn_snd_crtw   =       'none'                 ,    'no'    ,     ''      ,         ''           , 'U,V'
280   sn_snd_ifrac  =       'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
281   sn_snd_wlev   =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
282! receive
283   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
284   sn_rcv_taumod =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
285   sn_rcv_tau    =   'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V'
286   sn_rcv_dqnsdt =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
287   sn_rcv_qsr    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
288   sn_rcv_qns    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
289   sn_rcv_emp    =   'conservative'         ,    'no'    ,     ''      ,         ''          ,   ''
290   sn_rcv_rnf    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
291   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
292   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
293   sn_rcv_hsig   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
294   sn_rcv_iceflx =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
295   sn_rcv_mslp   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
296   sn_rcv_phioc  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
297   sn_rcv_sdrfx  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
298   sn_rcv_sdrfy  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
299   sn_rcv_wper   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
300   sn_rcv_wnum   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
301   sn_rcv_wstrf  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
302   sn_rcv_wdrag  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
303!
304   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentialy sending/receiving data
305   ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models
306   !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
307/
308!-----------------------------------------------------------------------
309&namsbc_sas    !   Stand Alone Surface boundary condition
310!-----------------------------------------------------------------------
311!              !  file name  ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
312!              !             !  (if <0  months)  !   name    !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
313   l_sasread   = .TRUE.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE.
314   sn_usp      = 'sas_grid_U',     120           , 'vozocrtx',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
315   sn_vsp      = 'sas_grid_V',     120           , 'vomecrty',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
316   sn_tem      = 'sas_grid_T',     120           , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
317   sn_sal      = 'sas_grid_T',     120           , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
318   sn_ssh      = 'sas_grid_T',     120           , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
319   sn_e3t      = 'sas_grid_T',     120           , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
320   sn_frq      = 'sas_grid_T',     120           , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
321
322   ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field
323   ln_read_frq = .false.   !  specify whether we must read frq or not
324   cn_dir      = './'      !  root directory for the location of the bulk files are
325/
326!-----------------------------------------------------------------------
327&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr=T)
328!-----------------------------------------------------------------------
329!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
330!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
331   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , ''
332
333   cn_dir      = './'      !  root directory for the location of the runoff files
334   ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration
335   ln_qsr_2bd  = .false.   !  2 bands              light penetration
336   ln_qsr_bio  = .false.   !  bio-model light penetration
337   nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0)
338   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
339   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
340   rn_si1      =   23.0    !  2 bands: longest depth of extinction
341   ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3
342/
343!-----------------------------------------------------------------------
344&namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf=T)
345!-----------------------------------------------------------------------
346!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
347!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
348   sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
349   sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
350   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
351   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
352   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , ''
353
354   cn_dir      = './'      !  root directory for the location of the runoff files
355   ln_rnf_mouth= .true.    !  specific treatment at rivers mouths
356      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T)
357      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
358   rn_rfact    =   1.e0    !  multiplicative factor for runoff
359   ln_rnf_depth= .false.   !  read in depth information for runoff
360   ln_rnf_tem  = .false.   !  read in temperature information for runoff
361   ln_rnf_sal  = .false.   !  read in salinity information for runoff
362   ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file
363      rn_rnf_max  = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
364      rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
365      nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0)
366/
367!-----------------------------------------------------------------------
368&namsbc_isf    !  Top boundary layer (ISF)                              (nn_isf >0)
369!-----------------------------------------------------------------------
370!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
371!              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
372! nn_isf == 4
373   sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
374! nn_isf == 3
375   sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
376! nn_isf == 2 and 3
377   sn_depmax_isf='rnfisf'  ,         -12       ,'sozisfmax',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
378   sn_depmin_isf='rnfisf'  ,         -12       ,'sozisfmin',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
379! nn_isf == 2
380   sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
381!
382! for all case
383   nn_isf      = 1         !  ice shelf melting/freezing
384                           !  1 = presence of ISF    2 = bg03 parametrisation
385                           !  3 = rnf file for isf   4 = ISF fwf specified
386                           !  option 1 and 4 need ln_isfcav = .true. (domzgr)
387! only for nn_isf = 1 or 2
388   rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula
389   rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula
390! only for nn_isf = 1 or 4
391   rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008)
392   !                       ! 0 => thickness of the tbl = thickness of the first wet cell
393! only for nn_isf = 1
394   nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006)
395   !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015)
396   nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s)
397   !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010)
398   !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999)
399/
400!-----------------------------------------------------------------------
401&namsbc_iscpl  !   land ice / ocean coupling option                     
402!-----------------------------------------------------------------------
403   nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells)
404   ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl)
405   nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency)
406/
407!-----------------------------------------------------------------------
408&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T)
409!-----------------------------------------------------------------------
410!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
411!              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
412   sn_apr      = 'patm'    ,         -1        ,'somslpre',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,      ''
413
414   cn_dir      = './'      !  root directory for the location of the bulk files
415   rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/
416   ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F)
417   ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data
418/
419!-----------------------------------------------------------------------
420&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr=T)
421!-----------------------------------------------------------------------
422!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
423!              !           !  (if <0  months)  !   name   !   (logical) !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
424   sn_sst      = 'sst_data',        24         ,  'sst'   ,    .false.  , .false., 'yearly'  ,    ''    ,    ''    ,     ''
425   sn_sss      = 'sss_data',        -1         ,  'sss'   ,    .true.   , .true. , 'yearly'  ,    ''    ,    ''    ,     ''
426
427   cn_dir      = './'      !  root directory for the location of the runoff files
428   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
429   nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2)
430                           !  or to SSS only (=1) or no damping term (=0)
431   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K]
432   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
433   ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2)
434   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
435/
436!-----------------------------------------------------------------------
437&namsbc_alb    !   albedo parameters
438!-----------------------------------------------------------------------
439   nn_ice_alb   =    1   !  parameterization of ice/snow albedo
440                         !     0: Shine & Henderson-Sellers (JGR 1985), giving clear-sky albedo
441                         !     1: "home made" based on Brandt et al. (JClim 2005) and Grenfell & Perovich (JGR 2004),
442                         !        giving cloud-sky albedo
443   rn_alb_sdry  =  0.85  !  dry snow albedo         : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky)
444   rn_alb_smlt  =  0.75  !  melting snow albedo     : 0.65 ( '' )          ; 0.75 ( '' )          ; obs 0.72-0.82 ( '' )
445   rn_alb_idry  =  0.60  !  dry ice albedo          : 0.72 ( '' )          ; 0.60 ( '' )          ; obs 0.54-0.65 ( '' )
446   rn_alb_imlt  =  0.50  !  bare puddled ice albedo : 0.53 ( '' )          ; 0.50 ( '' )          ; obs 0.49-0.58 ( '' )
447/
448!-----------------------------------------------------------------------
449&namsbc_wave   ! External fields from wave model                        (ln_wave=T)
450!-----------------------------------------------------------------------
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
462/
463!-----------------------------------------------------------------------
464&namberg       !   iceberg parameters                                   (default: No iceberg)
465!-----------------------------------------------------------------------
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 = './'
496/
497
498!!======================================================================
499!!               ***  Lateral boundary condition  ***
500!!======================================================================
501!!   namlbc        lateral momentum boundary condition
502!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")
503!!   nam_tide      Tidal forcing
504!!   nambdy        Unstructured open boundaries                         
505!!   nambdy_dta    Unstructured open boundaries - external data         
506!!   nambdy_tide   tidal forcing at open boundaries                     
507!!======================================================================
508!
509!-----------------------------------------------------------------------
510&namlbc        !   lateral momentum boundary condition
511!-----------------------------------------------------------------------
512   !                       !  free slip  !   partial slip  !   no slip   ! strong slip
513   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
514   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs.
515/
516!-----------------------------------------------------------------------
517&namagrif      !  AGRIF zoom                                            ("key_agrif")
518!-----------------------------------------------------------------------
519   nn_cln_update =    3    !  baroclinic update frequency
520   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
521   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
522   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
523   ln_chk_bathy  = .FALSE. !
524/
525!-----------------------------------------------------------------------
526&nam_tide      !   tide parameters
527!-----------------------------------------------------------------------
528   ln_tide     = .false.
529   ln_tide_pot = .true.    !  use tidal potential forcing
530   ln_tide_ramp= .false.   !
531   rdttideramp =    0.     !
532   clname(1)   = 'DUMMY'   !  name of constituent - all tidal components must be set in namelist_cfg
533/
534!-----------------------------------------------------------------------
535&nambdy        !  unstructured open boundaries                         
536!-----------------------------------------------------------------------
537    ln_bdy         = .false.              !  Use unstructured open boundaries
538    nb_bdy         = 0                    !  number of open boundary sets
539    ln_coords_file = .true.               !  =T : read bdy coordinates from file
540    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files
541    ln_mask_file   = .false.              !  =T : read mask from file
542    cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.)
543    cn_dyn2d       = 'none'               !
544    nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state
545                                          !  = 1, bdy data are read in 'bdydata   .nc' files
546                                          !  = 2, use tidal harmonic forcing data from files
547                                          !  = 3, use external data AND tidal harmonic forcing
548    cn_dyn3d      =  'none'               !
549    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state
550                                          !  = 1, bdy data are read in 'bdydata   .nc' files
551    cn_tra        =  'none'               !
552    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state
553                                          !  = 1, bdy data are read in 'bdydata   .nc' files
554    cn_ice_lim      =  'none'             !
555    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state
556                                          !  = 1, bdy data are read in 'bdydata   .nc' files
557    rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice
558    rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           --
559    rn_ice_age      = 30.                 !  lim3 only:      --   age                --
560
561    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers
562    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities
563    rn_time_dmp   =  1.                   ! Damping time scale in days
564    rn_time_dmp_out =  1.                 ! Outflow damping time scale
565    nn_rimwidth   = 10                    !  width of the relaxation zone
566    ln_vol        = .false.               !  total volume correction (see nn_volctl parameter)
567    nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero
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                       
572!-----------------------------------------------------------------------
573!              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
574!              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
575   bn_ssh      = 'amm12_bdyT_u2d',         24        , 'sossheig',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
576   bn_u2d      = 'amm12_bdyU_u2d',         24        , 'vobtcrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
577   bn_v2d      = 'amm12_bdyV_u2d',         24        , 'vobtcrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
578   bn_u3d      = 'amm12_bdyU_u3d',         24        , 'vozocrtx',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
579   bn_v3d      = 'amm12_bdyV_u3d',         24        , 'vomecrty',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
580   bn_tem      = 'amm12_bdyT_tra',         24        , 'votemper',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
581   bn_sal      = 'amm12_bdyT_tra',         24        , 'vosaline',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
582! for lim2
583!   bn_frld    = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
584!   bn_hicif   = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
585!   bn_hsnif   = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
586! for lim3
587!   bn_a_i     = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
588!   bn_ht_i    = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
589!   bn_ht_s    = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
590
591   cn_dir      = 'bdydta/' !  root directory for the location of the bulk files
592   ln_full_vel = .false.   
593/
594!-----------------------------------------------------------------------
595&nambdy_tide   !  tidal forcing at open boundaries
596!-----------------------------------------------------------------------
597   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files
598   ln_bdytide_2ddta = .false.                   !
599   ln_bdytide_conj  = .false.                   !
600/
601
602!!======================================================================
603!!                 ***  Bottom boundary condition  ***
604!!======================================================================
605!!   nambfr        bottom friction
606!!   nambbc        bottom temperature boundary condition
607!!   nambbl        bottom boundary layer scheme                         ("key_trabbl")
608!!======================================================================
609!
610!-----------------------------------------------------------------------
611&nambfr        !   bottom friction                                      (default: linear)
612!-----------------------------------------------------------------------
613   nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
614                           !                              = 2 : nonlinear friction
615   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case)
616   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
617   rn_bfri2_max=    1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T)
618   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2)
619   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T
620   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file )
621   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T)
622   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case)
623   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
624   rn_tfri2_max=    1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T)
625   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2)
626   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T
627   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file )
628   rn_tfrien   =   50.     !  local multiplying factor of tfr (ln_tfr2d=T)
629
630   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true)
631   ln_loglayer = .false.   !  logarithmic formulation (non linear case)
632/
633!-----------------------------------------------------------------------
634&nambbc        !   bottom temperature boundary condition                (default: NO)
635!-----------------------------------------------------------------------
636!              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
637!              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
638   sn_qgh      ='geothermal_heating.nc',  -12.       , 'heatflow',   .false.   , .true.  , 'yearly'  ,   ''     ,   ''     ,   ''
639   !
640   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
641   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux
642                           !     = 1 constant flux
643                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
644   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
645   cn_dir      = './'      !  root directory for the location of the runoff files
646/
647!-----------------------------------------------------------------------
648&nambbl        !   bottom boundary layer scheme                         ("key_trabbl")
649!-----------------------------------------------------------------------
650   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0)
651   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0)
652   rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s]
653   rn_gambbl   =  10.      !  advective bbl coefficient                 [s]
654/
655
656!!======================================================================
657!!                        Tracer (T & S ) namelists
658!!======================================================================
659!!   nameos           equation of state
660!!   namtra_adv       advection scheme
661!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.)
662!!   namtra_ldf       lateral diffusion scheme
663!!   namtra_ldfeiv    eddy induced velocity param.
664!!   namtra_dmp       T & S newtonian damping
665!!======================================================================
666!
667!-----------------------------------------------------------------------
668&nameos        !   ocean physical parameters
669!-----------------------------------------------------------------------
670   ln_teos10   = .false.         !  = Use TEOS-10 equation of state
671   ln_eos80    = .false.         !  = Use EOS80 equation of state
672   ln_seos     = .false.         !  = Use simplified equation of state (S-EOS)
673                                 !
674   !                     ! S-EOS coefficients (ln_seos=T):
675   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
676   rn_a0       =  1.6550e-1      !  thermal expension coefficient
677   rn_b0       =  7.6554e-1      !  saline  expension coefficient
678   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
679   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
680   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
681   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
682   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
683/
684!-----------------------------------------------------------------------
685&namtra_adv    !   advection scheme for tracer                          (default: NO advection)
686!-----------------------------------------------------------------------
687   ln_traadv_cen = .false. !  2nd order centered scheme
688      nn_cen_h   =  4            !  =2/4, horizontal 2nd order CEN / 4th order CEN
689      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT
690   ln_traadv_fct = .false. !  FCT scheme
691      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order
692      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order
693      nn_fct_zts =  0            !  >=1,  2nd order FCT scheme with vertical sub-timestepping
694      !                          !        (number of sub-timestep = nn_fct_zts)
695   ln_traadv_mus = .false. !  MUSCL scheme
696      ln_mus_ups = .false.       !  use upstream scheme near river mouths
697   ln_traadv_ubs = .false. !  UBS scheme
698      nn_ubs_v   =  2            !  =2  , vertical 2nd order FCT / COMPACT 4th order
699   ln_traadv_qck = .false. !  QUICKEST scheme
700/
701!-----------------------------------------------------------------------
702&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) (default: NO)
703!-----------------------------------------------------------------------
704   ln_mle      = .false.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation
705   rn_ce       = 0.06      ! magnitude of the MLE (typical value: 0.06 to 0.08)
706   nn_mle      = 1         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
707   rn_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_mle=0)
708   rn_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0)
709   rn_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
710   nn_mld_uv   = 0         ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
711   nn_conv     = 0         ! =1 no MLE in case of convection ; =0 always MLE
712   rn_rho_c_mle= 0.01      ! delta rho criterion used to calculate MLD for FK
713/
714!-----------------------------------------------------------------------
715&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO diffusion)
716!-----------------------------------------------------------------------
717   !                       !  Operator type:
718   !                           !  no diffusion: set ln_traldf_lap=..._blp=F
719   ln_traldf_lap   =  .false.  !    laplacian operator
720   ln_traldf_blp   =  .false.  !  bilaplacian operator
721   !
722   !                       !  Direction of action:
723   ln_traldf_lev   =  .false.  !  iso-level
724   ln_traldf_hor   =  .false.  !  horizontal (geopotential)
725   ln_traldf_iso   =  .false.  !  iso-neutral (standard operator)
726   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator)
727   !
728   !                       !  iso-neutral options:       
729   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators)
730   rn_slpmax       =   0.01    !  slope limit                      (both operators)
731   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only)
732   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only)
733   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only)
734   !
735   !                       !  Coefficients:
736   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef
737   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file
738   !                                !   =  0           constant
739   !                                !   = 10 F(k)      =ldf_c1d
740   !                                !   = 20 F(i,j)    =ldf_c2d
741   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
742   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
743   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
744   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s]
745   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s]
746/
747!-----------------------------------------------------------------------
748&namtra_ldfeiv !   eddy induced velocity param.                         (default: NO)
749!-----------------------------------------------------------------------
750   ln_ldfeiv     =.false.  ! use eddy induced velocity parameterization
751   ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities
752   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s]
753   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient
754   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
755   !                                !   =  0           constant
756   !                                !   = 10 F(k)      =ldf_c1d
757   !                                !   = 20 F(i,j)    =ldf_c2d
758   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
759   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d
760/
761!-----------------------------------------------------------------------
762&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO)
763!-----------------------------------------------------------------------
764   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F)
765   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column
766                           !                   =1 no damping in the mixing layer (kz  criteria)
767                           !                   =2 no damping in the mixed  layer (rho crieria)
768   cn_resto    ='resto.nc' !  Name of file containing restoration coeff. field (use dmp_tools to create this)
769/
770
771!!======================================================================
772!!                      ***  Dynamics namelists  ***
773!!======================================================================
774!!   namdyn_adv    formulation of the momentum advection
775!!   namdyn_vor    advection scheme
776!!   namdyn_hpg    hydrostatic pressure gradient
777!!   namdyn_spg    surface pressure gradient
778!!   namdyn_ldf    lateral diffusion scheme
779!!======================================================================
780!
781!-----------------------------------------------------------------------
782&namdyn_adv    !   formulation of the momentum advection                (default: vector form)
783!-----------------------------------------------------------------------
784   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)
785   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction
786   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
787   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
788   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection
789/
790!-----------------------------------------------------------------------
791&nam_vvl    !   vertical coordinate options                             (default: zstar)
792!-----------------------------------------------------------------------
793   ln_vvl_zstar  = .true.           !  zstar vertical coordinate
794   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations
795   ln_vvl_layer  = .false.          !  full layer vertical coordinate
796   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
797   ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator
798   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient
799   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days]
800   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days]
801   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation
802   ln_vvl_dbg    = .true.           !  debug prints    (T/F)
803/
804!-----------------------------------------------------------------------
805&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO)
806!-----------------------------------------------------------------------
807   ln_dynvor_ene = .false. !  enstrophy conserving scheme
808   ln_dynvor_ens = .false. !  energy conserving scheme
809   ln_dynvor_mix = .false. !  mixed scheme
810   ln_dynvor_een = .false. !  energy & enstrophy scheme
811      nn_een_e3f = 1          ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)
812   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes)  ! PLEASE DO NOT ACTIVATE
813/
814!-----------------------------------------------------------------------
815&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: zps)
816!-----------------------------------------------------------------------
817   ln_hpg_zco  = .false.   !  z-coordinate - full steps
818   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
819   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
820   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
821   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
822   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
823/
824!-----------------------------------------------------------------------
825&namdyn_spg    !   surface pressure gradient                            (default: NO)
826!-----------------------------------------------------------------------
827   ln_dynspg_exp  = .false.   ! explicit free surface
828   ln_dynspg_ts   = .false.   ! split-explicit free surface
829      ln_bt_fw      = .true.     ! Forward integration of barotropic Eqs.
830      ln_bt_av      = .true.     ! Time filtering of barotropic variables
831         nn_bt_flt     = 1          ! Time filter choice  = 0 None
832         !                          !                     = 1 Boxcar over   nn_baro sub-steps
833         !                          !                     = 2 Boxcar over 2*nn_baro  "    "
834      ln_bt_auto    = .true.     ! Number of sub-step defined from:
835         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed
836         nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds
837      rn_bt_alpha   = 0.1        ! Temporal diffusion parameter (if ln_bt_av=F)
838/
839!-----------------------------------------------------------------------
840&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO)
841!-----------------------------------------------------------------------
842   !                       !  Type of the operator :
843   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F
844   ln_dynldf_lap =  .false.    !    laplacian operator
845   ln_dynldf_blp =  .false.    !  bilaplacian operator
846   !                       !  Direction of action  :
847   ln_dynldf_lev =  .false.    !  iso-level
848   ln_dynldf_hor =  .false.    !  horizontal (geopotential)
849   ln_dynldf_iso =  .false.    !  iso-neutral
850   !                       !  Coefficient
851   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef
852   !                                !  =-30  read in eddy_viscosity_3D.nc file
853   !                                !  =-20  read in eddy_viscosity_2D.nc file
854   !                                !  =  0  constant
855   !                                !  = 10  F(k)=c1d
856   !                                !  = 20  F(i,j)=F(grid spacing)=c2d
857   !                                !  = 30  F(i,j,k)=c2d*c1d
858   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity)
859   !                                !  = 32  F(i,j,k)=F(local gridscale and deformation rate)
860   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km)
861   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s]
862   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s]
863   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s]
864   !                       !  Smagorinsky settings (nn_ahm_ijk_t  = 32) :
865   rn_csmc       = 3.5         !  Smagorinsky constant of proportionality
866   rn_minfac     = 1.0         !  multiplier of theorectical lower limit
867   rn_maxfac     = 1.0         !  multiplier of theorectical upper limit
868/
869
870!!======================================================================
871!!             Tracers & Dynamics vertical physics namelists
872!!======================================================================
873!!    namzdf        vertical physics
874!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric")
875!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke")
876!!    namzdf_gls    GLS vertical mixing                                 ("key_zdfgls")
877!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm")
878!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx")
879!!======================================================================
880!
881!-----------------------------------------------------------------------
882&namzdf        !   vertical physics
883!-----------------------------------------------------------------------
884   rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
885   rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
886   nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0)
887   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0)
888   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F)
889      nn_evdm     =    0        ! evd apply on tracer (=0) or on tracer and momentum (=1)
890      rn_avevd    =  100.       !  evd mixing coefficient [m2/s]
891   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F)
892      nn_npc      =    1        ! frequency of application of npc
893      nn_npcp     =  365        ! npc control print frequency
894   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
895      nn_zdfexp   =    3        ! number of sub-timestep for ln_zdfexp=T
896   ln_zdfqiao  = .false.   !  Enhanced wave vertical mixing Qiao (2010) (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
897/
898!-----------------------------------------------------------------------
899&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
900!-----------------------------------------------------------------------
901   rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity
902   rn_alp      =    5.     !  coefficient of the parameterization
903   nn_ric      =    2      !  coefficient of the parameterization
904   rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation
905   rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m)
906   rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m)
907   rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
908   rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
909   ln_mldw     =  .true.   !  Flag to use or not the mixed layer depth param.
910/
911!-----------------------------------------------------------------------
912&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
913!-----------------------------------------------------------------------
914   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
915   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
916   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
917   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
918   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
919   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
920   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
921                           !                 = 1 bounded by the local vertical scale factor
922                           !                 = 2 first vertical derivative of mixing length bounded by 1
923                           !                 = 3 as =2 with distinct disspipative an mixing length scale
924   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
925   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
926   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
927   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
928   rn_lc       =   0.15    !  coef. associated to Langmuir cells
929   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to near intertial waves
930                           !        = 0 no penetration
931                           !        = 1 add a tke source below the ML
932                           !        = 2 add a tke source just at the base of the ML
933                           !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T)
934   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
935   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
936                           !        = 0  constant 10 m length scale
937                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees
938/
939!-----------------------------------------------------------------------
940&namzdf_gls    !   GLS vertical diffusion                               ("key_zdfgls")
941!-----------------------------------------------------------------------
942   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
943   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
944   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
945   rn_clim_galp  = 0.267   !  galperin limit
946   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
947   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
948   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
949   rn_hsro       =  0.02   !  Minimum surface roughness
950   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2)
951   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3)
952   !                             ! =3 requires ln_wave=T
953   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
954   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
955   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
956   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
957/
958!-----------------------------------------------------------------------
959&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
960!-----------------------------------------------------------------------
961   rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
962   rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
963/
964!-----------------------------------------------------------------------
965&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
966!-----------------------------------------------------------------------
967   rn_htmx     = 500.      !  vertical decay scale for turbulence (meters)
968   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1)
969   rn_tfe      = 0.333     !  tidal dissipation efficiency
970   rn_me       = 0.2       !  mixing efficiency
971   ln_tmx_itf  = .true.    !  ITF specific parameterisation
972   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency
973/
974!-----------------------------------------------------------------------
975&namzdf_tmx_new !   internal wave-driven mixing parameterization        ("key_zdftmx_new" & "key_zdfddm")
976!-----------------------------------------------------------------------
977   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
978   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency
979   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F)
980/
981!!======================================================================
982!!                  ***  Miscellaneous namelists  ***
983!!======================================================================
984!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi)
985!!   namctl            Control prints
986!!   namsto            Stochastic parametrization of EOS
987!!======================================================================
988!
989!-----------------------------------------------------------------------
990&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
991!-----------------------------------------------------------------------
992   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
993                           !  buffer blocking send or immediate non-blocking sends, resp.
994   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
995   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold
996   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1)
997   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1)
998   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1)
999/
1000!-----------------------------------------------------------------------
1001&namctl        !   Control prints
1002!-----------------------------------------------------------------------
1003   ln_ctl      = .false.   !  trends control print (expensive!)
1004   nn_print    =    0      !  level of print (0 no extra print)
1005   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
1006   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
1007   nn_jctls    =    0      !  start j indice of control               over a subdomain)
1008   nn_jctle    =    0      !  end   j indice of control
1009   nn_isplt    =    1      !  number of processors in i-direction
1010   nn_jsplt    =    1      !  number of processors in j-direction
1011   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0)
1012   nn_diacfl   =    0      !  Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0)
1013/
1014!-----------------------------------------------------------------------
1015&namsto        ! Stochastic parametrization of EOS                      (default: NO)
1016!-----------------------------------------------------------------------
1017   ln_sto_eos  = .false.   ! stochastic equation of state
1018   nn_sto_eos  = 1         ! number of independent random walks
1019   rn_eos_stdxy= 1.4       ! random walk horz. standard deviation (in grid points)
1020   rn_eos_stdz = 0.7       ! random walk vert. standard deviation (in grid points)
1021   rn_eos_tcor = 1440.     ! random walk time correlation (in timesteps)
1022   nn_eos_ord  = 1         ! order of autoregressive processes
1023   nn_eos_flt  = 0         ! passes of Laplacian filter
1024   rn_eos_lim  = 2.0       ! limitation factor (default = 3.0)
1025   ln_rststo   = .false.   ! start from mean parameter (F) or from restart file (T)
1026   ln_rstseed  = .true.    ! read seed of RNG from restart file
1027   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
1028   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)
1029/
1030
1031!!======================================================================
1032!!                  ***  Diagnostics namelists  ***
1033!!======================================================================
1034!!   namtrd       dynamics and/or tracer trends                         (default F)
1035!!   namptr       Poleward Transport Diagnostics                        (default F)
1036!!   namhsb       Heat and salt budgets                                 (default F)
1037!!   namdiu       Cool skin and warm layer models                       (default F)
1038!!   namdiu       Cool skin and warm layer models                       (default F)
1039!!   namflo       float parameters                                      ("key_float")
1040!!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm")
1041!!   namdct       transports through some sections                      ("key_diadct")
1042!!   nam_diatmb   Top Middle Bottom Output                              (default F)
1043!!   nam_dia25h   25h Mean Output                                       (default F)
1044!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1045!!======================================================================
1046!
1047!-----------------------------------------------------------------------
1048&namtrd        !   trend diagnostics                                    (default F)
1049!-----------------------------------------------------------------------
1050   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1051   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1052   ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1053   ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet)
1054   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1055   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1056   ln_tra_trd  = .false.   ! (T) 3D tracer trend output
1057   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1058   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1059/
1060!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1061!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1062!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1063!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1064!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1065!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1066!!gm
1067!-----------------------------------------------------------------------
1068&namptr        !   Poleward Transport Diagnostic                         (default F)
1069!-----------------------------------------------------------------------
1070   ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F)
1071   ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not
1072/
1073!-----------------------------------------------------------------------
1074&namhsb        !  Heat and salt budgets                                  (default F)
1075!-----------------------------------------------------------------------
1076   ln_diahsb   = .false.   !  check the heat and salt budgets (T) or not (F)
1077/
1078!-----------------------------------------------------------------------
1079&namdiu        !   Cool skin and warm layer models                       (default F)
1080!-----------------------------------------------------------------------
1081   ln_diurnal      = .false.   !
1082   ln_diurnal_only = .false.   !
1083/
1084!-----------------------------------------------------------------------
1085&namflo        !   float parameters                                      ("key_float")
1086!-----------------------------------------------------------------------
1087   jpnfl       = 1         !  total number of floats during the run
1088   jpnnewflo   = 0         !  number of floats for the restart
1089   ln_rstflo   = .false.   !  float restart (T) or not (F)
1090   nn_writefl  =      75   !  frequency of writing in float output file
1091   nn_stockfl  =    5475   !  frequency of creation of the float restart file
1092   ln_argo     = .false.   !  Argo type floats (stay at the surface each 10 days)
1093   ln_flork4   = .false.   !  trajectories computed with a 4th order Runge-Kutta (T)
1094   !                       !  or computed with Blanke' scheme (F)
1095   ln_ariane   = .true.    !  Input with Ariane tool convention(T)
1096   ln_flo_ascii= .true.    !  Output with Ariane tool netcdf convention(F) or ascii file (T)
1097/
1098!-----------------------------------------------------------------------
1099&nam_diaharm   !   Harmonic analysis of tidal constituents               ("key_diaharm")
1100!-----------------------------------------------------------------------
1101    nit000_han = 1         ! First time step used for harmonic analysis
1102    nitend_han = 75        ! Last time step used for harmonic analysis
1103    nstep_han  = 15        ! Time step frequency for harmonic analysis
1104    tname(1)   = 'M2'      ! Name of tidal constituents
1105    tname(2)   = 'K1'
1106/
1107!-----------------------------------------------------------------------
1108&namdct        ! transports through some sections                        ("key_diadct")
1109!-----------------------------------------------------------------------
1110    nn_dct     = 15        !  time step frequency for transports computing
1111    nn_dctwri  = 15        !  time step frequency for transports writing
1112    nn_secdebug= 112       !      0 : no section to debug
1113    !                      !     -1 : debug all section
1114    !                      !  0 < n : debug section number n
1115/
1116!-----------------------------------------------------------------------
1117&nam_diatmb    !  Top Middle Bottom Output                               (default F)
1118!-----------------------------------------------------------------------
1119   ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not
1120/
1121!-----------------------------------------------------------------------
1122&nam_dia25h    !  25h Mean Output                                        (default F)
1123!-----------------------------------------------------------------------
1124   ln_dia25h   = .false.   ! Choose 25h mean output or not
1125/
1126!-----------------------------------------------------------------------
1127&namnc4        !   netcdf4 chunking and compression settings             ("key_netcdf4")
1128!-----------------------------------------------------------------------
1129   nn_nchunks_i=   4       !  number of chunks in i-dimension
1130   nn_nchunks_j=   4       !  number of chunks in j-dimension
1131   nn_nchunks_k=   31      !  number of chunks in k-dimension
1132   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1133   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs
1134   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1135   !                       !  (F) ignore chunking information and produce netcdf3-compatible files
1136/
1137
1138!!======================================================================
1139!!               ***  Observation & Assimilation  ***
1140!!======================================================================
1141!!   namobs       observation and model comparison
1142!!   nam_asminc   assimilation increments                               ('key_asminc')
1143!!======================================================================
1144!
1145!-----------------------------------------------------------------------
1146&namobs        !  observation usage switch
1147!-----------------------------------------------------------------------
1148   ln_diaobs   = .false.             ! Logical switch for the observation operator
1149   ln_t3d      = .false.             ! Logical switch for T profile observations
1150   ln_s3d      = .false.             ! Logical switch for S profile observations
1151   ln_sla      = .false.             ! Logical switch for SLA observations
1152   ln_sst      = .false.             ! Logical switch for SST observations
1153   ln_sic      = .false.             ! Logical switch for Sea Ice observations
1154   ln_vel3d    = .false.             ! Logical switch for velocity observations
1155   ln_altbias  = .false.             ! Logical switch for altimeter bias correction
1156   ln_nea      = .false.             ! Logical switch for rejection of observations near land
1157   ln_grid_global = .true.           ! Logical switch for global distribution of observations
1158   ln_grid_search_lookup = .false.   ! Logical switch for obs grid search w/lookup table
1159   ln_ignmis   = .true.              ! Logical switch for ignoring missing files
1160   ln_s_at_t   = .false.             ! Logical switch for computing model S at T obs if not there
1161   ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs
1162! All of the *files* variables below are arrays. Use namelist_cfg to add more files
1163   cn_profbfiles = 'profiles_01.nc'  ! Profile feedback input observation file names
1164   cn_slafbfiles = 'sla_01.nc'       ! SLA feedback input observation file names
1165   cn_sstfbfiles = 'sst_01.nc'       ! SST feedback input observation file names
1166   cn_sicfbfiles = 'sic_01.nc'       ! SIC feedback input observation file names
1167   cn_velfbfiles = 'vel_01.nc'       ! Velocity feedback input observation file names
1168   cn_altbiasfile = 'altbias.nc'     ! Altimeter bias input file name
1169   cn_gridsearchfile='gridsearch.nc' ! Grid search file name
1170   rn_gridsearchres = 0.5            ! Grid search resolution
1171   rn_dobsini  = 00010101.000000     ! Initial date in window YYYYMMDD.HHMMSS
1172   rn_dobsend  = 00010102.000000     ! Final date in window YYYYMMDD.HHMMSS
1173   nn_1dint    = 0                   ! Type of vertical interpolation method
1174   nn_2dint    = 0                   ! Type of horizontal interpolation method
1175   nn_msshc    = 0                   ! MSSH correction scheme
1176   rn_mdtcorr  = 1.61                ! MDT  correction
1177   rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction
1178   nn_profdavtypes = -1              ! Profile daily average types - array
1179   ln_sstbias  = .false.             !
1180   cn_sstbias_files = 'sstbias.nc'   !
1181/
1182!-----------------------------------------------------------------------
1183&nam_asminc    !   assimilation increments                              ('key_asminc')
1184!-----------------------------------------------------------------------
1185    ln_bkgwri  = .false.   !  Logical switch for writing out background state
1186    ln_trainc  = .false.   !  Logical switch for applying tracer increments
1187    ln_dyninc  = .false.   !  Logical switch for applying velocity increments
1188    ln_sshinc  = .false.   !  Logical switch for applying SSH increments
1189    ln_asmdin  = .false.   !  Logical switch for Direct Initialization (DI)
1190    ln_asmiau  = .false.   !  Logical switch for Incremental Analysis Updating (IAU)
1191    nitbkg     = 0         !  Timestep of background in [0,nitend-nit000-1]
1192    nitdin     = 0         !  Timestep of background for DI in [0,nitend-nit000-1]
1193    nitiaustr  = 1         !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1194    nitiaufin  = 15        !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1195    niaufn     = 0         !  Type of IAU weighting function
1196    ln_salfix  = .false.   !  Logical switch for ensuring that the sa > salfixmin
1197    salfixmin  = -9999     !  Minimum salinity after applying the increments
1198    nn_divdmp  = 0         !  Number of iterations of divergence damping operator
1199/
Note: See TracBrowser for help on using the repository browser.