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

source: trunk/CONFIG/ORCA2_LIM/EXP00/namelist @ 1055

Last change on this file since 1055 was 1055, checked in by rblod, 16 years ago

Desactived control print (ln_nicep) for LIM3, correct forcing field for tair, supress u_ice and v_ice as argument for buls clior

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