source: NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/cfgs/SHARED/namelist_ref @ 11518

Last change on this file since 11518 was 11518, checked in by clem, 14 months ago

add the final touch to the famous gaston's branch. More precisely, add the possibility to have melt ponds as input file when using bdy

  • Property svn:keywords set to Id
  • Property svn:mime-type set to text/x-fortran
File size: 104.8 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OCE :   Reference namelist_ref                                !!
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namtsd, namcrs, namc1d, namc1d_uvd)
5!! namelists    2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl,
6!!                                    namsbc_sas, namtra_qsr, namsbc_rnf,
7!!                                    namsbc_isf, namsbc_iscpl, namsbc_apr,
8!!                                    namsbc_ssr, namsbc_wave, namberg)
9!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide)
10!!              4 - top/bot boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl)
11!!              5 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_eiv, namtra_dmp)
12!!              6 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
13!!              7 - Vertical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm)
14!!              8 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb)
15!!              9 - Obs & Assim      (namobs, nam_asminc)
16!!             10 - miscellaneous    (nammpp, namctl, namsto)
17!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
18
19!!======================================================================
20!!              ***  Domain & Run management namelists  ***           !!
21!!                                                                    !!
22!!   namrun       parameters of the run
23!!   namdom       space and time domain
24!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
25!!   namwad       Wetting and drying                                    (default: OFF)
26!!   namtsd       data: temperature & salinity                          (default: OFF)
27!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
28!!   namc1d       1D configuration options                              ("key_c1d")
29!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
30!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
31!!======================================================================
32!
33!-----------------------------------------------------------------------
34&namrun        !   parameters of the run
35!-----------------------------------------------------------------------
36   nn_no       =       0   !  Assimilation cycle index
37   cn_exp      =  "ORCA2"  !  experience name
38   nn_it000    =       1   !  first time step
39   nn_itend    =    5840   !  last  time step (std 5840)
40   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
41   nn_time0    =       0   !  initial time of day in hhmm
42   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
43   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
44      nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T
45      nn_rstctl   =    0      !  restart control ==> activated only if ln_rstart=T
46      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
47      !                          !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
48      !                          !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart
49      cn_ocerst_in    = "restart"   !  suffix of ocean restart name (input)
50      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts
51      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output)
52      cn_ocerst_outdir = "."        !  directory in which to write output ocean restarts
53   ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model
54   nn_istate   =       0   !  output the initial state (1) or not (0)
55   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
56   nn_stock    =       0   !  used only if ln_rst_list = F: output restart freqeuncy (modulo referenced to 1)
57      !                          !    =  0 force to write restart files only at the end of the run
58      !                          !    = -1 do not do any restart
59   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
60   nn_write    =       0   !  used only if key_iomput is not defined: output frequency (modulo referenced to nn_it000)
61      !                          !    =  0 force to write output files only at the end of the run
62      !                          !    = -1 do not do any output file
63   ln_mskland  = .false.   !  mask land points in NetCDF outputs
64   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard
65   ln_clobber  = .true.    !  clobber (overwrite) an existing file
66   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
67   ln_xios_read = .FALSE.  !  use XIOS to read restart file (only for a single file restart)
68   nn_wxios = 0      !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output
69/
70!-----------------------------------------------------------------------
71&namdom        !   time and space domain
72!-----------------------------------------------------------------------
73   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time
74   rn_isfhmin  =    1.00   !  treshold [m] to discriminate grounding ice from floating ice
75   !
76   rn_rdt      = 5400.     !  time step for the dynamics and tracer
77   rn_atfp     =    0.1    !  asselin time filter parameter
78   !
79   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs)
80   !
81   ln_meshmask = .false.   !  =T create a mesh file
82/
83!-----------------------------------------------------------------------
84&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
85!-----------------------------------------------------------------------
86   ln_read_cfg = .false.   !  (=T) read the domain configuration file
87      !                    !  (=F) user defined configuration           (F => create/check namusr_def)
88      cn_domcfg = "domain_cfg"  ! domain configuration filename
89      !
90      ln_closea    = .false.    !  T => keep closed seas (defined by closea_mask field) in the 
91      !                         !       domain and apply special treatment of freshwater fluxes.
92      !                         !  F => suppress closed seas (defined by closea_mask field)
93      !                         !       from the bathymetry at runtime.
94      !                         !  If closea_mask field doesn't exist in the domain_cfg file
95      !                         !       then this logical does nothing.
96   ln_write_cfg = .false.   !  (=T) create the domain configuration file
97      cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename
98      !
99   ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present
100   !                       !  in netcdf input files, as the start j-row for reading
101/
102!-----------------------------------------------------------------------
103&namtsd        !    Temperature & Salinity Data  (init/dmp)             (default: OFF)
104!-----------------------------------------------------------------------
105   !                       ! =T  read T-S fields for:
106   ln_tsd_init = .false.         !  ocean initialisation
107   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp)
108   
109   cn_dir      = './'      !  root directory for the T-S data location
110   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
111   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
112   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
113   sn_tem = 'data_1m_potential_temperature_nomask',  -1.     , 'votemper',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
114   sn_sal = 'data_1m_salinity_nomask'             ,  -1.     , 'vosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
115/
116!-----------------------------------------------------------------------
117&namwad        !   Wetting and Drying (WaD)                             (default: OFF)
118!-----------------------------------------------------------------------
119   ln_wd_il    = .false.   !  T/F activation of iterative   limiter
120   ln_wd_dl    = .false.   !  T/F activation of directional limiter
121   ln_wd_dl_bc = .false.   !  T/F Directional limiteer Baroclinic option
122   ln_wd_dl_rmp = .false.  !  T/F Turn on directional limiter ramp
123   rn_wdmin0   =  0.30     !  depth at which WaD starts
124   rn_wdmin1   =  0.2      !  Minimum wet depth on dried cells
125   rn_wdmin2   =  0.0001   !  Tolerance of min wet depth on dried cells
126   rn_wdld     =  2.5      !  Land elevation below which WaD is allowed
127   nn_wdit     =   20      !  Max iterations for WaD limiter
128   rn_wd_sbcdep =  5.0     !  Depth at which to taper sbc fluxes
129   rn_wd_sbcfra =  0.999   !  Fraction of SBC fluxes at taper depth (Must be <1)
130/
131!-----------------------------------------------------------------------
132&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T)
133!-----------------------------------------------------------------------
134   nn_factx    = 3         !  Reduction factor of x-direction
135   nn_facty    = 3         !  Reduction factor of y-direction
136   nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT
137      !                    !  0, coarse grid is binned with preferential treatment of the north fold
138      !                    !  1, coarse grid is binned with centering at the equator
139      !                    !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.
140   ln_msh_crs  = .false.   ! =T create a mesh & mask file
141   nn_crs_kz   = 0         ! 0, MEAN of volume boxes
142      !                    ! 1, MAX of boxes
143      !                    ! 2, MIN of boxes
144   ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F )
145/
146!-----------------------------------------------------------------------
147&namc1d        !   1D configuration options                             ("key_c1d" default: PAPA station)
148!-----------------------------------------------------------------------
149   rn_lat1d    =      50   !  Column latitude
150   rn_lon1d    =    -145   !  Column longitude
151   ln_c1d_locpt =  .true.   ! Localization of 1D config in a grid (T) or independant point (F)
152/
153!-----------------------------------------------------------------------
154&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
155!-----------------------------------------------------------------------
156   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
157/
158!-----------------------------------------------------------------------
159&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
160!-----------------------------------------------------------------------
161   !                       !  =T read U-V fields for:
162   ln_uvd_init   = .false.       !  ocean initialisation
163   ln_uvd_dyndmp = .false.       !  U-V restoring
164
165   cn_dir      = './'      !  root directory for the U-V data location
166   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
167   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
168   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
169   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
170   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
171/
172
173!!======================================================================
174!!            ***  Surface Boundary Condition namelists  ***          !!
175!!                                                                    !!
176!!   namsbc          surface boundary condition manager                 (default: NO selection)
177!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
178!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
179!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
180!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
181!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
182!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
183!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
184!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
185!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
186!!   namsbc_isf      ice shelf melting/freezing                         (ln_isfcav  =T : read (ln_read_cfg=T) or set or usr_def_zgr )
187!!   namsbc_iscpl    coupling option between land ice model and ocean   (ln_isfcav  =T)
188!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
189!!   namberg         iceberg floats                                     (ln_icebergs=T)
190!!======================================================================
191!
192!-----------------------------------------------------------------------
193&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
194!-----------------------------------------------------------------------
195   nn_fsbc     = 2         !  frequency of SBC module call
196      !                    !  (control sea-ice & iceberg model call)
197                     ! Type of air-sea fluxes
198   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc)
199   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx )
200   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk )
201      !              ! Type of coupling (Ocean/Ice/Atmosphere) :
202   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 )
203   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 )
204   nn_components = 0       !  configuration of the opa-sas OASIS coupling
205      !                    !  =0 no opa-sas OASIS coupling: default single executable config.
206      !                    !  =1 opa-sas OASIS coupling: multi executable config., OPA component
207      !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component
208                     ! Sea-ice :
209   nn_ice      = 0         !  =0 no ice boundary condition   
210      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
211      !                    !  =2 or 3 automatically for SI3 or CICE    ("key_si3" or "key_cice")
212      !                    !          except in AGRIF zoom where it has to be specified
213   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges)
214      !                    !  =F levitating ice (no pressure, mass and salt exchanges)
215                     ! Misc. options of sbc :
216   ln_traqsr   = .false.   !  Light penetration in the ocean            (T => fill namtra_qsr)
217   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave
218   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
219   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked
220      !                    !     =1 global mean of e-p-r set to zero at each time step
221      !                    !     =2 annual global mean of e-p-r set to zero
222   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf)
223   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr )
224   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf & namsbc_iscpl)
225   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave)
226   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)
227   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)
228   nn_sdrift   =  0        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
229      !                    !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
230      !                    !   = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
231      !                    !   = 2 Phillips as (1) but using the wave frequency from a wave model
232   ln_tauwoc   = .false.   !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)
233   ln_tauw     = .false.   !  Activate ocean stress components from wave model
234   ln_stcor    = .false.   !  Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
235   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
236                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
237/
238!-----------------------------------------------------------------------
239&namsbc_flx    !   surface boundary condition : flux formulation        (ln_flx =T)
240!-----------------------------------------------------------------------
241   cn_dir      = './'      !  root directory for the fluxes data location
242   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
243   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
244   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
245   sn_utau     = 'utau'                  ,        24.        , 'utau'    , .false.     , .false., 'yearly'  , ''               , ''       , ''
246   sn_vtau     = 'vtau'                  ,        24.        , 'vtau'    , .false.     , .false., 'yearly'  , ''               , ''       , ''
247   sn_qtot     = 'qtot'                  ,        24.        , 'qtot'    , .false.     , .false., 'yearly'  , ''               , ''       , ''
248   sn_qsr      = 'qsr'                   ,        24.        , 'qsr'     , .false.     , .false., 'yearly'  , ''               , ''       , ''
249   sn_emp      = 'emp'                   ,        24.        , 'emp'     , .false.     , .false., 'yearly'  , ''               , ''       , ''
250/
251!-----------------------------------------------------------------------
252&namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T)
253!-----------------------------------------------------------------------
254   !                    !  bulk algorithm :
255   ln_NCAR     = .false.   ! "NCAR"      algorithm   (Large and Yeager 2008)
256   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003)
257   ln_COARE_3p5 = .false.   ! "COARE 3.5" algorithm   (Edson et al. 2013)
258   ln_ECMWF    = .false.   ! "ECMWF"     algorithm   (IFS cycle 31)
259      !
260      rn_zqt      = 10.       !  Air temperature & humidity reference height (m)
261      rn_zu       = 10.       !  Wind vector reference height (m)
262      ln_Cd_L12   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2012)
263      ln_Cd_L15   = .false.   !  air-ice drags = F(ice concentration) (Lupkes et al. 2015)
264      ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data
265      rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow)
266      rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.)
267      rn_vfac     = 0.        !  multiplicative factor for ocean & ice velocity used to
268      !                       !  calculate the wind stress (0.=absolute or 1.=relative winds)
269
270   cn_dir      = './'      !  root directory for the bulk data location
271   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
272   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
273   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! pairing  !    filename   !
274   sn_wndi     = 'u_10.15JUNE2009_fill'       ,    6.        , 'U_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Uwnd'   , ''
275   sn_wndj     = 'v_10.15JUNE2009_fill'       ,    6.        , 'V_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'  , 'Vwnd'   , ''
276   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'SWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
277   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,   24.        , 'LWDN_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
278   sn_tair     = 't_10.15JUNE2009_fill'       ,    6.        , 'T_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
279   sn_humi     = 'q_10.15JUNE2009_fill'       ,    6.        , 'Q_10_MOD',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
280   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
281   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
282   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
283   sn_tdif     = 'taudif_core'                ,   24         , 'taudif'  ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
284/
285!-----------------------------------------------------------------------
286&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
287!-----------------------------------------------------------------------
288   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data
289   ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models
290   !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
291   nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1)
292
293   !_____________!__________________________!____________!_____________!______________________!________!
294   !             !        description       !  multiple  !    vector   !       vector         ! vector !
295   !             !                          ! categories !  reference  !     orientation      ! grids  !
296!***   send    ***
297   sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   ''
298   sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
299   sn_snd_thick  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
300   sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T'
301   sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
302   sn_snd_crtw   =   'none'                 ,    'no'    ,     ''      ,         ''           , 'U,V'
303   sn_snd_ifrac  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
304   sn_snd_wlev   =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
305   sn_snd_cond   =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
306   sn_snd_thick1 =   'ice and snow'         ,    'no'    ,     ''      ,         ''           ,   ''
307   sn_snd_mpnd   =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
308   sn_snd_sstfrz =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
309   sn_snd_ttilyr =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
310!***  receive  ***
311   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
312   sn_rcv_taumod =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
313   sn_rcv_tau    =   'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward' ,  'U,V'
314   sn_rcv_dqnsdt =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
315   sn_rcv_qsr    =   'oce and ice'          ,    'no'    ,     ''      ,         ''           ,   ''
316   sn_rcv_qns    =   'oce and ice'          ,    'no'    ,     ''      ,         ''           ,   ''
317   sn_rcv_emp    =   'conservative'         ,    'no'    ,     ''      ,         ''           ,   ''
318   sn_rcv_rnf    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
319   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
320   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
321   sn_rcv_hsig   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
322   sn_rcv_iceflx =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
323   sn_rcv_mslp   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
324   sn_rcv_phioc  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
325   sn_rcv_sdrfx  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
326   sn_rcv_sdrfy  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
327   sn_rcv_wper   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
328   sn_rcv_wnum   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
329   sn_rcv_wstrf  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
330   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
331   sn_rcv_ts_ice =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
332   sn_rcv_isf    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
333   sn_rcv_icb    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
334   sn_rcv_tauwoc =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
335   sn_rcv_tauw   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
336   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
337/
338!-----------------------------------------------------------------------
339&namsbc_sas    !   Stand-Alone Surface module: ocean data               (SAS_SRC  only)
340!-----------------------------------------------------------------------
341   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
342      ln_3d_uve   = .false.   !  specify whether we are supplying a 3D u,v and e3 field
343      ln_read_frq = .false.   !  specify whether we must read frq or not
344
345   cn_dir      = './'      !  root directory for the ocean data location
346   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
347   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
348   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
349   sn_usp      = 'sas_grid_U'            ,       120.        , 'uos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
350   sn_vsp      = 'sas_grid_V'            ,       120.        , 'vos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
351   sn_tem      = 'sas_grid_T'            ,       120.        , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
352   sn_sal      = 'sas_grid_T'            ,       120.        , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
353   sn_ssh      = 'sas_grid_T'            ,       120.        , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
354   sn_e3t      = 'sas_grid_T'            ,       120.        , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
355   sn_frq      = 'sas_grid_T'            ,       120.        , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
356/
357!-----------------------------------------------------------------------
358&namsbc_iif    !   Ice-IF : use observed ice cover                      (nn_ice = 1)
359!-----------------------------------------------------------------------
360   cn_dir      = './'      !  root directory for the ice cover data location
361   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
362   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
363   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
364   sn_ice      ='ice_cover_clim.nc'      ,        -12.       ,'ice_cover',   .true.    , .true.  , 'yearly'  ,   ''            ,   ''     ,   ''
365/
366!-----------------------------------------------------------------------
367&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T)
368!-----------------------------------------------------------------------
369   !                       !  type of penetration                        (default: NO selection)
370   ln_qsr_rgb  = .false.      !  RGB light penetration (Red-Green-Blue)
371   ln_qsr_2bd  = .false.      !  2BD light penetration (two bands)
372   ln_qsr_bio  = .false.      !  bio-model light penetration
373   !                       !  RGB & 2BD choices:
374   rn_abs      =   0.58       !  RGB & 2BD: fraction absorbed in the very near surface
375   rn_si0      =   0.35       !  RGB & 2BD: shortess depth of extinction
376   nn_chldta   =      0       !  RGB : Chl data (=1) or cst value (=0)
377   rn_si1      =   23.0       !  2BD : longest depth of extinction
378   
379   cn_dir      = './'      !  root directory for the chlorophyl data location
380   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
381   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
382   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
383   sn_chl      ='chlorophyll'            ,        -1.        , 'CHLA'    ,   .true.    , .true. , 'yearly'  , ''               , ''       , ''
384/
385!-----------------------------------------------------------------------
386&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T)
387!-----------------------------------------------------------------------
388   nn_sstr     =     0     !  add a retroaction term to the surface heat flux (=1) or not (=0)
389      rn_dqdt     = -40.      !  magnitude of the retroaction on temperature   [W/m2/K]
390   nn_sssr     =     0     !  add a damping term to the surface freshwater flux (=2)
391      !                    !  or to SSS only (=1) or no damping term (=0)
392      rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
393      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2)
394      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
395
396   cn_dir      = './'      !  root directory for the SST/SSS data location
397   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
398   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
399   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
400   sn_sst      = 'sst_data'              ,        24.        ,  'sst'    ,    .false.  , .false., 'yearly'  ,    ''            ,    ''    ,     ''
401   sn_sss      = 'sss_data'              ,        -1.        ,  'sss'    ,    .true.   , .true. , 'yearly'  ,    ''            ,    ''    ,     ''
402/
403!-----------------------------------------------------------------------
404&namsbc_rnf    !   runoffs                                              (ln_rnf =T)
405!-----------------------------------------------------------------------
406   ln_rnf_mouth = .false.   !  specific treatment at rivers mouths
407      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T)
408      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
409   rn_rfact    =   1.e0    !  multiplicative factor for runoff
410   ln_rnf_depth = .false.   !  read in depth information for runoff
411   ln_rnf_tem  = .false.   !  read in temperature information for runoff
412   ln_rnf_sal  = .false.   !  read in salinity information for runoff
413   ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file
414      rn_rnf_max  = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
415      rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
416      nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0)
417
418   cn_dir      = './'      !  root directory for the runoff data location
419   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
420   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
421   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
422   sn_rnf      = 'runoff_core_monthly'   ,        -1.        , 'sorunoff',   .true.    , .true. , 'yearly'  , ''               , ''       , ''
423   sn_cnf      = 'runoff_core_monthly'   ,         0.        , 'socoefr0',   .false.   , .true. , 'yearly'  , ''               , ''       , ''
424   sn_s_rnf    = 'runoffs'               ,        24.        , 'rosaline',   .true.    , .true. , 'yearly'  , ''               , ''       , ''
425   sn_t_rnf    = 'runoffs'               ,        24.        , 'rotemper',   .true.    , .true. , 'yearly'  , ''               , ''       , ''
426   sn_dep_rnf  = 'runoffs'               ,         0.        , 'rodepth' ,   .false.   , .true. , 'yearly'  , ''               , ''       , ''
427/
428!-----------------------------------------------------------------------
429&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T)
430!-----------------------------------------------------------------------
431   rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/
432   ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F)
433   ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data
434
435   cn_dir = './'        !  root directory for the Patm data location
436   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
437   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
438   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
439   sn_apr      = 'patm'                  ,         -1.       ,'somslpre' ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,      ''
440/
441!-----------------------------------------------------------------------
442&namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T)
443!-----------------------------------------------------------------------             or set or usr_def_zgr )
444   !                 ! type of top boundary layer
445   nn_isf      = 1         !  ice shelf melting/freezing
446                           !  1 = presence of ISF   ;  2 = bg03 parametrisation
447                           !  3 = rnf file for ISF  ;  4 = ISF specified freshwater flux
448                           !  options 1 and 4 need ln_isfcav = .true. (domzgr)
449      !              !  nn_isf = 1 or 2 cases:
450      rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula
451      rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula
452      !              !  nn_isf = 1 or 4 cases:
453      rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008)
454      !                       ! 0 => thickness of the tbl = thickness of the first wet cell
455      !              ! nn_isf = 1 case
456      nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006)
457      !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015)
458      nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s)
459      !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010)
460      !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999)
461
462   !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________!
463   !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
464   !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
465!* nn_isf = 4 case
466   sn_fwfisf   = 'rnfisf'    ,         -12.      ,'sowflisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
467!* nn_isf = 3 case
468   sn_rnfisf   = 'rnfisf'    ,         -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
469!* nn_isf = 2 and 3 cases
470   sn_depmax_isf ='rnfisf'   ,         -12.      ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
471   sn_depmin_isf ='rnfisf'   ,         -12.      ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
472!* nn_isf = 2 case
473   sn_Leff_isf = 'rnfisf'    ,         -12.      ,'Leff'     ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
474/
475!-----------------------------------------------------------------------
476&namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T)
477!-----------------------------------------------------------------------             or set or usr_def_zgr )
478   nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells)
479   ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl)
480   nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency)
481/
482!-----------------------------------------------------------------------
483&namsbc_wave   ! External fields from wave model                        (ln_wave=T)
484!-----------------------------------------------------------------------
485   cn_dir      = './'      !  root directory for the waves data location
486   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
487   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
488   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
489   sn_cdg      =  'sdw_ecwaves_orca2'    ,        6.         , 'drag_coeff' ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
490   sn_usd      =  'sdw_ecwaves_orca2'    ,        6.         , 'u_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
491   sn_vsd      =  'sdw_ecwaves_orca2'    ,        6.         , 'v_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
492   sn_hsw      =  'sdw_ecwaves_orca2'    ,        6.         , 'hs'         ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
493   sn_wmp      =  'sdw_ecwaves_orca2'    ,        6.         , 'wmp'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
494   sn_wfr      =  'sdw_ecwaves_orca2'    ,        6.         , 'wfr'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
495   sn_wnum     =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_num'   ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
496   sn_tauwoc   =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
497   sn_tauwx    =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
498   sn_tauwy    =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
499/
500!-----------------------------------------------------------------------
501&namberg       !   iceberg parameters                                   (default: OFF)
502!-----------------------------------------------------------------------
503   ln_icebergs = .false.      ! activate iceberg floats (force =F with "key_agrif")
504   !
505   !                          ! diagnostics:
506   ln_bergdia        = .true.        ! Calculate budgets
507   nn_verbose_level  = 0             ! Turn on more verbose output if level > 0
508   nn_verbose_write  = 15            ! Timesteps between verbose messages
509   nn_sample_rate    = 1             ! Timesteps between sampling for trajectory storage
510   !
511   !                          ! iceberg setting:
512   !                                 ! Initial mass required for an iceberg of each class
513   rn_initial_mass   = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
514   !                                 ! Proportion of calving mass to apportion to each class
515   rn_distribution   = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
516   !                                 ! Ratio between effective and real iceberg mass (non-dim)
517   !                                 ! i.e. number of icebergs represented at a point
518   rn_mass_scaling   = 2000., 200., 50., 20., 10., 5., 2., 1., 1., 1.
519                                     ! thickness of newly calved bergs (m)
520   rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
521   !
522   rn_rho_bergs            = 850.    ! Density of icebergs
523   rn_LoW_ratio            = 1.5     ! Initial ratio L/W for newly calved icebergs
524   ln_operator_splitting   = .true.  ! Use first order operator splitting for thermodynamics
525   rn_bits_erosion_fraction = 0.     ! Fraction of erosion melt flux to divert to bergy bits
526   rn_sicn_shift           = 0.      ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
527   ln_passive_mode         = .false. ! iceberg - ocean decoupling
528   nn_test_icebergs        =  10     ! Create test icebergs of this class (-1 = no)
529   !                                 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
530   rn_test_box             = 108.0,  116.0, -66.0, -58.0
531   ln_use_calving          = .false. ! Use calving data even when nn_test_icebergs > 0
532   rn_speed_limit          = 0.      ! CFL speed limit for a berg
533
534   cn_dir      = './'      !  root directory for the calving data location
535   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
536   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
537   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
538   sn_icb     =  'calving'              ,         -1.        ,'calvingmask',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
539/
540
541!!======================================================================
542!!               ***  Lateral boundary condition  ***                 !!
543!!                                                                    !!
544!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
545!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
546!!   nam_tide      Tidal forcing                                        (default: OFF)
547!!   nambdy        Unstructured open boundaries                         (default: OFF)
548!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
549!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
550!!======================================================================
551!
552!-----------------------------------------------------------------------
553&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
554!-----------------------------------------------------------------------
555   !                       !  free slip  !   partial slip  !   no slip   ! strong slip
556   rn_shlat    =  -9999.   !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
557   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs.
558/
559!-----------------------------------------------------------------------
560&namagrif      !  AGRIF zoom                                            ("key_agrif")
561!-----------------------------------------------------------------------
562   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
563   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
564   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
565   ln_chk_bathy  = .false. !  =T  check the parent bathymetry
566/
567!-----------------------------------------------------------------------
568&nam_tide      !   tide parameters                                      (default: OFF)
569!-----------------------------------------------------------------------
570   ln_tide     = .false.      ! Activate tides
571      ln_tide_pot   = .true.                !  use tidal potential forcing
572         ln_scal_load  = .false.               ! Use scalar approximation for
573            rn_scal_load = 0.094               !     load potential
574         ln_read_load  = .false.               ! Or read load potential from file
575            cn_tide_load = 'tide_LOAD_grid_T.nc'  ! filename for load potential
576            !     
577      ln_tide_ramp  = .false.               !  Use linear ramp for tides at startup
578         rdttideramp   =    0.                 !  ramp duration in days
579      clname(1)     = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg
580/
581!-----------------------------------------------------------------------
582&nambdy        !  unstructured open boundaries                          (default: OFF)
583!-----------------------------------------------------------------------
584   ln_bdy         = .false.   !  Use unstructured open boundaries
585   nb_bdy         = 0         !  number of open boundary sets
586   ln_coords_file = .true.    !  =T : read bdy coordinates from file
587      cn_coords_file = 'coordinates.bdy.nc'  !  bdy coordinates files
588   ln_mask_file   = .false.   !  =T : read mask from file
589      cn_mask_file = ''        !  name of mask file (if ln_mask_file=.TRUE.)
590   cn_dyn2d    = 'none'       !
591   nn_dyn2d_dta   =  0        !  = 0, bdy data are equal to the initial state
592      !                       !  = 1, bdy data are read in 'bdydata   .nc' files
593      !                       !  = 2, use tidal harmonic forcing data from files
594      !                       !  = 3, use external data AND tidal harmonic forcing
595   cn_dyn3d      =  'none'    !
596   nn_dyn3d_dta  =  0         !  = 0, bdy data are equal to the initial state
597   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
598   cn_tra        =  'none'    !
599   nn_tra_dta    =  0         !  = 0, bdy data are equal to the initial state
600   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
601   cn_ice        =  'none'    !
602   nn_ice_dta    =  0         !  = 0, bdy data are equal to the initial state
603   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
604   !
605   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers
606   ln_dyn3d_dmp  =.false.     !  open boundary condition for baroclinic velocities
607   rn_time_dmp   =  1.        !  Damping time scale in days
608   rn_time_dmp_out = 1.       !  Outflow damping time scale
609   nn_rimwidth   = 10         !  width of the relaxation zone
610   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter)
611   nn_volctl     =  1         !  = 0, the total water flux across open boundaries is zero
612/
613!-----------------------------------------------------------------------
614&nambdy_dta    !  open boundaries - external data                       (see nam_bdy)
615!-----------------------------------------------------------------------
616   ln_zinterp  = .false.      !  T if a vertical interpolation is required. Variables gdep[tuv] and e3[tuv] must exist in the file
617   !                          !  automatically defined to T if the number of vertical levels in bdy dta /= jpk
618   ln_full_vel = .false.      !  T if [uv]3d are "full" velocities and not only its baroclinic components
619   !                          !  in this case, baroclinic and barotropic velocities will be recomputed -> [uv]2d not needed
620   !
621   cn_dir      = 'bdydta/'    !  root directory for the BDY data location
622   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
623   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
624   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
625   bn_ssh      = 'amm12_bdyT_u2d'        ,         24.       , 'sossheig',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
626   bn_u2d      = 'amm12_bdyU_u2d'        ,         24.       , 'vobtcrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
627   bn_v2d      = 'amm12_bdyV_u2d'        ,         24.       , 'vobtcrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
628   bn_u3d      = 'amm12_bdyU_u3d'        ,         24.       , 'vozocrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
629   bn_v3d      = 'amm12_bdyV_u3d'        ,         24.       , 'vomecrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
630   bn_tem      = 'amm12_bdyT_tra'        ,         24.       , 'votemper',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
631   bn_sal      = 'amm12_bdyT_tra'        ,         24.       , 'vosaline',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
632!* for si3
633   bn_a_i      = 'amm12_bdyT_ice'        ,         24.       , 'siconc'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
634   bn_h_i      = 'amm12_bdyT_ice'        ,         24.       , 'sithic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
635   bn_h_s      = 'amm12_bdyT_ice'        ,         24.       , 'snthic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
636   bn_t_i      = 'NOT USED'              ,         24.       , 'sitemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
637   bn_t_s      = 'NOT USED'              ,         24.       , 'sntemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
638   bn_tsu      = 'NOT USED'              ,         24.       , 'sittop'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
639   bn_s_i      = 'NOT USED'              ,         24.       , 'sisalt'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
640   ! melt ponds (be careful, bn_aip is the pond concentration (not fraction), so it differs from rn_iceapnd)
641   bn_aip      = 'NOT USED'              ,         24.       , 'siapnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
642   bn_hip      = 'NOT USED'              ,         24.       , 'sihpnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
643   ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds
644   rn_ice_tem  = 270.         !  arbitrary temperature               of incoming sea ice
645   rn_ice_sal  = 10.          !       --   salinity                            --
646   rn_ice_age  = 30.          !       --   age                                 --
647   rn_ice_apnd = 0.2          !       --   pond fraction = a_ip/a_i            --
648   rn_ice_hpnd = 0.05         !       --   pond depth                          --
649/
650!-----------------------------------------------------------------------
651&nambdy_tide   !  tidal forcing at open boundaries                      (default: OFF)
652!-----------------------------------------------------------------------
653   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files
654   ln_bdytide_2ddta = .false.                   !
655   ln_bdytide_conj  = .false.                   !
656/
657
658!!======================================================================
659!!                ***  Top/Bottom boundary condition  ***             !!
660!!                                                                    !!
661!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
662!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
663!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
664!!   nambbc        bottom temperature boundary condition                (default: OFF)
665!!   nambbl        bottom boundary layer scheme                         (default: OFF)
666!!======================================================================
667!
668!-----------------------------------------------------------------------
669&namdrg        !   top/bottom drag coefficient                          (default: NO selection)
670!-----------------------------------------------------------------------
671   ln_OFF      = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot
672   ln_lin      = .false.   !      linear  drag: Cd = Cd0 Uc0                   &   namdrg_top)
673   ln_non_lin  = .false.   !  non-linear  drag: Cd = Cd0 |U|
674   ln_loglayer = .false.   !  logarithmic drag: Cd = vkarmn/log(z/z0) |U|
675   !
676   ln_drgimp   = .true.    !  implicit top/bottom friction flag
677/
678!-----------------------------------------------------------------------
679&namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T)
680!-----------------------------------------------------------------------
681   rn_Cd0      =  1.e-3    !  drag coefficient [-]
682   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)
683   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag)
684   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases)
685   rn_z0       =  3.0e-3   !  roughness [m] (ln_loglayer=T)
686   ln_boost    = .false.   !  =T regional boost of Cd0 ; =F constant
687      rn_boost =  50.         !  local boost factor  [-]
688/
689!-----------------------------------------------------------------------
690&namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F)
691!-----------------------------------------------------------------------
692   rn_Cd0      =  1.e-3    !  drag coefficient [-]
693   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)
694   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag)
695   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases)
696   rn_z0       =  3.e-3    !  roughness [m] (ln_loglayer=T)
697   ln_boost    = .false.   !  =T regional boost of Cd0 ; =F constant
698      rn_boost =  50.         !  local boost factor  [-]
699/
700!-----------------------------------------------------------------------
701&nambbc        !   bottom temperature boundary condition                (default: OFF)
702!-----------------------------------------------------------------------
703   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
704      nn_geoflx     = 2       !  geothermal heat flux: = 1 constant flux
705      !                       !                        = 2 read variable flux [mW/m2]
706      rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux       [mW/m2]
707
708   cn_dir      = './'      !  root directory for the geothermal data location
709   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
710   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
711   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
712   sn_qgh      ='geothermal_heating.nc'  ,       -12.        , 'heatflow',   .false.   , .true. , 'yearly'  ,   ''             ,   ''     ,   ''
713/
714!-----------------------------------------------------------------------
715&nambbl        !   bottom boundary layer scheme                         (default: OFF)
716!-----------------------------------------------------------------------
717   ln_trabbl   = .false.   !  Bottom Boundary Layer parameterisation flag
718      nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0)
719      nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0)
720      rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s]
721      rn_gambbl   =  10.      !  advective bbl coefficient                 [s]
722/
723
724!!======================================================================
725!!                        Tracer (T-S) namelists                      !!
726!!                                                                    !!
727!!   nameos        equation of state                                    (default: NO selection)
728!!   namtra_adv    advection scheme                                     (default: NO selection)
729!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
730!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
731!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
732!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
733!!======================================================================
734!
735!-----------------------------------------------------------------------
736&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
737!-----------------------------------------------------------------------
738   ln_teos10   = .false.         !  = Use TEOS-10
739   ln_eos80    = .false.         !  = Use EOS80
740   ln_seos     = .false.         !  = Use S-EOS (simplified Eq.)
741                                 !
742   !                     ! S-EOS coefficients (ln_seos=T):
743   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
744   rn_a0       =  1.6550e-1      !  thermal expension coefficient
745   rn_b0       =  7.6554e-1      !  saline  expension coefficient
746   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
747   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
748   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
749   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
750   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
751/
752!-----------------------------------------------------------------------
753&namtra_adv    !   advection scheme for tracer                          (default: NO selection)
754!-----------------------------------------------------------------------
755   ln_traadv_OFF = .false. !  No tracer advection
756   ln_traadv_cen = .false. !  2nd order centered scheme
757      nn_cen_h   =  4            !  =2/4, horizontal 2nd order CEN / 4th order CEN
758      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT
759   ln_traadv_fct = .false. !  FCT scheme
760      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order
761      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order
762   ln_traadv_mus = .false. !  MUSCL scheme
763      ln_mus_ups = .false.       !  use upstream scheme near river mouths
764   ln_traadv_ubs = .false. !  UBS scheme
765      nn_ubs_v   =  2            !  =2  , vertical 2nd order FCT / COMPACT 4th order
766   ln_traadv_qck = .false. !  QUICKEST scheme
767/
768!-----------------------------------------------------------------------
769&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection)
770!-----------------------------------------------------------------------
771   !                       !  Operator type:
772   ln_traldf_OFF   = .false.   !  No explicit diffusion
773   ln_traldf_lap   = .false.   !    laplacian operator
774   ln_traldf_blp   = .false.   !  bilaplacian operator
775   !
776   !                       !  Direction of action:
777   ln_traldf_lev   = .false.   !  iso-level
778   ln_traldf_hor   = .false.   !  horizontal  (geopotential)
779   ln_traldf_iso   = .false.   !  iso-neutral (standard operator)
780   ln_traldf_triad = .false.   !  iso-neutral (triad    operator)
781   !
782   !                       !  iso-neutral options:       
783   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators)
784   rn_slpmax       =  0.01     !  slope limit                           (both operators)
785   ln_triad_iso    = .false.   !  pure horizontal mixing in ML              (triad only)
786   rn_sw_triad     = 1         !  =1 switching triad ; =0 all 4 triads used (triad only)
787   ln_botmix_triad = .false.   !  lateral mixing on bottom                  (triad only)
788   !
789   !                       !  Coefficients:
790   nn_aht_ijk_t    = 0         !  space/time variation of eddy coefficient:
791      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file
792      !                             !   =  0           constant
793      !                             !   = 10 F(k)      =ldf_c1d
794      !                             !   = 20 F(i,j)    =ldf_c2d
795      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
796      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
797      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
798      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)
799      !                             !                           or   = 1/12 Ud*Ld^3 (blp case)
800      rn_Ud        = 0.01           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
801      rn_Ld        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10)
802/
803!-----------------------------------------------------------------------
804&namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF)
805!-----------------------------------------------------------------------
806   ln_mle      = .false.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation
807   rn_ce       = 0.06      ! magnitude of the MLE (typical value: 0.06 to 0.08)
808   nn_mle      = 1         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
809   rn_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_mle=0)
810   rn_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0)
811   rn_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
812   nn_mld_uv   = 0         ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
813   nn_conv     = 0         ! =1 no MLE in case of convection ; =0 always MLE
814   rn_rho_c_mle = 0.01      ! delta rho criterion used to calculate MLD for FK
815/
816!-----------------------------------------------------------------------
817&namtra_eiv    !   eddy induced velocity param.                         (default: OFF)
818!-----------------------------------------------------------------------
819   ln_ldfeiv   = .false.   ! use eddy induced velocity parameterization
820      !
821      !                        !  Coefficients:
822      nn_aei_ijk_t    = 0           !  space/time variation of eddy coefficient:
823      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
824      !                             !   =  0           constant
825      !                             !   = 10 F(k)      =ldf_c1d
826      !                             !   = 20 F(i,j)    =ldf_c2d
827      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
828      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
829      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le
830      rn_Ue        = 0.02           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
831      rn_Le        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10)
832      !
833      ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities
834/
835!-----------------------------------------------------------------------
836&namtra_dmp    !   tracer: T & S newtonian damping                      (default: OFF)
837!-----------------------------------------------------------------------
838   ln_tradmp   =  .false.  !  add a damping term (using resto.nc coef.)
839      nn_zdmp  =    0         !  vertical shape =0    damping throughout the water column
840      !                       !                 =1 no damping in the mixing layer (kz  criteria)
841      !                       !                 =2 no damping in the mixed  layer (rho crieria)
842      cn_resto = 'resto.nc'   !  Name of file containing restoration coeff. field (use dmp_tools to create this)
843/
844
845!!======================================================================
846!!                      ***  Dynamics namelists  ***                  !!
847!!                                                                    !!
848!!   nam_vvl       vertical coordinate options                          (default: z-star)
849!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
850!!   namdyn_vor    advection scheme                                     (default: NO selection)
851!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
852!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
853!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
854!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
855!!======================================================================
856!
857!-----------------------------------------------------------------------
858&nam_vvl       !   vertical coordinate options                          (default: z-star)
859!-----------------------------------------------------------------------
860   ln_vvl_zstar  = .true.           !  z-star vertical coordinate
861   ln_vvl_ztilde = .false.          !  z-tilde vertical coordinate: only high frequency variations
862   ln_vvl_layer  = .false.          !  full layer vertical coordinate
863   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
864   ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator
865   rn_ahe3       =  0.0             !  thickness diffusion coefficient
866   rn_rst_e3t    = 30.0             !  ztilde to zstar restoration timescale [days]
867   rn_lf_cutoff  =  5.0             !  cutoff frequency for low-pass filter  [days]
868   rn_zdef_max   =  0.9             !  maximum fractional e3t deformation
869   ln_vvl_dbg    = .true.           !  debug prints    (T/F)
870/
871!-----------------------------------------------------------------------
872&namdyn_adv    !   formulation of the momentum advection                (default: NO selection)
873!-----------------------------------------------------------------------
874   ln_dynadv_OFF = .false. !  linear dynamics (no momentum advection)
875   ln_dynadv_vec = .false. !  vector form - 2nd centered scheme
876     nn_dynkeg     = 0        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction
877   ln_dynadv_cen2 = .false. !  flux form - 2nd order centered scheme
878   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
879/
880!-----------------------------------------------------------------------
881&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection)
882!-----------------------------------------------------------------------
883   ln_dynvor_ene = .false. !  energy    conserving scheme
884   ln_dynvor_ens = .false. !  enstrophy conserving scheme
885   ln_dynvor_mix = .false. !  mixed scheme
886   ln_dynvor_enT = .false. !  energy conserving scheme (T-point)
887   ln_dynvor_eeT = .false. !  energy conserving scheme (een using e3t)
888   ln_dynvor_een = .false. !  energy & enstrophy scheme
889      nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4
890      !                       ! =1  e3f = mi(mj(e3t))/mi(mj( tmask))
891   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T)        ==>>> PLEASE DO NOT ACTIVATE
892      !                    !  (f-point vorticity schemes only)
893/
894!-----------------------------------------------------------------------
895&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: NO selection)
896!-----------------------------------------------------------------------
897   ln_hpg_zco  = .false.   !  z-coordinate - full steps
898   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
899   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
900   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
901   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
902   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
903/
904!-----------------------------------------------------------------------
905&namdyn_spg    !   surface pressure gradient                            (default: NO selection)
906!-----------------------------------------------------------------------
907   ln_dynspg_exp  = .false.   ! explicit free surface
908   ln_dynspg_ts   = .false.   ! split-explicit free surface
909      ln_bt_fw      = .true.     ! Forward integration of barotropic Eqs.
910      ln_bt_av      = .true.     ! Time filtering of barotropic variables
911         nn_bt_flt     = 1          ! Time filter choice  = 0 None
912         !                          !                     = 1 Boxcar over   nn_baro sub-steps
913         !                          !                     = 2 Boxcar over 2*nn_baro  "    "
914      ln_bt_auto    = .true.     ! Number of sub-step defined from:
915         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed
916         nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds
917      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F)
918/
919!-----------------------------------------------------------------------
920&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO selection)
921!-----------------------------------------------------------------------
922   !                       !  Type of the operator :
923   ln_dynldf_OFF = .false.     !  No operator (i.e. no explicit diffusion)
924   ln_dynldf_lap = .false.     !    laplacian operator
925   ln_dynldf_blp = .false.     !  bilaplacian operator
926   !                       !  Direction of action  :
927   ln_dynldf_lev = .false.     !  iso-level
928   ln_dynldf_hor = .false.     !  horizontal  (geopotential)
929   ln_dynldf_iso = .false.     !  iso-neutral (lap only)
930   !                       !  Coefficient
931   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coefficient :
932      !                             !  =-30  read in eddy_viscosity_3D.nc file
933      !                             !  =-20  read in eddy_viscosity_2D.nc file
934      !                             !  =  0  constant
935      !                             !  = 10  F(k)=c1d
936      !                             !  = 20  F(i,j)=F(grid spacing)=c2d
937      !                             !  = 30  F(i,j,k)=c2d*c1d
938      !                             !  = 31  F(i,j,k)=F(grid spacing and local velocity)
939      !                             !  = 32  F(i,j,k)=F(local gridscale and deformation rate)
940      !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case)
941      !                             !                            or  = 1/12 Uv*Lv^3 (blp case)
942      rn_Uv      = 0.1              !  lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30)
943      rn_Lv      = 10.e+3           !  lateral viscous length   [m]   (nn_ahm_ijk_t= 0, 10)
944      !                       !  Smagorinsky settings  (nn_ahm_ijk_t= 32) :
945      rn_csmc       = 3.5         !  Smagorinsky constant of proportionality
946      rn_minfac     = 1.0         !  multiplier of theorectical lower limit
947      rn_maxfac     = 1.0         !  multiplier of theorectical upper limit
948      !                       !  iso-neutral laplacian operator (ln_dynldf_iso=T) :
949      rn_ahm_b      = 0.0         !  background eddy viscosity  [m2/s]
950/
951!-----------------------------------------------------------------------
952&namdta_dyn    !   offline ocean input files                            (OFF_SRC only)
953!-----------------------------------------------------------------------
954   ln_dynrnf       =  .false.    !  runoffs option enabled (T) or not (F)
955   ln_dynrnf_depth =  .false.    !  runoffs is spread in vertical (T) or not (F)
956!   fwbcorr        = 3.786e-06   !  annual global mean of empmr for ssh correction
957
958   cn_dir      = './'      !  root directory for the ocean data location
959   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
960   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
961   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
962   sn_tem      = 'dyna_grid_T'           ,       120.        , 'votemper'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
963   sn_sal      = 'dyna_grid_T'           ,       120.        , 'vosaline'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
964   sn_mld      = 'dyna_grid_T'           ,       120.        , 'somixhgt'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
965   sn_emp      = 'dyna_grid_T'           ,       120.        , 'sowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
966   sn_fmf      = 'dyna_grid_T'           ,       120.        , 'iowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
967   sn_ice      = 'dyna_grid_T'           ,       120.        , 'soicecov'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
968   sn_qsr      = 'dyna_grid_T'           ,       120.        , 'soshfldo'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
969   sn_wnd      = 'dyna_grid_T'           ,       120.        , 'sowindsp'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
970   sn_uwd      = 'dyna_grid_U'           ,       120.        , 'uocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
971   sn_vwd      = 'dyna_grid_V'           ,       120.        , 'vocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
972   sn_wwd      = 'dyna_grid_W'           ,       120.        , 'wocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
973   sn_avt      = 'dyna_grid_W'           ,       120.        , 'voddmavs'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
974   sn_ubl      = 'dyna_grid_U'           ,       120.        , 'sobblcox'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
975   sn_vbl      = 'dyna_grid_V'           ,       120.        , 'sobblcoy'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
976/
977
978!!======================================================================
979!!                     vertical physics namelists                     !!
980!!                                                                    !!
981!!    namzdf        vertical physics manager                            (default: NO selection)
982!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
983!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
984!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
985!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
986!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
987!!======================================================================
988!
989!-----------------------------------------------------------------------
990&namzdf        !   vertical physics manager                             (default: NO selection)
991!-----------------------------------------------------------------------
992   !                       ! adaptive-implicit vertical advection
993   ln_zad_Aimp = .false.      !  Courant number dependent scheme (Shchepetkin 2015)
994   !
995   !                       ! type of vertical closure (required)
996   ln_zdfcst   = .false.      !  constant mixing
997   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric)
998   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke)
999   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls)
1000   ln_zdfosm   = .false.      !  OSMOSIS BL closure                     (T =>   fill namzdf_osm)
1001   !
1002   !                       ! convection
1003   ln_zdfevd   = .false.      !  enhanced vertical diffusion
1004      nn_evdm     =    0         ! apply on tracer (=0) or on tracer and momentum (=1)
1005      rn_evd      =  100.        ! mixing coefficient [m2/s]
1006   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm
1007      nn_npc      =    1         ! frequency of application of npc
1008      nn_npcp     =  365         ! npc control print frequency
1009   !
1010   ln_zdfddm   = .false.   ! double diffusive mixing
1011      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity)
1012      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio
1013   !
1014   !                       ! gravity wave-driven vertical mixing
1015   ln_zdfiwm   = .false.      ! internal wave-induced mixing            (T =>   fill namzdf_iwm)
1016   ln_zdfswm   = .false.      ! surface  wave-induced mixing            (T => ln_wave=ln_sdw=T )
1017   !
1018   !                       ! coefficients
1019   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F)
1020   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F)
1021   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0)
1022   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0)
1023/
1024!-----------------------------------------------------------------------
1025&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T)
1026!-----------------------------------------------------------------------
1027   rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity
1028   rn_alp      =    5.     !  coefficient of the parameterization
1029   nn_ric      =    2      !  coefficient of the parameterization
1030   ln_mldw     =  .false.  !  enhanced mixing in the Ekman layer
1031      rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation
1032      rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m)
1033      rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m)
1034      rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
1035      rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
1036/
1037!-----------------------------------------------------------------------
1038&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T)
1039!-----------------------------------------------------------------------
1040   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
1041   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
1042   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
1043   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
1044   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
1045   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
1046   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
1047   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
1048   !                       !                 = 1 bounded by the local vertical scale factor
1049   !                       !                 = 2 first vertical derivative of mixing length bounded by 1
1050   !                       !                 = 3 as =2 with distinct dissipative an mixing length scale
1051   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
1052   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
1053   ln_drg      = .false.   !  top/bottom friction added as boundary condition of TKE
1054   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
1055      rn_lc       =   0.15    !  coef. associated to Langmuir cells
1056   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs
1057                              !        = 0 none ; = 1 add a tke source below the ML
1058                              !        = 2 add a tke source just at the base of the ML
1059                              !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T)
1060      rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
1061      nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
1062                              !        = 0  constant 10 m length scale
1063                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees
1064      rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4   
1065/
1066!-----------------------------------------------------------------------
1067&namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls =T)
1068!-----------------------------------------------------------------------
1069   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
1070   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
1071   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
1072   rn_clim_galp  = 0.267   !  galperin limit
1073   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
1074   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
1075   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
1076   rn_hsro       =  0.02   !  Minimum surface roughness
1077   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met>1)
1078   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3)
1079   !                             ! =3 requires ln_wave=T
1080   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
1081   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
1082   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
1083   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
1084/
1085!-----------------------------------------------------------------------
1086&namzdf_osm    !   OSM vertical diffusion                               (ln_zdfosm =T)
1087!-----------------------------------------------------------------------
1088   ln_use_osm_la = .false.      !  Use namelist  rn_osm_la
1089   rn_osm_la     = 0.3         !  Turbulent Langmuir number
1090   rn_osm_dstokes     = 5.     !  Depth scale of Stokes drift (m)
1091   nn_ave = 0                  ! choice of horizontal averaging on avt, avmu, avmv
1092   ln_dia_osm = .true.         ! output OSMOSIS-OBL variables
1093   rn_osm_hbl0 = 10.           ! initial hbl value
1094   ln_kpprimix = .true.        ! Use KPP-style Ri# mixing below BL
1095   rn_riinfty  = 0.7           ! Highest local Ri_g permitting shear instability
1096   rn_difri  =  0.005          ! max Ri# diffusivity at Ri_g = 0 (m^2/s)
1097   ln_convmix  = .true.        ! Use convective instability mixing below BL
1098   rn_difconv = 1.             ! diffusivity when unstable below BL  (m2/s)
1099   nn_osm_wave = 0             ! Method used to calculate Stokes drift
1100      !                        !  = 2: Use ECMWF wave fields
1101      !                        !  = 1: Pierson Moskowitz wave spectrum
1102      !                        !  = 0: Constant La# = 0.3
1103/
1104!-----------------------------------------------------------------------
1105&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T)
1106!-----------------------------------------------------------------------
1107   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
1108   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency
1109   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F)
1110/
1111
1112!!======================================================================
1113!!                  ***  Diagnostics namelists  ***                   !!
1114!!                                                                    !!
1115!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
1116!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
1117!!   namhsb       Heat and salt budgets                                 (default: OFF)
1118!!   namdiu       Cool skin and warm layer models                       (default: OFF)
1119!!   namdiu       Cool skin and warm layer models                       (default: OFF)
1120!!   namflo       float parameters                                      (default: OFF)
1121!!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF)
1122!!   nam_diadct   transports through some sections                      (default: OFF)
1123!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
1124!!   nam_dia25h   25h Mean Output                                       (default: OFF)
1125!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1126!!======================================================================
1127!
1128!-----------------------------------------------------------------------
1129&namtrd        !   trend diagnostics                                    (default: OFF)
1130!-----------------------------------------------------------------------
1131   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1132   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1133   ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1134   ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet)
1135   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1136   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1137   ln_tra_trd  = .false.   ! (T) 3D tracer trend output
1138   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1139   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1140/
1141!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1142!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1143!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1144!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1145!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1146!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1147!!gm
1148!-----------------------------------------------------------------------
1149&namptr        !   Poleward Transport Diagnostic                        (default: OFF)
1150!-----------------------------------------------------------------------
1151   ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F)
1152   ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not
1153/
1154!-----------------------------------------------------------------------
1155&namhsb        !  Heat and salt budgets                                 (default: OFF)
1156!-----------------------------------------------------------------------
1157   ln_diahsb   = .false.   !  output the heat and salt budgets (T) or not (F)
1158/
1159!-----------------------------------------------------------------------
1160&namdiu        !   Cool skin and warm layer models                      (default: OFF)
1161!-----------------------------------------------------------------------
1162   ln_diurnal      = .false.   !
1163   ln_diurnal_only = .false.   !
1164/
1165!-----------------------------------------------------------------------
1166&namflo        !   float parameters                                     (default: OFF)
1167!-----------------------------------------------------------------------
1168   ln_floats   = .false.      ! activate floats or not
1169      jpnfl       = 1         !    total number of floats during the run
1170      jpnnewflo   = 0         !    number of floats for the restart
1171      ln_rstflo   = .false.   !    float restart (T) or not (F)
1172      nn_writefl  =      75   !    frequency of writing in float output file
1173      nn_stockfl  =    5475   !    frequency of creation of the float restart file
1174      ln_argo     = .false.   !    Argo type floats (stay at the surface each 10 days)
1175      ln_flork4   = .false.   !    trajectories computed with a 4th order Runge-Kutta (T)
1176      !                       !    or computed with Blanke' scheme (F)
1177      ln_ariane   = .true.    !    Input with Ariane tool convention(T)
1178      ln_flo_ascii= .true.    !    Output with Ariane tool netcdf convention(F) or ascii file (T)
1179/
1180!-----------------------------------------------------------------------
1181&nam_diaharm   !   Harmonic analysis of tidal constituents              (default: OFF)
1182!-----------------------------------------------------------------------
1183    ln_diaharm = .false.   ! Choose tidal harmonic output or not
1184       nit000_han = 1      !    First time step used for harmonic analysis
1185       nitend_han = 75     !    Last time step used for harmonic analysis
1186       nstep_han  = 15     !    Time step frequency for harmonic analysis
1187       tname(1)   = 'M2'   !    Name of tidal constituents
1188       tname(2)   = 'K1'   !              ---
1189/
1190!-----------------------------------------------------------------------
1191&nam_diadct    !   transports through some sections                     (default: OFF)
1192!-----------------------------------------------------------------------
1193    ln_diadct  = .false.   ! Calculate transport thru sections or not
1194       nn_dct     = 15     !  time step frequency for transports computing
1195       nn_dctwri  = 15     !  time step frequency for transports writing
1196       nn_secdebug = 112   !      0 : no section to debug
1197       !                   !     -1 : debug all section
1198       !                   !  0 < n : debug section number n
1199/
1200!-----------------------------------------------------------------------
1201&nam_diatmb    !  Top Middle Bottom Output                              (default: OFF)
1202!-----------------------------------------------------------------------
1203   ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not
1204/
1205!-----------------------------------------------------------------------
1206&nam_dia25h    !  25h Mean Output                                       (default: OFF)
1207!-----------------------------------------------------------------------
1208   ln_dia25h   = .false.   ! Choose 25h mean output or not
1209/
1210!-----------------------------------------------------------------------
1211&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
1212!-----------------------------------------------------------------------
1213   nn_nchunks_i =   4       !  number of chunks in i-dimension
1214   nn_nchunks_j =   4       !  number of chunks in j-dimension
1215   nn_nchunks_k =   31      !  number of chunks in k-dimension
1216   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1217   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs
1218   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1219   !                       !  (F) ignore chunking information and produce netcdf3-compatible files
1220/
1221
1222!!======================================================================
1223!!               ***  Observation & Assimilation  ***                 !!
1224!!                                                                    !!
1225!!   namobs       observation and model comparison                      (default: OFF)
1226!!   nam_asminc   assimilation increments                               ('key_asminc')
1227!!======================================================================
1228!
1229!-----------------------------------------------------------------------
1230&namobs        !  observation usage switch                              (default: OFF)
1231!-----------------------------------------------------------------------
1232   ln_diaobs   = .false.             ! Logical switch for the observation operator
1233   !
1234   ln_t3d      = .false.             ! Logical switch for T profile observations
1235   ln_s3d      = .false.             ! Logical switch for S profile observations
1236   ln_sla      = .false.             ! Logical switch for SLA observations
1237   ln_sst      = .false.             ! Logical switch for SST observations
1238   ln_sss      = .false.             ! Logical swithc for SSS observations
1239   ln_sic      = .false.             ! Logical switch for Sea Ice observations
1240   ln_vel3d    = .false.             ! Logical switch for velocity observations
1241   ln_altbias  = .false.             ! Logical switch for altimeter bias correction
1242   ln_sstbias  = .false.             ! Logical switch for SST bias correction
1243   ln_nea      = .false.             ! Logical switch for rejection of observations near land
1244   ln_grid_global = .true.           ! Logical switch for global distribution of observations
1245   ln_grid_search_lookup = .false.   ! Logical switch for obs grid search w/lookup table
1246   ln_ignmis   = .true.              ! Logical switch for ignoring missing files
1247   ln_s_at_t   = .false.             ! Logical switch for computing model S at T obs if not there
1248   ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs
1249   ln_bound_reject  = .false.        ! Logical to remove obs near boundaries in LAMs.
1250   ln_sla_fp_indegs = .true.         ! Logical for SLA: T=> averaging footprint is in degrees, F=> in metres
1251   ln_sst_fp_indegs = .true.         ! Logical for SST: T=> averaging footprint is in degrees, F=> in metres
1252   ln_sss_fp_indegs = .true.         ! Logical for SSS: T=> averaging footprint is in degrees, F=> in metres
1253   ln_sic_fp_indegs = .true.         ! Logical for SIC: T=> averaging footprint is in degrees, F=> in metres
1254! All of the *files* variables below are arrays. Use namelist_cfg to add more files
1255   cn_profbfiles = 'profiles_01.nc'  ! Profile feedback input observation file names
1256   cn_slafbfiles = 'sla_01.nc'       ! SLA feedback input observation file names
1257   cn_sstfbfiles = 'sst_01.nc'       ! SST feedback input observation file names
1258   cn_sssfbfiles = 'sss_01.nc'       ! SSS feedback input observation file names
1259   cn_sicfbfiles = 'sic_01.nc'       ! SIC feedback input observation file names
1260   cn_velfbfiles = 'vel_01.nc'       ! Velocity feedback input observation file names
1261   cn_altbiasfile = 'altbias.nc'     ! Altimeter bias input file name
1262   cn_sstbiasfiles = 'sstbias.nc'    ! SST bias input file name
1263   cn_gridsearchfile ='gridsearch.nc' ! Grid search file name
1264   rn_gridsearchres = 0.5            ! Grid search resolution
1265   rn_mdtcorr  = 1.61                ! MDT  correction
1266   rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction
1267   rn_dobsini  = 00010101.000000     ! Initial date in window YYYYMMDD.HHMMSS
1268   rn_dobsend  = 00010102.000000     ! Final date in window YYYYMMDD.HHMMSS
1269   rn_sla_avglamscl = 0.             ! E/W diameter of SLA observation footprint (metres/degrees)
1270   rn_sla_avgphiscl = 0.             ! N/S diameter of SLA observation footprint (metres/degrees)
1271   rn_sst_avglamscl = 0.             ! E/W diameter of SST observation footprint (metres/degrees)
1272   rn_sst_avgphiscl = 0.             ! N/S diameter of SST observation footprint (metres/degrees)
1273   rn_sss_avglamscl = 0.             ! E/W diameter of SSS observation footprint (metres/degrees)
1274   rn_sss_avgphiscl = 0.             ! N/S diameter of SSS observation footprint (metres/degrees)
1275   rn_sic_avglamscl = 0.             ! E/W diameter of SIC observation footprint (metres/degrees)
1276   rn_sic_avgphiscl = 0.             ! N/S diameter of SIC observation footprint (metres/degrees)
1277   nn_1dint = 0                      ! Type of vertical interpolation method
1278   nn_2dint = 0                      ! Default horizontal interpolation method
1279   nn_2dint_sla = 0                  ! Horizontal interpolation method for SLA
1280   nn_2dint_sst = 0                  ! Horizontal interpolation method for SST
1281   nn_2dint_sss = 0                  ! Horizontal interpolation method for SSS
1282   nn_2dint_sic = 0                  ! Horizontal interpolation method for SIC
1283   nn_msshc     = 0                  ! MSSH correction scheme
1284   nn_profdavtypes = -1              ! Profile daily average types - array
1285/
1286!-----------------------------------------------------------------------
1287&nam_asminc    !   assimilation increments                              ('key_asminc')
1288!-----------------------------------------------------------------------
1289    ln_bkgwri  = .false.   !  Logical switch for writing out background state
1290    ln_trainc  = .false.   !  Logical switch for applying tracer increments
1291    ln_dyninc  = .false.   !  Logical switch for applying velocity increments
1292    ln_sshinc  = .false.   !  Logical switch for applying SSH increments
1293    ln_asmdin  = .false.   !  Logical switch for Direct Initialization (DI)
1294    ln_asmiau  = .false.   !  Logical switch for Incremental Analysis Updating (IAU)
1295    nitbkg     = 0         !  Timestep of background in [0,nitend-nit000-1]
1296    nitdin     = 0         !  Timestep of background for DI in [0,nitend-nit000-1]
1297    nitiaustr  = 1         !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1298    nitiaufin  = 15        !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1299    niaufn     = 0         !  Type of IAU weighting function
1300    ln_salfix  = .false.   !  Logical switch for ensuring that the sa > salfixmin
1301    salfixmin  = -9999     !  Minimum salinity after applying the increments
1302    nn_divdmp  = 0         !  Number of iterations of divergence damping operator
1303/
1304
1305!!======================================================================
1306!!                  ***  Miscellaneous namelists  ***                 !!
1307!!                                                                    !!
1308!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
1309!!   namctl            Control prints                                   (default: OFF)
1310!!   namsto            Stochastic parametrization of EOS                (default: OFF)
1311!!======================================================================
1312!
1313!-----------------------------------------------------------------------
1314&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
1315!-----------------------------------------------------------------------
1316   ln_listonly =  .false.  !  do nothing else than listing the best domain decompositions (with land domains suppression)
1317   !                       !  if T: the largest number of cores tested is defined by max(mppsize, jpni*jpnj)
1318   ln_nnogather =  .true.  !  activate code to avoid mpi_allgather use at the northfold
1319   jpni        =   0       !  number of processors following i (set automatically if < 1), see also ln_listonly = T
1320   jpnj        =   0       !  number of processors following j (set automatically if < 1), see also ln_listonly = T
1321/
1322!-----------------------------------------------------------------------
1323&namctl        !   Control prints                                       (default: OFF)
1324!-----------------------------------------------------------------------
1325   ln_ctl = .FALSE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T
1326     sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following
1327       sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings.
1328       sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure
1329       sn_cfctl%l_oceout  = .FALSE. ! that  all areas report.
1330       sn_cfctl%l_layout  = .FALSE. !
1331       sn_cfctl%l_mppout  = .FALSE. !
1332       sn_cfctl%l_mpptop  = .FALSE. !
1333       sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0]
1334       sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000]
1335       sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1]
1336       sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info
1337   nn_print    =    0      !  level of print (0 no extra print)
1338   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
1339   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
1340   nn_jctls    =    0      !  start j indice of control               over a subdomain)
1341   nn_jctle    =    0      !  end   j indice of control
1342   nn_isplt    =    1      !  number of processors in i-direction
1343   nn_jsplt    =    1      !  number of processors in j-direction
1344   ln_timing   = .false.   !  timing by routine write out in timing.output file
1345   ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii
1346/
1347!-----------------------------------------------------------------------
1348&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
1349!-----------------------------------------------------------------------
1350   ln_sto_eos  = .false.   ! stochastic equation of state
1351   nn_sto_eos  = 1         ! number of independent random walks
1352   rn_eos_stdxy = 1.4       ! random walk horz. standard deviation (in grid points)
1353   rn_eos_stdz = 0.7       ! random walk vert. standard deviation (in grid points)
1354   rn_eos_tcor = 1440.     ! random walk time correlation (in timesteps)
1355   nn_eos_ord  = 1         ! order of autoregressive processes
1356   nn_eos_flt  = 0         ! passes of Laplacian filter
1357   rn_eos_lim  = 2.0       ! limitation factor (default = 3.0)
1358   ln_rststo   = .false.   ! start from mean parameter (F) or from restart file (T)
1359   ln_rstseed  = .true.    ! read seed of RNG from restart file
1360   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
1361   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)
1362/
Note: See TracBrowser for help on using the repository browser.