source: branches/UKMO/dev_r5518_GO6_package/NEMOGCM/CONFIG/SHARED/namelist_ref @ 11101

Last change on this file since 11101 was 11101, checked in by frrh, 18 months ago

Merge changes from Met Office GMED ticket 450 to reduce unnecessary
text output from NEMO.
This output, which is typically not switchable, is rarely of interest
in normal (non-debugging) runs and simply redunantley consumes extra
file space.
Further, the presence of this text output has been shown to
significantly degrade performance of models which are run during
Met Office HPC RAID (disk) checks.
The new code introduces switches which are configurable via the
changes made in the associated Met Office MOCI ticket 399.

File size: 100.9 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  :  1 - run manager      (namrun)
3!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd)
4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas
5!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,
6!!                                    namsbc_apr, namsbc_ssr, namsbc_alb)
7!!              4 - lateral boundary (namlbc, namcla, namagrif, nambdy, nambdy_tide)
8!!              5 - bottom  boundary (nambfr, nambbc, nambbl)
9!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp)
10!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
11!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx, namzdf_tmx_new)
12!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto)
13!!             10 - miscellaneous    (namsol, nammpp, namctl)
14!!             11 - Obs & Assim      (namobs, nam_asminc)
15!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
16
17!!======================================================================
18!!                   ***  Run management namelists  ***
19!!======================================================================
20!!   namrun       parameters of the run
21!!======================================================================
22!
23!-----------------------------------------------------------------------
24&namrun        !   parameters of the run
25!-----------------------------------------------------------------------
26   nn_no       =       0   !  job number (no more used...)
27   cn_exp      =  "ORCA2"  !  experience name
28   nn_it000    =       1   !  first time step
29   nn_itend    =    5475   !  last  time step (std 5475)
30   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
31   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
32   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
33   nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=T
34   nn_rstctl   =       0   !  restart control ==> activated only if ln_rstart=T
35                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
36                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
37                           !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart
38   cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input)
39   cn_ocerst_indir = "."       !  directory from which to read input ocean restarts
40   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output)
41   cn_ocerst_outdir = "."      !  directory in which to write output ocean restarts
42   nn_istate   =       0   !  output the initial state (1) or not (0)
43   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
44   nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1)
45   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
46   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000)
47   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T)
48   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%)
49   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard
50   ln_clobber  = .false.   !  clobber (overwrite) an existing file
51   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
52/
53!
54!!======================================================================
55!!                      ***  Domain namelists  ***
56!!======================================================================
57!!   namcfg       parameters of the configuration
58!!   namzgr       vertical coordinate
59!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate
60!!   namdom       space and time domain (bathymetry, mesh, timestep)
61!!   namtsd       data: temperature & salinity
62!!======================================================================
63!
64!-----------------------------------------------------------------------
65&namcfg     !   parameters of the configuration
66!-----------------------------------------------------------------------
67   cp_cfg      =  "default"            !  name of the configuration
68   cp_cfz      =  "no zoom"            !  name of the zoom of configuration
69   jp_cfg      =       0               !  resolution of the configuration
70   jpidta      =      10               !  1st lateral dimension ( >= jpi )
71   jpjdta      =      12               !  2nd    "         "    ( >= jpj )
72   jpkdta      =      31               !  number of levels      ( >= jpk )
73   jpiglo      =      10               !  1st dimension of global domain --> i =jpidta
74   jpjglo      =      12               !  2nd    -                  -    --> j =jpjdta
75   jpizoom     =       1               !  left bottom (i,j) indices of the zoom
76   jpjzoom     =       1               !  in data domain indices
77   jperio      =       0               !  lateral cond. type (between 0 and 6)
78                                       !  = 0 closed                 ;   = 1 cyclic East-West
79                                       !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot
80                                       !  = 4 cyclic East-West AND North fold T-point pivot
81                                       !  = 5 North fold F-point pivot
82                                       !  = 6 cyclic East-West AND North fold F-point pivot
83   ln_use_jattr = .false.              !  use (T) the file attribute: open_ocean_jstart, if present
84                                       !  in netcdf input files, as the start j-row for reading
85/
86!-----------------------------------------------------------------------
87&namzgr        !   vertical coordinate
88!-----------------------------------------------------------------------
89   ln_zco      = .false.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined)
90   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F)
91   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F)
92   ln_isfcav   = .false.   !  ice shelf cavity               (T/F)
93/
94!-----------------------------------------------------------------------
95&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate
96!-----------------------------------------------------------------------
97   ln_s_sh94   = .true.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)|
98   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied
99   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch
100                           !  stretching coefficients for all functions
101   rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m)
102   rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m)
103   rn_hc       =  150.0    !  critical depth for transition to stretched coordinates
104                        !!!!!!!  Envelop bathymetry
105   rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1)
106                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.)
107   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20)
108   rn_bb       =    0.8    !  stretching with SH94 s-sigma
109                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.)
110   rn_alpha    =    4.4    !  stretching with SF12 s-sigma
111   rn_efold    =    0.0    !  efold length scale for transition to stretched coord
112   rn_zs       =    1.0    !  depth of surface grid box
113                           !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b
114   rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb
115   rn_zb_b     =   -0.2    !  offset for calculating Zb
116                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above]
117   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)
118/
119!-----------------------------------------------------------------------
120&namdom        !   space and time domain (bathymetry, mesh, timestep)
121!-----------------------------------------------------------------------
122   nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file
123   rn_bathy    =    0.     !  value of the bathymetry. if (=0) bottom flat at jpkm1
124   nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA)
125   nn_msh      =    1      !  create (=1) a mesh file or not (=0)
126   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0)
127   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of
128   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1
129                           !
130   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0)
131   rn_atfp     =    0.1    !  asselin time filter parameter
132   nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k)
133                                 !                          =0, not used, rdt = rdttra
134   rn_rdtmin   = 28800.          !  minimum time step on tracers (used if nn_acc=1)
135   rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nn_acc=1)
136   rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1)
137   ln_crs      = .false.      !  Logical switch for coarsening module
138   jphgr_msh   =       0               !  type of horizontal mesh
139                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc
140                                       !  = 1 geographical mesh on the sphere with regular grid-spacing
141                                       !  = 2 f-plane with regular grid-spacing
142                                       !  = 3 beta-plane with regular grid-spacing
143                                       !  = 4 Mercator grid with T/U point at the equator
144   ppglam0     =       0.0             !  longitude of first raw and column T-point (jphgr_msh = 1)
145   ppgphi0     =     -35.0             ! latitude  of first raw and column T-point (jphgr_msh = 1)
146   ppe1_deg    =       1.0             !  zonal      grid-spacing (degrees)
147   ppe2_deg    =       0.5             !  meridional grid-spacing (degrees)
148   ppe1_m      =    5000.0             !  zonal      grid-spacing (degrees)
149   ppe2_m      =    5000.0             !  meridional grid-spacing (degrees)
150   ppsur       =    -4762.96143546300  !  ORCA r4, r2 and r05 coefficients
151   ppa0        =      255.58049070440  ! (default coefficients)
152   ppa1        =      245.58132232490  !
153   ppkth       =       21.43336197938  !
154   ppacr       =        3.0            !
155   ppdzmin     =       10.             !  Minimum vertical spacing
156   pphmax      =     5000.             !  Maximum depth
157   ldbletanh   =    .TRUE.             !  Use/do not use double tanf function for vertical coordinates
158   ppa2        =      100.760928500000 !  Double tanh function parameters
159   ppkth2      =       48.029893720000 !
160   ppacr2      =       13.000000000000 !
161/
162!-----------------------------------------------------------------------
163&namsplit      !   time splitting parameters                            ("key_dynspg_ts")
164!-----------------------------------------------------------------------
165   ln_bt_fw      =    .TRUE.           !  Forward integration of barotropic equations
166   ln_bt_av      =    .TRUE.           !  Time filtering of barotropic variables
167   ln_bt_nn_auto =    .TRUE.           !  Set nn_baro automatically to be just below
168                                       !  a user defined maximum courant number (rn_bt_cmax)
169   nn_baro       =    30               !  Number of iterations of barotropic mode
170                                       !  during rn_rdt seconds. Only used if ln_bt_nn_auto=F
171   rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T
172   nn_bt_flt     =    1                !  Time filter choice
173                                       !  = 0 None
174                                       !  = 1 Boxcar over   nn_baro barotropic steps
175                                       !  = 2 Boxcar over 2*nn_baro     "        "
176/
177!-----------------------------------------------------------------------
178&namcrs        !   Grid coarsening for dynamics output and/or
179               !   passive tracer coarsened online simulations
180!-----------------------------------------------------------------------
181   nn_factx    = 3         !  Reduction factor of x-direction
182   nn_facty    = 3         !  Reduction factor of y-direction
183   nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT
184                           !  0, coarse grid is binned with preferential treatment of the north fold
185                           !  1, coarse grid is binned with centering at the equator
186                           !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.
187   nn_msh_crs  = 1         !  create (=1) a mesh file or not (=0)
188   nn_crs_kz   = 0         ! 0, MEAN of volume boxes
189                           ! 1, MAX of boxes
190                           ! 2, MIN of boxes
191   ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F )
192/
193!-----------------------------------------------------------------------
194&namc1d        !   1D configuration options                             ("key_c1d")
195!-----------------------------------------------------------------------
196   rn_lat1d    =      50   !  Column latitude (default at PAPA station)
197   rn_lon1d    =    -145   !  Column longitude (default at PAPA station)
198   ln_c1d_locpt=  .true.   ! Localization of 1D config in a grid (T) or independant point (F)
199/
200!-----------------------------------------------------------------------
201&namtsd    !   data : Temperature  & Salinity
202!-----------------------------------------------------------------------
203!-----------------------------------------------------------------------
204!          !  file name                            ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
205!          !                                       !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
206   sn_tem  = 'data_1m_potential_temperature_nomask',         -1        ,'votemper' ,    .true.    , .true. , 'yearly'   , ''       ,   ''    ,    ''
207   sn_sal  = 'data_1m_salinity_nomask'             ,         -1        ,'vosaline' ,    .true.    , .true. , 'yearly'   , ''       ,   ''    ,    ''
208   !
209   cn_dir        = './'     !  root directory for the location of the runoff files
210   ln_tsd_init   = .true.   !  Initialisation of ocean T & S with T &S input data (T) or not (F)
211   ln_tsd_tradmp = .true.   !  damping of ocean T & S toward T &S input data (T) or not (F)
212/
213!!======================================================================
214!!            ***  Surface Boundary Condition namelists  ***
215!!======================================================================
216!!   namsbc          surface boundary condition
217!!   namsbc_ana      analytical         formulation
218!!   namsbc_flx      flux               formulation
219!!   namsbc_clio     CLIO bulk formulae formulation
220!!   namsbc_core     CORE bulk formulae formulation
221!!   namsbc_mfs      MFS  bulk formulae formulation
222!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3")
223!!   namsbc_sas      StAndalone Surface module
224!!   namtra_qsr      penetrative solar radiation
225!!   namsbc_rnf      river runoffs
226!!   namsbc_isf      ice shelf melting/freezing
227!!   namsbc_apr      Atmospheric Pressure
228!!   namsbc_ssr      sea surface restoring term (for T and/or S)
229!!   namsbc_alb      albedo parameters
230!!======================================================================
231!
232!-----------------------------------------------------------------------
233&namsbc        !   Surface Boundary Condition (surface module)
234!-----------------------------------------------------------------------
235   nn_fsbc     = 5         !  frequency of surface boundary condition computation
236                           !     (also = the frequency of sea-ice model call)
237   ln_ana      = .false.   !  analytical formulation                    (T => fill namsbc_ana )
238   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx )
239   ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio)
240   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core)
241   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs )
242   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 )
243   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 )
244   nn_components = 0       !  configuration of the opa-sas OASIS coupling
245                           !  =0 no opa-sas OASIS coupling: default single executable configuration
246                           !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component
247                           !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component
248   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr )
249   nn_ice      = 2         !  =0 no ice boundary condition   ,
250                           !  =1 use observed ice-cover      ,
251                           !  =2 ice-model used                         ("key_lim3" or "key_lim2")
252   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect)
253                           !  =1 levitating ice with mass and salt exchange but no presure effect
254                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure)
255   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave
256   ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf)
257   nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf)
258                           !  0 =no isf                  1 = presence of ISF
259                           !  2 = bg03 parametrisation   3 = rnf file for isf
260                           !  4 = ISF fwf specified
261                           !  option 1 and 4 need ln_isfcav = .true. (domzgr)
262   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
263   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked
264                           !     =1 global mean of e-p-r set to zero at each time step
265                           !     =2 annual global mean of e-p-r set to zero
266   ln_wave = .false.       !  Activate coupling with wave (either Stokes Drift or Drag coefficient, or both)  (T => fill namsbc_wave)
267   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => fill namsbc_wave)
268   ln_sdw  = .false.       !  Computation of 3D stokes drift                (T => fill namsbc_wave)
269   nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
270                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
271   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used)
272                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled
273                           !  = 0  Average per-category fluxes (forced and coupled mode)
274                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled
275                           !  = 2  Redistribute a single flux over categories (coupled mode only)
276/
277!-----------------------------------------------------------------------
278&namsbc_ana    !   analytical surface boundary condition
279!-----------------------------------------------------------------------
280   nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps
281   rn_utau0    =   0.5     !  uniform value for the i-stress
282   rn_vtau0    =   0.e0    !  uniform value for the j-stress
283   rn_qns0     =   0.e0    !  uniform value for the total heat flux
284   rn_qsr0     =   0.e0    !  uniform value for the solar radiation
285   rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P)
286/
287!-----------------------------------------------------------------------
288&namsbc_flx    !   surface boundary condition : flux formulation
289!-----------------------------------------------------------------------
290!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
291!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
292   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
293   sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
294   sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
295   sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
296   sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
297
298   cn_dir      = './'      !  root directory for the location of the flux files
299/
300!-----------------------------------------------------------------------
301&namsbc_clio   !   namsbc_clio  CLIO bulk formulae
302!-----------------------------------------------------------------------
303!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
304!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
305   sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
306   sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
307   sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
308   sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
309   sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
310   sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
311   sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
312
313   cn_dir      = './'      !  root directory for the location of the bulk files are
314/
315!-----------------------------------------------------------------------
316&namsbc_core   !   namsbc_core  CORE bulk formulae
317!-----------------------------------------------------------------------
318!              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask !
319!              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      !
320   sn_wndi     = 'u_10.15JUNE2009_fill'        ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , ''
321   sn_wndj     = 'v_10.15JUNE2009_fill'        ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , ''
322   sn_qsr      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
323   sn_qlw      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
324   sn_tair     = 't_10.15JUNE2009_fill'        ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
325   sn_humi     = 'q_10.15JUNE2009_fill'        ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
326   sn_prec     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
327   sn_snow     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
328   sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
329
330   cn_dir      = './'      !  root directory for the location of the bulk files
331   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data
332   rn_zqt      = 10.        !  Air temperature and humidity reference height (m)
333   rn_zu       = 10.        !  Wind vector reference height (m)
334   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow)
335   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.)
336   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity
337                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)
338/
339!-----------------------------------------------------------------------
340&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae
341!-----------------------------------------------------------------------
342!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask !
343!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      !
344   sn_wndi     =   'ecmwf'   ,        6          , 'u10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
345   sn_wndj     =   'ecmwf'   ,        6          , 'v10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
346   sn_clc      =   'ecmwf'   ,        6          , 'clc'     ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , ''
347   sn_msl      =   'ecmwf'   ,        6          , 'msl'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
348   sn_tair     =   'ecmwf'   ,        6          , 't2'      ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
349   sn_rhm      =   'ecmwf'   ,        6          , 'rh'      ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , ''
350   sn_prec     =   'ecmwf'   ,        6          , 'precip'  ,    .true.    , .true.  , 'daily'  ,'bicubic.nc' , ''       , ''
351
352   cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files
353/
354!-----------------------------------------------------------------------
355&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
356!-----------------------------------------------------------------------
357!                    !     description       !  multiple  !    vector   !      vector          ! vector !
358!                    !                       ! categories !  reference  !    orientation       ! grids  !
359! send
360   sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   ''
361   sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
362   sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   ''
363   sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T'
364   sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
365   sn_snd_bio_co2 =      'medusa'               ,    'no'    ,     ''      ,         ''           ,   ''
366   sn_snd_bio_dms =      'medusa'               ,    'no'    ,     ''      ,         ''           ,   ''
367   sn_snd_bio_chloro =   'medusa'               ,    'no'    ,     ''      ,         ''           ,   ''
368   sn_snd_cond   =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
369   sn_snd_mpnd   =       'ice only'             ,    'no'    ,     ''      ,         ''           ,   ''
370   sn_snd_sstfrz =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
371   sn_snd_thick1 =       'ice and snow'         ,    'no'    ,     ''      ,         ''           ,   ''
372! receive
373   sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
374   sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
375   sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V'
376   sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
377   sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
378   sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
379   sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   ''
380   sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
381   sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
382   sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
383   sn_rcv_antm   =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
384   sn_rcv_grnm   =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
385   sn_rcv_iceflx =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
386   sn_rcv_ts_ice =       'ice'                  ,    'no'    ,     ''      ,         ''          ,   ''
387   sn_rcv_atm_dust =     'medusa'               ,    'no'    ,     ''      ,         ''          ,   ''
388   sn_rcv_atm_pco2 =     'medusa'               ,    'no'    ,     ''      ,         ''          ,   ''
389
390!
391   nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data
392   ln_usecplmask = .false.   !  use a coupling mask file to merge data received from several models
393                             !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
394   nn_coupled_iceshelf_fluxes = 0 ! =0 : total freshwater input from iceberg calving and ice shelf basal melting
395                                  ! taken from climatologies used (no action in coupling routines).
396                                  ! =1 :  use rate of change of mass of Greenland and Antarctic icesheets to set the
397                                  ! combined magnitude of the iceberg calving and iceshelf melting freshwater fluxes.
398                                  ! =2 :  specify constant freshwater inputs in this namelist to set the combined
399                                  ! magnitude of iceberg calving and iceshelf melting freshwater fluxes.
400   ln_iceshelf_init_atmos     = .true.  ! If true force ocean to initialise icesheet masses from atmospheric values rather than
401                                        ! from values in ocean restart file.
402   rn_greenland_total_fw_flux   = 0.0  ! Constant total rate of freshwater input (kg/s) for Greenland (if nn_coupled_iceshelf_fluxes=2)
403   rn_greenland_calving_fraction = 0.5  ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting.
404   rn_antarctica_total_fw_flux  = 0.0  ! Constant total rate of freshwater input (kg/s) for Antarctica (if nn_coupled_iceshelf_fluxes=2)
405   rn_antarctica_calving_fraction = 0.5 ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting.
406   rn_iceshelf_fluxes_tolerance = 1e-6  ! Fractional threshold for detecting differences in icesheet masses (must be positive definite).
407/
408!-----------------------------------------------------------------------
409&namsbc_sas    !   analytical surface boundary condition
410!-----------------------------------------------------------------------
411!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
412!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
413   sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
414   sn_vsp      = 'sas_grid_V' ,    120           , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
415   sn_tem      = 'sas_grid_T' ,    120           , 'sosstsst' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
416   sn_sal      = 'sas_grid_T' ,    120           , 'sosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
417   sn_ssh      = 'sas_grid_T' ,    120           , 'sossheig' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
418   sn_e3t      = 'sas_grid_T' ,    120           , 'e3t_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
419   sn_frq      = 'sas_grid_T' ,    120           , 'frq_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
420
421   ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field
422   ln_read_frq = .false.    !  specify whether we must read frq or not
423   cn_dir      = './'      !  root directory for the location of the bulk files are
424/
425!-----------------------------------------------------------------------
426&namtra_qsr    !   penetrative solar radiation
427!-----------------------------------------------------------------------
428!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
429!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
430   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , ''
431
432   cn_dir      = './'      !  root directory for the location of the runoff files
433   ln_traqsr   = .true.    !  Light penetration (T) or not (F)
434   ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration
435   ln_qsr_2bd  = .false.   !  2 bands              light penetration
436   ln_qsr_bio  = .false.   !  bio-model light penetration
437   nn_chldta   =      1    !  RGB : 2D Chl data (=1), 3D Chl data (=2) or cst value (=0)
438   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
439   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
440   rn_si1      =   23.0    !  2 bands: longest depth of extinction
441   ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3
442/
443!-----------------------------------------------------------------------
444&namsbc_rnf    !   runoffs namelist surface boundary condition
445!-----------------------------------------------------------------------
446!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
447!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
448   sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
449   sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
450   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
451   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
452   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , ''
453
454   cn_dir       = './'      !  root directory for the location of the runoff files
455   ln_rnf_mouth = .true.    !  specific treatment at rivers mouths
456   rn_hrnf      =  15.e0    !  depth over which enhanced vertical mixing is used
457   rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s]
458   rn_rfact     =   1.e0    !  multiplicative factor for runoff
459   ln_rnf_depth = .false.   !  read in depth information for runoff
460   ln_rnf_tem   = .false.   !  read in temperature information for runoff
461   ln_rnf_sal   = .false.   !  read in salinity information for runoff
462   ln_rnf_depth_ini = .false.  ! compute depth at initialisation from runoff file
463   rn_rnf_max   = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
464   rn_dep_max   = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
465   nn_rnf_depth_file = 0    !  create (=1) a runoff depth file or not (=0)
466/
467!-----------------------------------------------------------------------
468&namsbc_isf    !  Top boundary layer (ISF)
469!-----------------------------------------------------------------------
470!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
471!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
472! nn_isf == 4
473   sn_qisf      = 'rnfisf' ,         -12      ,'sohflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
474   sn_fwfisf    = 'rnfisf' ,         -12      ,'sowflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
475! nn_isf == 3
476   sn_rnfisf    = 'runoffs' ,         -12      ,'sofwfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
477! nn_isf == 2 and 3
478   sn_depmax_isf = 'runoffs' ,       -12        ,'sozisfmax' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
479   sn_depmin_isf = 'runoffs' ,       -12        ,'sozisfmin' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
480! nn_isf == 2
481   sn_Leff_isf = 'rnfisf' ,       0          ,'Leff'         ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
482! for all case
483   ln_divisf   = .true.  ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?)
484! only for nn_isf = 1 or 2
485   rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula
486   rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula
487! only for nn_isf = 1
488   nn_isfblk   =  1       ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??)
489   rn_hisf_tbl =  30.      ! thickness of the top boundary layer           (Losh et al. 2008)
490                          ! 0 => thickness of the tbl = thickness of the first wet cell
491   ln_conserve = .true.   ! conservative case (take into account meltwater advection)
492   nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s)
493                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010)
494                          !     if you want to keep the cd as in global config, adjust rn_gammat0 to compensate
495                          ! 2 = velocity and stability dependent Gamma    Holland et al. 1999
496/
497!-----------------------------------------------------------------------
498&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk
499!-----------------------------------------------------------------------
500!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
501!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
502   sn_apr      = 'patm'      ,         -1        ,'somslpre',    .true.     , .true. , 'yearly'  ,  ''      ,   ''     , ''
503
504   cn_dir      = './'       !  root directory for the location of the bulk files
505   rn_pref     = 101000.    !  reference atmospheric pressure   [N/m2]/
506   ln_ref_apr  = .false.    !  ref. pressure: global mean Patm (T) or a constant (F)
507   ln_apr_obc  = .false.    !  inverse barometer added to OBC ssh data
508/
509!-----------------------------------------------------------------------
510&namsbc_ssr    !   surface boundary condition : sea surface restoring
511!-----------------------------------------------------------------------
512!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
513!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
514   sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , ''       , ''
515   sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , ''       , ''
516
517   cn_dir      = './'      !  root directory for the location of the runoff files
518   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
519   nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2)
520                           !  or to SSS only (=1) or no damping term (=0)
521   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K]
522   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
523   ln_sssr_bnd =   .true.  !  flag to bound erp term (associated with nn_sssr=2)
524   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
525/
526!-----------------------------------------------------------------------
527&namsbc_alb    !   albedo parameters
528!-----------------------------------------------------------------------
529   nn_ice_alb  =    0   !  parameterization of ice/snow albedo
530                        !     0: Shine & Henderson-Sellers (JGR 1985)
531                        !     1: "home made" based on Brandt et al. (J. Climate 2005)
532                        !                         and Grenfell & Perovich (JGR 2004)
533   rn_albice   =  0.53  !  albedo of bare puddled ice (values from 0.49 to 0.58)
534                        !     0.53 (default) => if nn_ice_alb=0
535                        !     0.50 (default) => if nn_ice_alb=1
536/
537!-----------------------------------------------------------------------
538&namberg       !   iceberg parameters
539!-----------------------------------------------------------------------
540      ln_icebergs              = .false.
541      ln_bergdia               = .true.               ! Calculate budgets
542      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0
543      nn_verbose_write         = 15                   ! Timesteps between verbose messages
544      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage
545                                                      ! Initial mass required for an iceberg of each class
546      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
547                                                      ! Proportion of calving mass to apportion to each class
548      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
549                                                      ! Ratio between effective and real iceberg mass (non-dim)
550                                                      ! i.e. number of icebergs represented at a point
551      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1
552                                                      ! thickness of newly calved bergs (m)
553      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
554      rn_rho_bergs             = 850.                 ! Density of icebergs
555      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs
556      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics
557      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits
558      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
559      ln_passive_mode          = .false.              ! iceberg - ocean decoupling
560      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no)
561                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
562      rn_test_box              = 108.0,  116.0, -66.0, -58.0
563      rn_speed_limit           = 0.                   ! CFL speed limit for a berg
564
565!              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
566!              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
567      sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , ''
568
569      cn_dir = './'
570/
571
572!!======================================================================
573!!               ***  Lateral boundary condition  ***
574!!======================================================================
575!!   namlbc        lateral momentum boundary condition
576!!   namcla        cross land advection
577!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")
578!!   nambdy        Unstructured open boundaries                         ("key_bdy")
579!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides")
580!!======================================================================
581!
582!-----------------------------------------------------------------------
583&namlbc        !   lateral momentum boundary condition
584!-----------------------------------------------------------------------
585   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
586                           !  free slip  !   partial slip  !   no slip   ! strong slip
587   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs.
588/
589!-----------------------------------------------------------------------
590&namcla        !   cross land advection
591!-----------------------------------------------------------------------
592   nn_cla      =    0      !  advection between 2 ocean pts separates by land
593/
594!-----------------------------------------------------------------------
595&namagrif      !  AGRIF zoom                                            ("key_agrif")
596!-----------------------------------------------------------------------
597   nn_cln_update =    3    !  baroclinic update frequency
598   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
599   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
600   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
601/
602!-----------------------------------------------------------------------
603&nam_tide      !   tide parameters (#ifdef key_tide)
604!-----------------------------------------------------------------------
605   ln_tide_pot   = .true.   !  use tidal potential forcing
606   ln_tide_ramp  = .false.  !
607   rdttideramp   =    0.    !
608   clname(1)     = 'DUMMY'  !  name of constituent - all tidal components must be set in namelist_cfg
609/
610!-----------------------------------------------------------------------
611&nambdy        !  unstructured open boundaries                          ("key_bdy")
612!-----------------------------------------------------------------------
613    nb_bdy         = 0                    !  number of open boundary sets
614    ln_coords_file = .true.               !  =T : read bdy coordinates from file
615    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files
616    ln_mask_file   = .false.              !  =T : read mask from file
617    cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.)
618    cn_dyn2d       = 'none'               !
619    nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state
620                                          !  = 1, bdy data are read in 'bdydata   .nc' files
621                                          !  = 2, use tidal harmonic forcing data from files
622                                          !  = 3, use external data AND tidal harmonic forcing
623    cn_dyn3d      =  'none'               !
624    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state
625                                          !  = 1, bdy data are read in 'bdydata   .nc' files
626    cn_tra        =  'none'               !
627    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state
628                                          !  = 1, bdy data are read in 'bdydata   .nc' files
629    cn_ice_lim      =  'none'             !
630    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state
631                                          !  = 1, bdy data are read in 'bdydata   .nc' files
632    rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice
633    rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           --
634    rn_ice_age      = 30.                 !  lim3 only:      --   age                --
635
636    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers
637    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities
638    rn_time_dmp   =  1.                   ! Damping time scale in days
639    rn_time_dmp_out =  1.                 ! Outflow damping time scale
640    nn_rimwidth   = 10                    !  width of the relaxation zone
641    ln_vol        = .false.               !  total volume correction (see nn_volctl parameter)
642    nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero
643/
644!-----------------------------------------------------------------------
645&nambdy_dta      !  open boundaries - external data           ("key_bdy")
646!-----------------------------------------------------------------------
647!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
648!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
649   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
650   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
651   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
652   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
653   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
654   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
655   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
656! for lim2
657!   bn_frld  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
658!   bn_hicif =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
659!   bn_hsnif =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
660! for lim3
661!   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
662!   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
663!   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
664   cn_dir  =    'bdydta/'
665   ln_full_vel = .false.
666/
667!-----------------------------------------------------------------------
668&nambdy_tide     ! tidal forcing at open boundaries
669!-----------------------------------------------------------------------
670   filtide          = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files
671   ln_bdytide_2ddta = .false.
672   ln_bdytide_conj  = .false.
673/
674!!======================================================================
675!!                 ***  Bottom boundary condition  ***
676!!======================================================================
677!!   nambfr        bottom friction
678!!   nambbc        bottom temperature boundary condition
679!!   nambbl        bottom boundary layer scheme                         ("key_trabbl")
680!!======================================================================
681!
682!-----------------------------------------------------------------------
683&nambfr        !   bottom friction
684!-----------------------------------------------------------------------
685   nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
686                           !                              = 2 : nonlinear friction
687   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case)
688   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
689   rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T)
690   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2)
691   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T
692   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file )
693   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T)
694   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case)
695   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
696   rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T)
697   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2)
698   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T
699   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file )
700   rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T)
701
702   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true)
703   ln_loglayer = .false.   !  logarithmic formulation (non linear case)
704/
705!-----------------------------------------------------------------------
706&nambbc        !   bottom temperature boundary condition
707!-----------------------------------------------------------------------
708!              !                              !  (if <0  months)  ! 
709!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
710!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
711   sn_qgh      ='geothermal_heating.nc',  -12.  , 'heatflow'      ,   .false.      , .true.  , 'yearly'  , ''       , ''       , ''
712   !
713   cn_dir      = './'      !  root directory for the location of the runoff files
714   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom
715   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux
716                           !     = 1 constant flux
717                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
718   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
719
720/
721!-----------------------------------------------------------------------
722&nambbl        !   bottom boundary layer scheme
723!-----------------------------------------------------------------------
724   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0)
725   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0)
726   rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s]
727   rn_gambbl   =  10.    !  advective bbl coefficient                 [s]
728/
729
730!!======================================================================
731!!                        Tracer (T & S ) namelists
732!!======================================================================
733!!   nameos        equation of state
734!!   namtra_adv    advection scheme
735!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.)
736!!   namtra_ldf    lateral diffusion scheme
737!!   namtra_dmp    T & S newtonian damping
738!!======================================================================
739!
740!-----------------------------------------------------------------------
741&nameos        !   ocean physical parameters
742!-----------------------------------------------------------------------
743   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency
744                                 !  =-1, TEOS-10
745                                 !  = 0, EOS-80
746                                 !  = 1, S-EOS   (simplified eos)
747   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm
748   !                             !
749   !                     ! S-EOS coefficients :
750   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
751   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1)
752   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1)
753   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
754   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
755   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
756   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
757   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
758/
759!-----------------------------------------------------------------------
760&namtra_adv    !   advection scheme for tracer
761!-----------------------------------------------------------------------
762   ln_traadv_cen2   =  .false.   !  2nd order centered scheme
763   ln_traadv_tvd    =  .true.    !  TVD scheme
764   ln_traadv_muscl  =  .false.   !  MUSCL scheme
765   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries
766   ln_traadv_ubs    =  .false.   !  UBS scheme
767   ln_traadv_qck    =  .false.   !  QUICKEST scheme
768   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl
769   ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection
770/
771!-----------------------------------------------------------------------
772&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param)
773!-----------------------------------------------------------------------
774   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation
775   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08)
776   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
777   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0)
778   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0)
779   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
780   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
781   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE
782   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK
783/
784!----------------------------------------------------------------------------------
785&namtra_ldf    !   lateral diffusion scheme for tracers
786!----------------------------------------------------------------------------------
787   !                       !  Operator type:
788   ln_traldf_lap    =  .true.   !  laplacian operator
789   ln_traldf_bilap  =  .false.  !  bilaplacian operator
790   !                       !  Direction of action:
791   ln_traldf_level  =  .false.  !  iso-level
792   ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T)
793   ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp")
794   !                 !  Griffies parameters              (all need "key_ldfslp")
795   ln_traldf_grif   =  .false.  !  use griffies triads
796   ln_traldf_gdia   =  .false.  !  output griffies eddy velocities
797   ln_triad_iso     =  .false.  !  pure lateral mixing in ML
798   ln_botmix_grif   =  .false.  !  lateral mixing on bottom
799   !                       !  Coefficients
800   ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv"
801   ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv
802   !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05)
803   rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s]
804   rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s]
805   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s]
806   !                                           (normally=0; not used with Griffies)
807   rn_slpmax        =     0.01  !  slope limit
808   rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity
809   rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer
810   rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s)
811/
812!-----------------------------------------------------------------------
813&namtra_dmp    !   tracer: T & S newtonian damping
814!-----------------------------------------------------------------------
815   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F)
816   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column
817                           !                   =1 no damping in the mixing layer (kz  criteria)
818                           !                   =2 no damping in the mixed  layer (rho crieria)
819   cn_resto    = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this)
820/
821
822!!======================================================================
823!!                      ***  Dynamics namelists  ***
824!!======================================================================
825!!   namdyn_adv    formulation of the momentum advection
826!!   namdyn_vor    advection scheme
827!!   namdyn_hpg    hydrostatic pressure gradient
828!!   namdyn_spg    surface pressure gradient                            (CPP key only)
829!!   namdyn_ldf    lateral diffusion scheme
830!!======================================================================
831!
832!-----------------------------------------------------------------------
833&namdyn_adv    !   formulation of the momentum advection
834!-----------------------------------------------------------------------
835   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)
836   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction
837   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
838   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
839   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection
840/
841!-----------------------------------------------------------------------
842&nam_vvl    !   vertical coordinate options
843!-----------------------------------------------------------------------
844   ln_vvl_zstar  = .true.           !  zstar vertical coordinate
845   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations
846   ln_vvl_layer  = .false.          !  full layer vertical coordinate
847   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
848   ln_vvl_zstar_at_eqtor = .false.  !  ztilde near the equator
849   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient
850   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days]
851   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days]
852   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation
853   ln_vvl_dbg    = .true.           !  debug prints    (T/F)
854/
855!-----------------------------------------------------------------------
856&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
857!-----------------------------------------------------------------------
858   ln_dynvor_ene = .false. !  enstrophy conserving scheme
859   ln_dynvor_ens = .false. !  energy conserving scheme
860   ln_dynvor_mix = .false. !  mixed scheme
861   ln_dynvor_een = .true.  !  energy & enstrophy scheme
862   ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation
863/
864!-----------------------------------------------------------------------
865&namdyn_hpg    !   Hydrostatic pressure gradient option
866!-----------------------------------------------------------------------
867   ln_hpg_zco  = .false.   !  z-coordinate - full steps
868   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation)
869   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
870   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
871   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
872   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
873   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T)
874                                 !           centered      time scheme  (F)
875/
876!-----------------------------------------------------------------------
877!namdyn_spg    !   surface pressure gradient   (CPP key only)
878!-----------------------------------------------------------------------
879!                          !  explicit free surface                     ("key_dynspg_exp")
880!                          !  filtered free surface                     ("key_dynspg_flt")
881!                          !  split-explicit free surface               ("key_dynspg_ts")
882
883!-----------------------------------------------------------------------
884&namdyn_ldf    !   lateral diffusion on momentum
885!-----------------------------------------------------------------------
886   !                       !  Type of the operator :
887   ln_dynldf_lap    =  .true.   !  laplacian operator
888   ln_dynldf_bilap  =  .false.  !  bilaplacian operator
889   !                       !  Direction of action  :
890   ln_dynldf_level  =  .false.  !  iso-level
891   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.)
892   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp")
893   !                       !  Coefficient
894   rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s]
895   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s]
896   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s]
897   rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity
898   rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity
899   rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity
900   rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp)
901   rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap)
902/
903
904!!======================================================================
905!!             Tracers & Dynamics vertical physics namelists
906!!======================================================================
907!!    namzdf        vertical physics
908!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric")
909!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke")
910!!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp")
911!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm")
912!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx")
913!!    namzdf_tmx_new    new tidal mixing parameterization               ("key_zdftmx_new")
914!!    namzdf_mldzint vertically-interpolated mixed-layer depth parameters
915!!======================================================================
916!
917!-----------------------------------------------------------------------
918&namzdf        !   vertical physics
919!-----------------------------------------------------------------------
920   rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
921   rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
922   nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0)
923   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0)
924   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F)
925   nn_evdm     =    0      !  evd apply on tracer (=0) or on tracer and momentum (=1)
926   rn_avevd    =  100.     !  evd mixing coefficient [m2/s]
927   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F)
928   nn_npc      =    1            !  frequency of application of npc
929   nn_npcp     =  365            !  npc control print frequency
930   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
931   nn_zdfexp   =    3            !  number of sub-timestep for ln_zdfexp=T
932/
933!-----------------------------------------------------------------------
934&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
935!-----------------------------------------------------------------------
936   rn_avmri    = 100.e-4   !  maximum value of the vertical viscosity
937   rn_alp      =   5.      !  coefficient of the parameterization
938   nn_ric      =   2       !  coefficient of the parameterization
939   rn_ekmfc    =   0.7     !  Factor in the Ekman depth Equation
940   rn_mldmin   =   1.0     !  minimum allowable mixed-layer depth estimate (m)
941   rn_mldmax   =1000.0     !  maximum allowable mixed-layer depth estimate (m)
942   rn_wtmix    =  10.0     !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
943   rn_wvmix    =  10.0     !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
944   ln_mldw     = .true.    !  Flag to use or not the mized layer depth param.
945/
946!-----------------------------------------------------------------------
947&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
948!-----------------------------------------------------------------------
949   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
950   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
951   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
952   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
953   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
954   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
955   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
956                           !                 = 1 bounded by the local vertical scale factor
957                           !                 = 2 first vertical derivative of mixing length bounded by 1
958                           !                 = 3 as =2 with distinct disspipative an mixing length scale
959   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
960   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
961   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
962   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
963   rn_lc       =   0.15    !  coef. associated to Langmuir cells
964   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves
965                           !        = 0 no penetration
966                           !        = 1 add a tke source below the ML
967                           !        = 2 add a tke source just at the base of the ML
968                           !        = 3 as = 1 applied on HF part of the stress    ("key_oasis3")
969   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
970   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
971                           !        = 0  constant 10 m length scale
972                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees
973   rn_c        =   0.8     !  Default value only used when nn_htau = 2 (typically never!)
974/
975!------------------------------------------------------------------------
976&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally:
977!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
978   ln_kpprimix = .true.    !  shear instability mixing
979   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s]
980   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s]
981   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability
982   rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s]
983   rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]
984   rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s]
985   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv
986   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt
987/
988!-----------------------------------------------------------------------
989&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls")
990!-----------------------------------------------------------------------
991   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
992   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
993   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
994   rn_clim_galp  = 0.267   !  galperin limit
995   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
996   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
997   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
998   rn_hsro       =  0.02   !  Minimum surface roughness
999   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2)
1000   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2)
1001   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
1002   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
1003   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
1004   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
1005/
1006!-----------------------------------------------------------------------
1007&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
1008!-----------------------------------------------------------------------
1009   rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
1010   rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
1011/
1012!-----------------------------------------------------------------------
1013&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
1014!-----------------------------------------------------------------------
1015   rn_htmx     = 500.      !  vertical decay scale for turbulence (meters)
1016   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1)
1017   rn_tfe      = 0.333     !  tidal dissipation efficiency
1018   rn_me       = 0.2       !  mixing efficiency
1019   ln_tmx_itf  = .true.    !  ITF specific parameterisation
1020   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency
1021/
1022!-----------------------------------------------------------------------
1023&namzdf_tmx_new    !   new tidal mixing parameterization                ("key_zdftmx_new")
1024!-----------------------------------------------------------------------
1025   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
1026   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency
1027   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F)
1028/
1029!------------------------------------------------------------------------------------------
1030&namzdf_mldzint    !   Parameters for vertically-interpolated mixed-layer depth diagnostic
1031!------------------------------------------------------------------------------------------
1032   nn_mld_diag = 0         !  Number of MLD diagnostics to use from below
1033
1034!              ! MLD criterion ! Reference ! Finite difference ! Gradient layer !
1035!              ! type          ! depth     ! criterion         ! criterion      !
1036   sn_mld1     =       1       ,    10.0   ,        0.2        ,       0.1
1037   sn_mld2     =       0      ,      0.0   ,        0.0        ,       0.0 
1038   sn_mld3     =       0      ,      0.0   ,        0.0        ,       0.0 
1039   sn_mld4     =       0      ,      0.0   ,        0.0        ,       0.0 
1040   sn_mld5     =       0      ,      0.0   ,        0.0        ,       0.0 
1041/
1042
1043!!======================================================================
1044!!                  ***  Miscellaneous namelists  ***
1045!!======================================================================
1046!!   namsol            elliptic solver / island / free surface
1047!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi)
1048!!   namctl            Control prints & Benchmark
1049!!   namc1d            1D configuration options                         ("key_c1d")
1050!!   namc1d_uvd        data: U & V currents                             ("key_c1d")
1051!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d")
1052!!   namsto            Stochastic parametrization of EOS
1053!!======================================================================
1054!
1055!-----------------------------------------------------------------------
1056&namsol        !   elliptic solver / island / free surface
1057!-----------------------------------------------------------------------
1058   nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
1059                           !                   =2 successive-over-relaxation (sor)
1060   nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test
1061   rn_eps      =  1.e-6    !  absolute precision of the solver
1062   nn_nmin     =    300    !  minimum of iterations for the SOR solver
1063   nn_nmax     =    800    !  maximum of iterations for the SOR solver
1064   nn_nmod     =     10    !  frequency of test for the SOR solver
1065   rn_resmax   =  1.e-10   !  absolute precision for the SOR solver
1066   rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain)
1067/
1068!-----------------------------------------------------------------------
1069&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
1070!-----------------------------------------------------------------------
1071   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
1072                           !  buffer blocking send or immediate non-blocking sends, resp.
1073   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
1074   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold
1075   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1)
1076   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1)
1077   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1)
1078/
1079!-----------------------------------------------------------------------
1080&namctl        !   Control prints & Benchmark
1081!-----------------------------------------------------------------------
1082   ln_ctl      = .false.        !  trends control print (expensive!)
1083   sn_cfctl%l_config = .TRUE.   ! IF .true. then control which reports are written with the following
1084   sn_cfctl%l_runstat = .TRUE.  ! switches and which areas produce reports with the proc integer settings.
1085   sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure
1086   sn_cfctl%l_oceout  = .FALSE. ! that  all areas report.
1087   sn_cfctl%l_layout  = .FALSE. !
1088   sn_cfctl%l_mppout  = .FALSE. !
1089   sn_cfctl%l_mpptop  = .FALSE. !
1090   sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0]
1091   sn_cfctl%procmax   = 0       ! Maximum area number for reporting [default:0]
1092   sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1]
1093   sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info/solver.stat file
1094   nn_print    =    0       !  level of print (0 no extra print)
1095   ln_flush    =    .FALSE. !  flush numout after each print. NOT recommended for production run
1096   nn_ictls    =    0       !  start i indice of control sum (use to compare mono versus
1097   nn_ictle    =    0       !  end   i indice of control sum        multi processor runs
1098   nn_jctls    =    0       !  start j indice of control               over a subdomain)
1099   nn_jctle    =    0       !  end   j indice of control
1100   nn_isplt    =    1       !  number of processors in i-direction
1101   nn_jsplt    =    1       !  number of processors in j-direction
1102   nn_bench    =    0       !  Bench mode (1/0): CAUTION use zero except for bench
1103                            !     (no physical validity of the results)
1104   nn_timing   =    0       !  timing by routine activated (=1) creates timing.output file, or not (=0)
1105/
1106!-----------------------------------------------------------------------
1107&namc1d_uvd    !   data: U & V currents                                 ("key_c1d")
1108!-----------------------------------------------------------------------
1109!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1110!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
1111   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , ''
1112   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , ''
1113!
1114   cn_dir        = './'    !  root directory for the location of the files
1115   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F)
1116   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F)
1117/
1118!-----------------------------------------------------------------------
1119&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d")
1120!-----------------------------------------------------------------------
1121   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
1122/
1123!-----------------------------------------------------------------------
1124&namsto       ! Stochastic parametrization of EOS
1125!-----------------------------------------------------------------------
1126   ln_rststo = .false.           ! start from mean parameter (F) or from restart file (T)
1127   ln_rstseed = .true.           ! read seed of RNG from restart file
1128   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
1129   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)
1130
1131   ln_sto_eos = .false.          ! stochastic equation of state
1132   nn_sto_eos = 1                ! number of independent random walks
1133   rn_eos_stdxy = 1.4            ! random walk horz. standard deviation (in grid points)
1134   rn_eos_stdz  = 0.7            ! random walk vert. standard deviation (in grid points)
1135   rn_eos_tcor  = 1440.0         ! random walk time correlation (in timesteps)
1136   nn_eos_ord  = 1               ! order of autoregressive processes
1137   nn_eos_flt  = 0               ! passes of Laplacian filter
1138   rn_eos_lim  = 2.0             ! limitation factor (default = 3.0)
1139/
1140
1141!!======================================================================
1142!!                  ***  Diagnostics namelists  ***
1143!!======================================================================
1144!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1145!!   namtrd       dynamics and/or tracer trends
1146!!   namptr       Poleward Transport Diagnostics
1147!!   namflo       float parameters                                      ("key_float")
1148!!   namhsb       Heat and salt budgets
1149!!======================================================================
1150!
1151!-----------------------------------------------------------------------
1152&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
1153!-----------------------------------------------------------------------
1154   nn_nchunks_i=   4       !  number of chunks in i-dimension
1155   nn_nchunks_j=   4       !  number of chunks in j-dimension
1156   nn_nchunks_k=   31      !  number of chunks in k-dimension
1157                           !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1158                           !  is optimal for postprocessing which works exclusively with horizontal slabs
1159   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1160                           !  (F) ignore chunking information and produce netcdf3-compatible files
1161/
1162!-----------------------------------------------------------------------
1163&namtrd        !   diagnostics on dynamics and/or tracer trends
1164!              !       and/or mixed-layer trends and/or barotropic vorticity
1165!-----------------------------------------------------------------------
1166   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1167   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1168   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1169   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet)
1170   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1171   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1172   ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output
1173   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1174   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1175/
1176!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1177!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1178!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1179!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1180!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1181!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1182!!gm
1183!-----------------------------------------------------------------------
1184&namflo       !   float parameters                                      ("key_float")
1185!-----------------------------------------------------------------------
1186   jpnfl         = 1          !  total number of floats during the run
1187   jpnnewflo     = 0          !  number of floats for the restart
1188   ln_rstflo     = .false.    !  float restart (T) or not (F)
1189   nn_writefl    =      75    !  frequency of writing in float output file
1190   nn_stockfl    =    5475    !  frequency of creation of the float restart file
1191   ln_argo       = .false.    !  Argo type floats (stay at the surface each 10 days)
1192   ln_flork4     = .false.    !  trajectories computed with a 4th order Runge-Kutta (T)
1193                              !  or computed with Blanke' scheme (F)
1194   ln_ariane     = .true.     !  Input with Ariane tool convention(T)
1195   ln_flo_ascii  = .true.     !  Output with Ariane tool netcdf convention(F) or ascii file (T)
1196/
1197!-----------------------------------------------------------------------
1198&namptr       !   Poleward Transport Diagnostic
1199!-----------------------------------------------------------------------
1200   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F)
1201   ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not
1202/
1203!-----------------------------------------------------------------------
1204&namhsb       !  Heat and salt budgets
1205!-----------------------------------------------------------------------
1206   ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F)
1207/
1208!-----------------------------------------------------------------------
1209&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm')
1210!-----------------------------------------------------------------------
1211    nit000_han = 1         ! First time step used for harmonic analysis
1212    nitend_han = 75        ! Last time step used for harmonic analysis
1213    nstep_han  = 15        ! Time step frequency for harmonic analysis
1214    tname(1)   = 'M2'      ! Name of tidal constituents
1215    tname(2)   = 'K1'
1216/
1217!-----------------------------------------------------------------------
1218&namdct        ! transports through sections
1219!-----------------------------------------------------------------------
1220    nn_dct      = 15       !  time step frequency for transports computing
1221    nn_dctwri   = 15       !  time step frequency for transports writing
1222    nn_secdebug = 112      !      0 : no section to debug
1223                           !     -1 : debug all section
1224                           !  0 < n : debug section number n
1225/
1226
1227!!======================================================================
1228!!            ***  Observation & Assimilation namelists ***
1229!!======================================================================
1230!!   namobs       observation and model comparison                      ('key_diaobs')
1231!!   nam_asminc   assimilation increments                               ('key_asminc')
1232!!======================================================================
1233!
1234!-----------------------------------------------------------------------
1235&namobs       !  observation usage switch                               ('key_diaobs')
1236!-----------------------------------------------------------------------
1237   ln_t3d     = .false.    ! Logical switch for T profile observations
1238   ln_s3d     = .false.    ! Logical switch for S profile observations
1239   ln_ena     = .false.    ! Logical switch for ENACT insitu data set
1240   !                       !     ln_cor                  Logical switch for Coriolis insitu data set
1241   ln_profb   = .false.    ! Logical switch for feedback insitu data set
1242   ln_sla     = .false.    ! Logical switch for SLA observations
1243
1244   ln_sladt   = .false.    ! Logical switch for AVISO SLA data
1245
1246   ln_slafb   = .false.    ! Logical switch for feedback SLA data
1247                           !     ln_ssh                  Logical switch for SSH observations
1248
1249   ln_sst     = .false.     ! Logical switch for SST observations
1250   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations
1251   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations
1252
1253   ln_sstfb   = .false.    ! Logical switch for feedback SST data
1254                           !     ln_sss                  Logical switch for SSS observations
1255   ln_seaice  = .false.    ! Logical switch for Sea Ice observations
1256                           !     ln_vel3d                Logical switch for velocity observations
1257                           !     ln_velavcur             Logical switch for velocity daily av. cur.
1258                           !     ln_velhrcur             Logical switch for velocity high freq. cur.
1259                           !     ln_velavadcp            Logical switch for velocity daily av. ADCP
1260                           !     ln_velhradcp            Logical switch for velocity high freq. ADCP
1261                           !     ln_velfb                Logical switch for feedback velocity data
1262                           !     ln_grid_global          Global distribtion of observations
1263                           !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table
1264                           !     grid_search_file        Grid search lookup file header
1265                           !     enactfiles              ENACT input observation file names
1266                           !     coriofiles              Coriolis input observation file name
1267   !                       ! profbfiles: Profile feedback input observation file name
1268   profbfiles = 'profiles_01.nc'
1269                           !     ln_profb_enatim         Enact feedback input time setting switch
1270                           !     slafilesact             Active SLA input observation file name
1271                           !     slafilespas             Passive SLA input observation file name
1272   !                       ! slafbfiles: Feedback SLA input observation file name
1273   slafbfiles = 'sla_01.nc'
1274                           !     sstfiles                GHRSST input observation file name
1275   !                       ! sstfbfiles: Feedback SST input observation file name
1276   sstfbfiles = 'sst_01.nc'
1277                           !     seaicefiles             Sea Ice input observation file names
1278   seaicefiles = 'seaice_01.nc'
1279                           !     velavcurfiles           Vel. cur. daily av. input file name
1280                           !     velhvcurfiles           Vel. cur. high freq. input file name
1281                           !     velavadcpfiles          Vel. ADCP daily av. input file name
1282                           !     velhvadcpfiles          Vel. ADCP high freq. input file name
1283                           !     velfbfiles              Vel. feedback input observation file name
1284                           !     dobsini                 Initial date in window YYYYMMDD.HHMMSS
1285                           !     dobsend                 Final date in window YYYYMMDD.HHMMSS
1286                           !     n1dint                  Type of vertical interpolation method
1287                           !     n2dint                  Type of horizontal interpolation method
1288                           !     ln_nea                  Rejection of observations near land switch
1289   nmsshc     = 0          ! MSSH correction scheme
1290                           !     mdtcorr                 MDT  correction
1291                           !     mdtcutoff               MDT cutoff for computed correction
1292   ln_altbias = .false.    ! Logical switch for alt bias
1293   ln_ignmis  = .true.     ! Logical switch for ignoring missing files
1294                           !     endailyavtypes   ENACT daily average types
1295   ln_grid_global = .true.
1296   ln_grid_search_lookup = .false.
1297/
1298!-----------------------------------------------------------------------
1299&nam_asminc   !   assimilation increments                               ('key_asminc')
1300!-----------------------------------------------------------------------
1301    ln_bkgwri      = .false. !  Logical switch for writing out background state
1302    ln_balwri      = .false. !  Logical switch for writing out balancing increments
1303    ln_trainc      = .false. !  Logical switch for applying tracer increments
1304    ln_dyninc      = .false. !  Logical switch for applying velocity increments
1305    ln_sshinc      = .false. !  Logical switch for applying SSH increments
1306    ln_asmdin      = .false. !  Logical switch for Direct Initialization (DI)
1307    ln_asmiau      = .false. !  Logical switch for Incremental Analysis Updating (IAU)
1308    ln_seaiceinc   = .false. !  Logical switch for applying sea ice increments
1309    ln_phytobal    = .false. !  Logical switch for phytoplankton multivariate balancing
1310    ln_slchltotinc = .false. !  Logical switch for applying slchltot increments
1311    ln_slchldiainc = .false. !  Logical switch for applying slchldia increments
1312    ln_slchlnoninc = .false. !  Logical switch for applying slchlnon increments
1313    ln_schltotinc  = .false. !  Logical switch for applying schltot increments
1314    ln_slphytotinc = .false. !  Logical switch for applying slphytot increments
1315    ln_slphydiainc = .false. !  Logical switch for applying slphydia increments
1316    ln_slphynoninc = .false. !  Logical switch for applying slphynon increments
1317    ln_sfco2inc    = .false. !  Logical switch for applying sfCO2 increments
1318    ln_spco2inc    = .false. !  Logical switch for applying spCO2 increments
1319    ln_plchltotinc = .false. !  Logical switch for applying plchltot increments
1320    ln_pchltotinc  = .false. !  Logical switch for applying pchltot increments
1321    ln_pno3inc     = .false. !  Logical switch for applying pno3 increments
1322    ln_psi4inc     = .false. !  Logical switch for applying psi4 increments
1323    ln_pdicinc     = .false. !  Logical switch for applying pdic increments
1324    ln_palkinc     = .false. !  Logical switch for applying palk increments
1325    ln_pphinc      = .false. !  Logical switch for applying pph increments
1326    ln_po2inc      = .false. !  Logical switch for applying po2 increments
1327    ln_temnofreeze = .false. !  Logical to not add increments if temperature would fall below freezing
1328    nitbkg         = 0       !  Timestep of background in [0,nitend-nit000-1]
1329    nitdin         = 0       !  Timestep of background for DI in [0,nitend-nit000-1]
1330    nitiaustr      = 1       !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1331    nitiaufin      = 15      !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1332    niaufn         = 0       !  Type of IAU weighting function
1333    ln_salfix      = .false. !  Logical switch for ensuring that the sa > salfixmin
1334    salfixmin      = -9999   !  Minimum salinity after applying the increments
1335    nn_divdmp      = 0       !  Number of iterations of divergence damping operator
1336    mld_choice_bgc = 1       !  MLD criterion to use for biogeochemistry assimilation
1337    rn_maxchlinc   = -999.0  !  maximum absolute non-log chlorophyll increment from ocean colour assimilation
1338                             !  <= 0 implies no maximum applied (switch turned off)
1339                             !   > 0 implies maximum absolute chl increment capped at this value
1340/
1341!-----------------------------------------------------------------------
1342&namsbc_wave   ! External fields from wave model
1343!-----------------------------------------------------------------------
1344!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1345!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
1346   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1347   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1348   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1349   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1350!
1351   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files
1352/
1353!-----------------------------------------------------------------------
1354&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
1355!-----------------------------------------------------------------------
1356   ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model
1357   ln_neptsimp       = .false.  ! yes/no use simplified neptune
1358
1359   ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep
1360   rn_tslse          =  1.2e4   ! value of lengthscale L at the equator
1361   rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole
1362   ! Specify whether to ramp down the Neptune velocity in shallow
1363   ! water, and if so the depth range controlling such ramping down
1364   ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water
1365   rn_htrmin         =  100.0   ! min. depth of transition range
1366   rn_htrmax         =  200.0   ! max. depth of transition range
1367/
1368!-----------------------------------------------------------------------
1369&nambias   ! Bias pressure correctiom
1370!-----------------------------------------------------------------------
1371   ln_bias        = .false.
1372   ln_bias_asm    = .false.
1373   ln_bias_rlx    = .false.
1374   ln_bias_ofl    = .false.
1375   ln_bias_ts_app = .false.
1376   ln_bias_pc_app = .false.       
1377   fb_t_asm       = 0.0
1378   fb_t_rlx       = 0.0
1379   fb_t_ofl       = 1.0
1380   fb_p_asm       = 1.0
1381   fb_p_rlx       = 1.0
1382   fb_p_ofl       = 0.0
1383   eft_rlx        = 365.0
1384   eft_asm        = 365.0
1385   t_rlx_upd      = 0.1
1386   t_asm_upd      = 0.1
1387   nn_lat_ramp    = 0         
1388   bias_time_unit_asm = 86400.0
1389   bias_time_unit_rlx = 1.0
1390   bias_time_unit_ofl = 1.0 
1391   cn_bias_tot    = "bias_tot.nc" 
1392   cn_bias_asm    = "bias_asm.nc"
1393   cn_dir         = './' 
1394   ln_bsyncro     = .FALSE. 
1395   fctamp         = 1.
1396   rn_maxlat_bias = 23.0     
1397   rn_minlat_bias = 10.0
1398   nn_bias_itwrt  = 15
1399   ln_itdecay     = .FALSE.
1400   ln_incpc       = .FALSE.
1401/
Note: See TracBrowser for help on using the repository browser.