source: branches/UKMO/AMM15_v3_6_STABLE_package_collate/NEMOGCM/CONFIG/SHARED/namelist_ref @ 10268

Last change on this file since 10268 was 10268, checked in by kingr, 21 months ago

Updating ref namelist namobs to match OBS branch. Also added code to apply ssh increment in dynspg_ts which didnt make it into merged branch.

File size: 93.9 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  :  1 - run manager      (namrun)
3!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd)
4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas
5!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,
6!!                                    namsbc_apr, namsbc_ssr, namsbc_alb)
7!!              4 - lateral boundary (namlbc, namcla, 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_t3d     = .false.    ! Logical switch for T profile observations
1196   ln_s3d     = .false.    ! Logical switch for S profile observations
1197   ln_ena     = .false.    ! Logical switch for ENACT insitu data set
1198   !                       !     ln_cor                  Logical switch for Coriolis insitu data set
1199   ln_profb   = .false.    ! Logical switch for feedback insitu data set
1200   ln_sla     = .false.    ! Logical switch for SLA observations
1201
1202   ln_sladt   = .false.    ! Logical switch for AVISO SLA data
1203
1204   ln_slafb   = .false.    ! Logical switch for feedback SLA data
1205                           !     ln_ssh                  Logical switch for SSH observations
1206
1207   ln_sst     = .false.     ! Logical switch for SST observations
1208   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations
1209   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations
1210
1211   ln_sstfb   = .false.    ! Logical switch for feedback SST data
1212                           !     ln_sss                  Logical switch for SSS observations
1213   ln_seaice  = .false.    ! Logical switch for Sea Ice observations
1214                           !     ln_vel3d                Logical switch for velocity observations
1215                           !     ln_velavcur             Logical switch for velocity daily av. cur.
1216                           !     ln_velhrcur             Logical switch for velocity high freq. cur.
1217                           !     ln_velavadcp            Logical switch for velocity daily av. ADCP
1218                           !     ln_velhradcp            Logical switch for velocity high freq. ADCP
1219                           !     ln_velfb                Logical switch for feedback velocity data
1220                           !     ln_grid_global          Global distribtion of observations
1221                           !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table
1222                           !     grid_search_file        Grid search lookup file header
1223                           !     enactfiles              ENACT input observation file names
1224                           !     coriofiles              Coriolis input observation file name
1225   !                       ! profbfiles: Profile feedback input observation file name
1226   profbfiles = 'profiles_01.nc'
1227                           !     ln_profb_enatim         Enact feedback input time setting switch
1228                           !     slafilesact             Active SLA input observation file name
1229                           !     slafilespas             Passive SLA input observation file name
1230   !                       ! slafbfiles: Feedback SLA input observation file name
1231   slafbfiles = 'sla_01.nc'
1232                           !     sstfiles                GHRSST input observation file name
1233   !                       ! sstfbfiles: Feedback SST input observation file name
1234   sstfbfiles = 'sst_01.nc'
1235                           !     seaicefiles             Sea Ice input observation file names
1236   seaicefiles = 'seaice_01.nc'
1237                           !     velavcurfiles           Vel. cur. daily av. input file name
1238                           !     velhvcurfiles           Vel. cur. high freq. input file name
1239                           !     velavadcpfiles          Vel. ADCP daily av. input file name
1240                           !     velhvadcpfiles          Vel. ADCP high freq. input file name
1241                           !     velfbfiles              Vel. feedback input observation file name
1242                           !     dobsini                 Initial date in window YYYYMMDD.HHMMSS
1243                           !     dobsend                 Final date in window YYYYMMDD.HHMMSS
1244                           !     n1dint                  Type of vertical interpolation method
1245                           !     n2dint                  Type of horizontal interpolation method
1246                           !     ln_nea                  Rejection of observations near land switch
1247   nmsshc     = 0          ! MSSH correction scheme
1248                           !     mdtcorr                 MDT  correction
1249                           !     mdtcutoff               MDT cutoff for computed correction
1250   ln_altbias = .false.    ! Logical switch for alt bias
1251   ln_ignmis  = .true.     ! Logical switch for ignoring missing files
1252                           !     endailyavtypes   ENACT daily average types
1253   ln_grid_global = .true.
1254   ln_grid_search_lookup = .false.
1255/
1256!-----------------------------------------------------------------------
1257&nam_asminc   !   assimilation increments                               ('key_asminc')
1258!-----------------------------------------------------------------------
1259    ln_bkgwri = .false.    !  Logical switch for writing out background state
1260    ln_trainc = .false.    !  Logical switch for applying tracer increments
1261    ln_dyninc = .false.    !  Logical switch for applying velocity increments
1262    ln_sshinc = .false.    !  Logical switch for applying SSH increments
1263    ln_asmdin = .false.    !  Logical switch for Direct Initialization (DI)
1264    ln_asmiau = .false.    !  Logical switch for Incremental Analysis Updating (IAU)
1265    nitbkg    = 0          !  Timestep of background in [0,nitend-nit000-1]
1266    nitdin    = 0          !  Timestep of background for DI in [0,nitend-nit000-1]
1267    nitiaustr = 1          !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1268    nitiaufin = 15         !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1269    niaufn    = 0          !  Type of IAU weighting function
1270    ln_salfix = .false.    !  Logical switch for ensuring that the sa > salfixmin
1271    salfixmin = -9999      !  Minimum salinity after applying the increments
1272    nn_divdmp = 0          !  Number of iterations of divergence damping operator
1273/
1274!-----------------------------------------------------------------------
1275&namsbc_wave   ! External fields from wave model
1276!-----------------------------------------------------------------------
1277!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1278!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
1279   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1280   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1281   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1282   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1283!
1284   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files
1285/
1286!-----------------------------------------------------------------------
1287&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
1288!-----------------------------------------------------------------------
1289   ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model
1290   ln_neptsimp       = .false.  ! yes/no use simplified neptune
1291
1292   ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep
1293   rn_tslse          =  1.2e4   ! value of lengthscale L at the equator
1294   rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole
1295   ! Specify whether to ramp down the Neptune velocity in shallow
1296   ! water, and if so the depth range controlling such ramping down
1297   ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water
1298   rn_htrmin         =  100.0   ! min. depth of transition range
1299   rn_htrmax         =  200.0   ! max. depth of transition range
1300/
Note: See TracBrowser for help on using the repository browser.