source: branches/2016/dev_CNRS_2016/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/1_namelist_ref @ 7283

Last change on this file since 7283 was 7283, checked in by flavoni, 5 years ago

add default ORCA2_LIM3_PISCES, & remove old useless config

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