New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
namelist_ref in branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED – NEMO

source: branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED/namelist_ref @ 9169

Last change on this file since 9169 was 9169, checked in by gm, 6 years ago

dev_merge_2017: all SRC: finalize the removal of useless warning when reading namelist_cfg + remove all nn_closea + nn_msh replaced by a logical

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