source: trunk/CONFIG/GYRE_LOBSTER/EXP00/namelist @ 907

Last change on this file since 907 was 907, checked in by ctlod, 13 years ago

Re-organisation of ocean namelists to improve readability, see ticket: #125

File size: 42.0 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  :       1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl)
3!! namelists         2 - Domain           (nam_zgr, nam_zgr_sco, namdom)
4!!                   3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core
5!!                                         namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb)
6!!                   4 - lateral boundary (namlbc, namcla, namobc)
7!!                   5 - bottom  boundary (nambfr, nambbc, nambbl)
8!!                   6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp)
9!!                   7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf)
10!!                   8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm)
11!!                   9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr)
12!!                   9 - miscellaneous    (namsol)
13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
14!  CAUTION: some scripts does not support CAPITALs for logical definition
15!  *******       use .true. or .false. and NOT .TRUE. or .FALSE.
16
17!!======================================================================
18!!                   ***  Run management namelists  ***
19!!======================================================================
20!!   namrun            parameters of the run
21!!   nam_mpp           Massively Parallel Processing
22!!   nam_mpp_dyndist   
23!!   nam_ctl           Control prints & Benchmark
24!!======================================================================
25
26!-----------------------------------------------------------------------
27!       namrun   parameters of the run
28!-----------------------------------------------------------------------
29!  no         job number
30!  cexper     experience name for vairmer format
31!  ln_rstart  boolean term for restart (true or false)
32!  nrstdt     restart control = 0 restart, do not control nit000 in the restart file.
33!             !               = 1 restart, control nit000 in the restart file. Do not
34!             !                   use the date in the restart file (use ndate0 in namelist)
35!             !               = 2 restart, control nit000 in the restart file, use the date
36!             !                   in the restart file. ndate0 in the namelist is ignored.
37!  nit000     number of the first time step
38!  nitend     number of the last time step
39!  ndate0     initial calendar date aammjj
40!  nleapy     Leap year calendar (0/1)
41!  ninist     initial state output flag (0/1)
42!  nstock     frequency of restart file
43!  nwrite     frequency of OUTPUT file
44!  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc)
45&namrun
46   no         =       0
47   cexper     =  "GYRE"
48   ln_rstart  = .false.
49   nrstdt     =       0
50   nit000     =       1
51   nitend     =    4320
52   ndate0     =  010101
53   nleapy     =      30
54   ninist     =       0
55   nstock     =    4320 
56   nwrite     =     360 
57   ln_dimgnnn = .false.
58/
59!-----------------------------------------------------------------------
60!       nam_mpp      Massively Parallel Processing
61!-----------------------------------------------------------------------
62!  c_mpi_send         mpi send/recieve type
63!                      = 'S'  : standard blocking send
64!                      = 'B'  : buffer blocking send
65!                      = 'I'  : immediate non-blocking send
66&nam_mpp
67   c_mpi_send =  'S'
68/
69!-----------------------------------------------------------------------
70!       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist")
71!-----------------------------------------------------------------------
72!  jpni   number of processors following i
73!  jpnj   number of processors following j
74!  jpnij  number of local domains
75&nam_mpp_dyndist
76   jpni  = 1
77   jpnj  = 1
78   jpnij = 1
79/
80!-----------------------------------------------------------------------
81!       nam_ctl      Control prints & Benchmark
82!-----------------------------------------------------------------------
83!  ln_ctl     trends control print (expensive!)
84!  nprint     level of print (0 no print)
85!  nictls     start i indice to make the control SUM (very usefull to compare mono-
86!  nictle     end   i indice to make the control SUM (-versus multi processor runs)
87!  njctls     start j indice to make the control SUM (very usefull to compare mono-
88!  njctle     end   j indice to make the control SUM (-versus multi processor runs)
89!  nisplt     number of processors following i
90!  njsplt     number of processors following j
91!  nbench     Bench parameter (0/1): CAUTION it must be zero except for bench
92!             for which we don't care about physical meaning of the results
93!  nbit_cmp   bit comparison mode parameter (0/1): enables bit comparison between
94!             single and multiple processor runs.
95&namctl
96   ln_ctl =  .false.
97   nprint =       0
98   nictls =       0
99   nictle =       0
100   njctls =       0
101   njctle =       0
102   isplt  =       1
103   jsplt  =       1
104   nbench =       0
105   nbit_cmp =     0
106/
107
108!!======================================================================
109!!                      ***  Domain namelists  ***
110!!======================================================================
111!!   nam_zgr       vertical coordinate
112!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
113!!   namdom        space and time domain (bathymetry, mesh, timestep)
114!!======================================================================
115
116!-----------------------------------------------------------------------
117!       nam_zgr       vertical coordinate
118!-----------------------------------------------------------------------
119!  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined)
120!  ln_zps     z-coordinate - partial steps   (T/F)
121!  ln_sco     s- or hybrid z-s-coordinate    (T/F)
122&nam_zgr
123   ln_zco   =  .true.
124   ln_zps   =  .false.
125   ln_sco   =  .false.
126/
127!-----------------------------------------------------------------------
128!       nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
129!-----------------------------------------------------------------------
130!  sbot_min   minimum depth of s-bottom surface (>0) (m)
131!  sbot_max   maximum depth of s-bottom surface (= ocean depth) (>0) (m)
132!  theta      surface control parameter (0<=theta<=20)
133!  thetb      bottom control parameter  (0<=thetb<= 1)
134!  r_max      maximum cut-off r-value allowed (0<r_max<1)
135&nam_zgr_sco
136   sbot_min =  300.
137   sbot_max = 5250.
138   theta    =    6.0
139   thetb    =    0.75
140   r_max    =    0.15
141/
142!-----------------------------------------------------------------------
143!       namdom   space and time domain (bathymetry, mesh, timestep)
144!-----------------------------------------------------------------------
145!  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels)
146!  e3zps_min  the thickness of the partial step is set larger than the
147!  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1)
148!  nmsh       =1 create a mesh file (coordinates, scale factors, masks)
149!  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra
150!                                                  = 1,     used, rdt < rdttra(k)
151!  atfp       asselin time filter parameter
152!  rdt        time step for the dynamics (and tracer if nacc=0)
153!  rdtmin     minimum time step on tracers
154!  rdtmax     maximum time step on tracers
155!  rdth       depth variation of tracer time step
156!  rdtbt      barotropic time step (for the time splitting algorithm)
157!  nclosea    =0 remove the closed sea from the global domain (orca configuration)
158!             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)
159&namdom
160   ntopo     =     0
161   e3zps_min =     5.
162   e3zps_rat =     0.1
163   nmsh      =     0
164   nacc      =     0
165   atfp      =     0.1
166   rdt       =  7200.
167   rdtmin    =  7200.
168   rdtmax    =  7200.
169   rdth      =   800.
170   rdtbt     =   120.
171   nclosea   =     0
172/
173
174!!======================================================================
175!!            ***  Surface Boundary Condition namelists  ***
176!!======================================================================
177!!   namsbc        surface boundary condition
178!!   namsbc_ana    analytical         formulation
179!!   namsbc_flx    flux               formulation
180!!   namsbc_clio   CLIO bulk formulea formulation
181!!   namsbc_core   CORE bulk formulea formulation
182!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled")
183!!   namqsr        penetrative solar radiation
184!!   namsbc_rnf    river runoffs
185!!   namsbc_ssr    sea surface restoring term (for T and/or S)
186!!   namalb        albedo parameters
187!!======================================================================
188
189!-----------------------------------------------------------------------
190!       namsbc       Surface Boundary Condition (surface module)
191!-----------------------------------------------------------------------
192!  nn_fsbc       frequency of surface boundary condition computation
193!                (= the frequency of sea-ice model call)
194!  ln_ana        analytical formulation (fill namsbc_ana )
195!  ln_flx        flux formulation       (fill namsbc_flx )
196!  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core)
197!  ln_blk_core   CORE bulk formulation  (fill namsbc_clio)
198!  ln_cpl        Coupled formulation    (fill namsbc_cpl )
199!  nn_ice        =0 no ice boundary condition   ,
200!                =1 observed ice-cover          ,
201!                =2 ice-model used                                      ("key_lim3" or "key_lim2)
202!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr)
203!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr)
204!  nn_fwb       FreshWater Budget: =0 no control                              ,
205!                                  =1 annual global mean of e-p-r set to zero   ,
206!                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step
207&namsbc
208   nn_fsbc     = 1
209   ln_ana      = .true.
210   ln_flx      = .false.
211   ln_blk_clio = .false.
212   ln_blk_core = .false.
213   ln_cpl      = .false.
214   nn_ice      = 0
215   ln_dm2dc    = .false.
216   ln_rnf      = .false.
217   ln_ssr      = .false.
218   nn_fwb      = 0
219/
220!-----------------------------------------------------------------------
221!       namsbc_ana   analytical surface boundary condition
222!-----------------------------------------------------------------------
223!  nn_tau000   gently increase the stress over the first ntau_rst time-steps
224!  rn_utau0    uniform value for the i-stress
225!  rn_vtau0    uniform value for the j-stress
226!  rn_q0       uniform value for the total heat flux
227!  rn_qsr0     uniform value for the solar radiation
228!  rn_emp0     uniform value for the freswater budget (E-P)
229&namsbc_ana
230   nn_tau000  =   100
231   rn_utau0   =   0.1e0
232   rn_vtau0   =   0.e0
233   rn_q0      =   0.e0
234   rn_qsr0    =   0.e0
235   rn_emp0    =   0.e0
236/
237!-----------------------------------------------------------------------
238!       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux")
239!-----------------------------------------------------------------------
240!  cn_dir   directory in which the model is executed
241!
242!  THE ORDER OF THE FILES MATTER:
243!               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting !
244!               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  !
245&namsbc_flx
246   cn_dir       = './'
247   sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0
248   sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0
249   sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0
250   sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0
251   sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0
252/     
253!-----------------------------------------------------------------------
254!       namsbc_clio  CLIO bulk formulea
255!-----------------------------------------------------------------------
256!  cn_dir   directory in which the model is executed
257!
258!  THE ORDER OF THE FILES MATTER:
259!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting !
260!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  !
261&namsbc_clio
262   cn_dir       = './'
263   sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0   
264   sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0   
265   sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0     
266   sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0     
267   sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0   
268   sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0     
269   sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0       
270/
271!-----------------------------------------------------------------------
272!       namsbc_core  CORE bulk formulea
273!-----------------------------------------------------------------------
274!  cn_dir         directory in which the model is executed
275!  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m
276!  alpha_precip   multiplication factor for precipitation (total & snow)
277!
278!  THE ORDER OF THE FILES MATTER:
279!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting !
280!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  !
281&namsbc_core
282   cn_dir       = './'
283   ln_2m        = .false.
284   alpha_precip = 1.
285   sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0   
286   sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0   
287   sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0     
288   sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0     
289   sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0     
290   sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0   
291   sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0     
292   sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0       
293/
294!-----------------------------------------------------------------------
295!       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled")
296!-----------------------------------------------------------------------
297&namsbc_cpl
298/
299!-----------------------------------------------------------------------
300!       namqsr   penetrative solar radiation
301!-----------------------------------------------------------------------
302!  ln_traqsr : penetrative solar radiation (T) or not (F)     (Default=T)
303!  rabs       fraction of qsr associated with xsi1
304!  xsi1       first depth of extinction
305!  xsi2       second depth of extinction
306&namqsr
307   ln_traqsr = .true.
308   rabs     =   0.58
309   xsi1     =   0.35
310   xsi2     =   23.0
311/
312!-----------------------------------------------------------------------
313!       namsbc_rnf   runoffs namelist surface boundary condition
314!-----------------------------------------------------------------------
315!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv
316!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used
317!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s]
318!cn_dir      Root directory for location of ssr files
319!sn_rnf      informations about the runoff file to be read
320!sn_cnf      informations about the runoff mouth file to be read
321!  THE ORDER OF THE FILES MATTER:
322!               !   file name        ! frequency (hours) ! variable   ! time interpolation !
323!               !                    !   (if <0  months) !   name     !     (logical)      !
324&namsbc_rnf
325   nn_runoff  = 0
326   rn_hrnf    = 0.e0
327   rn_avt_rnf = 1.e-3
328   cn_dir     = './'
329   sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0     
330   sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0             
331/
332!-----------------------------------------------------------------------
333!       namsbc_ssr   surface boundary condition : sea surface restoring
334!-----------------------------------------------------------------------
335!  nn_sstr   =0/1   add a retroaction term in the surface heat flux
336!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux
337!             (=1, salt flux, concentration/dillution emps only)
338!             (=2, volume flux, both emp and emps are updated)
339!  dqdt       magnitude of the retroaction on temperature [W/m2/K]
340!  deds       magnitude of the damping on salinity [mm/day]
341!  THE ORDER OF THE FILES MATTER:
342!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting !
343!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  !
344&namsbc_ssr
345   nn_sstr      = 0
346   nn_sssr      = 0
347   dqdt         = -40.
348   deds         =  27.7
349   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0
350   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0
351/     
352!-----------------------------------------------------------------------
353!       namalb   albedo parameters
354!-----------------------------------------------------------------------
355!  cgren    correction of the snow or ice albedo to take into account
356!  albice   albedo of melting ice in the arctic and antarctic
357!  alphd    coefficients for linear interpolation used to compute albedo
358!           between two extremes values (Pyane, 1972)
359!  alphc     "                                         "
360!  alphdi    "                                         "
361&namalb
362   cgren    =      0.06
363   albice   =      0.5
364   alphd    =      0.80
365   alphc    =      0.65
366   alphdi   =      0.72
367/
368
369!!======================================================================
370!!               ***  Lateral boundary condition  ***
371!!======================================================================
372!!   namlbc        lateral momentum boundary condition
373!!   namcla        cross land advection
374!!   namobc        open boundaries parameters                           ("key_obc")
375!!======================================================================
376
377!-----------------------------------------------------------------------
378!       namlbc   lateral momentum boundary condition
379!-----------------------------------------------------------------------
380!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip
381!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip
382&namlbc
383   shlat  =      0.
384/
385!-----------------------------------------------------------------------
386!       namcla   cross land advection
387!-----------------------------------------------------------------------
388!  n_cla   advection between 2 ocean pts separates by land
389&namcla
390   n_cla   = 0
391/
392!-----------------------------------------------------------------------
393!       namobc    open boundaries parameters (#ifdef key_obc)
394!-----------------------------------------------------------------------
395!  nobc_dta   = 0 the obc data are equal to the initial state
396!             = 1 the obc data are read in 'obc   .dta' files
397!  rdpeob  time relaxation (days) for the east open boundary
398!  rdpwob  time relaxation (days) for the west open boundary
399!  rdpnob  time relaxation (days) for the north open boundary
400!  rdpsob  time relaxation (days) for the south open boundary
401!  zbsic1  barotropic stream function on isolated coastline 1
402!  zbsic2  barotropic stream function on isolated coastline 2
403!  zbsic3  barotropic stream function on isolated coastline 3
404!  ln_obc_clim  climatological obc data files (default T)
405!  ln_vol_cst   total volume conserved
406&namobc
407    nobc_dta =    0
408    rdpein   =    1.
409    rdpwin   =    1.
410    rdpnin   =   30.
411    rdpsin   =    1.
412    rdpeob   = 1500.
413    rdpwob   =   15.
414    rdpnob   =  150.
415    rdpsob   =   15.
416    zbsic1   =  140.e+6
417    zbsic2   =    1.e+6
418    zbsic3   =    0.
419    ln_obc_clim = .true.
420    ln_vol_cst  = .false.
421/
422
423!!======================================================================
424!!                 ***  Bottom boundary condition  ***
425!!======================================================================
426!!   nambfr        bottom friction
427!!   nambbc        bottom temperature boundary condition                ("key_trabbc")
428!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv")
429!!======================================================================
430!-----------------------------------------------------------------------
431!       nambfr   bottom friction
432!-----------------------------------------------------------------------
433!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction
434!                                     = 3 : free slip  ;  = 2 : non linear friction
435!  bfri1    bottom drag coefficient (linear case)
436!  bfri2    bottom drag coefficient (non linear case)
437!  bfeb2    bottom turbulent kinetic energy  (m2/s2)
438&nambfr
439   nbotfr =       2
440   bfri1  =   4.e-4
441   bfri2  =   1.e-3
442   bfeb2  =  2.5e-3
443/
444!-----------------------------------------------------------------------
445!       nambbc   bottom temperature boundary condition
446!-----------------------------------------------------------------------
447!  ngeo_flux  = 0 no geothermal heat flux
448!             = 1 constant geothermal heat flux
449!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2)
450!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)
451&nambbc
452   ngeo_flux =  0
453   ngeo_flux_const = 86.4e-3
454/
455!-----------------------------------------------------------------------
456!       nambbl   bottom boundary layer scheme
457!-----------------------------------------------------------------------
458!                                !  diffusive bbl                       ("key_trabbl")
459!                                !  advective bbl                       ("key_trabbl_adv")
460!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s)
461&nambbl
462   atrbbl = 10000.
463/
464
465!!======================================================================
466!!                        Tracer (T & S ) namelists
467!!======================================================================
468!!   nameos        equation of state
469!!   nam_traadv    advection scheme
470!!   nam_traldf    lateral diffusion scheme
471!!   namtdp        tracer newtonian damping                             ("key_tradmp")
472!!======================================================================
473
474!-----------------------------------------------------------------------
475!       nameos   ocean physical parameters
476!-----------------------------------------------------------------------
477!  neos    type of equation of state and Brunt-Vaisala frequency
478!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) )
479!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T )
480!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T )
481!  ralpha  thermal expension coefficient (neos= 1 or 2)
482!  rbeta   saline  expension coefficient (neos= 2)
483&nameos
484   neos   =      2
485   ralpha =  2.e-4
486   rbeta  =  7.7e-4
487/
488!-----------------------------------------------------------------------
489!       nam_traadv   advection scheme for tracer
490!-----------------------------------------------------------------------
491!  ln_traadv_cen2     2nd order centered scheme          (default T)
492!  ln_traadv_tvd      TVD scheme                         (default F)
493!  ln_traadv_muscl    MUSCL scheme                       (default F)
494!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F)
495!  ln_traadv_ubs      UBS scheme                         (default F)
496!  ln_traadv_qck      QUICKEST scheme                    (default F)
497&nam_traadv
498   ln_traadv_cen2   =  .false.
499   ln_traadv_tvd    =  .true.
500   ln_traadv_muscl  =  .false.
501   ln_traadv_muscl2 =  .false.
502   ln_traadv_ubs    =  .false.
503   ln_traadv_qck    =  .false.
504/
505!-----------------------------------------------------------------------
506!       nam_traldf   lateral diffusion scheme for tracer
507!-----------------------------------------------------------------------
508!  Type of the operator :
509!     ln_traldf_lap    laplacian operator          (default T)
510!     ln_traldf_bilap  bilaplacian operator        (default F)
511!  Direction of action  :
512!     ln_traldf_level  iso-level                   (default F)
513!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T)
514!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
515!  Coefficient
516!     aht0    horizontal eddy diffusivity for tracers (m2/s)
517!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s)
518!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv")
519&nam_traldf
520   ln_traldf_lap    =  .true.
521   ln_traldf_bilap  =  .false.
522   ln_traldf_level  =  .false.
523   ln_traldf_hor    =  .false.
524   ln_traldf_iso    =  .true.
525   aht0    =  1000.
526   ahtb0   =     0.
527   aeiv0   =     0.
528/
529!-----------------------------------------------------------------------
530!       namtdp   tracer newtonian damping ('key_tradmp')
531!-----------------------------------------------------------------------
532!  ndmp    type of damping in temperature and salinity
533!          !   ='latitude', damping poleward of 'ndmp' degrees and function
534!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1
535!          !   =-1          damping only in Med and Red Seas
536!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array)
537!  nmldmp  type of damping in the mixed layer
538!          !   =0 damping throughout the water column
539!          !   =1 no damping in the mixed layer defined by avt >5cm2/s
540!          !   =2 no damping in the mixed layer defined rho<rho(surf)+.01
541!  sdmp    surface time scale for internal damping (days)
542!  bdmp    bottom time scale for internal damping (days)
543!  hdmp    depth of transition between sdmp and bdmp (meters)
544&namtdp
545   ndmp   =   -1
546   ndmpf  =    1
547   nmldmp =    1
548   sdmp   =  50.
549   bdmp   = 360.
550   hdmp   = 800.
551/
552
553!!======================================================================
554!!                      ***  Dynamics namelists  ***
555!!======================================================================
556!!   nam_dynadv    formulation of the momentum advection
557!!   nam_dynvor    advection scheme
558!!   nam_dynhpg    hydrostatic pressure gradient
559!!   namflg        hydrostatic pressure gradient time stepping
560!!   nam_dynspg    surface pressure gradient                            (CPP key only)
561!!   nam_dynldf    lateral diffusion scheme
562!!======================================================================
563
564!-----------------------------------------------------------------------
565!       nam_dynadv   formulation of the momentum advection
566!-----------------------------------------------------------------------
567!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T)
568!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T)
569!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F)
570&nam_dynadv   
571   ln_dynadv_vec   = .true.
572   ln_dynadv_cen2  = .false.
573   ln_dynadv_ubs   = .false.
574/ 
575!-----------------------------------------------------------------------
576!       nam_dynvor   option of physics/algorithm (not control by CPP keys)
577!-----------------------------------------------------------------------
578!  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T)
579!  ln_dynvor_ene      "         "  : energy conserving scheme    (default F)
580!  ln_dynvor_mix      "         "  : mixed scheme                (default F)
581!  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F)
582&nam_dynvor
583   ln_dynvor_ene = .true.
584   ln_dynvor_ens = .false.
585   ln_dynvor_mix = .false.
586   ln_dynvor_een = .false.
587/
588!-----------------------------------------------------------------------
589!       nam_dynhpg   Hydrostatic pressure gradient option
590!-----------------------------------------------------------------------
591!  type of pressure gradient scheme (choose one only!)
592!     ln_hpg_zco    z-coordinate - full steps                   (default T)
593!     ln_hpg_zps    z-coordinate - partial steps (interpolation)
594!     ln_hpg_sco    s-coordinate (standard jacobian formulation)
595!     ln_hpg_hel    s-coordinate (helsinki modification)
596!     ln_hpg_wdj    s-coordinate (weighted density jacobian)
597!     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial)
598!     ln_hpg_rot    s-coordinate (ROTated axes scheme)
599!  parameters
600!     gamm          weighting coefficient (wdj scheme)
601&nam_dynhpg
602   ln_hpg_zco = .true.
603   ln_hpg_zps = .false.
604   ln_hpg_sco = .false.
605   ln_hpg_hel = .false.
606   ln_hpg_wdj = .false.
607   ln_hpg_djc = .false.
608   ln_hpg_rot = .false.
609   gamm       = 0.e0
610/
611!-----------------------------------------------------------------------
612!       namflg   algorithm flags (algorithm not control by CPP keys)
613!-----------------------------------------------------------------------
614!  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T)
615!                                                  centered     time scheme  (F)
616!   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0)
617&namflg
618   ln_dynhpg_imp   =  .false.
619   nn_dynhpg_rst   =  0
620/
621!-----------------------------------------------------------------------
622!       nam_dynspg   surface pressure gradient   (CPP key only)
623!-----------------------------------------------------------------------
624!                                !  explicit free surface               ("key_dynspg_exp")
625!                                !  filtered free surface               ("key_dynspg_flt")
626!                                !  split-explicit free surface         ("key_dynspg_ts")
627!                                !  rigid-lid                           ("key_dynspg_rl")
628
629!-----------------------------------------------------------------------
630!       nam_dynldf   lateral diffusion on momentum
631!-----------------------------------------------------------------------
632!  Type of the operator :
633!     ln_dynldf_lap    laplacian operator          (default T)
634!     ln_dynldf_bilap  bilaplacian operator        (default F)
635!  Direction of action  :
636!     ln_dynldf_level  iso-level                   (default F)
637!     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T)
638!     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
639!  Coefficient
640!     ahm0    horizontal eddy viscosity for the dynamics (m2/s)
641!     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s)
642&nam_dynldf
643   ln_dynldf_lap    =  .true.
644   ln_dynldf_bilap  =  .false.
645   ln_dynldf_level  =  .false.
646   ln_dynldf_hor    =  .true.
647   ln_dynldf_iso    =  .false.
648   ahm0    = 100000.
649   ahmb0   =      0.
650/
651
652!!======================================================================
653!!             Tracers & Dynamics vertical physics namelists
654!!======================================================================
655!!       namzdf        vertical physics
656!!       namnpc        non penetrative convection                       
657!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      )
658!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      )
659!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      )
660!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      )
661!!======================================================================
662
663!-----------------------------------------------------------------------
664!       namzdf   vertical physics
665!-----------------------------------------------------------------------
666!  ln_zdfevd  enhanced vertical diffusion         (default T)
667!  ln_zdfnpc  Non-Penetrative Convection          (default T)
668!  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst")
669!  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst")
670!  avevd      vertical coefficient for enhanced diffusion scheme (m2/s)
671!  n_evdm     = 0  apply enhanced mixing on tracer only
672!             = 1  apply enhanced mixing on both tracer and momentum
673!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F)
674!                               (=F)  euler backward (F)
675!  n_zdfexp   number of sub-timestep for time splitting scheme
676&namzdf
677   ln_zdfevd = .true.
678   ln_zdfnpc = .false.
679   avm0      = 1.2e-4
680   avt0      = 1.2e-5
681   avevd     =   100.
682   n_evdm    =     1
683   ln_zdfexp =  .false.
684   n_zdfexp  =     3
685/
686!-----------------------------------------------------------------------
687!       namnpc   non penetrative convection
688!-----------------------------------------------------------------------
689!  nnpc1   non penetrative convective scheme frequency   
690!  nnpc2   non penetrative convective scheme print frequency
691&namnpc
692   nnpc1  =      1
693   nnpc2  =    365
694/
695!-----------------------------------------------------------------------
696!       namric   richardson number dependent vertical diffusion         ("key_zdfric" )
697!-----------------------------------------------------------------------
698!  avmri   maximum value of the vertical viscosity
699!  alp     coefficient of the parameterization
700!  nric    coefficient of the parameterization
701&namric
702   avmri = 100.e-4
703   alp   =      5.
704   nric  =       2
705/
706!-----------------------------------------------------------------------
707!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke")
708!-----------------------------------------------------------------------
709!  ln_rstke flag to restart with tke from a run without tke (default F)
710!  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) )
711!  ediss    coef. of the Kolmogoroff dissipation 
712!  ebb      coef. of the surface input of tke
713!  efave    coef. to applied to the tke diffusion ( avtke=efave*avm )
714!  emin     minimum value of tke (m^2/s^2)
715!  emin0    surface minimum value of tke (m^2/s^2)
716!  nitke    number of restart iterative loops
717!  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded)
718!  nmxl     length used   = 0 bounded by the distance to surface and bottom
719!           !             = 1 bounded by the local vertical scale factor
720!           !             = 2 first vertical derivative of mixing length bounded by 1
721!  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm)
722!           !                = 1 prandtl number function of richarson number (avt=pdl*avm)
723!           !                = 2 same as = 1 but a shapiro filter is applied on pdl
724!  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1)
725!  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb
726&namtke
727   ln_rstke = .false.
728   ediff =       0.1
729   ediss =       0.7
730   ebb   =      3.75
731   efave =        1.
732   emin  =     1.e-5
733   emin0 =     1.e-4
734   nitke =        50
735   nmxl  =         2
736   npdl  =         1
737   navb  =         0
738/
739!-----------------------------------------------------------------------
740!       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" )
741!-----------------------------------------------------------------------
742!  ln_kpprimix   shear instability mixing  (default T)
743!  difmiw        constant internal wave viscosity (m2/s)
744!  difsiw        constant internal wave diffusivity (m2/s)
745!  Riinfty       local Richardson Number limit for shear instability
746!  difri         maximum shear mixing at Rig = 0    (m2/s)
747!  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection
748!  difcon        maximum mixing in interior convection (m2/s)
749!  nave          = 0/1 flag for horizontal average on avt, avmu, avmv
750!  navb          = 0/1 flag for constant or profile background avt
751&namkpp
752   ln_kpprimix  = .true.
753   difmiw       =  1.e-04
754   difsiw       =  0.1e-04
755   Riinfty      =  0.8
756   difri        =  0.0050
757   bvsqcon      = -0.01e-07
758   difcon       =  1.
759   navb         =  0
760   nave         =  1
761/
762!-----------------------------------------------------------------------
763!       namddm   double diffusive mixing parameterization               ("key_zdfddm")
764!-----------------------------------------------------------------------
765!   avts    maximum avs for dd mixing
766!   hsbfr   heat/salt buoyancy flux ratio
767&namddm
768      avts  = 1.e-4
769      hsbfr = 1.6
770/
771
772!!======================================================================
773!!                  ***  Miscelaneous namelists  ***
774!!======================================================================
775!!   namsol      elliptic solver / island / free surface
776!!======================================================================
777
778!-----------------------------------------------------------------------
779!       namsol   elliptic solver / island / free surface
780!-----------------------------------------------------------------------
781!  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg)
782!                            (=2 successive-over-relaxation: sor)
783!                            (=3 FETI currently it does not work!       ("key_feti")
784!  nsol_arp  absolute/relative (0/1) precision convergence test
785!  nmin      minimum of iterations for the SOR solver
786!  nmax      maximum of iterations for the SOR solver
787!  nmod      frequency of test for the SOR solver
788!  eps       absolute precision of the solver
789!  resmax    absolute precision for the SOR solver
790!  sor       optimal coefficient for SOR solver
791!  epsisl    absolute precision on stream function solver
792!  nmisl     maximum pcg iterations for island
793!  rnu       strength of the additional force used in free surface b.c.
794&namsol
795   nsolv     =      2
796   nsol_arp  =      0
797   nmin      =    100
798   nmax      =    800
799   nmod      =     10
800   eps       =  1.E-6
801   resmax    = 1.E-10
802   sor       =   1.96
803   epsisl    = 1.e-10
804   nmisl     =   4000
805   rnu       =     1.
806/
807!!======================================================================
808!!                       ***  Diagnostics namelists  ***
809!!======================================================================
810!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld")
811!!   namgap       level mean model-data gap                             ("key_diagap")
812!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr")
813!!   namflo       float parameters                                      ("key_float")
814!!   namptr       Poleward Transport Diagnostics
815!!======================================================================
816
817!-----------------------------------------------------------------------
818!       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra")
819!                          or mixed-layer trends                        ('key_trdmld')
820!                          or barotropic vorticity                      ("key_trdvor")
821!-----------------------------------------------------------------------
822!  ntrd                time step frequency dynamics and tracers trends
823!  nctls               control surface type in mixed-layer trends (0,1 or n<jpk)
824!  ln_trdmld_restart   restart for ML diagnostics
825!  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day)
826!  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S
827&namtrd
828   ntrd              = 365
829   nctls             =   0
830   ln_trdmld_restart = .false.
831   ucf               =  1.
832   ln_trdmld_instant = .false.
833/
834!-----------------------------------------------------------------------
835!       namgap    level mean model-data gap ('key_diagap')
836!-----------------------------------------------------------------------
837!  ngap    time-step frequency of model-data gap computation
838!  nprg    time-step frequency of gap print in model output
839&namgap
840   ngap =  15
841   nprg =  10
842/
843!-----------------------------------------------------------------------
844!       namspr  surface pressure diagnostic
845!-----------------------------------------------------------------------
846!  nmaxp   maximum of iterations for the solver
847!  epsp    absolute precision of the solver
848!  niterp  number of iteration done by the solver
849&namspr
850   nmaxp   =   1000
851   epsp    =  1.e-3
852   niterp  =    400
853/
854!-----------------------------------------------------------------------
855!       namflo    float parameters                                      ("key_float")
856!-----------------------------------------------------------------------
857!  ln_rstflo   boolean term for float restart (true or false)
858!  nwritefl   frequency of float output file
859!  nstockfl   frequency of float restart file
860!  ln_argo    Argo type floats (stay at the surface each 10 days)
861!  ln_flork4  = T trajectories computed with a 4th order Runge-Kutta
862!             = F  (default)   computed with Blanke' scheme
863&namflo
864    ln_rstflo = .false.
865    nwritefl  =      75
866    nstockfl  =    5475
867    ln_argo   = .false.
868    ln_flork4 = .false.
869/
870!-----------------------------------------------------------------------
871!       namptr   Poleward Transport Diagnostic
872!-----------------------------------------------------------------------
873!  ln_diaptr  logical flag for Poleward transport computation
874!  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation
875!             need input basins mask file named "subbasins.nc"
876!  nf_ptr     Frequency of computation
877&namptr
878   ln_diaptr = .false.
879   ln_subbas = .false.
880   nf_ptr    =  15
881/
Note: See TracBrowser for help on using the repository browser.