New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
namelist_ref in branches/2017/dev_r8183_ICEMODEL/NEMOGCM/TOOLS/DOMAINcfg – NEMO

source: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/TOOLS/DOMAINcfg/namelist_ref @ 8306

Last change on this file since 8306 was 8306, checked in by clem, 7 years ago

step1: remove LIM2 from the code

File size: 92.5 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!!                            namelist_ref
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!! NEMO/OPA  :  1 - run manager      (namrun)
5!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd, namcrs, namc1d, namc1d_uvd)
6!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas
7!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,
8!!                                    namsbc_apr, namsbc_ssr, namsbc_alb, namsbc_wave)
9!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide)
10!!              5 - bottom  boundary (nambfr, nambbc, nambbl)
11!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp)
12!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
13!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx)
14!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto)
15!!             10 - miscellaneous    (nammpp, namctl)
16!!             11 - Obs & Assim      (namobs, nam_asminc)
17!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
18
19!!======================================================================
20!!                   ***  Run management namelists  ***
21!!======================================================================
22!!   namrun       parameters of the run
23!!======================================================================
24!
25!-----------------------------------------------------------------------
26&namrun        !   parameters of the run
27!-----------------------------------------------------------------------
28   nn_no       =       0   !  job number (no more used...)
29   cn_exp      =  "ORCA2"  !  experience name
30   nn_it000    =       1   !  first time step
31   nn_itend    =    5475   !  last  time step (std 5475)
32   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
33   nn_time0    =       0   !  initial time of day in hhmm
34   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
35   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
36      nn_euler    =    1            !  = 0 : start with forward time step if ln_rstart=T
37      nn_rstctl   =    0            !  restart control ==> activated only if ln_rstart=T
38      !                             !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
39      !                             !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
40      !                             !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart
41      cn_ocerst_in    = "restart"   !  suffix of ocean restart name (input)
42      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts
43      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output)
44      cn_ocerst_outdir= "."         !  directory in which to write output ocean restarts
45   ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model
46   nn_istate   =       0   !  output the initial state (1) or not (0)
47   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
48   nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1)
49   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
50   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000)
51   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%)
52   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard
53   ln_clobber  = .true.    !  clobber (overwrite) an existing file
54   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
55/
56!
57!!======================================================================
58!!                      ***  Domain namelists  ***
59!!======================================================================
60!!   namcfg       parameters of the configuration
61!!   namzgr       vertical coordinate                                   (default: NO selection)
62!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate
63!!   namdom       space and time domain (bathymetry, mesh, timestep)
64!!   namwad       Wetting and drying                                    (default F)
65!!   namtsd       data: temperature & salinity
66!!   namcrs       coarsened grid (for outputs and/or TOP)               ("key_crs")
67!!   namc1d       1D configuration options                              ("key_c1d")
68!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
69!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
70!!======================================================================
71!
72!-----------------------------------------------------------------------
73&namcfg        !   parameters of the configuration
74!-----------------------------------------------------------------------
75   !
76   ln_e3_dep   = .true.    ! =T : e3=dk[depth] in discret sens.
77   !                       !      ===>>> will become the only possibility in v4.0
78   !                       ! =F : e3 analytical derivative of depth function
79   !                       !      only there for backward compatibility test with v3.6
80   !
81   cp_cfg      = "default" !  name of the configuration
82   cp_cfz      = "no zoom" !  name of the zoom of configuration
83   jp_cfg      =      0    !  resolution of the configuration
84   jpidta      =     10    !  1st lateral dimension ( >= jpi )
85   jpjdta      =     12    !  2nd    "         "    ( >= jpj )
86   jpkdta      =     31    !  number of levels      ( >= jpk )
87   jpiglo      =     10    !  1st dimension of global domain --> i =jpidta
88   jpjglo      =     12    !  2nd    -                  -    --> j =jpjdta
89   jpizoom     =      1    !  left bottom (i,j) indices of the zoom
90   jpjzoom     =      1    !  in data domain indices
91   jperio      =      0    !  lateral cond. type (between 0 and 6)
92                                 !  = 0 closed                 ;   = 1 cyclic East-West
93                                 !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot
94                                 !  = 4 cyclic East-West AND North fold T-point pivot
95                                 !  = 5 North fold F-point pivot
96                                 !  = 6 cyclic East-West AND North fold F-point pivot
97   ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present
98                           !  in netcdf input files, as the start j-row for reading
99/
100!-----------------------------------------------------------------------
101&namzgr        !   vertical coordinate                                  (default: NO selection)
102!-----------------------------------------------------------------------
103   ln_zco      = .false.   !  z-coordinate - full    steps
104   ln_zps      = .false.   !  z-coordinate - partial steps
105   ln_sco      = .false.   !  s- or hybrid z-s-coordinate
106   ln_isfcav   = .false.   !  ice shelf cavity
107   ln_linssh   = .false.   !  linear free surface
108/
109!-----------------------------------------------------------------------
110&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default F)
111!-----------------------------------------------------------------------
112   ln_s_sh94   = .false.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)|
113   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied
114   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch
115                           !  stretching coefficients for all functions
116   rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m)
117   rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m)
118   rn_hc       =  150.0    !  critical depth for transition to stretched coordinates
119                        !!!!!!!  Envelop bathymetry
120   rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1)
121                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.)
122   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20)
123   rn_bb       =    0.8    !  stretching with SH94 s-sigma
124                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.)
125   rn_alpha    =    4.4    !  stretching with SF12 s-sigma
126   rn_efold    =    0.0    !  efold length scale for transition to stretched coord
127   rn_zs       =    1.0    !  depth of surface grid box
128                           !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b
129   rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb
130   rn_zb_b     =   -0.2    !  offset for calculating Zb
131                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above]
132   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)
133/
134!-----------------------------------------------------------------------
135&namdom        !   space and time domain (bathymetry, mesh, timestep)
136!-----------------------------------------------------------------------
137   nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file
138   rn_bathy    =    0.     !  value of the bathymetry. if (=0) bottom flat at jpkm1
139   nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA)
140   nn_msh      =    1      !  create (=1) a mesh file or not (=0)
141   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0)
142   rn_isfhmin  =    1.00   !  treshold (m) to discriminate grounding ice to floating ice
143   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of
144   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1
145                           !
146   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0)
147   rn_atfp     =    0.1    !  asselin time filter parameter
148   ln_crs      = .false.      !  Logical switch for coarsening module
149   jphgr_msh   =       0               !  type of horizontal mesh
150                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc
151                                       !  = 1 geographical mesh on the sphere with regular grid-spacing
152                                       !  = 2 f-plane with regular grid-spacing
153                                       !  = 3 beta-plane with regular grid-spacing
154                                       !  = 4 Mercator grid with T/U point at the equator
155   ppglam0     =       0.0             !  longitude of first raw and column T-point (jphgr_msh = 1)
156   ppgphi0     =     -35.0             ! latitude  of first raw and column T-point (jphgr_msh = 1)
157   ppe1_deg    =       1.0             !  zonal      grid-spacing (degrees)
158   ppe2_deg    =       0.5             !  meridional grid-spacing (degrees)
159   ppe1_m      =    5000.0             !  zonal      grid-spacing (degrees)
160   ppe2_m      =    5000.0             !  meridional grid-spacing (degrees)
161   ppsur       =    -4762.96143546300  !  ORCA r4, r2 and r05 coefficients
162   ppa0        =      255.58049070440  ! (default coefficients)
163   ppa1        =      245.58132232490  !
164   ppkth       =       21.43336197938  !
165   ppacr       =        3.0            !
166   ppdzmin     =       10.             !  Minimum vertical spacing
167   pphmax      =     5000.             !  Maximum depth
168   ldbletanh   =    .TRUE.             !  Use/do not use double tanf function for vertical coordinates
169   ppa2        =      100.760928500000 !  Double tanh function parameters
170   ppkth2      =       48.029893720000 !
171   ppacr2      =       13.000000000000 !
172/
173!-----------------------------------------------------------------------
174&namwad        !   Wetting and drying                                   (default F)
175!-----------------------------------------------------------------------
176   ln_wd       = .false.   !  T/F activation of wetting and drying
177   rn_wdmin1   =  0.1      !  Minimum wet depth on dried cells
178   rn_wdmin2   =  0.01     !  Tolerance of min wet depth on dried cells
179   rn_wdld     =  20.0     !  Land elevation below which wetting/drying is allowed
180   nn_wdit     =  10       !  Max iterations for W/D limiter
181/
182!-----------------------------------------------------------------------
183&namtsd        !   data : Temperature  & Salinity
184!-----------------------------------------------------------------------
185!              !  file name                 ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
186!              !                            !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
187   sn_tem = 'data_1m_potential_temperature_nomask',     -1      ,'votemper',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
188   sn_sal = 'data_1m_salinity_nomask'             ,     -1      ,'vosaline',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,    ''
189   !
190   cn_dir      = './'      !  root directory for the location of the runoff files
191   ln_tsd_init = .true.    !  Initialisation of ocean T & S with T & S input data (T) or not (F)
192   ln_tsd_tradmp = .true.  !  damping of ocean T & S toward T & S input data (T) or not (F)
193/
194!-----------------------------------------------------------------------
195&namcrs        !   coarsened grid (for outputs and/or TOP)              ("key_crs")
196!-----------------------------------------------------------------------
197   nn_factx    = 3         !  Reduction factor of x-direction
198   nn_facty    = 3         !  Reduction factor of y-direction
199   nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT
200                           !  0, coarse grid is binned with preferential treatment of the north fold
201                           !  1, coarse grid is binned with centering at the equator
202                           !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.
203   nn_msh_crs  = 1         !  create (=1) a mesh file or not (=0)
204   nn_crs_kz   = 0         ! 0, MEAN of volume boxes
205                           ! 1, MAX of boxes
206                           ! 2, MIN of boxes
207   ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F )
208/
209!-----------------------------------------------------------------------
210&namc1d        !   1D configuration options                             ("key_c1d")
211!-----------------------------------------------------------------------
212   rn_lat1d    =      50   !  Column latitude (default at PAPA station)
213   rn_lon1d    =    -145   !  Column longitude (default at PAPA station)
214   ln_c1d_locpt=  .true.   ! Localization of 1D config in a grid (T) or independant point (F)
215/
216!-----------------------------------------------------------------------
217&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d")
218!-----------------------------------------------------------------------
219   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
220/
221!-----------------------------------------------------------------------
222&namc1d_uvd    !   data: U & V currents                                 ("key_c1d")
223!-----------------------------------------------------------------------
224!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
225!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
226   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , ''
227   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , ''
228!
229   cn_dir        = './'    !  root directory for the location of the files
230   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F)
231   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F)
232/
233
234!!======================================================================
235!!            ***  Surface Boundary Condition namelists  ***
236!!======================================================================
237!!   namsbc          surface boundary condition
238!!   namsbc_ana      analytical         formulation                     (ln_ana     =T)
239!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
240!!   namsbc_clio     CLIO bulk formulae formulation                     (ln_blk_clio=T)
241!!   namsbc_core     CORE bulk formulae formulation                     (ln_blk_core=T)
242!!   namsbc_mfs      MFS  bulk formulae formulation                     (ln_blk_mfs =T)
243!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
244!!   namsbc_sas      StAndalone Surface module
245!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
246!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
247!!   namsbc_isf      ice shelf melting/freezing                         (nn_isf     >0)
248!!   namsbc_iscpl    coupling option between land ice model and ocean
249!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
250!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
251!!   namsbc_alb      albedo parameters
252!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
253!!   namberg         iceberg floats                                     (ln_icebergs=T)
254!!======================================================================
255!
256!-----------------------------------------------------------------------
257&namsbc        !   Surface Boundary Condition (surface module)
258!-----------------------------------------------------------------------
259   nn_fsbc     = 5         !  frequency of surface boundary condition computation
260                           !     (also = the frequency of sea-ice & iceberg model call)
261                     ! Type of air-sea fluxes
262   ln_ana      = .false.   !  analytical formulation                    (T => fill namsbc_ana )
263   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx )
264   ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio)
265   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core)
266   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs )
267                     ! Type of coupling (Ocean/Ice/Atmosphere) :
268   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 )
269   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 )
270   nn_components = 0       !  configuration of the opa-sas OASIS coupling
271                           !  =0 no opa-sas OASIS coupling: default single executable configuration
272                           !  =1 opa-sas OASIS coupling: multi executable configuration, OPA component
273                           !  =2 opa-sas OASIS coupling: multi executable configuration, SAS component
274   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used)
275                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled
276                           !  = 0  Average per-category fluxes (forced and coupled mode)
277                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled
278                           !  = 2  Redistribute a single flux over categories (coupled mode only)
279                     ! Sea-ice :
280   nn_ice      = 3         !  =0 no ice boundary condition   ,
281                           !  =1 use observed ice-cover      ,
282                           !  =3-4 ice-model used                         ("key_lim3", "key_cice")
283   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect)
284                           !  =1 levitating ice with mass and salt exchange but no presure effect
285                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure)
286                     ! Misc. options of sbc :
287   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr )
288   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave
289   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf)
290   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
291   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked
292                           !     =1 global mean of e-p-r set to zero at each time step
293                           !     =2 annual global mean of e-p-r set to zero
294   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr )
295   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf)
296   ln_wave     = .false.   !  coupling with surface wave                (T => fill namsbc_wave)
297   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
298                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
299/
300!-----------------------------------------------------------------------
301&namsbc_ana    !   analytical surface boundary condition
302!-----------------------------------------------------------------------
303   nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps
304   rn_utau0    =   0.5     !  uniform value for the i-stress
305   rn_vtau0    =   0.e0    !  uniform value for the j-stress
306   rn_qns0     =   0.e0    !  uniform value for the total heat flux
307   rn_qsr0     =   0.e0    !  uniform value for the solar radiation
308   rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P)
309/
310!-----------------------------------------------------------------------
311&namsbc_flx    !   surface boundary condition : flux formulation
312!-----------------------------------------------------------------------
313!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
314!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
315   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
316   sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
317   sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , ''       , ''
318   sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
319   sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , ''
320
321   cn_dir      = './'      !  root directory for the location of the flux files
322/
323!-----------------------------------------------------------------------
324&namsbc_clio   !   namsbc_clio  CLIO bulk formulae
325!-----------------------------------------------------------------------
326!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
327!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
328   sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
329   sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
330   sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
331   sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
332   sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
333   sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
334   sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
335
336   cn_dir      = './'      !  root directory for the location of the bulk files are
337/
338!-----------------------------------------------------------------------
339&namsbc_core   !   namsbc_core  CORE bulk formulae
340!-----------------------------------------------------------------------
341!              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask !
342!              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      !
343   sn_wndi     = 'u_10.15JUNE2009_fill'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , ''
344   sn_wndj     = 'v_10.15JUNE2009_fill'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , ''
345   sn_qsr      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
346   sn_qlw      = 'ncar_rad.15JUNE2009_fill'   ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
347   sn_tair     = 't_10.15JUNE2009_fill'       ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
348   sn_humi     = 'q_10.15JUNE2009_fill'       ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
349   sn_prec     = 'ncar_precip.15JUNE2009_fill',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
350   sn_snow     = 'ncar_precip.15JUNE2009_fill',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
351   sn_tdif     = 'taudif_core'                ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , ''
352
353   cn_dir      = './'      !  root directory for the location of the bulk files
354   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data
355   rn_zqt      = 10.       !  Air temperature and humidity reference height (m)
356   rn_zu       = 10.       !  Wind vector reference height (m)
357   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow)
358   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.)
359   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity
360                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)
361/
362!-----------------------------------------------------------------------
363&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae
364!-----------------------------------------------------------------------
365!              !  file name  ! frequency (hours) ! variable ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask !
366!              !             !  (if <0  months)  !   name   !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      !
367   sn_wndi     =   'ecmwf'   ,        6          , 'u10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   ''
368   sn_wndj     =   'ecmwf'   ,        6          , 'v10'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   ''
369   sn_clc      =   'ecmwf'   ,        6          , 'clc'    ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,   ''
370   sn_msl      =   'ecmwf'   ,        6          , 'msl'    ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   ''
371   sn_tair     =   'ecmwf'   ,        6          , 't2'     ,    .true.    , .false., 'daily'   ,'bicubic.nc' ,   ''     ,   ''
372   sn_rhm      =   'ecmwf'   ,        6          , 'rh'     ,    .true.    , .false., 'daily'   ,'bilinear.nc',   ''     ,   ''
373   sn_prec     =   'ecmwf'   ,        6          , 'precip' ,    .true.    , .true. , 'daily'   ,'bicubic.nc' ,   ''     ,   ''
374
375   cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files
376/
377!-----------------------------------------------------------------------
378&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
379!-----------------------------------------------------------------------
380!                    !     description      !  multiple  !    vector   !      vector          ! vector !
381!                    !                      ! categories !  reference  !    orientation       ! grids  !
382! send
383   sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   ''
384   sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
385   sn_snd_thick  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
386   sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T'
387   sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
388! receive
389   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
390   sn_rcv_taumod =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
391   sn_rcv_tau    =   'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V'
392   sn_rcv_dqnsdt =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
393   sn_rcv_qsr    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
394   sn_rcv_qns    =   'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''
395   sn_rcv_emp    =   'conservative'         ,    'no'    ,     ''      ,         ''          ,   ''
396   sn_rcv_rnf    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
397   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
398   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''
399!
400   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentialy sending/receiving data
401   ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models
402   !                       !   -> file cplmask.nc with the float variable called cplmask (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                          (ln_traqsr=T)
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_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration
430   ln_qsr_2bd  = .false.   !  2 bands              light penetration
431   ln_qsr_bio  = .false.   !  bio-model light penetration
432   nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0)
433   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
434   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
435   rn_si1      =   23.0    !  2 bands: longest depth of extinction
436   ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3
437/
438!-----------------------------------------------------------------------
439&namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf=T)
440!-----------------------------------------------------------------------
441!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
442!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
443   sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
444   sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , ''
445   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
446   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
447   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , ''
448
449   cn_dir      = './'      !  root directory for the location of the runoff files
450   ln_rnf_mouth= .true.    !  specific treatment at rivers mouths
451      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T)
452      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
453   rn_rfact    =   1.e0    !  multiplicative factor for runoff
454   ln_rnf_depth= .false.   !  read in depth information for runoff
455   ln_rnf_tem  = .false.   !  read in temperature information for runoff
456   ln_rnf_sal  = .false.   !  read in salinity information for runoff
457   ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file
458      rn_rnf_max  = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
459      rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
460      nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0)
461/
462!-----------------------------------------------------------------------
463&namsbc_isf    !  Top boundary layer (ISF)                              (nn_isf >0)
464!-----------------------------------------------------------------------
465!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
466!              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
467! nn_isf == 4
468   sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
469! nn_isf == 3
470   sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
471! nn_isf == 2 and 3
472   sn_depmax_isf='rnfisf'  ,         -12       ,'sozisfmax',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
473   sn_depmin_isf='rnfisf'  ,         -12       ,'sozisfmin',  .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
474! nn_isf == 2
475   sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,   .false.   , .true. , 'yearly'  ,    ''    ,   ''     ,    ''
476!
477! for all case
478   nn_isf      = 1         !  ice shelf melting/freezing
479                           !  1 = presence of ISF    2 = bg03 parametrisation
480                           !  3 = rnf file for isf   4 = ISF fwf specified
481                           !  option 1 and 4 need ln_isfcav = .true. (domzgr)
482! only for nn_isf = 1 or 2
483   rn_gammat0  = 1.e-4     ! gammat coefficient used in blk formula
484   rn_gammas0  = 1.e-4     ! gammas coefficient used in blk formula
485! only for nn_isf = 1 or 4
486   rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008)
487   !                       ! 0 => thickness of the tbl = thickness of the first wet cell
488! only for nn_isf = 1
489   nn_isfblk   = 1         ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006)
490   !                       ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015)
491   nn_gammablk = 1         ! 0 = cst Gammat (= gammat/s)
492   !                       ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010)
493   !                       ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999)
494/
495!-----------------------------------------------------------------------
496&namsbc_iscpl  !   land ice / ocean coupling option                     
497!-----------------------------------------------------------------------
498   nn_drown    = 10        ! number of iteration of the extrapolation loop (fill the new wet cells)
499   ln_hsb      = .false.   ! activate conservation module (conservation exact after a time of rn_fiscpl)
500   nn_fiscpl   = 43800     ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency)
501/
502!-----------------------------------------------------------------------
503&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk
504!-----------------------------------------------------------------------
505!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
506!              !           !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
507   sn_apr      = 'patm'    ,         -1        ,'somslpre',   .true.    , .true. , 'yearly'  ,    ''    ,    ''    ,      ''
508
509   cn_dir      = './'      !  root directory for the location of the bulk files
510   rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/
511   ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F)
512   ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data
513/
514!-----------------------------------------------------------------------
515&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr=T)
516!-----------------------------------------------------------------------
517!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
518!              !           !  (if <0  months)  !   name   !   (logical) !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
519   sn_sst      = 'sst_data',        24         ,  'sst'   ,    .false.  , .false., 'yearly'  ,    ''    ,    ''    ,     ''
520   sn_sss      = 'sss_data',        -1         ,  'sss'   ,    .true.   , .true. , 'yearly'  ,    ''    ,    ''    ,     ''
521
522   cn_dir      = './'      !  root directory for the location of the runoff files
523   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
524   nn_sssr     =     2     !  add a damping     term in the surface freshwater flux (=2)
525                           !  or to SSS only (=1) or no damping term (=0)
526   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K]
527   rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
528   ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2)
529   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
530/
531!-----------------------------------------------------------------------
532&namsbc_alb    !   albedo parameters
533!-----------------------------------------------------------------------
534   nn_ice_alb  =    0      !  parameterization of ice/snow albedo
535                           !     0: Shine & Henderson-Sellers (JGR 1985)
536                           !     1: "home made" based on Brandt et al. (J. Climate 2005)
537                           !                         and Grenfell & Perovich (JGR 2004)
538   rn_albice   =  0.53     !  albedo of bare puddled ice (values from 0.49 to 0.58)
539                           !     0.53 (default) => if nn_ice_alb=0
540                           !     0.50 (default) => if nn_ice_alb=1
541/
542!-----------------------------------------------------------------------
543&namsbc_wave   ! External fields from wave model                        (ln_wave=T)
544!-----------------------------------------------------------------------
545!              ! file name ! frequency (hours) ! variable    ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
546!              !           !  (if <0  months)  !   name      !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
547   sn_cdg      = 'cdg_wave',        1          , 'drag_coeff',   .true.    , .false., 'daily'   ,    ''      , ''       , ''
548   sn_usd      = 'sdw_wave',        1          , 'u_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , ''
549   sn_vsd      = 'sdw_wave',        1          , 'v_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , ''
550   sn_wn       = 'sdw_wave',        1          , 'wave_num'  ,   .true.    , .false., 'daily'   ,    ''      , ''       , ''
551!
552   cn_dir_cdg  = './'      !  root directory for the location of drag coefficient files
553   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model
554   ln_sdw      = .false.   !  Computation of 3D stokes drift               
555/
556!-----------------------------------------------------------------------
557&namberg       !   iceberg parameters                                   (default: No iceberg)
558!-----------------------------------------------------------------------
559      ln_icebergs              = .false.              ! iceberg floats or not
560      ln_bergdia               = .true.               ! Calculate budgets
561      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0
562      nn_verbose_write         = 15                   ! Timesteps between verbose messages
563      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage
564                                                      ! Initial mass required for an iceberg of each class
565      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
566                                                      ! Proportion of calving mass to apportion to each class
567      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
568                                                      ! Ratio between effective and real iceberg mass (non-dim)
569                                                      ! i.e. number of icebergs represented at a point
570      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1
571                                                      ! thickness of newly calved bergs (m)
572      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
573      rn_rho_bergs             = 850.                 ! Density of icebergs
574      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs
575      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics
576      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits
577      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
578      ln_passive_mode          = .false.              ! iceberg - ocean decoupling
579      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no)
580                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
581      rn_test_box              = 108.0,  116.0, -66.0, -58.0
582      rn_speed_limit           = 0.                   ! CFL speed limit for a berg
583
584!            ! file name ! frequency (hours) !   variable   ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
585!            !           !  (if <0  months)  !     name     !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
586      sn_icb =  'calving',       -1          , 'calvingmask',   .true.    , .true.  , 'yearly'  ,    ''    ,    ''    ,     ''
587
588      cn_dir = './'
589/
590
591!!======================================================================
592!!               ***  Lateral boundary condition  ***
593!!======================================================================
594!!   namlbc        lateral momentum boundary condition
595!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")
596!!   nam_tide      Tidal forcing
597!!   nambdy        Unstructured open boundaries                         ("key_bdy")
598!!   nambdy_dta    Unstructured open boundaries - external data         ("key_bdy")
599!!   nambdy_tide   tidal forcing at open boundaries                     ("key_bdy_tides")
600!!======================================================================
601!
602!-----------------------------------------------------------------------
603&namlbc        !   lateral momentum boundary condition
604!-----------------------------------------------------------------------
605   !                       !  free slip  !   partial slip  !   no slip   ! strong slip
606   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
607   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs.
608/
609!-----------------------------------------------------------------------
610&namagrif      !  AGRIF zoom                                            ("key_agrif")
611!-----------------------------------------------------------------------
612   nn_cln_update =    3    !  baroclinic update frequency
613   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
614   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
615   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
616   ln_chk_bathy  = .FALSE. !
617/
618!-----------------------------------------------------------------------
619&nam_tide      !   tide parameters                                      ("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 lim3
673!   bn_a_i     = 'amm12_bdyT_ice',         24        , 'ileadfra',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
674!   bn_ht_i    = 'amm12_bdyT_ice',         24        , 'iicethic',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
675!   bn_ht_s    = 'amm12_bdyT_ice',         24        , 'isnowthi',    .true.   , .false. ,  'daily'  ,    ''    ,   ''     ,     ''
676
677   cn_dir      = 'bdydta/' !  root directory for the location of the bulk files
678   ln_full_vel = .false.   ! 
679/
680!-----------------------------------------------------------------------
681&nambdy_tide   !  tidal forcing at open boundaries
682!-----------------------------------------------------------------------
683   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files
684   ln_bdytide_2ddta = .false.                   !
685   ln_bdytide_conj  = .false.                   !
686/
687
688!!======================================================================
689!!                 ***  Bottom boundary condition  ***
690!!======================================================================
691!!   nambfr        bottom friction
692!!   nambbc        bottom temperature boundary condition
693!!   nambbl        bottom boundary layer scheme                         ("key_trabbl")
694!!======================================================================
695!
696!-----------------------------------------------------------------------
697&nambfr        !   bottom friction                                      (default: linear)
698!-----------------------------------------------------------------------
699   nn_bfr      =    1      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
700                           !                              = 2 : nonlinear friction
701   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case)
702   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
703   rn_bfri2_max=    1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T)
704   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2)
705   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T
706   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file )
707   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T)
708   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case)
709   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T
710   rn_tfri2_max=    1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T)
711   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2)
712   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T
713   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file )
714   rn_tfrien   =   50.     !  local multiplying factor of tfr (ln_tfr2d=T)
715
716   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true)
717   ln_loglayer = .false.   !  logarithmic formulation (non linear case)
718/
719!-----------------------------------------------------------------------
720&nambbc        !   bottom temperature boundary condition                (default: NO)
721!-----------------------------------------------------------------------
722!              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
723!              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
724   sn_qgh      ='geothermal_heating.nc',  -12.       , 'heatflow',   .false.   , .true.  , 'yearly'  ,   ''     ,   ''     ,   ''
725   !
726   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
727   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux
728                           !     = 1 constant flux
729                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
730   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
731   cn_dir      = './'      !  root directory for the location of the runoff files
732/
733!-----------------------------------------------------------------------
734&nambbl        !   bottom boundary layer scheme                         ("key_trabbl")
735!-----------------------------------------------------------------------
736   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0)
737   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0)
738   rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s]
739   rn_gambbl   =  10.      !  advective bbl coefficient                 [s]
740/
741
742!!======================================================================
743!!                        Tracer (T & S ) namelists
744!!======================================================================
745!!   nameos           equation of state
746!!   namtra_adv       advection scheme
747!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.)
748!!   namtra_ldf       lateral diffusion scheme
749!!   namtra_ldfeiv    eddy induced velocity param.
750!!   namtra_dmp       T & S newtonian damping
751!!======================================================================
752!
753!-----------------------------------------------------------------------
754&nameos        !   ocean physical parameters
755!-----------------------------------------------------------------------
756   ln_teos10   = .false.         !  = Use TEOS-10 equation of state
757   ln_eos80    = .false.         !  = Use EOS80 equation of state
758   ln_seos     = .false.         !  = Use simplified equation of state (S-EOS)
759                                 !
760   !                     ! S-EOS coefficients (ln_seos=T):
761   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
762   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1)
763   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1)
764   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
765   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
766   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
767   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
768   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
769/
770!-----------------------------------------------------------------------
771&namtra_adv    !   advection scheme for tracer                          (default: NO advection)
772!-----------------------------------------------------------------------
773   ln_traadv_cen = .false. !  2nd order centered scheme
774      nn_cen_h   =  4            !  =2/4, horizontal 2nd order CEN / 4th order CEN
775      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT
776   ln_traadv_fct = .false. !  FCT scheme
777      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order
778      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order
779      nn_fct_zts =  0            !  >=1,  2nd order FCT scheme with vertical sub-timestepping
780      !                          !        (number of sub-timestep = nn_fct_zts)
781   ln_traadv_mus = .false. !  MUSCL scheme
782      ln_mus_ups = .false.       !  use upstream scheme near river mouths
783   ln_traadv_ubs = .false. !  UBS scheme
784      nn_ubs_v   =  2            !  =2  , vertical 2nd order FCT / COMPACT 4th order
785   ln_traadv_qck = .false. !  QUICKEST scheme
786/
787!-----------------------------------------------------------------------
788&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) (default: NO)
789!-----------------------------------------------------------------------
790   ln_mle      = .false.   ! (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                 (default: NO diffusion)
802!-----------------------------------------------------------------------
803   !                       !  Operator type:
804   !                           !  no diffusion: set ln_traldf_lap=..._blp=F
805   ln_traldf_lap   =  .false.  !    laplacian operator
806   ln_traldf_blp   =  .false.  !  bilaplacian operator
807   !
808   !                       !  Direction of action:
809   ln_traldf_lev   =  .false.  !  iso-level
810   ln_traldf_hor   =  .false.  !  horizontal (geopotential)
811   ln_traldf_iso   =  .false.  !  iso-neutral (standard operator)
812   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator)
813   !
814   !                       !  iso-neutral options:       
815   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators)
816   rn_slpmax       =   0.01    !  slope limit                      (both operators)
817   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only)
818   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only)
819   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only)
820   !
821   !                       !  Coefficients:
822   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef
823   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file
824   !                                !   =  0           constant
825   !                                !   = 10 F(k)      =ldf_c1d
826   !                                !   = 20 F(i,j)    =ldf_c2d
827   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
828   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
829   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
830   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s]
831   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s]
832/
833!-----------------------------------------------------------------------
834&namtra_ldfeiv !   eddy induced velocity param.                         (default: NO)
835!-----------------------------------------------------------------------
836   ln_ldfeiv     =.false.  ! use eddy induced velocity parameterization
837   ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities
838   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s]
839   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient
840   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
841   !                                !   =  0           constant
842   !                                !   = 10 F(k)      =ldf_c1d
843   !                                !   = 20 F(i,j)    =ldf_c2d
844   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
845   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d
846/
847!-----------------------------------------------------------------------
848&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO)
849!-----------------------------------------------------------------------
850   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F)
851   nn_zdmp     =    0      !  vertical   shape =0    damping throughout the water column
852                           !                   =1 no damping in the mixing layer (kz  criteria)
853                           !                   =2 no damping in the mixed  layer (rho crieria)
854   cn_resto    ='resto.nc' !  Name of file containing restoration coeff. field (use dmp_tools to create this)
855/
856
857!!======================================================================
858!!                      ***  Dynamics namelists  ***
859!!======================================================================
860!!   namdyn_adv    formulation of the momentum advection
861!!   namdyn_vor    advection scheme
862!!   namdyn_hpg    hydrostatic pressure gradient
863!!   namdyn_spg    surface pressure gradient
864!!   namdyn_ldf    lateral diffusion scheme
865!!======================================================================
866!
867!-----------------------------------------------------------------------
868&namdyn_adv    !   formulation of the momentum advection                (default: vector form)
869!-----------------------------------------------------------------------
870   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)
871   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction
872   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
873   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
874   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection
875/
876!-----------------------------------------------------------------------
877&nam_vvl    !   vertical coordinate options                             (default: zstar)
878!-----------------------------------------------------------------------
879   ln_vvl_zstar  = .true.           !  zstar vertical coordinate
880   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations
881   ln_vvl_layer  = .false.          !  full layer vertical coordinate
882   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
883   ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator
884   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient
885   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days]
886   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days]
887   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation
888   ln_vvl_dbg    = .true.           !  debug prints    (T/F)
889/
890!-----------------------------------------------------------------------
891&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO)
892!-----------------------------------------------------------------------
893   ln_dynvor_ene = .false. !  enstrophy conserving scheme
894   ln_dynvor_ens = .false. !  energy conserving scheme
895   ln_dynvor_mix = .false. !  mixed scheme
896   ln_dynvor_een = .false. !  energy & enstrophy scheme
897      nn_een_e3f = 1          ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)
898   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes)  ! PLEASE DO NOT ACTIVATE
899/
900!-----------------------------------------------------------------------
901&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: zps)
902!-----------------------------------------------------------------------
903   ln_hpg_zco  = .false.   !  z-coordinate - full steps
904   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation)
905   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
906   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
907   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
908   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
909/
910!-----------------------------------------------------------------------
911&namdyn_spg    !   surface pressure gradient                            (default: NO)
912!-----------------------------------------------------------------------
913   ln_dynspg_exp  = .false.   ! explicit free surface
914   ln_dynspg_ts   = .false.   ! split-explicit free surface
915      ln_bt_fw      = .true.     ! Forward integration of barotropic Eqs.
916      ln_bt_av      = .true.     ! Time filtering of barotropic variables
917         nn_bt_flt     = 1          ! Time filter choice  = 0 None
918         !                          !                     = 1 Boxcar over   nn_baro sub-steps
919         !                          !                     = 2 Boxcar over 2*nn_baro  "    "
920      ln_bt_auto    = .true.     ! Number of sub-step defined from:
921         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed
922         nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds
923/
924!-----------------------------------------------------------------------
925&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO)
926!-----------------------------------------------------------------------
927   !                       !  Type of the operator :
928   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F
929   ln_dynldf_lap =  .false.    !    laplacian operator
930   ln_dynldf_blp =  .false.    !  bilaplacian operator
931   !                       !  Direction of action  :
932   ln_dynldf_lev =  .false.    !  iso-level
933   ln_dynldf_hor =  .false.    !  horizontal (geopotential)
934   ln_dynldf_iso =  .false.    !  iso-neutral
935   !                       !  Coefficient
936   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef
937   !                                !  =-30  read in eddy_viscosity_3D.nc file
938   !                                !  =-20  read in eddy_viscosity_2D.nc file
939   !                                !  =  0  constant
940   !                                !  = 10  F(k)=c1d
941   !                                !  = 20  F(i,j)=F(grid spacing)=c2d
942   !                                !  = 30  F(i,j,k)=c2d*c1d
943   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity)
944   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s]
945   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s]
946   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s]
947   !
948   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km)
949/
950
951!!======================================================================
952!!             Tracers & Dynamics vertical physics namelists
953!!======================================================================
954!!    namzdf        vertical physics
955!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric")
956!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke")
957!!    namzdf_gls    GLS vertical mixing                                 ("key_zdfgls")
958!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm")
959!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx")
960!!======================================================================
961!
962!-----------------------------------------------------------------------
963&namzdf        !   vertical physics
964!-----------------------------------------------------------------------
965   rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
966   rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
967   nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0)
968   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0)
969   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F)
970      nn_evdm     =    0        ! evd apply on tracer (=0) or on tracer and momentum (=1)
971      rn_avevd    =  100.       !  evd mixing coefficient [m2/s]
972   ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F)
973      nn_npc      =    1        ! frequency of application of npc
974      nn_npcp     =  365        ! npc control print frequency
975   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
976      nn_zdfexp   =    3        ! number of sub-timestep for ln_zdfexp=T
977/
978!-----------------------------------------------------------------------
979&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
980!-----------------------------------------------------------------------
981   rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity
982   rn_alp      =    5.     !  coefficient of the parameterization
983   nn_ric      =    2      !  coefficient of the parameterization
984   rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation
985   rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m)
986   rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m)
987   rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
988   rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
989   ln_mldw     =  .true.   !  Flag to use or not the mixed layer depth param.
990/
991!-----------------------------------------------------------------------
992&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
993!-----------------------------------------------------------------------
994   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
995   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
996   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
997   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
998   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
999   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
1000   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
1001                           !                 = 1 bounded by the local vertical scale factor
1002                           !                 = 2 first vertical derivative of mixing length bounded by 1
1003                           !                 = 3 as =2 with distinct disspipative an mixing length scale
1004   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
1005   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
1006   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
1007   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
1008   rn_lc       =   0.15    !  coef. associated to Langmuir cells
1009   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to near intertial waves
1010                           !        = 0 no penetration
1011                           !        = 1 add a tke source below the ML
1012                           !        = 2 add a tke source just at the base of the ML
1013                           !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T)
1014   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
1015   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
1016                           !        = 0  constant 10 m length scale
1017                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees
1018/
1019!-----------------------------------------------------------------------
1020&namzdf_gls    !   GLS vertical diffusion                               ("key_zdfgls")
1021!-----------------------------------------------------------------------
1022   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
1023   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
1024   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
1025   rn_clim_galp  = 0.267   !  galperin limit
1026   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
1027   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
1028   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
1029   rn_hsro       =  0.02   !  Minimum surface roughness
1030   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met=2)
1031   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2)
1032   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
1033   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
1034   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
1035   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
1036/
1037!-----------------------------------------------------------------------
1038&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
1039!-----------------------------------------------------------------------
1040   rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
1041   rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
1042/
1043!-----------------------------------------------------------------------
1044&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
1045!-----------------------------------------------------------------------
1046   rn_htmx     = 500.      !  vertical decay scale for turbulence (meters)
1047   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1)
1048   rn_tfe      = 0.333     !  tidal dissipation efficiency
1049   rn_me       = 0.2       !  mixing efficiency
1050   ln_tmx_itf  = .true.    !  ITF specific parameterisation
1051   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency
1052/
1053!-----------------------------------------------------------------------
1054&namzdf_tmx_new !   internal wave-driven mixing parameterization        ("key_zdftmx_new" & "key_zdfddm")
1055!-----------------------------------------------------------------------
1056   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
1057   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency
1058   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F)
1059/
1060
1061
1062!!======================================================================
1063!!                  ***  Miscellaneous namelists  ***
1064!!======================================================================
1065!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi)
1066!!   namctl            Control prints & Benchmark
1067!!   namsto            Stochastic parametrization of EOS
1068!!======================================================================
1069!
1070!-----------------------------------------------------------------------
1071&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
1072!-----------------------------------------------------------------------
1073   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
1074                           !  buffer blocking send or immediate non-blocking sends, resp.
1075   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
1076   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold
1077   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1)
1078   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1)
1079   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1)
1080/
1081!-----------------------------------------------------------------------
1082&namctl        !   Control prints & Benchmark
1083!-----------------------------------------------------------------------
1084   ln_ctl      = .false.   !  trends control print (expensive!)
1085   nn_print    =    0      !  level of print (0 no extra print)
1086   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
1087   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
1088   nn_jctls    =    0      !  start j indice of control               over a subdomain)
1089   nn_jctle    =    0      !  end   j indice of control
1090   nn_isplt    =    1      !  number of processors in i-direction
1091   nn_jsplt    =    1      !  number of processors in j-direction
1092   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench
1093                           !     (no physical validity of the results)
1094   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0)
1095   nn_diacfl   =    0      !  Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0)
1096/
1097!-----------------------------------------------------------------------
1098&namsto        ! Stochastic parametrization of EOS                      (default: NO)
1099!-----------------------------------------------------------------------
1100   ln_sto_eos  = .false.   ! stochastic equation of state
1101   nn_sto_eos  = 1         ! number of independent random walks
1102   rn_eos_stdxy= 1.4       ! random walk horz. standard deviation (in grid points)
1103   rn_eos_stdz = 0.7       ! random walk vert. standard deviation (in grid points)
1104   rn_eos_tcor = 1440.     ! random walk time correlation (in timesteps)
1105   nn_eos_ord  = 1         ! order of autoregressive processes
1106   nn_eos_flt  = 0         ! passes of Laplacian filter
1107   rn_eos_lim  = 2.0       ! limitation factor (default = 3.0)
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
1114!!======================================================================
1115!!                  ***  Diagnostics namelists  ***
1116!!======================================================================
1117!!   namtrd       dynamics and/or tracer trends                         (default F)
1118!!   namptr       Poleward Transport Diagnostics                        (default F)
1119!!   namhsb       Heat and salt budgets                                 (default F)
1120!!   namdiu       Cool skin and warm layer models                       (default F)
1121!!   namflo       float parameters                                      ("key_float")
1122!!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm")
1123!!   namdct       transports through some sections                      ("key_diadct")
1124!!   nam_diatmb   Top Middle Bottom Output                              (default F)
1125!!   nam_dia25h   25h Mean Output                                       (default F)
1126!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1127!!======================================================================
1128!
1129!-----------------------------------------------------------------------
1130&namtrd        !   trend diagnostics                                    (default F)
1131!-----------------------------------------------------------------------
1132   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1133   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1134   ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1135   ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet)
1136   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1137   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1138   ln_tra_trd  = .false.   ! (T) 3D tracer trend output
1139   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1140   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1141/
1142!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1143!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1144!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1145!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1146!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1147!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1148!!gm
1149!-----------------------------------------------------------------------
1150&namptr        !   Poleward Transport Diagnostic                         (default F)
1151!-----------------------------------------------------------------------
1152   ln_diaptr   = .false.   !  Poleward heat and salt transport (T) or not (F)
1153   ln_subbas   = .false.   !  Atlantic/Pacific/Indian basins computation (T) or not
1154/
1155!-----------------------------------------------------------------------
1156&namhsb        !  Heat and salt budgets                                  (default F)
1157!-----------------------------------------------------------------------
1158   ln_diahsb   = .false.   !  check the heat and salt budgets (T) or not (F)
1159/
1160!-----------------------------------------------------------------------
1161&namdiu        !   Cool skin and warm layer models                       (default F)
1162!-----------------------------------------------------------------------
1163   ln_diurnal      = .false.   !
1164   ln_diurnal_only = .false.   !
1165/
1166!-----------------------------------------------------------------------
1167&namflo        !   float parameters                                      ("key_float")
1168!-----------------------------------------------------------------------
1169   jpnfl       = 1         !  total number of floats during the run
1170   jpnnewflo   = 0         !  number of floats for the restart
1171   ln_rstflo   = .false.   !  float restart (T) or not (F)
1172   nn_writefl  =      75   !  frequency of writing in float output file
1173   nn_stockfl  =    5475   !  frequency of creation of the float restart file
1174   ln_argo     = .false.   !  Argo type floats (stay at the surface each 10 days)
1175   ln_flork4   = .false.   !  trajectories computed with a 4th order Runge-Kutta (T)
1176   !                       !  or computed with Blanke' scheme (F)
1177   ln_ariane   = .true.    !  Input with Ariane tool convention(T)
1178   ln_flo_ascii= .true.    !  Output with Ariane tool netcdf convention(F) or ascii file (T)
1179/
1180!-----------------------------------------------------------------------
1181&nam_diaharm   !   Harmonic analysis of tidal constituents               ("key_diaharm")
1182!-----------------------------------------------------------------------
1183    nit000_han = 1         ! First time step used for harmonic analysis
1184    nitend_han = 75        ! Last time step used for harmonic analysis
1185    nstep_han  = 15        ! Time step frequency for harmonic analysis
1186    tname(1)   = 'M2'      ! Name of tidal constituents
1187    tname(2)   = 'K1'
1188/
1189!-----------------------------------------------------------------------
1190&namdct        ! transports through some sections                        ("key_diadct")
1191!-----------------------------------------------------------------------
1192    nn_dct     = 15        !  time step frequency for transports computing
1193    nn_dctwri  = 15        !  time step frequency for transports writing
1194    nn_secdebug= 112       !      0 : no section to debug
1195    !                      !     -1 : debug all section
1196    !                      !  0 < n : debug section number n
1197/
1198!-----------------------------------------------------------------------
1199&nam_diatmb    !  Top Middle Bottom Output                               (default F)
1200!-----------------------------------------------------------------------
1201   ln_diatmb   = .false.   !  Choose Top Middle and Bottom output or not
1202/
1203!-----------------------------------------------------------------------
1204&nam_dia25h    !  25h Mean Output                                        (default F)
1205!-----------------------------------------------------------------------
1206   ln_dia25h   = .false.   ! Choose 25h mean output or not
1207/
1208!-----------------------------------------------------------------------
1209&namnc4        !   netcdf4 chunking and compression settings             ("key_netcdf4")
1210!-----------------------------------------------------------------------
1211   nn_nchunks_i=   4       !  number of chunks in i-dimension
1212   nn_nchunks_j=   4       !  number of chunks in j-dimension
1213   nn_nchunks_k=   31      !  number of chunks in k-dimension
1214   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1215   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs
1216   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1217   !                       !  (F) ignore chunking information and produce netcdf3-compatible files
1218/
1219
1220!!======================================================================
1221!!               ***  Observation & Assimilation  ***
1222!!======================================================================
1223!!   namobs       observation and model comparison
1224!!   nam_asminc   assimilation increments                               ('key_asminc')
1225!!======================================================================
1226!
1227!-----------------------------------------------------------------------
1228&namobs        !  observation usage switch
1229!-----------------------------------------------------------------------
1230   ln_diaobs   = .false.             ! Logical switch for the observation operator
1231   ln_t3d      = .false.             ! Logical switch for T profile observations
1232   ln_s3d      = .false.             ! Logical switch for S profile observations
1233   ln_sla      = .false.             ! Logical switch for SLA observations
1234   ln_sst      = .false.             ! Logical switch for SST observations
1235   ln_sic      = .false.             ! Logical switch for Sea Ice observations
1236   ln_vel3d    = .false.             ! Logical switch for velocity observations
1237   ln_altbias  = .false.             ! Logical switch for altimeter bias correction
1238   ln_nea      = .false.             ! Logical switch for rejection of observations near land
1239   ln_grid_global = .true.           ! Logical switch for global distribution of observations
1240   ln_grid_search_lookup = .false.   ! Logical switch for obs grid search w/lookup table
1241   ln_ignmis   = .true.              ! Logical switch for ignoring missing files
1242   ln_s_at_t   = .false.             ! Logical switch for computing model S at T obs if not there
1243   ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs
1244! All of the *files* variables below are arrays. Use namelist_cfg to add more files
1245   cn_profbfiles = 'profiles_01.nc'  ! Profile feedback input observation file names
1246   cn_slafbfiles = 'sla_01.nc'       ! SLA feedback input observation file names
1247   cn_sstfbfiles = 'sst_01.nc'       ! SST feedback input observation file names
1248   cn_sicfbfiles = 'sic_01.nc'       ! SIC feedback input observation file names
1249   cn_velfbfiles = 'vel_01.nc'       ! Velocity feedback input observation file names
1250   cn_altbiasfile = 'altbias.nc'     ! Altimeter bias input file name
1251   cn_gridsearchfile='gridsearch.nc' ! Grid search file name
1252   rn_gridsearchres = 0.5            ! Grid search resolution
1253   rn_dobsini  = 00010101.000000     ! Initial date in window YYYYMMDD.HHMMSS
1254   rn_dobsend  = 00010102.000000     ! Final date in window YYYYMMDD.HHMMSS
1255   nn_1dint    = 0                   ! Type of vertical interpolation method
1256   nn_2dint    = 0                   ! Type of horizontal interpolation method
1257   nn_msshc    = 0                   ! MSSH correction scheme
1258   rn_mdtcorr  = 1.61                ! MDT  correction
1259   rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction
1260   nn_profdavtypes = -1              ! Profile daily average types - array
1261   ln_sstbias  = .false.             !
1262   cn_sstbias_files = 'sstbias.nc'   !
1263/
1264!-----------------------------------------------------------------------
1265&nam_asminc    !   assimilation increments                              ('key_asminc')
1266!-----------------------------------------------------------------------
1267    ln_bkgwri  = .false.   !  Logical switch for writing out background state
1268    ln_trainc  = .false.   !  Logical switch for applying tracer increments
1269    ln_dyninc  = .false.   !  Logical switch for applying velocity increments
1270    ln_sshinc  = .false.   !  Logical switch for applying SSH increments
1271    ln_asmdin  = .false.   !  Logical switch for Direct Initialization (DI)
1272    ln_asmiau  = .false.   !  Logical switch for Incremental Analysis Updating (IAU)
1273    nitbkg     = 0         !  Timestep of background in [0,nitend-nit000-1]
1274    nitdin     = 0         !  Timestep of background for DI in [0,nitend-nit000-1]
1275    nitiaustr  = 1         !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1276    nitiaufin  = 15        !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1277    niaufn     = 0         !  Type of IAU weighting function
1278    ln_salfix  = .false.   !  Logical switch for ensuring that the sa > salfixmin
1279    salfixmin  = -9999     !  Minimum salinity after applying the increments
1280    nn_divdmp  = 0         !  Number of iterations of divergence damping operator
1281/
Note: See TracBrowser for help on using the repository browser.