source: branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/CONFIG/SHARED/namelist_ref @ 11559

Last change on this file since 11559 was 11559, checked in by mattmartin, 13 months ago

Added code to include variable slip condition which is used by orca12 into the main GO6 FOAM branch so that the branch list doesn't need to be different for orca12 and orca025.

File size: 100.5 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!              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
586!              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
587   sn_shlat2d  =  'shlat_coef.nc',     -12     ,   'shlat2d'  ,    .false.     ,  .true. ,  'yearly' , ''       , ''       , ''
588   ln_shlat2d  = .false.   !  read lateral momentum boundary condition from a file (T)
589   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
590                           !  free slip  !   partial slip  !   no slip   ! strong slip
591   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs.
592/
593!-----------------------------------------------------------------------
594&namcla        !   cross land advection
595!-----------------------------------------------------------------------
596   nn_cla      =    0      !  advection between 2 ocean pts separates by land
597/
598!-----------------------------------------------------------------------
599&namagrif      !  AGRIF zoom                                            ("key_agrif")
600!-----------------------------------------------------------------------
601   nn_cln_update =    3    !  baroclinic update frequency
602   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
603   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
604   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
605/
606!-----------------------------------------------------------------------
607&nam_tide      !   tide parameters (#ifdef key_tide)
608!-----------------------------------------------------------------------
609   ln_tide_pot   = .true.   !  use tidal potential forcing
610   ln_tide_ramp  = .false.  !
611   rdttideramp   =    0.    !
612   clname(1)     = 'DUMMY'  !  name of constituent - all tidal components must be set in namelist_cfg
613/
614!-----------------------------------------------------------------------
615&nambdy        !  unstructured open boundaries                          ("key_bdy")
616!-----------------------------------------------------------------------
617    nb_bdy         = 0                    !  number of open boundary sets
618    ln_coords_file = .true.               !  =T : read bdy coordinates from file
619    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files
620    ln_mask_file   = .false.              !  =T : read mask from file
621    cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.)
622    cn_dyn2d       = 'none'               !
623    nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state
624                                          !  = 1, bdy data are read in 'bdydata   .nc' files
625                                          !  = 2, use tidal harmonic forcing data from files
626                                          !  = 3, use external data AND tidal harmonic forcing
627    cn_dyn3d      =  'none'               !
628    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state
629                                          !  = 1, bdy data are read in 'bdydata   .nc' files
630    cn_tra        =  'none'               !
631    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state
632                                          !  = 1, bdy data are read in 'bdydata   .nc' files
633    cn_ice_lim      =  'none'             !
634    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state
635                                          !  = 1, bdy data are read in 'bdydata   .nc' files
636    rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice
637    rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           --
638    rn_ice_age      = 30.                 !  lim3 only:      --   age                --
639
640    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers
641    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities
642    rn_time_dmp   =  1.                   ! Damping time scale in days
643    rn_time_dmp_out =  1.                 ! Outflow damping time scale
644    nn_rimwidth   = 10                    !  width of the relaxation zone
645    ln_vol        = .false.               !  total volume correction (see nn_volctl parameter)
646    nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero
647/
648!-----------------------------------------------------------------------
649&nambdy_dta      !  open boundaries - external data           ("key_bdy")
650!-----------------------------------------------------------------------
651!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
652!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
653   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
654   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
655   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
656   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
657   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
658   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
659   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
660! for lim2
661!   bn_frld  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
662!   bn_hicif =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
663!   bn_hsnif =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
664! for lim3
665!   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
666!   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
667!   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
668   cn_dir  =    'bdydta/'
669   ln_full_vel = .false.
670/
671!-----------------------------------------------------------------------
672&nambdy_tide     ! tidal forcing at open boundaries
673!-----------------------------------------------------------------------
674   filtide          = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files
675   ln_bdytide_2ddta = .false.
676   ln_bdytide_conj  = .false.
677/
678!!======================================================================
679!!                 ***  Bottom boundary condition  ***
680!!======================================================================
681!!   nambfr        bottom friction
682!!   nambbc        bottom temperature boundary condition
683!!   nambbl        bottom boundary layer scheme                         ("key_trabbl")
684!!======================================================================
685!
686!-----------------------------------------------------------------------
687&nambfr        !   bottom friction
688!-----------------------------------------------------------------------
689   nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
690                           !                              = 2 : nonlinear friction
691   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case)
692   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
693   rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T)
694   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2)
695   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T
696   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file )
697   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T)
698   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case)
699   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
700   rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T)
701   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2)
702   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T
703   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file )
704   rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T)
705
706   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true)
707   ln_loglayer = .false.   !  logarithmic formulation (non linear case)
708/
709!-----------------------------------------------------------------------
710&nambbc        !   bottom temperature boundary condition
711!-----------------------------------------------------------------------
712!              !                              !  (if <0  months)  ! 
713!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
714!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
715   sn_qgh      ='geothermal_heating.nc',  -12.  , 'heatflow'      ,   .false.      , .true.  , 'yearly'  , ''       , ''       , ''
716   !
717   cn_dir      = './'      !  root directory for the location of the runoff files
718   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom
719   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux
720                           !     = 1 constant flux
721                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
722   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
723
724/
725!-----------------------------------------------------------------------
726&nambbl        !   bottom boundary layer scheme
727!-----------------------------------------------------------------------
728   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0)
729   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0)
730   rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s]
731   rn_gambbl   =  10.    !  advective bbl coefficient                 [s]
732/
733
734!!======================================================================
735!!                        Tracer (T & S ) namelists
736!!======================================================================
737!!   nameos        equation of state
738!!   namtra_adv    advection scheme
739!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.)
740!!   namtra_ldf    lateral diffusion scheme
741!!   namtra_dmp    T & S newtonian damping
742!!======================================================================
743!
744!-----------------------------------------------------------------------
745&nameos        !   ocean physical parameters
746!-----------------------------------------------------------------------
747   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency
748                                 !  =-1, TEOS-10
749                                 !  = 0, EOS-80
750                                 !  = 1, S-EOS   (simplified eos)
751   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm
752   !                             !
753   !                     ! S-EOS coefficients :
754   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
755   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1)
756   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1)
757   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
758   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
759   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
760   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
761   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
762/
763!-----------------------------------------------------------------------
764&namtra_adv    !   advection scheme for tracer
765!-----------------------------------------------------------------------
766   ln_traadv_cen2   =  .false.   !  2nd order centered scheme
767   ln_traadv_tvd    =  .true.    !  TVD scheme
768   ln_traadv_muscl  =  .false.   !  MUSCL scheme
769   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries
770   ln_traadv_ubs    =  .false.   !  UBS scheme
771   ln_traadv_qck    =  .false.   !  QUICKEST scheme
772   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl
773   ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection
774/
775!-----------------------------------------------------------------------
776&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param)
777!-----------------------------------------------------------------------
778   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation
779   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08)
780   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
781   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0)
782   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0)
783   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
784   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
785   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE
786   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK
787/
788!----------------------------------------------------------------------------------
789&namtra_ldf    !   lateral diffusion scheme for tracers
790!----------------------------------------------------------------------------------
791   !                       !  Operator type:
792   ln_traldf_lap    =  .true.   !  laplacian operator
793   ln_traldf_bilap  =  .false.  !  bilaplacian operator
794   !                       !  Direction of action:
795   ln_traldf_level  =  .false.  !  iso-level
796   ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T)
797   ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp")
798   !                 !  Griffies parameters              (all need "key_ldfslp")
799   ln_traldf_grif   =  .false.  !  use griffies triads
800   ln_traldf_gdia   =  .false.  !  output griffies eddy velocities
801   ln_triad_iso     =  .false.  !  pure lateral mixing in ML
802   ln_botmix_grif   =  .false.  !  lateral mixing on bottom
803   !                       !  Coefficients
804   ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv"
805   ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv
806   !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05)
807   rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s]
808   rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s]
809   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s]
810   !                                           (normally=0; not used with Griffies)
811   rn_slpmax        =     0.01  !  slope limit
812   rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity
813   rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer
814   rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s)
815/
816!-----------------------------------------------------------------------
817&namtra_dmp    !   tracer: T & S newtonian damping
818!-----------------------------------------------------------------------
819   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F)
820   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column
821                           !                   =1 no damping in the mixing layer (kz  criteria)
822                           !                   =2 no damping in the mixed  layer (rho crieria)
823   cn_resto    = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this)
824/
825
826!!======================================================================
827!!                      ***  Dynamics namelists  ***
828!!======================================================================
829!!   namdyn_adv    formulation of the momentum advection
830!!   namdyn_vor    advection scheme
831!!   namdyn_hpg    hydrostatic pressure gradient
832!!   namdyn_spg    surface pressure gradient                            (CPP key only)
833!!   namdyn_ldf    lateral diffusion scheme
834!!======================================================================
835!
836!-----------------------------------------------------------------------
837&namdyn_adv    !   formulation of the momentum advection
838!-----------------------------------------------------------------------
839   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)
840   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction
841   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
842   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
843   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection
844/
845!-----------------------------------------------------------------------
846&nam_vvl    !   vertical coordinate options
847!-----------------------------------------------------------------------
848   ln_vvl_zstar  = .true.           !  zstar vertical coordinate
849   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations
850   ln_vvl_layer  = .false.          !  full layer vertical coordinate
851   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
852   ln_vvl_zstar_at_eqtor = .false.  !  ztilde near the equator
853   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient
854   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days]
855   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days]
856   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation
857   ln_vvl_dbg    = .true.           !  debug prints    (T/F)
858/
859!-----------------------------------------------------------------------
860&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
861!-----------------------------------------------------------------------
862   ln_dynvor_ene = .false. !  enstrophy conserving scheme
863   ln_dynvor_ens = .false. !  energy conserving scheme
864   ln_dynvor_mix = .false. !  mixed scheme
865   ln_dynvor_een = .true.  !  energy & enstrophy scheme
866   ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation
867/
868!-----------------------------------------------------------------------
869&namdyn_hpg    !   Hydrostatic pressure gradient option
870!-----------------------------------------------------------------------
871   ln_hpg_zco  = .false.   !  z-coordinate - full steps
872   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation)
873   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
874   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
875   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
876   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
877   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T)
878                                 !           centered      time scheme  (F)
879/
880!-----------------------------------------------------------------------
881!namdyn_spg    !   surface pressure gradient   (CPP key only)
882!-----------------------------------------------------------------------
883!                          !  explicit free surface                     ("key_dynspg_exp")
884!                          !  filtered free surface                     ("key_dynspg_flt")
885!                          !  split-explicit free surface               ("key_dynspg_ts")
886
887!-----------------------------------------------------------------------
888&namdyn_ldf    !   lateral diffusion on momentum
889!-----------------------------------------------------------------------
890   !                       !  Type of the operator :
891   ln_dynldf_lap    =  .true.   !  laplacian operator
892   ln_dynldf_bilap  =  .false.  !  bilaplacian operator
893   !                       !  Direction of action  :
894   ln_dynldf_level  =  .false.  !  iso-level
895   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.)
896   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp")
897   !                       !  Coefficient
898   rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s]
899   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s]
900   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s]
901   rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity
902   rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity
903   rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity
904   rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp)
905   rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap)
906/
907
908!!======================================================================
909!!             Tracers & Dynamics vertical physics namelists
910!!======================================================================
911!!    namzdf        vertical physics
912!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric")
913!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke")
914!!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp")
915!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm")
916!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx")
917!!    namzdf_tmx_new    new tidal mixing parameterization               ("key_zdftmx_new")
918!!    namzdf_mldzint vertically-interpolated mixed-layer depth parameters
919!!======================================================================
920!
921!-----------------------------------------------------------------------
922&namzdf        !   vertical physics
923!-----------------------------------------------------------------------
924   rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
925   rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
926   nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0)
927   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0)
928   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F)
929   nn_evdm     =    0      !  evd apply on tracer (=0) or on tracer and momentum (=1)
930   rn_avevd    =  100.     !  evd mixing coefficient [m2/s]
931   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F)
932   nn_npc      =    1            !  frequency of application of npc
933   nn_npcp     =  365            !  npc control print frequency
934   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
935   nn_zdfexp   =    3            !  number of sub-timestep for ln_zdfexp=T
936/
937!-----------------------------------------------------------------------
938&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
939!-----------------------------------------------------------------------
940   rn_avmri    = 100.e-4   !  maximum value of the vertical viscosity
941   rn_alp      =   5.      !  coefficient of the parameterization
942   nn_ric      =   2       !  coefficient of the parameterization
943   rn_ekmfc    =   0.7     !  Factor in the Ekman depth Equation
944   rn_mldmin   =   1.0     !  minimum allowable mixed-layer depth estimate (m)
945   rn_mldmax   =1000.0     !  maximum allowable mixed-layer depth estimate (m)
946   rn_wtmix    =  10.0     !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
947   rn_wvmix    =  10.0     !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
948   ln_mldw     = .true.    !  Flag to use or not the mized layer depth param.
949/
950!-----------------------------------------------------------------------
951&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
952!-----------------------------------------------------------------------
953   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
954   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
955   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
956   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
957   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
958   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
959   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
960                           !                 = 1 bounded by the local vertical scale factor
961                           !                 = 2 first vertical derivative of mixing length bounded by 1
962                           !                 = 3 as =2 with distinct disspipative an mixing length scale
963   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
964   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
965   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
966   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
967   rn_lc       =   0.15    !  coef. associated to Langmuir cells
968   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves
969                           !        = 0 no penetration
970                           !        = 1 add a tke source below the ML
971                           !        = 2 add a tke source just at the base of the ML
972                           !        = 3 as = 1 applied on HF part of the stress    ("key_oasis3")
973   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
974   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
975                           !        = 0  constant 10 m length scale
976                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees
977   rn_c        =   0.8     !  Default value only used when nn_htau = 2 (typically never!)
978/
979!------------------------------------------------------------------------
980&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally:
981!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
982   ln_kpprimix = .true.    !  shear instability mixing
983   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s]
984   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s]
985   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability
986   rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s]
987   rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]
988   rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s]
989   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv
990   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt
991/
992!-----------------------------------------------------------------------
993&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls")
994!-----------------------------------------------------------------------
995   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
996   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
997   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
998   rn_clim_galp  = 0.267   !  galperin limit
999   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
1000   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
1001   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
1002   rn_hsro       =  0.02   !  Minimum surface roughness
1003   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2)
1004   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2)
1005   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
1006   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
1007   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
1008   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
1009/
1010!-----------------------------------------------------------------------
1011&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
1012!-----------------------------------------------------------------------
1013   rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
1014   rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
1015/
1016!-----------------------------------------------------------------------
1017&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
1018!-----------------------------------------------------------------------
1019   rn_htmx     = 500.      !  vertical decay scale for turbulence (meters)
1020   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1)
1021   rn_tfe      = 0.333     !  tidal dissipation efficiency
1022   rn_me       = 0.2       !  mixing efficiency
1023   ln_tmx_itf  = .true.    !  ITF specific parameterisation
1024   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency
1025/
1026!-----------------------------------------------------------------------
1027&namzdf_tmx_new    !   new tidal mixing parameterization                ("key_zdftmx_new")
1028!-----------------------------------------------------------------------
1029   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
1030   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency
1031   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F)
1032/
1033!------------------------------------------------------------------------------------------
1034&namzdf_mldzint    !   Parameters for vertically-interpolated mixed-layer depth diagnostic
1035!------------------------------------------------------------------------------------------
1036   nn_mld_diag = 0         !  Number of MLD diagnostics to use from below
1037
1038!              ! MLD criterion ! Reference ! Finite difference ! Gradient layer !
1039!              ! type          ! depth     ! criterion         ! criterion      !
1040   sn_mld1     =       1       ,    10.0   ,        0.2        ,       0.1
1041   sn_mld2     =       0      ,      0.0   ,        0.0        ,       0.0 
1042   sn_mld3     =       0      ,      0.0   ,        0.0        ,       0.0 
1043   sn_mld4     =       0      ,      0.0   ,        0.0        ,       0.0 
1044   sn_mld5     =       0      ,      0.0   ,        0.0        ,       0.0 
1045/
1046
1047!!======================================================================
1048!!                  ***  Miscellaneous namelists  ***
1049!!======================================================================
1050!!   namsol            elliptic solver / island / free surface
1051!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi)
1052!!   namctl            Control prints & Benchmark
1053!!   namc1d            1D configuration options                         ("key_c1d")
1054!!   namc1d_uvd        data: U & V currents                             ("key_c1d")
1055!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d")
1056!!   namsto            Stochastic parametrization of EOS
1057!!======================================================================
1058!
1059!-----------------------------------------------------------------------
1060&namsol        !   elliptic solver / island / free surface
1061!-----------------------------------------------------------------------
1062   nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
1063                           !                   =2 successive-over-relaxation (sor)
1064   nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test
1065   rn_eps      =  1.e-6    !  absolute precision of the solver
1066   nn_nmin     =    300    !  minimum of iterations for the SOR solver
1067   nn_nmax     =    800    !  maximum of iterations for the SOR solver
1068   nn_nmod     =     10    !  frequency of test for the SOR solver
1069   rn_resmax   =  1.e-10   !  absolute precision for the SOR solver
1070   rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain)
1071/
1072!-----------------------------------------------------------------------
1073&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
1074!-----------------------------------------------------------------------
1075   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
1076                           !  buffer blocking send or immediate non-blocking sends, resp.
1077   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
1078   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold
1079   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1)
1080   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1)
1081   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1)
1082/
1083!-----------------------------------------------------------------------
1084&namctl        !   Control prints & Benchmark
1085!-----------------------------------------------------------------------
1086   ln_ctl      = .false.   !  trends control print (expensive!)
1087   nn_print    =    0      !  level of print (0 no extra print)
1088   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
1089   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
1090   nn_jctls    =    0      !  start j indice of control               over a subdomain)
1091   nn_jctle    =    0      !  end   j indice of control
1092   nn_isplt    =    1      !  number of processors in i-direction
1093   nn_jsplt    =    1      !  number of processors in j-direction
1094   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench
1095                           !     (no physical validity of the results)
1096   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0)
1097/
1098!-----------------------------------------------------------------------
1099&namc1d_uvd    !   data: U & V currents                                 ("key_c1d")
1100!-----------------------------------------------------------------------
1101!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1102!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
1103   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , ''
1104   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , ''
1105!
1106   cn_dir        = './'    !  root directory for the location of the files
1107   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F)
1108   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F)
1109/
1110!-----------------------------------------------------------------------
1111&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d")
1112!-----------------------------------------------------------------------
1113   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
1114/
1115!-----------------------------------------------------------------------
1116&namsto       ! Stochastic parametrization of EOS
1117!-----------------------------------------------------------------------
1118   ln_rststo = .false.           ! start from mean parameter (F) or from restart file (T)
1119   ln_rstseed = .true.           ! read seed of RNG from restart file
1120   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
1121   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)
1122
1123   ln_sto_eos = .false.          ! stochastic equation of state
1124   nn_sto_eos = 1                ! number of independent random walks
1125   rn_eos_stdxy = 1.4            ! random walk horz. standard deviation (in grid points)
1126   rn_eos_stdz  = 0.7            ! random walk vert. standard deviation (in grid points)
1127   rn_eos_tcor  = 1440.0         ! random walk time correlation (in timesteps)
1128   nn_eos_ord  = 1               ! order of autoregressive processes
1129   nn_eos_flt  = 0               ! passes of Laplacian filter
1130   rn_eos_lim  = 2.0             ! limitation factor (default = 3.0)
1131/
1132
1133!!======================================================================
1134!!                  ***  Diagnostics namelists  ***
1135!!======================================================================
1136!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1137!!   namtrd       dynamics and/or tracer trends
1138!!   namptr       Poleward Transport Diagnostics
1139!!   namflo       float parameters                                      ("key_float")
1140!!   namhsb       Heat and salt budgets
1141!!======================================================================
1142!
1143!-----------------------------------------------------------------------
1144&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
1145!-----------------------------------------------------------------------
1146   nn_nchunks_i=   4       !  number of chunks in i-dimension
1147   nn_nchunks_j=   4       !  number of chunks in j-dimension
1148   nn_nchunks_k=   31      !  number of chunks in k-dimension
1149                           !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1150                           !  is optimal for postprocessing which works exclusively with horizontal slabs
1151   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1152                           !  (F) ignore chunking information and produce netcdf3-compatible files
1153/
1154!-----------------------------------------------------------------------
1155&namtrd        !   diagnostics on dynamics and/or tracer trends
1156!              !       and/or mixed-layer trends and/or barotropic vorticity
1157!-----------------------------------------------------------------------
1158   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1159   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1160   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1161   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet)
1162   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1163   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1164   ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output
1165   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1166   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1167/
1168!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1169!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1170!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1171!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1172!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1173!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1174!!gm
1175!-----------------------------------------------------------------------
1176&namflo       !   float parameters                                      ("key_float")
1177!-----------------------------------------------------------------------
1178   jpnfl         = 1          !  total number of floats during the run
1179   jpnnewflo     = 0          !  number of floats for the restart
1180   ln_rstflo     = .false.    !  float restart (T) or not (F)
1181   nn_writefl    =      75    !  frequency of writing in float output file
1182   nn_stockfl    =    5475    !  frequency of creation of the float restart file
1183   ln_argo       = .false.    !  Argo type floats (stay at the surface each 10 days)
1184   ln_flork4     = .false.    !  trajectories computed with a 4th order Runge-Kutta (T)
1185                              !  or computed with Blanke' scheme (F)
1186   ln_ariane     = .true.     !  Input with Ariane tool convention(T)
1187   ln_flo_ascii  = .true.     !  Output with Ariane tool netcdf convention(F) or ascii file (T)
1188/
1189!-----------------------------------------------------------------------
1190&namptr       !   Poleward Transport Diagnostic
1191!-----------------------------------------------------------------------
1192   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F)
1193   ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not
1194/
1195!-----------------------------------------------------------------------
1196&namhsb       !  Heat and salt budgets
1197!-----------------------------------------------------------------------
1198   ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F)
1199/
1200!-----------------------------------------------------------------------
1201&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm')
1202!-----------------------------------------------------------------------
1203    nit000_han = 1         ! First time step used for harmonic analysis
1204    nitend_han = 75        ! Last time step used for harmonic analysis
1205    nstep_han  = 15        ! Time step frequency for harmonic analysis
1206    tname(1)   = 'M2'      ! Name of tidal constituents
1207    tname(2)   = 'K1'
1208/
1209!-----------------------------------------------------------------------
1210&namdct        ! transports through sections
1211!-----------------------------------------------------------------------
1212    nn_dct      = 15       !  time step frequency for transports computing
1213    nn_dctwri   = 15       !  time step frequency for transports writing
1214    nn_secdebug = 112      !      0 : no section to debug
1215                           !     -1 : debug all section
1216                           !  0 < n : debug section number n
1217/
1218
1219!!======================================================================
1220!!            ***  Observation & Assimilation namelists ***
1221!!======================================================================
1222!!   namobs       observation and model comparison                      ('key_diaobs')
1223!!   nam_asminc   assimilation increments                               ('key_asminc')
1224!!======================================================================
1225!
1226!-----------------------------------------------------------------------
1227&namobs       !  observation usage switch                               ('key_diaobs')
1228!-----------------------------------------------------------------------
1229   ln_t3d     = .false.    ! Logical switch for T profile observations
1230   ln_s3d     = .false.    ! Logical switch for S profile observations
1231   ln_ena     = .false.    ! Logical switch for ENACT insitu data set
1232   !                       !     ln_cor                  Logical switch for Coriolis insitu data set
1233   ln_profb   = .false.    ! Logical switch for feedback insitu data set
1234   ln_sla     = .false.    ! Logical switch for SLA observations
1235
1236   ln_sladt   = .false.    ! Logical switch for AVISO SLA data
1237
1238   ln_slafb   = .false.    ! Logical switch for feedback SLA data
1239                           !     ln_ssh                  Logical switch for SSH observations
1240
1241   ln_sst     = .false.     ! Logical switch for SST observations
1242   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations
1243   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations
1244
1245   ln_sstfb   = .false.    ! Logical switch for feedback SST data
1246                           !     ln_sss                  Logical switch for SSS observations
1247   ln_seaice  = .false.    ! Logical switch for Sea Ice observations
1248                           !     ln_vel3d                Logical switch for velocity observations
1249                           !     ln_velavcur             Logical switch for velocity daily av. cur.
1250                           !     ln_velhrcur             Logical switch for velocity high freq. cur.
1251                           !     ln_velavadcp            Logical switch for velocity daily av. ADCP
1252                           !     ln_velhradcp            Logical switch for velocity high freq. ADCP
1253                           !     ln_velfb                Logical switch for feedback velocity data
1254                           !     ln_grid_global          Global distribtion of observations
1255                           !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table
1256                           !     grid_search_file        Grid search lookup file header
1257                           !     enactfiles              ENACT input observation file names
1258                           !     coriofiles              Coriolis input observation file name
1259   !                       ! profbfiles: Profile feedback input observation file name
1260   profbfiles = 'profiles_01.nc'
1261                           !     ln_profb_enatim         Enact feedback input time setting switch
1262                           !     slafilesact             Active SLA input observation file name
1263                           !     slafilespas             Passive SLA input observation file name
1264   !                       ! slafbfiles: Feedback SLA input observation file name
1265   slafbfiles = 'sla_01.nc'
1266                           !     sstfiles                GHRSST input observation file name
1267   !                       ! sstfbfiles: Feedback SST input observation file name
1268   sstfbfiles = 'sst_01.nc'
1269                           !     seaicefiles             Sea Ice input observation file names
1270   seaicefiles = 'seaice_01.nc'
1271                           !     velavcurfiles           Vel. cur. daily av. input file name
1272                           !     velhvcurfiles           Vel. cur. high freq. input file name
1273                           !     velavadcpfiles          Vel. ADCP daily av. input file name
1274                           !     velhvadcpfiles          Vel. ADCP high freq. input file name
1275                           !     velfbfiles              Vel. feedback input observation file name
1276                           !     dobsini                 Initial date in window YYYYMMDD.HHMMSS
1277                           !     dobsend                 Final date in window YYYYMMDD.HHMMSS
1278                           !     n1dint                  Type of vertical interpolation method
1279                           !     n2dint                  Type of horizontal interpolation method
1280                           !     ln_nea                  Rejection of observations near land switch
1281   nmsshc     = 0          ! MSSH correction scheme
1282                           !     mdtcorr                 MDT  correction
1283                           !     mdtcutoff               MDT cutoff for computed correction
1284   ln_altbias = .false.    ! Logical switch for alt bias
1285   ln_ignmis  = .true.     ! Logical switch for ignoring missing files
1286                           !     endailyavtypes   ENACT daily average types
1287   ln_grid_global = .true.
1288   ln_grid_search_lookup = .false.
1289/
1290!-----------------------------------------------------------------------
1291&nam_asminc   !   assimilation increments                               ('key_asminc')
1292!-----------------------------------------------------------------------
1293    ln_bkgwri      = .false. !  Logical switch for writing out background state
1294    ln_balwri      = .false. !  Logical switch for writing out balancing increments
1295    ln_trainc      = .false. !  Logical switch for applying tracer increments
1296    ln_dyninc      = .false. !  Logical switch for applying velocity increments
1297    ln_sshinc      = .false. !  Logical switch for applying SSH increments
1298    ln_asmdin      = .false. !  Logical switch for Direct Initialization (DI)
1299    ln_asmiau      = .false. !  Logical switch for Incremental Analysis Updating (IAU)
1300    ln_seaiceinc   = .false. !  Logical switch for applying sea ice increments
1301    ln_phytobal    = .false. !  Logical switch for phytoplankton multivariate balancing
1302    ln_slchltotinc = .false. !  Logical switch for applying slchltot increments
1303    ln_slchldiainc = .false. !  Logical switch for applying slchldia increments
1304    ln_slchlnoninc = .false. !  Logical switch for applying slchlnon increments
1305    ln_schltotinc  = .false. !  Logical switch for applying schltot increments
1306    ln_slphytotinc = .false. !  Logical switch for applying slphytot increments
1307    ln_slphydiainc = .false. !  Logical switch for applying slphydia increments
1308    ln_slphynoninc = .false. !  Logical switch for applying slphynon increments
1309    ln_sfco2inc    = .false. !  Logical switch for applying sfCO2 increments
1310    ln_spco2inc    = .false. !  Logical switch for applying spCO2 increments
1311    ln_plchltotinc = .false. !  Logical switch for applying plchltot increments
1312    ln_pchltotinc  = .false. !  Logical switch for applying pchltot increments
1313    ln_pno3inc     = .false. !  Logical switch for applying pno3 increments
1314    ln_psi4inc     = .false. !  Logical switch for applying psi4 increments
1315    ln_pdicinc     = .false. !  Logical switch for applying pdic increments
1316    ln_palkinc     = .false. !  Logical switch for applying palk increments
1317    ln_pphinc      = .false. !  Logical switch for applying pph increments
1318    ln_po2inc      = .false. !  Logical switch for applying po2 increments
1319    ln_temnofreeze = .false. !  Logical to not add increments if temperature would fall below freezing
1320    nitbkg         = 0       !  Timestep of background in [0,nitend-nit000-1]
1321    nitdin         = 0       !  Timestep of background for DI in [0,nitend-nit000-1]
1322    nitiaustr      = 1       !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1323    nitiaufin      = 15      !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1324    niaufn         = 0       !  Type of IAU weighting function
1325    ln_salfix      = .false. !  Logical switch for ensuring that the sa > salfixmin
1326    salfixmin      = -9999   !  Minimum salinity after applying the increments
1327    nn_divdmp      = 0       !  Number of iterations of divergence damping operator
1328    mld_choice_bgc = 1       !  MLD criterion to use for biogeochemistry assimilation
1329    rn_maxchlinc   = -999.0  !  maximum absolute non-log chlorophyll increment from ocean colour assimilation
1330                             !  <= 0 implies no maximum applied (switch turned off)
1331                             !   > 0 implies maximum absolute chl increment capped at this value
1332/
1333!-----------------------------------------------------------------------
1334&namsbc_wave   ! External fields from wave model
1335!-----------------------------------------------------------------------
1336!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1337!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
1338   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1339   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1340   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1341   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1342!
1343   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files
1344/
1345!-----------------------------------------------------------------------
1346&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
1347!-----------------------------------------------------------------------
1348   ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model
1349   ln_neptsimp       = .false.  ! yes/no use simplified neptune
1350
1351   ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep
1352   rn_tslse          =  1.2e4   ! value of lengthscale L at the equator
1353   rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole
1354   ! Specify whether to ramp down the Neptune velocity in shallow
1355   ! water, and if so the depth range controlling such ramping down
1356   ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water
1357   rn_htrmin         =  100.0   ! min. depth of transition range
1358   rn_htrmax         =  200.0   ! max. depth of transition range
1359/
1360!-----------------------------------------------------------------------
1361&nambias   ! Bias pressure correctiom
1362!-----------------------------------------------------------------------
1363   ln_bias        = .false.
1364   ln_bias_asm    = .false.
1365   ln_bias_rlx    = .false.
1366   ln_bias_ofl    = .false.
1367   ln_bias_ts_app = .false.
1368   ln_bias_pc_app = .false.       
1369   fb_t_asm       = 0.0
1370   fb_t_rlx       = 0.0
1371   fb_t_ofl       = 1.0
1372   fb_p_asm       = 1.0
1373   fb_p_rlx       = 1.0
1374   fb_p_ofl       = 0.0
1375   eft_rlx        = 365.0
1376   eft_asm        = 365.0
1377   t_rlx_upd      = 0.1
1378   t_asm_upd      = 0.1
1379   nn_lat_ramp    = 0         
1380   bias_time_unit_asm = 86400.0
1381   bias_time_unit_rlx = 1.0
1382   bias_time_unit_ofl = 1.0 
1383   cn_bias_tot    = "bias_tot.nc" 
1384   cn_bias_asm    = "bias_asm.nc"
1385   cn_dir         = './' 
1386   ln_bsyncro     = .FALSE. 
1387   fctamp         = 1.
1388   rn_maxlat_bias = 23.0     
1389   rn_minlat_bias = 10.0
1390   nn_bias_itwrt  = 15
1391   ln_itdecay     = .FALSE.
1392   ln_incpc       = .FALSE.
1393/
Note: See TracBrowser for help on using the repository browser.