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.
1_namelist in trunk/CONFIG/ORCA2_LIM/EXP00 – NEMO

source: trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist @ 908

Last change on this file since 908 was 908, checked in by ctlod, 16 years ago

add cn_dir parameter in namsbc_rnf & namsbc_ssr namelists, see ticket: #125

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 42.1 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     =  "Agulhas"
48   ln_rstart  = .false.
49   nrstdt     =       0
50   nit000     =       1
51   nitend     =   10950
52   ndate0     =  010101
53   nleapy     =       0
54   ninist     =       0
55   nstock     =   10950
56   nwrite     =   10950
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   =  .false.
124   ln_zps   =  .true.
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     =     1
161   e3zps_min =     5.
162   e3zps_rat =     0.1
163   nmsh      =     0
164   nacc      =     0
165   atfp      =     0.1
166   rdt       =  2880.
167   rdtmin    =  2880.
168   rdtmax    =  2880.
169   rdth      =   800.
170   rdtbt     =    90.
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     = 5
209   ln_ana      = .false.
210   ln_flx      = .false.
211   ln_blk_clio = .true.
212   ln_blk_core = .false.
213   ln_cpl      = .false.
214   nn_ice      = 2
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  =   0
231   rn_utau0   =   0.5
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!cn_dir      Root directory for location of ssr files
316!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv
317!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used
318!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s]
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   cn_dir     = './'
326   nn_runoff  = 2
327   rn_hrnf    = 0.e0
328   rn_avt_rnf = 1.e-3
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!  cn_dir    directory in which the model is executed
336!  nn_sstr   =0/1   add a retroaction term in the surface heat flux
337!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux
338!             (=1, salt flux, concentration/dillution emps only)
339!             (=2, volume flux, both emp and emps are updated)
340!  dqdt       magnitude of the retroaction on temperature [W/m2/K]
341!  deds       magnitude of the damping on salinity [mm/day]
342!  THE ORDER OF THE FILES MATTER:
343!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting !
344!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  !
345&namsbc_ssr
346   cn_dir       = './'
347   nn_sstr      = 0
348   nn_sssr      = 0
349   dqdt         = -40.
350   deds         =  27.7
351   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0
352   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0
353/     
354!-----------------------------------------------------------------------
355!       namalb   albedo parameters
356!-----------------------------------------------------------------------
357!  cgren    correction of the snow or ice albedo to take into account
358!  albice   albedo of melting ice in the arctic and antarctic
359!  alphd    coefficients for linear interpolation used to compute albedo
360!           between two extremes values (Pyane, 1972)
361!  alphc     "                                         "
362!  alphdi    "                                         "
363&namalb
364   cgren    =      0.06
365   albice   =      0.5
366   alphd    =      0.80
367   alphc    =      0.65
368   alphdi   =      0.72
369/
370
371!!======================================================================
372!!               ***  Lateral boundary condition  ***
373!!======================================================================
374!!   namlbc        lateral momentum boundary condition
375!!   namcla        cross land advection
376!!   namobc        open boundaries parameters                           ("key_obc")
377!!======================================================================
378
379!-----------------------------------------------------------------------
380!       namlbc   lateral momentum boundary condition
381!-----------------------------------------------------------------------
382!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip
383!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip
384&namlbc
385   shlat  =      2.
386/
387!-----------------------------------------------------------------------
388!       namcla   cross land advection
389!-----------------------------------------------------------------------
390!  n_cla   advection between 2 ocean pts separates by land
391&namcla
392   n_cla   = 0
393/
394!-----------------------------------------------------------------------
395!       namobc    open boundaries parameters (#ifdef key_obc)
396!-----------------------------------------------------------------------
397!  nobc_dta   = 0 the obc data are equal to the initial state
398!             = 1 the obc data are read in 'obc   .dta' files
399!  rdpeob  time relaxation (days) for the east open boundary
400!  rdpwob  time relaxation (days) for the west open boundary
401!  rdpnob  time relaxation (days) for the north open boundary
402!  rdpsob  time relaxation (days) for the south open boundary
403!  zbsic1  barotropic stream function on isolated coastline 1
404!  zbsic2  barotropic stream function on isolated coastline 2
405!  zbsic3  barotropic stream function on isolated coastline 3
406!  ln_obc_clim  climatological obc data files (default T)
407!  ln_vol_cst   total volume conserved
408&namobc
409    nobc_dta =    0
410    rdpein   =    1.
411    rdpwin   =    1.
412    rdpnin   =   30.
413    rdpsin   =    1.
414    rdpeob   = 1500.
415    rdpwob   =   15.
416    rdpnob   =  150.
417    rdpsob   =   15.
418    zbsic1   =  140.e+6
419    zbsic2   =    1.e+6
420    zbsic3   =    0.
421    ln_obc_clim = .true.
422    ln_vol_cst  = .false.
423/
424
425!!======================================================================
426!!                 ***  Bottom boundary condition  ***
427!!======================================================================
428!!   nambfr        bottom friction
429!!   nambbc        bottom temperature boundary condition                ("key_trabbc")
430!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv")
431!!======================================================================
432!-----------------------------------------------------------------------
433!       nambfr   bottom friction
434!-----------------------------------------------------------------------
435!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction
436!                                     = 3 : free slip  ;  = 2 : non linear friction
437!  bfri1    bottom drag coefficient (linear case)
438!  bfri2    bottom drag coefficient (non linear case)
439!  bfeb2    bottom turbulent kinetic energy  (m2/s2)
440&nambfr
441   nbotfr =       1
442   bfri1  =   4.e-4
443   bfri2  =   1.e-3
444   bfeb2  =  2.5e-3
445/
446!-----------------------------------------------------------------------
447!       nambbc   bottom temperature boundary condition
448!-----------------------------------------------------------------------
449!  ngeo_flux  = 0 no geothermal heat flux
450!             = 1 constant geothermal heat flux
451!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2)
452!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)
453&nambbc
454   ngeo_flux =  2
455   ngeo_flux_const = 86.4e-3
456/
457!-----------------------------------------------------------------------
458!       nambbl   bottom boundary layer scheme
459!-----------------------------------------------------------------------
460!                                !  diffusive bbl                       ("key_trabbl")
461!                                !  advective bbl                       ("key_trabbl_adv")
462!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s)
463&nambbl
464   atrbbl = 10000.
465/
466
467!!======================================================================
468!!                        Tracer (T & S ) namelists
469!!======================================================================
470!!   nameos        equation of state
471!!   nam_traadv    advection scheme
472!!   nam_traldf    lateral diffusion scheme
473!!   namtdp        tracer newtonian damping                             ("key_tradmp")
474!!======================================================================
475
476!-----------------------------------------------------------------------
477!       nameos   ocean physical parameters
478!-----------------------------------------------------------------------
479!  neos    type of equation of state and Brunt-Vaisala frequency
480!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) )
481!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T )
482!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T )
483!  ralpha  thermal expension coefficient (neos= 1 or 2)
484!  rbeta   saline  expension coefficient (neos= 2)
485&nameos
486   neos   =      0
487   ralpha =  2.e-4
488   rbeta  =  0.001
489/
490!-----------------------------------------------------------------------
491!       nam_traadv   advection scheme for tracer
492!-----------------------------------------------------------------------
493!  ln_traadv_cen2     2nd order centered scheme          (default T)
494!  ln_traadv_tvd      TVD scheme                         (default F)
495!  ln_traadv_muscl    MUSCL scheme                       (default F)
496!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F)
497!  ln_traadv_ubs      UBS scheme                         (default F)
498!  ln_traadv_qck      QUICKEST scheme                    (default F)
499&nam_traadv
500   ln_traadv_cen2   =  .false.
501   ln_traadv_tvd    =  .true.
502   ln_traadv_muscl  =  .false.
503   ln_traadv_muscl2 =  .false.
504   ln_traadv_ubs    =  .false.
505   ln_traadv_qck    =  .false.
506/
507!-----------------------------------------------------------------------
508!       nam_traldf   lateral diffusion scheme for tracer
509!-----------------------------------------------------------------------
510!  Type of the operator :
511!     ln_traldf_lap    laplacian operator          (default T)
512!     ln_traldf_bilap  bilaplacian operator        (default F)
513!  Direction of action  :
514!     ln_traldf_level  iso-level                   (default F)
515!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T)
516!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
517!  Coefficient
518!     aht0    horizontal eddy diffusivity for tracers (m2/s)
519!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s)
520!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv")
521&nam_traldf
522   ln_traldf_lap    =  .true.
523   ln_traldf_bilap  =  .false.
524   ln_traldf_level  =  .false.
525   ln_traldf_hor    =  .false.
526   ln_traldf_iso    =  .true.
527   aht0    =  1000.
528   ahtb0   =     0.
529   aeiv0   =     0.
530/
531!-----------------------------------------------------------------------
532!       namtdp   tracer newtonian damping ('key_tradmp')
533!-----------------------------------------------------------------------
534!  ndmp    type of damping in temperature and salinity
535!          !   ='latitude', damping poleward of 'ndmp' degrees and function
536!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1
537!          !   =-1          damping only in Med and Red Seas
538!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array)
539!  nmldmp  type of damping in the mixed layer
540!          !   =0 damping throughout the water column
541!          !   =1 no damping in the mixed layer defined by avt >5cm2/s
542!          !   =2 no damping in the mixed layer defined rho<rho(surf)+.01
543!  sdmp    surface time scale for internal damping (days)
544!  bdmp    bottom time scale for internal damping (days)
545!  hdmp    depth of transition between sdmp and bdmp (meters)
546&namtdp
547   ndmp   =   90
548   ndmpf  =    1
549   nmldmp =    1
550   sdmp   =  50.
551   bdmp   = 360.
552   hdmp   = 800.
553/
554
555!!======================================================================
556!!                      ***  Dynamics namelists  ***
557!!======================================================================
558!!   nam_dynadv    formulation of the momentum advection
559!!   nam_dynvor    advection scheme
560!!   nam_dynhpg    hydrostatic pressure gradient
561!!   namflg        hydrostatic pressure gradient time stepping
562!!   nam_dynspg    surface pressure gradient                            (CPP key only)
563!!   nam_dynldf    lateral diffusion scheme
564!!======================================================================
565
566!-----------------------------------------------------------------------
567!       nam_dynadv   formulation of the momentum advection
568!-----------------------------------------------------------------------
569!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T)
570!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T)
571!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F)
572&nam_dynadv   
573   ln_dynadv_vec   = .true.
574   ln_dynadv_cen2  = .false.
575   ln_dynadv_ubs   = .false.
576
577!-----------------------------------------------------------------------
578!       nam_dynvor   option of physics/algorithm (not control by CPP keys)
579!-----------------------------------------------------------------------
580!  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T)
581!  ln_dynvor_ene      "         "  : energy conserving scheme    (default F)
582!  ln_dynvor_mix      "         "  : mixed scheme                (default F)
583!  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F)
584&nam_dynvor
585   ln_dynvor_ene = .false.
586   ln_dynvor_ens = .false.
587   ln_dynvor_mix = .false.
588   ln_dynvor_een = .true.
589/
590!-----------------------------------------------------------------------
591!       nam_dynhpg   Hydrostatic pressure gradient option
592!-----------------------------------------------------------------------
593!  type of pressure gradient scheme (choose one only!)
594!     ln_hpg_zco    z-coordinate - full steps                   (default T)
595!     ln_hpg_zps    z-coordinate - partial steps (interpolation)
596!     ln_hpg_sco    s-coordinate (standard jacobian formulation)
597!     ln_hpg_hel    s-coordinate (helsinki modification)
598!     ln_hpg_wdj    s-coordinate (weighted density jacobian)
599!     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial)
600!     ln_hpg_rot    s-coordinate (ROTated axes scheme)
601!  parameters
602!     gamm          weighting coefficient (wdj scheme)
603&nam_dynhpg
604   ln_hpg_zco = .false.
605   ln_hpg_zps = .true.
606   ln_hpg_sco = .false.
607   ln_hpg_hel = .false.
608   ln_hpg_wdj = .false.
609   ln_hpg_djc = .false.
610   ln_hpg_rot = .false.
611   gamm       = 0.e0
612/
613!-----------------------------------------------------------------------
614!       namflg   algorithm flags (algorithm not control by CPP keys)
615!-----------------------------------------------------------------------
616!  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T)
617!                                                  centered     time scheme  (F)
618!   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0)
619&namflg
620   ln_dynhpg_imp   =  .false.
621   nn_dynhpg_rst   =  0
622/
623!-----------------------------------------------------------------------
624!       nam_dynspg   surface pressure gradient   (CPP key only)
625!-----------------------------------------------------------------------
626!                                !  explicit free surface               ("key_dynspg_exp")
627!                                !  filtered free surface               ("key_dynspg_flt")
628!                                !  split-explicit free surface         ("key_dynspg_ts")
629!                                !  rigid-lid                           ("key_dynspg_rl")
630
631!-----------------------------------------------------------------------
632!       nam_dynldf   lateral diffusion on momentum
633!-----------------------------------------------------------------------
634!  Type of the operator :
635!     ln_dynldf_lap    laplacian operator          (default T)
636!     ln_dynldf_bilap  bilaplacian operator        (default F)
637!  Direction of action  :
638!     ln_dynldf_level  iso-level                   (default F)
639!     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T)
640!     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
641!  Coefficient
642!     ahm0    horizontal eddy viscosity for the dynamics (m2/s)
643!     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s)
644&nam_dynldf
645   ln_dynldf_lap    =  .false.
646   ln_dynldf_bilap  =  .true.
647   ln_dynldf_level  =  .false.
648   ln_dynldf_hor    =  .true.
649   ln_dynldf_iso    =  .false.
650   ahm0    = -8.5e+11
651   ahmb0   =     0.
652/
653
654!!======================================================================
655!!             Tracers & Dynamics vertical physics namelists
656!!======================================================================
657!!       namzdf        vertical physics
658!!       namnpc        non penetrative convection                       
659!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      )
660!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      )
661!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      )
662!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      )
663!!======================================================================
664
665!-----------------------------------------------------------------------
666!       namzdf   vertical physics
667!-----------------------------------------------------------------------
668!  ln_zdfevd  enhanced vertical diffusion         (default T)
669!  ln_zdfnpc  Non-Penetrative Convection          (default T)
670!  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst")
671!  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst")
672!  avevd      vertical coefficient for enhanced diffusion scheme (m2/s)
673!  n_evdm     = 0  apply enhanced mixing on tracer only
674!             = 1  apply enhanced mixing on both tracer and momentum
675!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F)
676!                               (=F)  euler backward (F)
677!  n_zdfexp   number of sub-timestep for time splitting scheme
678&namzdf
679   ln_zdfevd = .true.
680   ln_zdfnpc = .false.
681   avm0      = 1.2e-4
682   avt0      = 1.2e-5
683   avevd     =   100.
684   n_evdm    =     1
685   ln_zdfexp =  .false.
686   n_zdfexp  =     3
687/
688!-----------------------------------------------------------------------
689!       namnpc   non penetrative convection
690!-----------------------------------------------------------------------
691!  nnpc1   non penetrative convective scheme frequency   
692!  nnpc2   non penetrative convective scheme print frequency
693&namnpc
694   nnpc1  =      1
695   nnpc2  =    365
696/
697!-----------------------------------------------------------------------
698!       namric   richardson number dependent vertical diffusion         ("key_zdfric" )
699!-----------------------------------------------------------------------
700!  avmri   maximum value of the vertical viscosity
701!  alp     coefficient of the parameterization
702!  nric    coefficient of the parameterization
703&namric
704   avmri = 100.e-4
705   alp   =      5.
706   nric  =       2
707/
708!-----------------------------------------------------------------------
709!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke")
710!-----------------------------------------------------------------------
711!  ln_rstke flag to restart with tke from a run without tke (default F)
712!  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) )
713!  ediss    coef. of the Kolmogoroff dissipation 
714!  ebb      coef. of the surface input of tke
715!  efave    coef. to applied to the tke diffusion ( avtke=efave*avm )
716!  emin     minimum value of tke (m^2/s^2)
717!  emin0    surface minimum value of tke (m^2/s^2)
718!  nitke    number of restart iterative loops
719!  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded)
720!  nmxl     length used   = 0 bounded by the distance to surface and bottom
721!           !             = 1 bounded by the local vertical scale factor
722!           !             = 2 first vertical derivative of mixing length bounded by 1
723!  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm)
724!           !                = 1 prandtl number function of richarson number (avt=pdl*avm)
725!           !                = 2 same as = 1 but a shapiro filter is applied on pdl
726!  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1)
727!  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb
728&namtke
729   ln_rstke = .false.
730   ediff =       0.1
731   ediss =       0.7
732   ebb   =       60.
733   efave =        1.
734   emin  =     1.e-6
735   emin0 =     1.e-4
736   nitke =        50
737   nmxl  =         2
738   npdl  =         1
739   navb  =         0
740/
741!-----------------------------------------------------------------------
742!       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" )
743!-----------------------------------------------------------------------
744!  ln_kpprimix   shear instability mixing  (default T)
745!  difmiw        constant internal wave viscosity (m2/s)
746!  difsiw        constant internal wave diffusivity (m2/s)
747!  Riinfty       local Richardson Number limit for shear instability
748!  difri         maximum shear mixing at Rig = 0    (m2/s)
749!  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection
750!  difcon        maximum mixing in interior convection (m2/s)
751!  nave          = 0/1 flag for horizontal average on avt, avmu, avmv
752!  navb          = 0/1 flag for constant or profile background avt
753&namkpp
754   ln_kpprimix  = .true.
755   difmiw       =  1.e-04
756   difsiw       =  0.1e-04
757   Riinfty      =  0.8
758   difri        =  0.0050
759   bvsqcon      = -0.01e-07
760   difcon       =  1.
761   navb         =  0
762   nave         =  1
763/
764!-----------------------------------------------------------------------
765!       namddm   double diffusive mixing parameterization               ("key_zdfddm")
766!-----------------------------------------------------------------------
767!   avts    maximum avs for dd mixing
768!   hsbfr   heat/salt buoyancy flux ratio
769&namddm
770      avts  = 1.e-4
771      hsbfr = 1.6
772/
773
774!!======================================================================
775!!                  ***  Miscelaneous namelists  ***
776!!======================================================================
777!!   namsol      elliptic solver / island / free surface
778!!======================================================================
779
780!-----------------------------------------------------------------------
781!       namsol   elliptic solver / island / free surface
782!-----------------------------------------------------------------------
783!  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg)
784!                            (=2 successive-over-relaxation: sor)
785!                            (=3 FETI currently it does not work!       ("key_feti")
786!  nsol_arp  absolute/relative (0/1) precision convergence test
787!  nmin      minimum of iterations for the SOR solver
788!  nmax      maximum of iterations for the SOR solver
789!  nmod      frequency of test for the SOR solver
790!  eps       absolute precision of the solver
791!  resmax    absolute precision for the SOR solver
792!  sor       optimal coefficient for SOR solver
793!  epsisl    absolute precision on stream function solver
794!  nmisl     maximum pcg iterations for island
795!  rnu       strength of the additional force used in free surface b.c.
796&namsol
797   nsolv     =      1
798   nsol_arp  =      0
799   nmin      =    300
800   nmax      =    800
801   nmod      =     10
802   eps       =  1.E-6
803   resmax    = 1.E-10
804   sor       =   1.92
805   epsisl    = 1.e-10
806   nmisl     =   4000
807   rnu       =     1.
808/
809!!======================================================================
810!!                       ***  Diagnostics namelists  ***
811!!======================================================================
812!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld")
813!!   namgap       level mean model-data gap                             ("key_diagap")
814!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr")
815!!   namflo       float parameters                                      ("key_float")
816!!   namptr       Poleward Transport Diagnostics
817!!======================================================================
818
819!-----------------------------------------------------------------------
820!       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra")
821!                          or mixed-layer trends                        ('key_trdmld')
822!                          or barotropic vorticity                      ("key_trdvor")
823!-----------------------------------------------------------------------
824!  ntrd                time step frequency dynamics and tracers trends
825!  nctls               control surface type in mixed-layer trends (0,1 or n<jpk)
826!  ln_trdmld_restart   restart for ML diagnostics
827!  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day)
828!  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S
829&namtrd
830   ntrd              = 365
831   nctls             =   0
832   ln_trdmld_restart = .false.
833   ucf               =  1.
834   ln_trdmld_instant = .false.
835/
836!-----------------------------------------------------------------------
837!       namgap    level mean model-data gap ('key_diagap')
838!-----------------------------------------------------------------------
839!  ngap    time-step frequency of model-data gap computation
840!  nprg    time-step frequency of gap print in model output
841&namgap
842   ngap =  15
843   nprg =  10
844/
845!-----------------------------------------------------------------------
846!       namspr  surface pressure diagnostic
847!-----------------------------------------------------------------------
848!  nmaxp   maximum of iterations for the solver
849!  epsp    absolute precision of the solver
850!  niterp  number of iteration done by the solver
851&namspr
852   nmaxp   =   1000
853   epsp    =  1.e-3
854   niterp  =    400
855/
856!-----------------------------------------------------------------------
857!       namflo    float parameters                                      ("key_float")
858!-----------------------------------------------------------------------
859!  ln_rstflo   boolean term for float restart (true or false)
860!  nwritefl   frequency of float output file
861!  nstockfl   frequency of float restart file
862!  ln_argo    Argo type floats (stay at the surface each 10 days)
863!  ln_flork4  = T trajectories computed with a 4th order Runge-Kutta
864!             = F  (default)   computed with Blanke' scheme
865&namflo
866    ln_rstflo = .false.
867    nwritefl  =      75
868    nstockfl  =    5475
869    ln_argo   = .false.
870    ln_flork4 = .false.
871/
872!-----------------------------------------------------------------------
873!       namptr   Poleward Transport Diagnostic
874!-----------------------------------------------------------------------
875!  ln_diaptr  logical flag for Poleward transport computation
876!  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation
877!             need input basins mask file named "subbasins.nc"
878!  nf_ptr     Frequency of computation
879&namptr
880   ln_diaptr = .false.
881   ln_subbas = .false.
882   nf_ptr    =  15
883/
Note: See TracBrowser for help on using the repository browser.