source: branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/CONFIG/SHARED/namelist_ref @ 12140

Last change on this file since 12140 was 12140, checked in by kingr, 20 months ago

Added option to remove tides from SLA bkg by taking average over 24h50m.

File size: 98.4 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, namobc, 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)
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! receive
366   sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
367   sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
368   sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V'
369   sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
370   sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
371   sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
372   sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   ''
373   sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
374   sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
375   sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
376!
377   nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data
378   ln_usecplmask = .false.   !  use a coupling mask file to merge data received from several models
379                             !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
380/
381!-----------------------------------------------------------------------
382&namsbc_sas    !   analytical surface boundary condition
383!-----------------------------------------------------------------------
384!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
385!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
386   sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
387   sn_vsp      = 'sas_grid_V' ,    120           , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
388   sn_tem      = 'sas_grid_T' ,    120           , 'sosstsst' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
389   sn_sal      = 'sas_grid_T' ,    120           , 'sosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
390   sn_ssh      = 'sas_grid_T' ,    120           , 'sossheig' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
391   sn_e3t      = 'sas_grid_T' ,    120           , 'e3t_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
392   sn_frq      = 'sas_grid_T' ,    120           , 'frq_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
393
394   ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field
395   ln_read_frq = .false.    !  specify whether we must read frq or not
396   cn_dir      = './'      !  root directory for the location of the bulk files are
397/
398!-----------------------------------------------------------------------
399&namtra_qsr    !   penetrative solar radiation
400!-----------------------------------------------------------------------
401!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
402!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
403   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , ''
404
405   cn_dir      = './'      !  root directory for the location of the runoff files
406   ln_traqsr   = .true.    !  Light penetration (T) or not (F)
407   ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration
408   ln_qsr_2bd  = .false.   !  2 bands              light penetration
409   ln_qsr_bio  = .false.   !  bio-model light penetration
410   nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0)
411   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
412   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
413   rn_si1      =   23.0    !  2 bands: longest depth of extinction
414   ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3
415/
416!-----------------------------------------------------------------------
417&namsbc_rnf    !   runoffs namelist surface boundary condition
418!-----------------------------------------------------------------------
419!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
420!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
421   sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
422   sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
423   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
424   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
425   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , ''
426
427   cn_dir       = './'      !  root directory for the location of the runoff files
428   ln_rnf_mouth = .true.    !  specific treatment at rivers mouths
429   rn_hrnf      =  15.e0    !  depth over which enhanced vertical mixing is used
430   rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s]
431   rn_rfact     =   1.e0    !  multiplicative factor for runoff
432   ln_rnf_depth = .false.   !  read in depth information for runoff
433   ln_rnf_tem   = .false.   !  read in temperature information for runoff
434   ln_rnf_sal   = .false.   !  read in salinity information for runoff
435   ln_rnf_depth_ini = .false.  ! compute depth at initialisation from runoff file
436   rn_rnf_max   = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
437   rn_dep_max   = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
438   nn_rnf_depth_file = 0    !  create (=1) a runoff depth file or not (=0)
439/
440!-----------------------------------------------------------------------
441&namsbc_isf    !  Top boundary layer (ISF)
442!-----------------------------------------------------------------------
443!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
444!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
445! nn_isf == 4
446   sn_qisf      = 'rnfisf' ,         -12      ,'sohflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
447   sn_fwfisf    = 'rnfisf' ,         -12      ,'sowflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
448! nn_isf == 3
449   sn_rnfisf    = 'runoffs' ,         -12      ,'sofwfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
450! nn_isf == 2 and 3
451   sn_depmax_isf = 'runoffs' ,       -12        ,'sozisfmax' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
452   sn_depmin_isf = 'runoffs' ,       -12        ,'sozisfmin' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
453! nn_isf == 2
454   sn_Leff_isf = 'rnfisf' ,       0          ,'Leff'         ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
455! for all case
456   ln_divisf   = .true.  ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?)
457! only for nn_isf = 1 or 2
458   rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula
459   rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula
460! only for nn_isf = 1
461   nn_isfblk   =  1       ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??)
462   rn_hisf_tbl =  30.      ! thickness of the top boundary layer           (Losh et al. 2008)
463                          ! 0 => thickness of the tbl = thickness of the first wet cell
464   ln_conserve = .true.   ! conservative case (take into account meltwater advection)
465   nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s)
466                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010)
467                          !     if you want to keep the cd as in global config, adjust rn_gammat0 to compensate
468                          ! 2 = velocity and stability dependent Gamma    Holland et al. 1999
469/
470!-----------------------------------------------------------------------
471&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk
472!-----------------------------------------------------------------------
473!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
474!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
475   sn_apr      = 'patm'      ,         -1        ,'somslpre',    .true.     , .true. , 'yearly'  ,  ''      ,   ''     , ''
476
477   cn_dir      = './'       !  root directory for the location of the bulk files
478   rn_pref     = 101000.    !  reference atmospheric pressure   [N/m2]/
479   ln_ref_apr  = .false.    !  ref. pressure: global mean Patm (T) or a constant (F)
480   ln_apr_obc  = .false.    !  inverse barometer added to OBC ssh data
481/
482!-----------------------------------------------------------------------
483&namsbc_ssr    !   surface boundary condition : sea surface restoring
484!-----------------------------------------------------------------------
485!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
486!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
487   sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , ''       , ''
488   sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , ''       , ''
489
490   cn_dir      = './'      !  root directory for the location of the runoff files
491   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
492   nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2)
493                           !  or to SSS only (=1) or no damping term (=0)
494   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K]
495   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
496   ln_sssr_bnd =   .true.  !  flag to bound erp term (associated with nn_sssr=2)
497   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
498/
499!-----------------------------------------------------------------------
500&namsbc_alb    !   albedo parameters
501!-----------------------------------------------------------------------
502   rn_cloud    =    0.06   !  cloud correction to snow and ice albedo
503   rn_albice   =    0.53   !  albedo of melting ice in the arctic and antarctic
504   rn_alphd    =    0.80   !  coefficients for linear interpolation used to
505   rn_alphc    =    0.65   !  compute albedo between two extremes values
506   rn_alphdi   =    0.72   !  (Pyane, 1972)
507/
508!-----------------------------------------------------------------------
509&namberg       !   iceberg parameters
510!-----------------------------------------------------------------------
511      ln_icebergs              = .false.
512      ln_bergdia               = .true.               ! Calculate budgets
513      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0
514      nn_verbose_write         = 15                   ! Timesteps between verbose messages
515      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage
516                                                      ! Initial mass required for an iceberg of each class
517      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
518                                                      ! Proportion of calving mass to apportion to each class
519      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
520                                                      ! Ratio between effective and real iceberg mass (non-dim)
521                                                      ! i.e. number of icebergs represented at a point
522      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1
523                                                      ! thickness of newly calved bergs (m)
524      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
525      rn_rho_bergs             = 850.                 ! Density of icebergs
526      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs
527      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics
528      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits
529      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
530      ln_passive_mode          = .false.              ! iceberg - ocean decoupling
531      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no)
532                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
533      rn_test_box              = 108.0,  116.0, -66.0, -58.0
534      rn_speed_limit           = 0.                   ! CFL speed limit for a berg
535
536!              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
537!              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
538      sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , ''
539
540      cn_dir = './'
541/
542
543!!======================================================================
544!!               ***  Lateral boundary condition  ***
545!!======================================================================
546!!   namlbc        lateral momentum boundary condition
547!!   namcla        cross land advection
548!!   namobc        open boundaries parameters                           ("key_obc")
549!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")
550!!   nambdy        Unstructured open boundaries                         ("key_bdy")
551!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides")
552!!======================================================================
553!
554!-----------------------------------------------------------------------
555&namlbc        !   lateral momentum boundary condition
556!-----------------------------------------------------------------------
557   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
558                           !  free slip  !   partial slip  !   no slip   ! strong slip
559   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs.
560/
561!-----------------------------------------------------------------------
562&namcla        !   cross land advection
563!-----------------------------------------------------------------------
564   nn_cla      =    0      !  advection between 2 ocean pts separates by land
565/
566!-----------------------------------------------------------------------
567&namobc        !   open boundaries parameters                           ("key_obc")
568!-----------------------------------------------------------------------
569   ln_obc_clim = .false.   !  climatological obc data files (T) or not (F)
570   ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F)
571   ln_obc_fla  = .false.   !  Flather open boundary condition
572   nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state
573                           !  = 1 the obc data are read in 'obc.dta' files
574   cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data
575                           !  set to monthly if obc datafile hold 1 month of data
576   rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary
577   rn_dpwin    =    1.     !     -           -         -       west    -      -
578   rn_dpnin    =    1.     !     -           -         -       north   -      -
579   rn_dpsin    =    1.     !     -           -         -       south   -      -
580   rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary
581   rn_dpwob    =   15.     !     -           -         -     west    -      -
582   rn_dpnob    = 3000.     !     -           -         -     north   -      -
583   rn_dpsob    =   15.     !     -           -         -     south   -      -
584   rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R)
585                           !  = 1 the total volume remains constant
586/
587!-----------------------------------------------------------------------
588&namagrif      !  AGRIF zoom                                            ("key_agrif")
589!-----------------------------------------------------------------------
590   nn_cln_update =    3    !  baroclinic update frequency
591   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
592   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
593   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
594/
595!-----------------------------------------------------------------------
596&nam_tide      !   tide parameters (#ifdef key_tide)
597!-----------------------------------------------------------------------
598   ln_tide_pot   = .true.   !  use tidal potential forcing
599   ln_tide_ramp  = .false.  !
600   rdttideramp   =    0.    !
601   clname(1)     = 'DUMMY'  !  name of constituent - all tidal components must be set in namelist_cfg
602/
603!-----------------------------------------------------------------------
604&nambdy        !  unstructured open boundaries                          ("key_bdy")
605!-----------------------------------------------------------------------
606    nb_bdy         = 0                    !  number of open boundary sets
607    ln_coords_file = .true.               !  =T : read bdy coordinates from file
608    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files
609    ln_mask_file   = .false.              !  =T : read mask from file
610    cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.)
611    cn_dyn2d       = 'none'               !
612    nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state
613                                          !  = 1, bdy data are read in 'bdydata   .nc' files
614                                          !  = 2, use tidal harmonic forcing data from files
615                                          !  = 3, use external data AND tidal harmonic forcing
616    cn_dyn3d      =  'none'               !
617    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state
618                                          !  = 1, bdy data are read in 'bdydata   .nc' files
619    cn_tra        =  'none'               !
620    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state
621                                          !  = 1, bdy data are read in 'bdydata   .nc' files
622    cn_ice_lim      =  'none'             !
623    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state
624                                          !  = 1, bdy data are read in 'bdydata   .nc' files
625    rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice
626    rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           --
627    rn_ice_age      = 30.                 !  lim3 only:      --   age                --
628
629    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers
630    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities
631    rn_time_dmp   =  1.                   ! Damping time scale in days
632    rn_time_dmp_out =  1.                 ! Outflow damping time scale
633    nn_rimwidth   = 10                    !  width of the relaxation zone
634    ln_vol        = .false.               !  total volume correction (see nn_volctl parameter)
635    nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero
636/
637!-----------------------------------------------------------------------
638&nambdy_dta      !  open boundaries - external data           ("key_bdy")
639!-----------------------------------------------------------------------
640!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
641!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
642   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
643   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
644   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
645   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
646   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
647   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
648   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
649! for lim2
650!   bn_frld  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
651!   bn_hicif =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
652!   bn_hsnif =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
653! for lim3
654!   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
655!   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
656!   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
657   cn_dir  =    'bdydta/'
658   ln_full_vel = .false.
659/
660!-----------------------------------------------------------------------
661&nambdy_tide     ! tidal forcing at open boundaries
662!-----------------------------------------------------------------------
663   filtide          = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files
664   ln_bdytide_2ddta = .false.
665   ln_bdytide_conj  = .false.
666/
667!!======================================================================
668!!                 ***  Bottom boundary condition  ***
669!!======================================================================
670!!   nambfr        bottom friction
671!!   nambbc        bottom temperature boundary condition
672!!   nambbl        bottom boundary layer scheme                         ("key_trabbl")
673!!======================================================================
674!
675!-----------------------------------------------------------------------
676&nambfr        !   bottom friction
677!-----------------------------------------------------------------------
678   nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
679                           !                              = 2 : nonlinear friction
680   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case)
681   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
682   rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T)
683   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2)
684   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T
685   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file )
686   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T)
687   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case)
688   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
689   rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T)
690   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2)
691   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T
692   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file )
693   rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T)
694
695   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true)
696   ln_loglayer = .false.   !  logarithmic formulation (non linear case)
697/
698!-----------------------------------------------------------------------
699&nambbc        !   bottom temperature boundary condition
700!-----------------------------------------------------------------------
701!              !                              !  (if <0  months)  ! 
702!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
703!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
704   sn_qgh      ='geothermal_heating.nc',  -12.  , 'heatflow'      ,   .false.      , .true.  , 'yearly'  , ''       , ''       , ''
705   !
706   cn_dir      = './'      !  root directory for the location of the runoff files
707   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom
708   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux
709                           !     = 1 constant flux
710                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
711   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
712
713/
714!-----------------------------------------------------------------------
715&nambbl        !   bottom boundary layer scheme
716!-----------------------------------------------------------------------
717   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0)
718   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0)
719   rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s]
720   rn_gambbl   =  10.    !  advective bbl coefficient                 [s]
721/
722
723!!======================================================================
724!!                        Tracer (T & S ) namelists
725!!======================================================================
726!!   nameos        equation of state
727!!   namtra_adv    advection scheme
728!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.)
729!!   namtra_ldf    lateral diffusion scheme
730!!   namtra_dmp    T & S newtonian damping
731!!======================================================================
732!
733!-----------------------------------------------------------------------
734&nameos        !   ocean physical parameters
735!-----------------------------------------------------------------------
736   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency
737                                 !  =-1, TEOS-10
738                                 !  = 0, EOS-80
739                                 !  = 1, S-EOS   (simplified eos)
740   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm
741   !                             !
742   !                     ! S-EOS coefficients :
743   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
744   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1)
745   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1)
746   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
747   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
748   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
749   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
750   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
751/
752!-----------------------------------------------------------------------
753&namtra_adv    !   advection scheme for tracer
754!-----------------------------------------------------------------------
755   ln_traadv_cen2   =  .false.   !  2nd order centered scheme
756   ln_traadv_tvd    =  .true.    !  TVD scheme
757   ln_traadv_muscl  =  .false.   !  MUSCL scheme
758   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries
759   ln_traadv_ubs    =  .false.   !  UBS scheme
760   ln_traadv_qck    =  .false.   !  QUICKEST scheme
761   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl
762   ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection
763/
764!-----------------------------------------------------------------------
765&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param)
766!-----------------------------------------------------------------------
767   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation
768   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08)
769   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
770   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0)
771   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0)
772   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
773   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
774   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE
775   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK
776/
777!----------------------------------------------------------------------------------
778&namtra_ldf    !   lateral diffusion scheme for tracers
779!----------------------------------------------------------------------------------
780   !                       !  Operator type:
781   ln_traldf_lap    =  .true.   !  laplacian operator
782   ln_traldf_bilap  =  .false.  !  bilaplacian operator
783   !                       !  Direction of action:
784   ln_traldf_level  =  .false.  !  iso-level
785   ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T)
786   ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp")
787   !                 !  Griffies parameters              (all need "key_ldfslp")
788   ln_traldf_grif   =  .false.  !  use griffies triads
789   ln_traldf_gdia   =  .false.  !  output griffies eddy velocities
790   ln_triad_iso     =  .false.  !  pure lateral mixing in ML
791   ln_botmix_grif   =  .false.  !  lateral mixing on bottom
792   !                       !  Coefficients
793   ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv"
794   ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv
795   !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05)
796   rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s]
797   rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s]
798   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s]
799   !                                           (normally=0; not used with Griffies)
800   rn_slpmax        =     0.01  !  slope limit
801   rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity
802   rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer
803   rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s)
804/
805!-----------------------------------------------------------------------
806&namtra_dmp    !   tracer: T & S newtonian damping
807!-----------------------------------------------------------------------
808   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F)
809   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column
810                           !                   =1 no damping in the mixing layer (kz  criteria)
811                           !                   =2 no damping in the mixed  layer (rho crieria)
812   cn_resto    = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this)
813/
814
815!!======================================================================
816!!                      ***  Dynamics namelists  ***
817!!======================================================================
818!!   namdyn_adv    formulation of the momentum advection
819!!   namdyn_vor    advection scheme
820!!   namdyn_hpg    hydrostatic pressure gradient
821!!   namdyn_spg    surface pressure gradient                            (CPP key only)
822!!   namdyn_ldf    lateral diffusion scheme
823!!======================================================================
824!
825!-----------------------------------------------------------------------
826&namdyn_adv    !   formulation of the momentum advection
827!-----------------------------------------------------------------------
828   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)
829   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction
830   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
831   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
832   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection
833/
834!-----------------------------------------------------------------------
835&nam_vvl    !   vertical coordinate options
836!-----------------------------------------------------------------------
837   ln_vvl_zstar  = .true.           !  zstar vertical coordinate
838   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations
839   ln_vvl_layer  = .false.          !  full layer vertical coordinate
840   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
841   ln_vvl_zstar_at_eqtor = .false.  !  ztilde near the equator
842   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient
843   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days]
844   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days]
845   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation
846   ln_vvl_dbg    = .true.           !  debug prints    (T/F)
847/
848!-----------------------------------------------------------------------
849&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
850!-----------------------------------------------------------------------
851   ln_dynvor_ene = .false. !  enstrophy conserving scheme
852   ln_dynvor_ens = .false. !  energy conserving scheme
853   ln_dynvor_mix = .false. !  mixed scheme
854   ln_dynvor_een = .true.  !  energy & enstrophy scheme
855   ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation
856/
857!-----------------------------------------------------------------------
858&namdyn_hpg    !   Hydrostatic pressure gradient option
859!-----------------------------------------------------------------------
860   ln_hpg_zco  = .false.   !  z-coordinate - full steps
861   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation)
862   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
863   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
864   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
865   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
866   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T)
867                                 !           centered      time scheme  (F)
868/
869!-----------------------------------------------------------------------
870!namdyn_spg    !   surface pressure gradient   (CPP key only)
871!-----------------------------------------------------------------------
872!                          !  explicit free surface                     ("key_dynspg_exp")
873!                          !  filtered free surface                     ("key_dynspg_flt")
874!                          !  split-explicit free surface               ("key_dynspg_ts")
875
876!-----------------------------------------------------------------------
877&namdyn_ldf    !   lateral diffusion on momentum
878!-----------------------------------------------------------------------
879   !                       !  Type of the operator :
880   ln_dynldf_lap    =  .true.   !  laplacian operator
881   ln_dynldf_bilap  =  .false.  !  bilaplacian operator
882   !                       !  Direction of action  :
883   ln_dynldf_level  =  .false.  !  iso-level
884   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.)
885   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp")
886   !                       !  Coefficient
887   rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s]
888   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s]
889   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s]
890   rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity
891   rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity
892   rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity
893   rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp)
894   rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap)
895/
896
897!!======================================================================
898!!             Tracers & Dynamics vertical physics namelists
899!!======================================================================
900!!    namzdf        vertical physics
901!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric")
902!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke")
903!!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp")
904!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm")
905!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx")
906!!======================================================================
907!
908!-----------------------------------------------------------------------
909&namzdf        !   vertical physics
910!-----------------------------------------------------------------------
911   rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
912   rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
913   nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0)
914   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0)
915   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F)
916   nn_evdm     =    0      !  evd apply on tracer (=0) or on tracer and momentum (=1)
917   rn_avevd    =  100.     !  evd mixing coefficient [m2/s]
918   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F)
919   nn_npc      =    1            !  frequency of application of npc
920   nn_npcp     =  365            !  npc control print frequency
921   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
922   nn_zdfexp   =    3            !  number of sub-timestep for ln_zdfexp=T
923/
924!-----------------------------------------------------------------------
925&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
926!-----------------------------------------------------------------------
927   rn_avmri    = 100.e-4   !  maximum value of the vertical viscosity
928   rn_alp      =   5.      !  coefficient of the parameterization
929   nn_ric      =   2       !  coefficient of the parameterization
930   rn_ekmfc    =   0.7     !  Factor in the Ekman depth Equation
931   rn_mldmin   =   1.0     !  minimum allowable mixed-layer depth estimate (m)
932   rn_mldmax   =1000.0     !  maximum allowable mixed-layer depth estimate (m)
933   rn_wtmix    =  10.0     !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
934   rn_wvmix    =  10.0     !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
935   ln_mldw     = .true.    !  Flag to use or not the mized layer depth param.
936/
937!-----------------------------------------------------------------------
938&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
939!-----------------------------------------------------------------------
940   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
941   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
942   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
943   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
944   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
945   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
946   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
947                           !                 = 1 bounded by the local vertical scale factor
948                           !                 = 2 first vertical derivative of mixing length bounded by 1
949                           !                 = 3 as =2 with distinct disspipative an mixing length scale
950   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
951   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
952   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
953   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
954   rn_lc       =   0.15    !  coef. associated to Langmuir cells
955   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves
956                           !        = 0 no penetration
957                           !        = 1 add a tke source below the ML
958                           !        = 2 add a tke source just at the base of the ML
959                           !        = 3 as = 1 applied on HF part of the stress    ("key_oasis3")
960   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
961   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
962                           !        = 0  constant 10 m length scale
963                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees
964/
965!------------------------------------------------------------------------
966&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally:
967!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
968   ln_kpprimix = .true.    !  shear instability mixing
969   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s]
970   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s]
971   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability
972   rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s]
973   rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]
974   rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s]
975   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv
976   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt
977/
978!-----------------------------------------------------------------------
979&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls")
980!-----------------------------------------------------------------------
981   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
982   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
983   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
984   rn_clim_galp  = 0.267   !  galperin limit
985   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
986   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
987   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
988   rn_hsro       =  0.02   !  Minimum surface roughness
989   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2)
990   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2)
991   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
992   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
993   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
994   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
995/
996!-----------------------------------------------------------------------
997&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
998!-----------------------------------------------------------------------
999   rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
1000   rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
1001/
1002!-----------------------------------------------------------------------
1003&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
1004!-----------------------------------------------------------------------
1005   rn_htmx     = 500.      !  vertical decay scale for turbulence (meters)
1006   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1)
1007   rn_tfe      = 0.333     !  tidal dissipation efficiency
1008   rn_me       = 0.2       !  mixing efficiency
1009   ln_tmx_itf  = .true.    !  ITF specific parameterisation
1010   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency
1011/
1012
1013!!======================================================================
1014!!                  ***  Miscellaneous namelists  ***
1015!!======================================================================
1016!!   namsol            elliptic solver / island / free surface
1017!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi)
1018!!   namctl            Control prints & Benchmark
1019!!   namc1d            1D configuration options                         ("key_c1d")
1020!!   namc1d_uvd        data: U & V currents                             ("key_c1d")
1021!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d")
1022!!   namsto            Stochastic parametrization of EOS
1023!!======================================================================
1024!
1025!-----------------------------------------------------------------------
1026&namsol        !   elliptic solver / island / free surface
1027!-----------------------------------------------------------------------
1028   nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
1029                           !                   =2 successive-over-relaxation (sor)
1030   nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test
1031   rn_eps      =  1.e-6    !  absolute precision of the solver
1032   nn_nmin     =    300    !  minimum of iterations for the SOR solver
1033   nn_nmax     =    800    !  maximum of iterations for the SOR solver
1034   nn_nmod     =     10    !  frequency of test for the SOR solver
1035   rn_resmax   =  1.e-10   !  absolute precision for the SOR solver
1036   rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain)
1037/
1038!-----------------------------------------------------------------------
1039&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
1040!-----------------------------------------------------------------------
1041   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
1042                           !  buffer blocking send or immediate non-blocking sends, resp.
1043   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
1044   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold
1045   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1)
1046   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1)
1047   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1)
1048/
1049!-----------------------------------------------------------------------
1050&namctl        !   Control prints & Benchmark
1051!-----------------------------------------------------------------------
1052   ln_ctl      = .false.   !  trends control print (expensive!)
1053   nn_print    =    0      !  level of print (0 no extra print)
1054   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
1055   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
1056   nn_jctls    =    0      !  start j indice of control               over a subdomain)
1057   nn_jctle    =    0      !  end   j indice of control
1058   nn_isplt    =    1      !  number of processors in i-direction
1059   nn_jsplt    =    1      !  number of processors in j-direction
1060   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench
1061                           !     (no physical validity of the results)
1062   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0)
1063/
1064!-----------------------------------------------------------------------
1065&namc1d_uvd    !   data: U & V currents                                 ("key_c1d")
1066!-----------------------------------------------------------------------
1067!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1068!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
1069   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , ''
1070   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , ''
1071!
1072   cn_dir        = './'    !  root directory for the location of the files
1073   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F)
1074   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F)
1075/
1076!-----------------------------------------------------------------------
1077&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d")
1078!-----------------------------------------------------------------------
1079   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
1080/
1081!-----------------------------------------------------------------------
1082&namsto       ! Stochastic parametrization of EOS
1083!-----------------------------------------------------------------------
1084   ln_rststo = .false.           ! start from mean parameter (F) or from restart file (T)
1085   ln_rstseed = .true.           ! read seed of RNG from restart file
1086   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
1087   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)
1088
1089   ln_sto_eos = .false.          ! stochastic equation of state
1090   nn_sto_eos = 1                ! number of independent random walks
1091   rn_eos_stdxy = 1.4            ! random walk horz. standard deviation (in grid points)
1092   rn_eos_stdz  = 0.7            ! random walk vert. standard deviation (in grid points)
1093   rn_eos_tcor  = 1440.0         ! random walk time correlation (in timesteps)
1094   nn_eos_ord  = 1               ! order of autoregressive processes
1095   nn_eos_flt  = 0               ! passes of Laplacian filter
1096   rn_eos_lim  = 2.0             ! limitation factor (default = 3.0)
1097/
1098
1099!!======================================================================
1100!!                  ***  Diagnostics namelists  ***
1101!!======================================================================
1102!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1103!!   namtrd       dynamics and/or tracer trends
1104!!   namptr       Poleward Transport Diagnostics
1105!!   namflo       float parameters                                      ("key_float")
1106!!   namhsb       Heat and salt budgets
1107!!======================================================================
1108!
1109!-----------------------------------------------------------------------
1110&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
1111!-----------------------------------------------------------------------
1112   nn_nchunks_i=   4       !  number of chunks in i-dimension
1113   nn_nchunks_j=   4       !  number of chunks in j-dimension
1114   nn_nchunks_k=   31      !  number of chunks in k-dimension
1115                           !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1116                           !  is optimal for postprocessing which works exclusively with horizontal slabs
1117   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1118                           !  (F) ignore chunking information and produce netcdf3-compatible files
1119/
1120!-----------------------------------------------------------------------
1121&namtrd        !   diagnostics on dynamics and/or tracer trends
1122!              !       and/or mixed-layer trends and/or barotropic vorticity
1123!-----------------------------------------------------------------------
1124   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1125   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1126   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1127   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet)
1128   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1129   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1130   ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output
1131   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1132   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1133/
1134!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1135!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1136!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1137!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1138!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1139!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1140!!gm
1141!-----------------------------------------------------------------------
1142&namflo       !   float parameters                                      ("key_float")
1143!-----------------------------------------------------------------------
1144   jpnfl         = 1          !  total number of floats during the run
1145   jpnnewflo     = 0          !  number of floats for the restart
1146   ln_rstflo     = .false.    !  float restart (T) or not (F)
1147   nn_writefl    =      75    !  frequency of writing in float output file
1148   nn_stockfl    =    5475    !  frequency of creation of the float restart file
1149   ln_argo       = .false.    !  Argo type floats (stay at the surface each 10 days)
1150   ln_flork4     = .false.    !  trajectories computed with a 4th order Runge-Kutta (T)
1151                              !  or computed with Blanke' scheme (F)
1152   ln_ariane     = .true.     !  Input with Ariane tool convention(T)
1153   ln_flo_ascii  = .true.     !  Output with Ariane tool netcdf convention(F) or ascii file (T)
1154/
1155!-----------------------------------------------------------------------
1156&namptr       !   Poleward Transport Diagnostic
1157!-----------------------------------------------------------------------
1158   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F)
1159   ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not
1160/
1161!-----------------------------------------------------------------------
1162&namhsb       !  Heat and salt budgets
1163!-----------------------------------------------------------------------
1164   ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F)
1165/
1166!-----------------------------------------------------------------------
1167&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm')
1168!-----------------------------------------------------------------------
1169    nit000_han = 1         ! First time step used for harmonic analysis
1170    nitend_han = 75        ! Last time step used for harmonic analysis
1171    nstep_han  = 15        ! Time step frequency for harmonic analysis
1172    tname(1)   = 'M2'      ! Name of tidal constituents
1173    tname(2)   = 'K1'
1174/
1175!-----------------------------------------------------------------------
1176&namdct        ! transports through sections
1177!-----------------------------------------------------------------------
1178    nn_dct      = 15       !  time step frequency for transports computing
1179    nn_dctwri   = 15       !  time step frequency for transports writing
1180    nn_secdebug = 112      !      0 : no section to debug
1181                           !     -1 : debug all section
1182                           !  0 < n : debug section number n
1183/
1184
1185!!======================================================================
1186!!            ***  Observation & Assimilation namelists ***
1187!!======================================================================
1188!!   namobs       observation and model comparison                      ('key_diaobs')
1189!!   nam_asminc   assimilation increments                               ('key_asminc')
1190!!======================================================================
1191!
1192!-----------------------------------------------------------------------
1193&namobs       !  observation usage switch                               ('key_diaobs')
1194!-----------------------------------------------------------------------
1195   ln_diaobs  = .false.             ! Logical switch for the observation operator
1196   ln_t3d     = .false.             ! Logical switch for T profile observations
1197   ln_s3d     = .false.             ! Logical switch for S profile observations
1198   ln_sla     = .false.             ! Logical switch for SLA observations
1199   ln_sst     = .false.             ! Logical switch for SST observations
1200   ln_sic     = .false.             ! Logical switch for Sea Ice observations
1201   ln_vel3d   = .false.             ! Logical switch for velocity observations
1202   ln_sss     = .false.             ! Logical swithc for SSS observations
1203   ln_slchltot = .false.            ! Logical switch for surface total              log10(chlorophyll) obs
1204   ln_slchldia = .false.            ! Logical switch for surface diatom             log10(chlorophyll) obs
1205   ln_slchlnon = .false.            ! Logical switch for surface non-diatom         log10(chlorophyll) obs
1206   ln_slchldin = .false.            ! Logical switch for surface dinoflagellate     log10(chlorophyll) obs
1207   ln_slchlmic = .false.            ! Logical switch for surface microphytoplankton log10(chlorophyll) obs
1208   ln_slchlnan = .false.            ! Logical switch for surface nanophytoplankton  log10(chlorophyll) obs
1209   ln_slchlpic = .false.            ! Logical switch for surface picophytoplankton  log10(chlorophyll) obs
1210   ln_schltot  = .false.            ! Logical switch for surface total              chlorophyll        obs
1211   ln_slphytot = .false.            ! Logical switch for surface total      log10(phytoplankton carbon) obs
1212   ln_slphydia = .false.            ! Logical switch for surface diatom     log10(phytoplankton carbon) obs
1213   ln_slphynon = .false.            ! Logical switch for surface non-diatom log10(phytoplankton carbon) obs
1214   ln_sspm     = .false.            ! Logical switch for surface suspended particulate matter obs
1215   ln_skd490   = .false.            ! Logical switch for surface attenuation coefficient of downwelling radiation at 490 nm
1216   ln_sfco2    = .false.            ! Logical switch for surface fugacity         of carbon dioxide obs
1217   ln_spco2    = .false.            ! Logical switch for surface partial pressure of carbon dioxide obs
1218   ln_plchltot = .false.            ! Logical switch for profile total log10(chlorophyll) obs
1219   ln_pchltot  = .false.            ! Logical switch for profile total chlorophyll obs
1220   ln_pno3     = .false.            ! Logical switch for profile nitrate obs
1221   ln_psi4     = .false.            ! Logical switch for profile silicate obs
1222   ln_ppo4     = .false.            ! Logical switch for profile phosphate obs
1223   ln_pdic     = .false.            ! Logical switch for profile dissolved inorganic carbon obs
1224   ln_palk     = .false.            ! Logical switch for profile alkalinity obs
1225   ln_pph      = .false.            ! Logical switch for profile pH obs
1226   ln_po2      = .false.            ! Logical switch for profile dissolved oxygen obs
1227   ln_altbias = .false.             ! Logical switch for altimeter bias correction
1228   ln_sstbias = .false.             ! Logical switch for SST bias correction
1229   ln_nea     = .false.             ! Logical switch for rejection of observations near land
1230   ln_grid_global = .true.          ! Logical switch for global distribution of observations
1231   ln_grid_search_lookup = .false.  ! Logical switch for obs grid search w/lookup table
1232   ln_ignmis  = .true.              ! Logical switch for ignoring missing files
1233   ln_s_at_t  = .false.             ! Logical switch for computing model S at T obs if not there
1234   ln_sstnight = .false.            ! Logical switch for calculating night-time average for SST obs
1235   ln_output_clim = .false.         ! Logical switch for writing climatological values to fdbk files
1236   ln_time_mean_sla_bkg = .false.   ! Logical switch for applying time mean of SLA background to remove tidal signal
1237   ln_default_fp_indegs = .true.    ! Logical: T=> averaging footprint is in degrees, F=> in metres
1238   ln_sla_fp_indegs = .true.        ! Logical: T=> averaging footprint is in degrees, F=> in metres
1239   ln_sst_fp_indegs = .true.
1240   ln_sss_fp_indegs = .true.
1241   ln_sic_fp_indegs = .true.
1242! All of the *files* variables below are arrays. Use namelist_cfg to add more files
1243   cn_profbfiles = 'profiles_01.nc'      ! Profile feedback input observation file names
1244   cn_slafbfiles = 'sla_01.nc'           ! SLA feedback input observation file names
1245   cn_sstfbfiles = 'sst_01.nc'           ! SST feedback input observation file names
1246   cn_sicfbfiles = 'sic_01.nc'           ! SIC feedback input observation file names
1247   cn_velfbfiles = 'vel_01.nc'           ! Velocity feedback input observation file names
1248   cn_sssfbfiles = 'sss_01.nc'           ! SSS feedback input observation file names
1249   cn_slchltotfbfiles = 'slchltot_01.nc' ! Surface total              log10(chlorophyll) obs file names
1250   cn_slchldiafbfiles = 'slchldia_01.nc' ! Surface diatom             log10(chlorophyll) obs file names
1251   cn_slchlnonfbfiles = 'slchlnon_01.nc' ! Surface non-diatom         log10(chlorophyll) obs file names
1252   cn_slchldinfbfiles = 'slchldin_01.nc' ! Surface dinoflagellate     log10(chlorophyll) obs file names
1253   cn_slchlmicfbfiles = 'slchlmic_01.nc' ! Surface microphytoplankton log10(chlorophyll) obs file names
1254   cn_slchlnanfbfiles = 'slchlnan_01.nc' ! Surface nanophytoplankton  log10(chlorophyll) obs file names
1255   cn_slchlpicfbfiles = 'slchlpic_01.nc' ! Surface picophytoplankton  log10(chlorophyll) obs file names
1256   cn_schltotfbfiles  = 'schltot_01.nc'  ! Surface total              chlorophyll        obs file names
1257   cn_slphytotfbfiles = 'slphytot_01.nc' ! Surface total      log10(phytoplankton carbon) obs file names
1258   cn_slphydiafbfiles = 'slphydia_01.nc' ! Surface diatom     log10(phytoplankton carbon) obs file names
1259   cn_slphynonfbfiles = 'slphynon_01.nc' ! Surface non-diatom log10(phytoplankton carbon) obs file names
1260   cn_sspmfbfiles     = 'sspm_01.nc'     ! Surface suspended particulate matter obs file names
1261   cn_skd490fbfiles   = 'skd490_01.nc'   ! Surface Kd490 obs file names
1262   cn_sfco2fbfiles    = 'sfco2_01.nc'    ! Surface fugacity         of carbon dioxide obs file names
1263   cn_spco2fbfiles    = 'spco2_01.nc'    ! Surface partial pressure of carbon dioxide obs file names
1264   cn_plchltotfbfiles = 'plchltot_01.nc' ! Profile total log10(chlorophyll) obs file names
1265   cn_pchltotfbfiles  = 'pchltot_01.nc'  ! Profile total chlorophyll obs file names
1266   cn_pno3fbfiles     = 'pno3_01.nc'     ! Profile nitrate obs file names
1267   cn_psi4fbfiles     = 'psi4_01.nc'     ! Profile silicate obs file names
1268   cn_ppo4fbfiles     = 'ppo4_01.nc'     ! Profile phosphate obs file names
1269   cn_pdicfbfiles     = 'pdic_01.nc'     ! Profile dissolved inorganic carbon obs file names
1270   cn_palkfbfiles     = 'palk_01.nc'     ! Profile alkalinity obs file names
1271   cn_pphfbfiles      = 'pph_01.nc'      ! Profile pH obs file names
1272   cn_po2fbfiles      = 'po2_01.nc'      ! Profile dissolved oxygen obs file names
1273   cn_altbiasfile = 'altbias.nc'         ! Altimeter bias input file name
1274   cn_sstbiasfiles = 'sstbias.nc'        ! SST bias input file names
1275   cn_gridsearchfile = 'gridsearch.nc'   ! Grid search file name
1276   rn_gridsearchres = 0.5                ! Grid search resolution
1277   rn_default_avglamscl = 0.             ! Default E/W diameter of observation footprint (metres/degrees)
1278   rn_default_avgphiscl = 0.             ! Default N/S diameter of observation footprint (metres/degrees)
1279   rn_sla_avglamscl = 0.                 ! E/W diameter of SLA observation footprint (metres/degrees)
1280   rn_sla_avgphiscl = 0.                 ! N/S diameter of SLA observation footprint (metres/degrees)
1281   rn_sst_avglamscl = 0.                 ! E/W diameter of SST observation footprint (metres/degrees)
1282   rn_sst_avgphiscl = 0.                 ! N/S diameter of SST observation footprint (metres/degrees)
1283   rn_sss_avglamscl = 0.                 ! E/W diameter of SSS observation footprint (metres/degrees)
1284   rn_sss_avgphiscl = 0.                 ! N/S diameter of SSS observation footprint (metres/degrees)
1285   rn_sic_avglamscl = 0.                 ! E/W diameter of SIC observation footprint (metres/degrees)
1286   rn_sic_avgphiscl = 0.                 ! N/S diameter of SIC observation footprint (metres/degrees)
1287   nn_1dint = 0                          ! Type of vertical interpolation method
1288   nn_2dint_default = 0                  ! Default horizontal interpolation method
1289   nn_2dint_sla = -1                     ! Horizontal interpolation method for SLA
1290   nn_2dint_sst = -1                     ! Horizontal interpolation method for SST
1291   nn_2dint_sss = -1                     ! Horizontal interpolation method for SSS
1292   nn_2dint_sic = -1                     ! Horizontal interpolation method for SIC
1293   nn_msshc = 0                          ! MSSH correction scheme
1294   rn_mdtcorr = 1.61                     ! MDT  correction
1295   rn_mdtcutoff = 65.0                   ! MDT cutoff for computed correction
1296   nn_profdavtypes = -1                  ! Profile daily average types - array
1297   
1298/
1299!-----------------------------------------------------------------------
1300&nam_asminc   !   assimilation increments                               ('key_asminc')
1301!-----------------------------------------------------------------------
1302    ln_bkgwri = .false.    !  Logical switch for writing out background state
1303    ln_trainc = .false.    !  Logical switch for applying tracer increments
1304    ln_dyninc = .false.    !  Logical switch for applying velocity increments
1305    ln_sshinc = .false.    !  Logical switch for applying SSH increments
1306    ln_asmdin = .false.    !  Logical switch for Direct Initialization (DI)
1307    ln_asmiau = .false.    !  Logical switch for Incremental Analysis Updating (IAU)
1308    nitbkg    = 0          !  Timestep of background in [0,nitend-nit000-1]
1309    nitdin    = 0          !  Timestep of background for DI in [0,nitend-nit000-1]
1310    nitiaustr = 1          !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1311    nitiaufin = 15         !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1312    niaufn    = 0          !  Type of IAU weighting function
1313    ln_salfix = .false.    !  Logical switch for ensuring that the sa > salfixmin
1314    salfixmin = -9999      !  Minimum salinity after applying the increments
1315    nn_divdmp = 0          !  Number of iterations of divergence damping operator
1316/
1317!-----------------------------------------------------------------------
1318&namsbc_wave   ! External fields from wave model
1319!-----------------------------------------------------------------------
1320!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1321!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
1322   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1323   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1324   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1325   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1326!
1327   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files
1328/
1329!-----------------------------------------------------------------------
1330&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
1331!-----------------------------------------------------------------------
1332   ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model
1333   ln_neptsimp       = .false.  ! yes/no use simplified neptune
1334
1335   ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep
1336   rn_tslse          =  1.2e4   ! value of lengthscale L at the equator
1337   rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole
1338   ! Specify whether to ramp down the Neptune velocity in shallow
1339   ! water, and if so the depth range controlling such ramping down
1340   ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water
1341   rn_htrmin         =  100.0   ! min. depth of transition range
1342   rn_htrmax         =  200.0   ! max. depth of transition range
1343/
Note: See TracBrowser for help on using the repository browser.