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

Last change on this file since 11277 was 11277, checked in by kingr, 17 months ago

Merged Juan's changes for running AMM15 woth wave coupling.
Corrected minor logic error to allow AMM7-uncoupled to reproduce earlier results.
Few line spacing changes to allow merging with OBS br and trunk rvn 5518.

File size: 100.5 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  :  1 - run manager      (namrun)
3!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd)
4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas
5!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,
6!!                                    namsbc_apr, namsbc_ssr, namsbc_alb)
7!!              4 - lateral boundary (namlbc, namcla, 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.   !  coupled   formulation                       ( requires key_oasis3 )
243   ln_mixcpl   = .false.   !  forced-coupled mixed atmosphere formulation ( requires key_oasis3 )
244   ln_wavcpl   = .false.   !  forced-coupled mixed wave       formulation ( requires key_oasis3 )
245   nn_components = 0       !  configuration of the opa-sas OASIS coupling
246                           !  =0 no opa-sas OASIS coupling: default single executable configuration
247                           !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component
248                           !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component
249   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr )
250   nn_ice      = 2         !  =0 no ice boundary condition   ,
251                           !  =1 use observed ice-cover      ,
252                           !  =2 ice-model used                         ("key_lim3" or "key_lim2")
253   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect)
254                           !  =1 levitating ice with mass and salt exchange but no presure effect
255                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure)
256   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave
257   ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf)
258   nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf)
259                           !  0 =no isf                  1 = presence of ISF
260                           !  2 = bg03 parametrisation   3 = rnf file for isf
261                           !  4 = ISF fwf specified
262                           !  option 1 and 4 need ln_isfcav = .true. (domzgr)
263   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
264   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked
265                           !     =1 global mean of e-p-r set to zero at each time step
266                           !     =2 annual global mean of e-p-r set to zero
267   ln_wave = .false.       !  Activate coupling with wave (T => fill namsbc_wave)   
268   nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
269                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
270   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used)
271                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled
272                           !  = 0  Average per-category fluxes (forced and coupled mode)
273                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled
274                           !  = 2  Redistribute a single flux over categories (coupled mode only)
275   nn_drag   = 0      !  formula to calculate momentum from the wind components
276                           ! = 0 UKMO SHELF formulation
277                           ! = 1 standard formulation with forced of coupled drag coefficient
278                           ! = 2 standard formulation with constant drag coefficient
279                           ! = 3 momentum calculated from core forcing fields
280/
281!-----------------------------------------------------------------------
282&namsbc_ana    !   analytical surface boundary condition
283!-----------------------------------------------------------------------
284   nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps
285   rn_utau0    =   0.5     !  uniform value for the i-stress
286   rn_vtau0    =   0.e0    !  uniform value for the j-stress
287   rn_qns0     =   0.e0    !  uniform value for the total heat flux
288   rn_qsr0     =   0.e0    !  uniform value for the solar radiation
289   rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P)
290/
291!-----------------------------------------------------------------------
292&namsbc_flx    !   surface boundary condition : flux formulation
293!-----------------------------------------------------------------------
294!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
295!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
296   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
297   sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
298   sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
299   sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
300   sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
301
302   cn_dir       = './'      !  root directory for the location of the flux files
303   ln_shelf_flx = .false.   !  UKMO SHELF specific flux flag - read from file wind components instead of momentum 
304   ln_rel_wind  = .false.   !  UKMO SHELF - calculate momentum from wind speed relative to currents 
305   rn_wfac      = 1.0       !  UKMO SHELF - multiplicative factor for ocean/wind velocity
306/
307!-----------------------------------------------------------------------
308&namsbc_clio   !   namsbc_clio  CLIO bulk formulae
309!-----------------------------------------------------------------------
310!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
311!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
312   sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
313   sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
314   sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
315   sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
316   sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
317   sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
318   sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
319
320   cn_dir      = './'      !  root directory for the location of the bulk files are
321/
322!-----------------------------------------------------------------------
323&namsbc_core   !   namsbc_core  CORE bulk formulae
324!-----------------------------------------------------------------------
325!              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask !
326!              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      !
327   sn_wndi     = 'u_10.15JUNE2009_fill'        ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , ''
328   sn_wndj     = 'v_10.15JUNE2009_fill'        ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , ''
329   sn_qsr      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
330   sn_qlw      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
331   sn_tair     = 't_10.15JUNE2009_fill'        ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
332   sn_humi     = 'q_10.15JUNE2009_fill'        ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
333   sn_prec     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
334   sn_snow     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
335   sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
336
337   cn_dir      = './'      !  root directory for the location of the bulk files
338   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data
339   rn_zqt      = 10.        !  Air temperature and humidity reference height (m)
340   rn_zu       = 10.        !  Wind vector reference height (m)
341   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow)
342   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.)
343   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity
344                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)
345/
346!-----------------------------------------------------------------------
347&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae
348!-----------------------------------------------------------------------
349!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask !
350!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      !
351   sn_wndi     =   'ecmwf'   ,        6          , 'u10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
352   sn_wndj     =   'ecmwf'   ,        6          , 'v10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
353   sn_clc      =   'ecmwf'   ,        6          , 'clc'     ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , ''
354   sn_msl      =   'ecmwf'   ,        6          , 'msl'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
355   sn_tair     =   'ecmwf'   ,        6          , 't2'      ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , ''
356   sn_rhm      =   'ecmwf'   ,        6          , 'rh'      ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , ''
357   sn_prec     =   'ecmwf'   ,        6          , 'precip'  ,    .true.    , .true.  , 'daily'  ,'bicubic.nc' , ''       , ''
358
359   cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files
360/
361!-----------------------------------------------------------------------
362&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
363!-----------------------------------------------------------------------
364!                    !     description       !  multiple  !    vector   !      vector          ! vector !
365!                    !                       ! categories !  reference  !    orientation       ! grids  !
366! send
367   sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   ''
368   sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
369   sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   ''
370   sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T'
371   sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
372   sn_snd_crtw   =       'none'                 ,    'no'    ,     ''      ,         ''           , 'U,V' 
373   sn_snd_ifrac  =       'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
374   sn_snd_wlev   =       'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
375! receive
376   sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
377   sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
378   sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V'
379   sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
380   sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
381   sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
382   sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   ''
383   sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
384   sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
385   sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
386   sn_rcv_hsig   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
387   sn_rcv_iceflx =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
388   sn_rcv_mslp   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
389   sn_rcv_phioc  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
390   sn_rcv_sdrft  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
391   sn_rcv_wper   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
392   sn_rcv_wfreq  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
393   sn_rcv_wnum   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
394   sn_rcv_tauoc  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
395   sn_rcv_tauw   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
396   sn_rcv_wdrag  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
397!
398   nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data
399   ln_usecplmask = .false.   !  use a coupling mask file to merge data received from several models
400                             !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
401   ln_usernfmask = .false.   !  use a runoff mask file to merge data received from several models
402                             !   -> file rnfmask.nc with the float variable called rnfmask (jpi,jpj,nn_cplmodel)
403/
404!-----------------------------------------------------------------------
405&namsbc_sas    !   analytical surface boundary condition
406!-----------------------------------------------------------------------
407!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
408!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
409   sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
410   sn_vsp      = 'sas_grid_V' ,    120           , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
411   sn_tem      = 'sas_grid_T' ,    120           , 'sosstsst' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
412   sn_sal      = 'sas_grid_T' ,    120           , 'sosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
413   sn_ssh      = 'sas_grid_T' ,    120           , 'sossheig' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
414   sn_e3t      = 'sas_grid_T' ,    120           , 'e3t_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
415   sn_frq      = 'sas_grid_T' ,    120           , 'frq_m'    ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''
416
417   ln_3d_uve   = .true.    !  specify whether we are supplying a 3D u,v and e3 field
418   ln_read_frq = .false.    !  specify whether we must read frq or not
419   cn_dir      = './'      !  root directory for the location of the bulk files are
420/
421!-----------------------------------------------------------------------
422&namtra_qsr    !   penetrative solar radiation
423!-----------------------------------------------------------------------
424!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
425!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
426   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , ''
427
428   cn_dir      = './'      !  root directory for the location of the runoff files
429   ln_traqsr   = .true.    !  Light penetration (T) or not (F)
430   ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration
431   ln_qsr_2bd  = .false.   !  2 bands              light penetration
432   ln_qsr_bio  = .false.   !  bio-model light penetration
433   nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0)
434   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
435   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
436   rn_si1      =   23.0    !  2 bands: longest depth of extinction
437   ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3
438/
439!-----------------------------------------------------------------------
440&namsbc_rnf    !   runoffs namelist surface boundary condition
441!-----------------------------------------------------------------------
442!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
443!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
444   sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
445   sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
446   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
447   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
448   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , ''
449
450   cn_dir       = './'      !  root directory for the location of the runoff files
451   ln_rnf_mouth = .true.    !  specific treatment at rivers mouths
452   rn_hrnf      =  15.e0    !  depth over which enhanced vertical mixing is used
453   rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s]
454   rn_rfact     =   1.e0    !  multiplicative factor for runoff
455   ln_rnf_depth = .false.   !  read in depth information for runoff
456   ln_rnf_tem   = .false.   !  read in temperature information for runoff
457   ln_rnf_sal   = .false.   !  read in salinity information for runoff
458   ln_rnf_depth_ini = .false.  ! compute depth at initialisation from runoff file
459   rn_rnf_max   = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
460   rn_dep_max   = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
461   nn_rnf_depth_file = 0    !  create (=1) a runoff depth file or not (=0)
462/
463!-----------------------------------------------------------------------
464&namsbc_isf    !  Top boundary layer (ISF)
465!-----------------------------------------------------------------------
466!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
467!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
468! nn_isf == 4
469   sn_qisf      = 'rnfisf' ,         -12      ,'sohflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
470   sn_fwfisf    = 'rnfisf' ,         -12      ,'sowflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
471! nn_isf == 3
472   sn_rnfisf    = 'runoffs' ,         -12      ,'sofwfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   ''
473! nn_isf == 2 and 3
474   sn_depmax_isf = 'runoffs' ,       -12        ,'sozisfmax' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
475   sn_depmin_isf = 'runoffs' ,       -12        ,'sozisfmin' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
476! nn_isf == 2
477   sn_Leff_isf = 'rnfisf' ,       0          ,'Leff'         ,   .false.  , .true.  , 'yearly'  ,  ''      ,   ''
478! for all case
479   ln_divisf   = .true.  ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?)
480! only for nn_isf = 1 or 2
481   rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula
482   rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula
483! only for nn_isf = 1
484   nn_isfblk   =  1       ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??)
485   rn_hisf_tbl =  30.      ! thickness of the top boundary layer           (Losh et al. 2008)
486                          ! 0 => thickness of the tbl = thickness of the first wet cell
487   ln_conserve = .true.   ! conservative case (take into account meltwater advection)
488   nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s)
489                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010)
490                          !     if you want to keep the cd as in global config, adjust rn_gammat0 to compensate
491                          ! 2 = velocity and stability dependent Gamma    Holland et al. 1999
492/
493!-----------------------------------------------------------------------
494&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk
495!-----------------------------------------------------------------------
496!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
497!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
498   sn_apr      = 'patm'      ,         -1        ,'somslpre',    .true.     , .true. , 'yearly'  ,  ''      ,   ''     , ''
499
500   cn_dir      = './'       !  root directory for the location of the bulk files
501   rn_pref     = 101000.    !  reference atmospheric pressure   [N/m2]/
502   ln_ref_apr  = .false.    !  ref. pressure: global mean Patm (T) or a constant (F)
503   ln_apr_obc  = .false.    !  inverse barometer added to OBC ssh data
504/
505!-----------------------------------------------------------------------
506&namsbc_ssr    !   surface boundary condition : sea surface restoring
507!-----------------------------------------------------------------------
508!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
509!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
510   sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , ''       , ''
511   sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , ''       , ''
512
513   cn_dir      = './'      !  root directory for the location of the runoff files
514   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
515   nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2)
516                           !  or to SSS only (=1) or no damping term (=0)
517   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K]
518   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
519   ln_sssr_bnd =   .true.  !  flag to bound erp term (associated with nn_sssr=2)
520   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
521/
522!-----------------------------------------------------------------------
523&namsbc_alb    !   albedo parameters
524!-----------------------------------------------------------------------
525   rn_cloud    =    0.06   !  cloud correction to snow and ice albedo
526   rn_albice   =    0.53   !  albedo of melting ice in the arctic and antarctic
527   rn_alphd    =    0.80   !  coefficients for linear interpolation used to
528   rn_alphc    =    0.65   !  compute albedo between two extremes values
529   rn_alphdi   =    0.72   !  (Pyane, 1972)
530/
531!-----------------------------------------------------------------------
532&namberg       !   iceberg parameters
533!-----------------------------------------------------------------------
534      ln_icebergs              = .false.
535      ln_bergdia               = .true.               ! Calculate budgets
536      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0
537      nn_verbose_write         = 15                   ! Timesteps between verbose messages
538      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage
539                                                      ! Initial mass required for an iceberg of each class
540      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
541                                                      ! Proportion of calving mass to apportion to each class
542      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
543                                                      ! Ratio between effective and real iceberg mass (non-dim)
544                                                      ! i.e. number of icebergs represented at a point
545      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1
546                                                      ! thickness of newly calved bergs (m)
547      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
548      rn_rho_bergs             = 850.                 ! Density of icebergs
549      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs
550      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics
551      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits
552      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
553      ln_passive_mode          = .false.              ! iceberg - ocean decoupling
554      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no)
555                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
556      rn_test_box              = 108.0,  116.0, -66.0, -58.0
557      rn_speed_limit           = 0.                   ! CFL speed limit for a berg
558
559!              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
560!              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
561      sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , ''
562
563      cn_dir = './'
564/
565
566!!======================================================================
567!!               ***  Lateral boundary condition  ***
568!!======================================================================
569!!   namlbc        lateral momentum boundary condition
570!!   namcla        cross land advection
571!!   namobc        open boundaries parameters                           ("key_obc")
572!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")
573!!   nambdy        Unstructured open boundaries                         ("key_bdy")
574!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides")
575!!======================================================================
576!
577!-----------------------------------------------------------------------
578&namlbc        !   lateral momentum boundary condition
579!-----------------------------------------------------------------------
580   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
581                           !  free slip  !   partial slip  !   no slip   ! strong slip
582   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs.
583/
584!-----------------------------------------------------------------------
585&namcla        !   cross land advection
586!-----------------------------------------------------------------------
587   nn_cla      =    0      !  advection between 2 ocean pts separates by land
588/
589!-----------------------------------------------------------------------
590&namobc        !   open boundaries parameters                           ("key_obc")
591!-----------------------------------------------------------------------
592   ln_obc_clim = .false.   !  climatological obc data files (T) or not (F)
593   ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F)
594   ln_obc_fla  = .false.   !  Flather open boundary condition
595   nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state
596                           !  = 1 the obc data are read in 'obc.dta' files
597   cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data
598                           !  set to monthly if obc datafile hold 1 month of data
599   rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary
600   rn_dpwin    =    1.     !     -           -         -       west    -      -
601   rn_dpnin    =    1.     !     -           -         -       north   -      -
602   rn_dpsin    =    1.     !     -           -         -       south   -      -
603   rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary
604   rn_dpwob    =   15.     !     -           -         -     west    -      -
605   rn_dpnob    = 3000.     !     -           -         -     north   -      -
606   rn_dpsob    =   15.     !     -           -         -     south   -      -
607   rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R)
608                           !  = 1 the total volume remains constant
609/
610!-----------------------------------------------------------------------
611&namagrif      !  AGRIF zoom                                            ("key_agrif")
612!-----------------------------------------------------------------------
613   nn_cln_update =    3    !  baroclinic update frequency
614   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
615   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
616   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
617/
618!-----------------------------------------------------------------------
619&nam_tide      !   tide parameters (#ifdef key_tide)
620!-----------------------------------------------------------------------
621   ln_tide_pot   = .true.   !  use tidal potential forcing
622   ln_tide_ramp  = .false.  !
623   rdttideramp   =    0.    !
624   clname(1)     = 'DUMMY'  !  name of constituent - all tidal components must be set in namelist_cfg
625/
626!-----------------------------------------------------------------------
627&nambdy        !  unstructured open boundaries                          ("key_bdy")
628!-----------------------------------------------------------------------
629    nb_bdy         = 0                    !  number of open boundary sets
630    ln_coords_file = .true.               !  =T : read bdy coordinates from file
631    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files
632    ln_mask_file   = .false.              !  =T : read mask from file
633    cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.)
634    cn_dyn2d       = 'none'               !
635    nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state
636                                          !  = 1, bdy data are read in 'bdydata   .nc' files
637                                          !  = 2, use tidal harmonic forcing data from files
638                                          !  = 3, use external data AND tidal harmonic forcing
639    cn_dyn3d      =  'none'               !
640    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state
641                                          !  = 1, bdy data are read in 'bdydata   .nc' files
642    cn_tra        =  'none'               !
643    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state
644                                          !  = 1, bdy data are read in 'bdydata   .nc' files
645    cn_ice_lim      =  'none'             !
646    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state
647                                          !  = 1, bdy data are read in 'bdydata   .nc' files
648    rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice
649    rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           --
650    rn_ice_age      = 30.                 !  lim3 only:      --   age                --
651
652    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers
653    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities
654    rn_time_dmp   =  1.                   ! Damping time scale in days
655    rn_time_dmp_out =  1.                 ! Outflow damping time scale
656    nn_rimwidth   = 10                    !  width of the relaxation zone
657    ln_vol        = .false.               !  total volume correction (see nn_volctl parameter)
658    nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero
659/
660!-----------------------------------------------------------------------
661&nambdy_dta      !  open boundaries - external data           ("key_bdy")
662!-----------------------------------------------------------------------
663!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
664!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
665   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
666   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
667   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
668   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
669   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
670   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
671   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
672! for lim2
673!   bn_frld  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
674!   bn_hicif =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
675!   bn_hsnif =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
676! for lim3
677!   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
678!   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
679!   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , ''
680   cn_dir  =    'bdydta/'
681   ln_full_vel = .false.
682/
683!-----------------------------------------------------------------------
684&nambdy_tide     ! tidal forcing at open boundaries
685!-----------------------------------------------------------------------
686   filtide          = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files
687   ln_bdytide_2ddta = .false.
688   ln_bdytide_conj  = .false.
689/
690!!======================================================================
691!!                 ***  Bottom boundary condition  ***
692!!======================================================================
693!!   nambfr        bottom friction
694!!   nambbc        bottom temperature boundary condition
695!!   nambbl        bottom boundary layer scheme                         ("key_trabbl")
696!!======================================================================
697!
698!-----------------------------------------------------------------------
699&nambfr        !   bottom friction
700!-----------------------------------------------------------------------
701   nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
702                           !                              = 2 : nonlinear friction
703   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case)
704   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
705   rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T)
706   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2)
707   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T
708   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file )
709   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T)
710   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case)
711   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
712   rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T)
713   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2)
714   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T
715   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file )
716   rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T)
717
718   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true)
719   ln_loglayer = .false.   !  logarithmic formulation (non linear case)
720/
721!-----------------------------------------------------------------------
722&nambbc        !   bottom temperature boundary condition
723!-----------------------------------------------------------------------
724!              !                              !  (if <0  months)  ! 
725!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
726!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
727   sn_qgh      ='geothermal_heating.nc',  -12.  , 'heatflow'      ,   .false.      , .true.  , 'yearly'  , ''       , ''       , ''
728   !
729   cn_dir      = './'      !  root directory for the location of the runoff files
730   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom
731   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux
732                           !     = 1 constant flux
733                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
734   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
735
736/
737!-----------------------------------------------------------------------
738&nambbl        !   bottom boundary layer scheme
739!-----------------------------------------------------------------------
740   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0)
741   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0)
742   rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s]
743   rn_gambbl   =  10.    !  advective bbl coefficient                 [s]
744/
745
746!!======================================================================
747!!                        Tracer (T & S ) namelists
748!!======================================================================
749!!   nameos        equation of state
750!!   namtra_adv    advection scheme
751!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.)
752!!   namtra_ldf    lateral diffusion scheme
753!!   namtra_dmp    T & S newtonian damping
754!!======================================================================
755!
756!-----------------------------------------------------------------------
757&nameos        !   ocean physical parameters
758!-----------------------------------------------------------------------
759   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency
760                                 !  =-1, TEOS-10
761                                 !  = 0, EOS-80
762                                 !  = 1, S-EOS   (simplified eos)
763   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm
764   !                             !
765   !                     ! S-EOS coefficients :
766   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
767   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1)
768   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1)
769   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
770   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
771   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
772   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
773   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
774/
775!-----------------------------------------------------------------------
776&namtra_adv    !   advection scheme for tracer
777!-----------------------------------------------------------------------
778   ln_traadv_cen2   =  .false.   !  2nd order centered scheme
779   ln_traadv_tvd    =  .true.    !  TVD scheme
780   ln_traadv_muscl  =  .false.   !  MUSCL scheme
781   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries
782   ln_traadv_ubs    =  .false.   !  UBS scheme
783   ln_traadv_qck    =  .false.   !  QUICKEST scheme
784   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl
785   ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection
786/
787!-----------------------------------------------------------------------
788&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param)
789!-----------------------------------------------------------------------
790   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation
791   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08)
792   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
793   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0)
794   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0)
795   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
796   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
797   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE
798   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK
799/
800!----------------------------------------------------------------------------------
801&namtra_ldf    !   lateral diffusion scheme for tracers
802!----------------------------------------------------------------------------------
803   !                       !  Operator type:
804   ln_traldf_lap    =  .true.   !  laplacian operator
805   ln_traldf_bilap  =  .false.  !  bilaplacian operator
806   !                       !  Direction of action:
807   ln_traldf_level  =  .false.  !  iso-level
808   ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T)
809   ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp")
810   !                 !  Griffies parameters              (all need "key_ldfslp")
811   ln_traldf_grif   =  .false.  !  use griffies triads
812   ln_traldf_gdia   =  .false.  !  output griffies eddy velocities
813   ln_triad_iso     =  .false.  !  pure lateral mixing in ML
814   ln_botmix_grif   =  .false.  !  lateral mixing on bottom
815   !                       !  Coefficients
816   ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv"
817   ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv
818   !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05)
819   rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s]
820   rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s]
821   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s]
822   !                                           (normally=0; not used with Griffies)
823   rn_slpmax        =     0.01  !  slope limit
824   rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity
825   rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer
826   rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s)
827/
828!-----------------------------------------------------------------------
829&namtra_dmp    !   tracer: T & S newtonian damping
830!-----------------------------------------------------------------------
831   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F)
832   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column
833                           !                   =1 no damping in the mixing layer (kz  criteria)
834                           !                   =2 no damping in the mixed  layer (rho crieria)
835   cn_resto    = 'resto.nc' ! Name of file containing restoration coefficient field (use dmp_tools to create this)
836/
837
838!!======================================================================
839!!                      ***  Dynamics namelists  ***
840!!======================================================================
841!!   namdyn_adv    formulation of the momentum advection
842!!   namdyn_vor    advection scheme
843!!   namdyn_hpg    hydrostatic pressure gradient
844!!   namdyn_spg    surface pressure gradient                            (CPP key only)
845!!   namdyn_ldf    lateral diffusion scheme
846!!======================================================================
847!
848!-----------------------------------------------------------------------
849&namdyn_adv    !   formulation of the momentum advection
850!-----------------------------------------------------------------------
851   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)
852   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction
853   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
854   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
855   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection
856/
857!-----------------------------------------------------------------------
858&nam_vvl    !   vertical coordinate options
859!-----------------------------------------------------------------------
860   ln_vvl_zstar  = .true.           !  zstar vertical coordinate
861   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations
862   ln_vvl_layer  = .false.          !  full layer vertical coordinate
863   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
864   ln_vvl_zstar_at_eqtor = .false.  !  ztilde near the equator
865   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient
866   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days]
867   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days]
868   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation
869   ln_vvl_dbg    = .true.           !  debug prints    (T/F)
870/
871!-----------------------------------------------------------------------
872&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
873!-----------------------------------------------------------------------
874   ln_dynvor_ene = .false. !  enstrophy conserving scheme
875   ln_dynvor_ens = .false. !  energy conserving scheme
876   ln_dynvor_mix = .false. !  mixed scheme
877   ln_dynvor_een = .true.  !  energy & enstrophy scheme
878   ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation
879/
880!-----------------------------------------------------------------------
881&namdyn_hpg    !   Hydrostatic pressure gradient option
882!-----------------------------------------------------------------------
883   ln_hpg_zco  = .false.   !  z-coordinate - full steps
884   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation)
885   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
886   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
887   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
888   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
889   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T)
890                                 !           centered      time scheme  (F)
891/
892!-----------------------------------------------------------------------
893!namdyn_spg    !   surface pressure gradient   (CPP key only)
894!-----------------------------------------------------------------------
895!                          !  explicit free surface                     ("key_dynspg_exp")
896!                          !  filtered free surface                     ("key_dynspg_flt")
897!                          !  split-explicit free surface               ("key_dynspg_ts")
898
899!-----------------------------------------------------------------------
900&namdyn_ldf    !   lateral diffusion on momentum
901!-----------------------------------------------------------------------
902   !                       !  Type of the operator :
903   ln_dynldf_lap    =  .true.   !  laplacian operator
904   ln_dynldf_bilap  =  .false.  !  bilaplacian operator
905   !                       !  Direction of action  :
906   ln_dynldf_level  =  .false.  !  iso-level
907   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.)
908   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp")
909   !                       !  Coefficient
910   rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s]
911   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s]
912   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s]
913   rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity
914   rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity
915   rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity
916   rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp)
917   rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap)
918/
919
920!!======================================================================
921!!             Tracers & Dynamics vertical physics namelists
922!!======================================================================
923!!    namzdf        vertical physics
924!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric")
925!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke")
926!!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp")
927!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm")
928!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx")
929!!======================================================================
930!
931!-----------------------------------------------------------------------
932&namzdf        !   vertical physics
933!-----------------------------------------------------------------------
934   rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
935   rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
936   nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0)
937   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0)
938   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F)
939   nn_evdm     =    0      !  evd apply on tracer (=0) or on tracer and momentum (=1)
940   rn_avevd    =  100.     !  evd mixing coefficient [m2/s]
941   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F)
942   nn_npc      =    1            !  frequency of application of npc
943   nn_npcp     =  365            !  npc control print frequency
944   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
945   nn_zdfexp   =    3            !  number of sub-timestep for ln_zdfexp=T
946/
947!-----------------------------------------------------------------------
948&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
949!-----------------------------------------------------------------------
950   rn_avmri    = 100.e-4   !  maximum value of the vertical viscosity
951   rn_alp      =   5.      !  coefficient of the parameterization
952   nn_ric      =   2       !  coefficient of the parameterization
953   rn_ekmfc    =   0.7     !  Factor in the Ekman depth Equation
954   rn_mldmin   =   1.0     !  minimum allowable mixed-layer depth estimate (m)
955   rn_mldmax   =1000.0     !  maximum allowable mixed-layer depth estimate (m)
956   rn_wtmix    =  10.0     !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
957   rn_wvmix    =  10.0     !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
958   ln_mldw     = .true.    !  Flag to use or not the mized layer depth param.
959/
960!-----------------------------------------------------------------------
961&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
962!-----------------------------------------------------------------------
963   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
964   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
965   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
966   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
967   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
968   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
969   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
970                           !                 = 1 bounded by the local vertical scale factor
971                           !                 = 2 first vertical derivative of mixing length bounded by 1
972                           !                 = 3 as =2 with distinct disspipative an mixing length scale
973   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
974   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
975   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
976   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
977   rn_lc       =   0.15    !  coef. associated to Langmuir cells
978   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves
979                           !        = 0 no penetration
980                           !        = 1 add a tke source below the ML
981                           !        = 2 add a tke source just at the base of the ML
982                           !        = 3 as = 1 applied on HF part of the stress    ("key_oasis3")
983   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
984   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
985                           !        = 0  constant 10 m length scale
986                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees
987/
988!------------------------------------------------------------------------
989&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally:
990!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
991   ln_kpprimix = .true.    !  shear instability mixing
992   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s]
993   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s]
994   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability
995   rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s]
996   rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]
997   rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s]
998   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv
999   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt
1000/
1001!-----------------------------------------------------------------------
1002&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls")
1003!-----------------------------------------------------------------------
1004   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
1005   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
1006   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
1007   rn_clim_galp  = 0.267   !  galperin limit
1008   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
1009   rn_crban_default = 100.    !  Craig and Banner 1994 constant for wb tke flux
1010   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
1011   rn_hsro       =  0.02   !  Minimum surface roughness
1012   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2)
1013   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3)
1014      !                             ! =3 requires ln_wave=T
1015   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
1016   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
1017   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
1018   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
1019/
1020!-----------------------------------------------------------------------
1021&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
1022!-----------------------------------------------------------------------
1023   rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
1024   rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
1025/
1026!-----------------------------------------------------------------------
1027&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
1028!-----------------------------------------------------------------------
1029   rn_htmx     = 500.      !  vertical decay scale for turbulence (meters)
1030   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1)
1031   rn_tfe      = 0.333     !  tidal dissipation efficiency
1032   rn_me       = 0.2       !  mixing efficiency
1033   ln_tmx_itf  = .true.    !  ITF specific parameterisation
1034   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency
1035/
1036
1037!!======================================================================
1038!!                  ***  Miscellaneous namelists  ***
1039!!======================================================================
1040!!   namsol            elliptic solver / island / free surface
1041!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi)
1042!!   namctl            Control prints & Benchmark
1043!!   namc1d            1D configuration options                         ("key_c1d")
1044!!   namc1d_uvd        data: U & V currents                             ("key_c1d")
1045!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d")
1046!!   namsto            Stochastic parametrization of EOS
1047!!======================================================================
1048!
1049!-----------------------------------------------------------------------
1050&namsol        !   elliptic solver / island / free surface
1051!-----------------------------------------------------------------------
1052   nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
1053                           !                   =2 successive-over-relaxation (sor)
1054   nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test
1055   rn_eps      =  1.e-6    !  absolute precision of the solver
1056   nn_nmin     =    300    !  minimum of iterations for the SOR solver
1057   nn_nmax     =    800    !  maximum of iterations for the SOR solver
1058   nn_nmod     =     10    !  frequency of test for the SOR solver
1059   rn_resmax   =  1.e-10   !  absolute precision for the SOR solver
1060   rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain)
1061/
1062!-----------------------------------------------------------------------
1063&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
1064!-----------------------------------------------------------------------
1065   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
1066                           !  buffer blocking send or immediate non-blocking sends, resp.
1067   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
1068   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold
1069   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1)
1070   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1)
1071   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1)
1072/
1073!-----------------------------------------------------------------------
1074&namctl        !   Control prints & Benchmark
1075!-----------------------------------------------------------------------
1076   ln_ctl      = .false.   !  trends control print (expensive!)
1077   nn_print    =    0      !  level of print (0 no extra print)
1078   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
1079   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
1080   nn_jctls    =    0      !  start j indice of control               over a subdomain)
1081   nn_jctle    =    0      !  end   j indice of control
1082   nn_isplt    =    1      !  number of processors in i-direction
1083   nn_jsplt    =    1      !  number of processors in j-direction
1084   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench
1085                           !     (no physical validity of the results)
1086   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0)
1087/
1088!-----------------------------------------------------------------------
1089&namc1d_uvd    !   data: U & V currents                                 ("key_c1d")
1090!-----------------------------------------------------------------------
1091!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1092!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
1093   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , ''
1094   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , ''
1095!
1096   cn_dir        = './'    !  root directory for the location of the files
1097   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F)
1098   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F)
1099/
1100!-----------------------------------------------------------------------
1101&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d")
1102!-----------------------------------------------------------------------
1103   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
1104/
1105!-----------------------------------------------------------------------
1106&namsto       ! Stochastic parametrization of EOS
1107!-----------------------------------------------------------------------
1108   ln_rststo = .false.           ! start from mean parameter (F) or from restart file (T)
1109   ln_rstseed = .true.           ! read seed of RNG from restart file
1110   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
1111   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)
1112
1113   ln_sto_eos = .false.          ! stochastic equation of state
1114   nn_sto_eos = 1                ! number of independent random walks
1115   rn_eos_stdxy = 1.4            ! random walk horz. standard deviation (in grid points)
1116   rn_eos_stdz  = 0.7            ! random walk vert. standard deviation (in grid points)
1117   rn_eos_tcor  = 1440.0         ! random walk time correlation (in timesteps)
1118   nn_eos_ord  = 1               ! order of autoregressive processes
1119   nn_eos_flt  = 0               ! passes of Laplacian filter
1120   rn_eos_lim  = 2.0             ! limitation factor (default = 3.0)
1121/
1122
1123!!======================================================================
1124!!                  ***  Diagnostics namelists  ***
1125!!======================================================================
1126!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1127!!   namtrd       dynamics and/or tracer trends
1128!!   namptr       Poleward Transport Diagnostics
1129!!   namflo       float parameters                                      ("key_float")
1130!!   namhsb       Heat and salt budgets
1131!!======================================================================
1132!
1133!-----------------------------------------------------------------------
1134&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
1135!-----------------------------------------------------------------------
1136   nn_nchunks_i=   4       !  number of chunks in i-dimension
1137   nn_nchunks_j=   4       !  number of chunks in j-dimension
1138   nn_nchunks_k=   31      !  number of chunks in k-dimension
1139                           !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1140                           !  is optimal for postprocessing which works exclusively with horizontal slabs
1141   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1142                           !  (F) ignore chunking information and produce netcdf3-compatible files
1143/
1144!-----------------------------------------------------------------------
1145&namtrd        !   diagnostics on dynamics and/or tracer trends
1146!              !       and/or mixed-layer trends and/or barotropic vorticity
1147!-----------------------------------------------------------------------
1148   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1149   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1150   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1151   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet)
1152   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1153   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1154   ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output
1155   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1156   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1157/
1158!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1159!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1160!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1161!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1162!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1163!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1164!!gm
1165!-----------------------------------------------------------------------
1166&namflo       !   float parameters                                      ("key_float")
1167!-----------------------------------------------------------------------
1168   jpnfl         = 1          !  total number of floats during the run
1169   jpnnewflo     = 0          !  number of floats for the restart
1170   ln_rstflo     = .false.    !  float restart (T) or not (F)
1171   nn_writefl    =      75    !  frequency of writing in float output file
1172   nn_stockfl    =    5475    !  frequency of creation of the float restart file
1173   ln_argo       = .false.    !  Argo type floats (stay at the surface each 10 days)
1174   ln_flork4     = .false.    !  trajectories computed with a 4th order Runge-Kutta (T)
1175                              !  or computed with Blanke' scheme (F)
1176   ln_ariane     = .true.     !  Input with Ariane tool convention(T)
1177   ln_flo_ascii  = .true.     !  Output with Ariane tool netcdf convention(F) or ascii file (T)
1178/
1179!-----------------------------------------------------------------------
1180&namptr       !   Poleward Transport Diagnostic
1181!-----------------------------------------------------------------------
1182   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F)
1183   ln_subbas  = .false.     !  Atlantic/Pacific/Indian basins computation (T) or not
1184/
1185!-----------------------------------------------------------------------
1186&namhsb       !  Heat and salt budgets
1187!-----------------------------------------------------------------------
1188   ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F)
1189/
1190!-----------------------------------------------------------------------
1191&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm')
1192!-----------------------------------------------------------------------
1193    nit000_han = 1         ! First time step used for harmonic analysis
1194    nitend_han = 75        ! Last time step used for harmonic analysis
1195    nstep_han  = 15        ! Time step frequency for harmonic analysis
1196    tname(1)   = 'M2'      ! Name of tidal constituents
1197    tname(2)   = 'K1'
1198/
1199!-----------------------------------------------------------------------
1200&namdct        ! transports through sections
1201!-----------------------------------------------------------------------
1202    nn_dct      = 15       !  time step frequency for transports computing
1203    nn_dctwri   = 15       !  time step frequency for transports writing
1204    nn_secdebug = 112      !      0 : no section to debug
1205                           !     -1 : debug all section
1206                           !  0 < n : debug section number n
1207/
1208
1209!!======================================================================
1210!!            ***  Observation & Assimilation namelists ***
1211!!======================================================================
1212!!   namobs       observation and model comparison                      ('key_diaobs')
1213!!   nam_asminc   assimilation increments                               ('key_asminc')
1214!!======================================================================
1215!
1216!-----------------------------------------------------------------------
1217&namobs       !  observation usage switch                               ('key_diaobs')
1218!-----------------------------------------------------------------------
1219   ln_t3d     = .false.    ! Logical switch for T profile observations
1220   ln_s3d     = .false.    ! Logical switch for S profile observations
1221   ln_ena     = .false.    ! Logical switch for ENACT insitu data set
1222   !                       !     ln_cor                  Logical switch for Coriolis insitu data set
1223   ln_profb   = .false.    ! Logical switch for feedback insitu data set
1224   ln_sla     = .false.    ! Logical switch for SLA observations
1225
1226   ln_sladt   = .false.    ! Logical switch for AVISO SLA data
1227
1228   ln_slafb   = .false.    ! Logical switch for feedback SLA data
1229                           !     ln_ssh                  Logical switch for SSH observations
1230
1231   ln_sst     = .false.     ! Logical switch for SST observations
1232   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations
1233   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations
1234
1235   ln_sstfb   = .false.    ! Logical switch for feedback SST data
1236                           !     ln_sss                  Logical switch for SSS observations
1237   ln_seaice  = .false.    ! Logical switch for Sea Ice observations
1238                           !     ln_vel3d                Logical switch for velocity observations
1239                           !     ln_velavcur             Logical switch for velocity daily av. cur.
1240                           !     ln_velhrcur             Logical switch for velocity high freq. cur.
1241                           !     ln_velavadcp            Logical switch for velocity daily av. ADCP
1242                           !     ln_velhradcp            Logical switch for velocity high freq. ADCP
1243                           !     ln_velfb                Logical switch for feedback velocity data
1244                           !     ln_grid_global          Global distribtion of observations
1245                           !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table
1246                           !     grid_search_file        Grid search lookup file header
1247                           !     enactfiles              ENACT input observation file names
1248                           !     coriofiles              Coriolis input observation file name
1249   !                       ! profbfiles: Profile feedback input observation file name
1250   profbfiles = 'profiles_01.nc'
1251                           !     ln_profb_enatim         Enact feedback input time setting switch
1252                           !     slafilesact             Active SLA input observation file name
1253                           !     slafilespas             Passive SLA input observation file name
1254   !                       ! slafbfiles: Feedback SLA input observation file name
1255   slafbfiles = 'sla_01.nc'
1256                           !     sstfiles                GHRSST input observation file name
1257   !                       ! sstfbfiles: Feedback SST input observation file name
1258   sstfbfiles = 'sst_01.nc'
1259                           !     seaicefiles             Sea Ice input observation file names
1260   seaicefiles = 'seaice_01.nc'
1261                           !     velavcurfiles           Vel. cur. daily av. input file name
1262                           !     velhvcurfiles           Vel. cur. high freq. input file name
1263                           !     velavadcpfiles          Vel. ADCP daily av. input file name
1264                           !     velhvadcpfiles          Vel. ADCP high freq. input file name
1265                           !     velfbfiles              Vel. feedback input observation file name
1266                           !     dobsini                 Initial date in window YYYYMMDD.HHMMSS
1267                           !     dobsend                 Final date in window YYYYMMDD.HHMMSS
1268                           !     n1dint                  Type of vertical interpolation method
1269                           !     n2dint                  Type of horizontal interpolation method
1270                           !     ln_nea                  Rejection of observations near land switch
1271   nmsshc     = 0          ! MSSH correction scheme
1272                           !     mdtcorr                 MDT  correction
1273                           !     mdtcutoff               MDT cutoff for computed correction
1274   ln_altbias = .false.    ! Logical switch for alt bias
1275   ln_ignmis  = .true.     ! Logical switch for ignoring missing files
1276                           !     endailyavtypes   ENACT daily average types
1277   ln_grid_global = .true.
1278   ln_grid_search_lookup = .false.
1279/
1280!-----------------------------------------------------------------------
1281&nam_asminc   !   assimilation increments                               ('key_asminc')
1282!-----------------------------------------------------------------------
1283    ln_bkgwri      = .false. !  Logical switch for writing out background state
1284    ln_balwri      = .false. !  Logical switch for writing out balancing increments
1285    ln_trainc      = .false. !  Logical switch for applying tracer increments
1286    ln_dyninc      = .false. !  Logical switch for applying velocity increments
1287    ln_sshinc      = .false. !  Logical switch for applying SSH increments
1288    ln_asmdin      = .false. !  Logical switch for Direct Initialization (DI)
1289    ln_asmiau      = .false. !  Logical switch for Incremental Analysis Updating (IAU)
1290    ln_phytobal    = .false. !  Logical switch for phytoplankton multivariate balancing
1291    ln_slchltotinc = .false. !  Logical switch for applying slchltot increments
1292    ln_slchldiainc = .false. !  Logical switch for applying slchldia increments
1293    ln_slchlnoninc = .false. !  Logical switch for applying slchlnon increments
1294    ln_slchlnaninc = .false. !  Logical switch for applying slchlnan increments
1295    ln_slchlpicinc = .false. !  Logical switch for applying slchlpic increments
1296    ln_slchldininc = .false. !  Logical switch for applying slchldin increments
1297    ln_schltotinc  = .false. !  Logical switch for applying schltot increments
1298    ln_slphytotinc = .false. !  Logical switch for applying slphytot increments
1299    ln_slphydiainc = .false. !  Logical switch for applying slphydia increments
1300    ln_slphynoninc = .false. !  Logical switch for applying slphynon increments
1301    ln_sfco2inc    = .false. !  Logical switch for applying sfCO2 increments
1302    ln_spco2inc    = .false. !  Logical switch for applying spCO2 increments
1303    ln_plchltotinc = .false. !  Logical switch for applying plchltot increments
1304    ln_pchltotinc  = .false. !  Logical switch for applying pchltot increments
1305    ln_pno3inc     = .false. !  Logical switch for applying pno3 increments
1306    ln_psi4inc     = .false. !  Logical switch for applying psi4 increments
1307    ln_ppo4inc     = .false. !  Logical switch for applying ppo4 increments
1308    ln_pdicinc     = .false. !  Logical switch for applying pdic increments
1309    ln_palkinc     = .false. !  Logical switch for applying palk increments
1310    ln_pphinc      = .false. !  Logical switch for applying pph increments
1311    ln_po2inc      = .false. !  Logical switch for applying po2 increments
1312    nitbkg         = 0       !  Timestep of background in [0,nitend-nit000-1]
1313    nitdin         = 0       !  Timestep of background for DI in [0,nitend-nit000-1]
1314    nitiaustr      = 1       !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1315    nitiaufin      = 15      !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1316    niaufn         = 0       !  Type of IAU weighting function
1317    ln_salfix      = .false. !  Logical switch for ensuring that the sa > salfixmin
1318    salfixmin      = -9999   !  Minimum salinity after applying the increments
1319    nn_divdmp      = 0       !  Number of iterations of divergence damping operator
1320    mld_choice_bgc = 1       !  MLD criterion to use for biogeochemistry assimilation
1321    rn_maxchlinc   = -999.0  !  maximum absolute non-log chlorophyll increment from ocean colour assimilation
1322                             !  <= 0 implies no maximum applied (switch turned off)
1323                             !   > 0 implies maximum absolute chl increment capped at this value
1324/
1325!-----------------------------------------------------------------------
1326&namsbc_wave   ! External fields from wave model                        (ln_wave=T)
1327!-----------------------------------------------------------------------
1328!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
1329!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
1330   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model
1331   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift
1332   ln_stcor    = .false.   !  Activate Stokes Coriolis term
1333   ln_tauoc    = .false.   !  Activate ocean stress modified by external wave induced stress
1334   ln_tauw     = .false.   !  Activate ocean stress components from wave model
1335   ln_phioc    = .false.   !  Activate wave to ocean energy
1336   ln_rough    = .false.   !  Wave roughness equals the significant wave height
1337   ln_zdfqiao  = .false.   !  Enhanced wave vertical mixing Qiao (2010)
1338   nn_sdrift   =  0   !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
1339                           ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
1340                           ! = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
1341   nn_wmix     =  0   !  type of wave breaking mixing
1342                           ! = 0 Craig and Banner formulation (original NEMO formulation) 
1343                           ! = 1 Janssen formulation (no assumption of direct energy conversion)
1344   sn_cdg      =  'cdw_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1345   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1346   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1347   sn_hsw      =  'sdw_wave' ,        1          , 'hs'         ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
1348   sn_wmp      =  'sdw_wave' ,        1          , 'wmp'        ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
1349   sn_wfr      =  'sdw_wave' ,        1          , 'wave_freq'  ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
1350   sn_wnum     =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
1351   sn_tauoc    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
1352   sn_tauwx    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
1353   sn_tauwy    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
1354   sn_phioc    =  'sdw_wave' ,        1          , 'wave_energy',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
1355   cn_dir      = './'  !  root directory for the location of wave forcing files
1356
1357/
1358!-----------------------------------------------------------------------
1359&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
1360!-----------------------------------------------------------------------
1361   ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model
1362   ln_neptsimp       = .false.  ! yes/no use simplified neptune
1363
1364   ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep
1365   rn_tslse          =  1.2e4   ! value of lengthscale L at the equator
1366   rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole
1367   ! Specify whether to ramp down the Neptune velocity in shallow
1368   ! water, and if so the depth range controlling such ramping down
1369   ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water
1370   rn_htrmin         =  100.0   ! min. depth of transition range
1371   rn_htrmax         =  200.0   ! max. depth of transition range
1372/
Note: See TracBrowser for help on using the repository browser.