source: CONFIG/trunk/IPSLCM5_WORK/EXP00/PARAM/namelist_ORCA2 @ 396

Last change on this file since 396 was 396, checked in by aclsce, 14 years ago

Creation of IPSLCM5_WORK coupled configuration as working configuration.

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