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_cfg in branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/GYRE/EXP00 – NEMO

source: branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg @ 5760

Last change on this file since 5760 was 5760, checked in by gm, 9 years ago

#1593: LDF-ADV, step II.1: phasing the improvements/simplifications of diffusive trend: suppression of KPP

File size: 23.9 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  : GYRE_PISCES Configuration namelist to overwrite reference dynamical namelist
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!-----------------------------------------------------------------------
5&namrun        !   parameters of the run
6!-----------------------------------------------------------------------
7   cn_exp      =  "GYRE"   !  experience name
8   nn_it000    =       1   !  first time step
9   nn_itend    =    4320   !  last  time step
10   nn_leapy    =      30   !  Leap year calendar (1) or not (0)
11   nn_stock    =    4320   !  frequency of creation of a restart file (modulo referenced to 1)
12   nn_write    =      60   !  frequency of write in the output file   (modulo referenced to nn_it000)
13
14   ln_clobber  = .true.    !  clobber (overwrite) an existing file
15
16/
17!-----------------------------------------------------------------------
18&namcfg     !   parameters of the configuration   
19!-----------------------------------------------------------------------
20   cp_cfg      =  "gyre"                 !  name of the configuration
21   jp_cfg      =       1                 !  resolution of the configuration
22   jpidta      =      32                 !  1st lateral dimension ( >= jpi ) = 30*jp_cfg+2
23   jpjdta      =      22                 !  2nd    "         "    ( >= jpj ) = 20*jp_cfg+2
24   jpkdta      =      31                 !  number of levels      ( >= jpk )
25   jpiglo      =      32                 !  1st dimension of global domain --> i  = jpidta
26   jpjglo      =      22                 !  2nd    -                  -    --> j  = jpjdta
27   jpizoom     =       1                 !  left bottom (i,j) indices of the zoom
28   jpjzoom     =       1                 !  in data domain indices
29   jperio      =       0                 !  lateral cond. type (between 0 and 6)
30/
31!-----------------------------------------------------------------------
32&namzgr        !   vertical coordinate
33!-----------------------------------------------------------------------
34   ln_zco      = .true.    !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined)
35   ln_zps      = .false.   !  z-coordinate - partial steps   (T/F)
36/
37!-----------------------------------------------------------------------
38&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate
39!-----------------------------------------------------------------------
40/
41!-----------------------------------------------------------------------
42&namdom        !   space and time domain (bathymetry, mesh, timestep)
43!-----------------------------------------------------------------------
44   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file
45   rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0)
46   rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1)
47   rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1)
48   jphgr_msh   =       5                 !  type of horizontal mesh
49   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1)
50   ppgphi0     =      29.0               ! latitude  of first raw and column T-point (jphgr_msh = 1)
51   ppe1_deg    =  999999.0               !  zonal      grid-spacing (degrees)
52   ppe2_deg    =  999999.0               !  meridional grid-spacing (degrees)
53   ppe1_m      =  999999.0               !  zonal      grid-spacing (degrees)
54   ppe2_m      =  999999.0               !  meridional grid-spacing (degrees)
55   ppsur       =   -2033.194295283385    !  ORCA r4, r2 and r05 coefficients
56   ppa0        =     155.8325369664153   ! (default coefficients)
57   ppa1        =     146.3615918601890   !
58   ppkth       =      17.28520372419791  !
59   ppacr       =       5.0               !
60   ppdzmin     =  999999.0               !  Minimum vertical spacing
61   pphmax      =  999999.0               !  Maximum depth
62   ldbletanh   =  .FALSE.                !  Use/do not use double tanf function for vertical coordinates
63   ppa2        =  999999.0               !  Double tanh function parameters
64   ppkth2      =  999999.0               !
65   ppacr2      =  999999.0               !
66/
67!-----------------------------------------------------------------------
68&namsplit      !   time splitting parameters                            ("key_dynspg_ts")
69!-----------------------------------------------------------------------
70/
71!-----------------------------------------------------------------------
72&namcrs        !   Grid coarsening for dynamics output and/or
73               !   passive tracer coarsened online simulations
74!-----------------------------------------------------------------------
75/
76!-----------------------------------------------------------------------
77&namtsd    !   data : Temperature  & Salinity
78!-----------------------------------------------------------------------
79   cn_dir        = './'      !  root directory for the location of the runoff files
80   ln_tsd_init   = .false.   !  Initialisation of ocean T & S with T &S input data (T) or not (F)
81   ln_tsd_tradmp = .false.   !  damping of ocean T & S toward T &S input data (T) or not (F)
82/
83!-----------------------------------------------------------------------
84&namsbc        !   Surface Boundary Condition (surface module)
85!-----------------------------------------------------------------------
86   nn_fsbc     = 1         !  frequency of surface boundary condition computation
87   !                       !     (also = the frequency of sea-ice model call)
88   ln_ana      = .true.    !  analytical formulation                    (T => fill namsbc_ana )
89   ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core)
90   nn_ice      = 0         !  =0 no ice boundary condition   ,
91   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf)
92   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
93   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked
94/
95!-----------------------------------------------------------------------
96&namsbc_ana    !   analytical surface boundary condition
97!-----------------------------------------------------------------------
98   nn_tau000   =   100     !  gently increase the stress over the first ntau_rst time-steps
99   rn_utau0    =   0.1e0   !  uniform value for the i-stress
100/
101!-----------------------------------------------------------------------
102&namsbc_flx    !   surface boundary condition : flux formulation
103!-----------------------------------------------------------------------
104/
105!-----------------------------------------------------------------------
106&namsbc_clio   !   namsbc_clio  CLIO bulk formulae
107!-----------------------------------------------------------------------
108/
109!-----------------------------------------------------------------------
110&namsbc_core   !   namsbc_core  CORE bulk formulae
111!-----------------------------------------------------------------------
112/
113!-----------------------------------------------------------------------
114&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae
115!-----------------------------------------------------------------------
116/
117!-----------------------------------------------------------------------
118&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
119!-----------------------------------------------------------------------
120/
121!-----------------------------------------------------------------------
122&namtra_qsr    !   penetrative solar radiation
123!-----------------------------------------------------------------------
124   ln_qsr_rgb  = .false.   !  RGB (Red-Green-Blue) light penetration
125   ln_qsr_2bd  = .true.    !  2 bands              light penetration
126   nn_chldta   =      0    !  RGB : Chl data (=1) or cst value (=0)
127/
128!-----------------------------------------------------------------------
129&namsbc_rnf    !   runoffs namelist surface boundary condition
130!-----------------------------------------------------------------------
131   ln_rnf_mouth = .false.   !  specific treatment at rivers mouths
132/
133!-----------------------------------------------------------------------
134&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk
135!-----------------------------------------------------------------------
136/
137!-----------------------------------------------------------------------
138&namsbc_ssr    !   surface boundary condition : sea surface restoring
139!-----------------------------------------------------------------------
140   nn_sssr     =     0     !  add a damping     term in the surface freshwater flux (=2)
141   rn_deds     =   -27.7   !  magnitude of the damping on salinity   [mm/day]
142   ln_sssr_bnd =   .false. !  flag to bound erp term (associated with nn_sssr=2)
143/
144!-----------------------------------------------------------------------
145&namsbc_alb    !   albedo parameters
146!-----------------------------------------------------------------------
147/
148!-----------------------------------------------------------------------
149&namberg       !   iceberg parameters
150!-----------------------------------------------------------------------
151/
152!-----------------------------------------------------------------------
153&namlbc        !   lateral momentum boundary condition
154!-----------------------------------------------------------------------
155   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
156/
157!-----------------------------------------------------------------------
158&namcla        !   cross land advection
159!-----------------------------------------------------------------------
160/
161!-----------------------------------------------------------------------
162&namobc        !   open boundaries parameters                           ("key_obc")
163!-----------------------------------------------------------------------
164/
165!-----------------------------------------------------------------------
166&namagrif      !  AGRIF zoom                                            ("key_agrif")
167!-----------------------------------------------------------------------
168/
169!-----------------------------------------------------------------------
170&nam_tide      !    tide parameters (#ifdef key_tide)
171!-----------------------------------------------------------------------
172/
173!-----------------------------------------------------------------------
174&nambdy        !  unstructured open boundaries                          ("key_bdy")
175!-----------------------------------------------------------------------
176/
177!-----------------------------------------------------------------------
178&nambdy_dta      !  open boundaries - external data           ("key_bdy")
179!-----------------------------------------------------------------------
180/
181!-----------------------------------------------------------------------
182&nambdy_tide     ! tidal forcing at open boundaries
183!-----------------------------------------------------------------------
184/
185!-----------------------------------------------------------------------
186&nambfr        !   bottom friction
187!-----------------------------------------------------------------------
188   nn_bfr      =    2      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
189/
190!-----------------------------------------------------------------------
191&nambbc        !   bottom temperature boundary condition
192!-----------------------------------------------------------------------
193   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
194   nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux
195/
196!-----------------------------------------------------------------------
197&nambbl        !   bottom boundary layer scheme
198!-----------------------------------------------------------------------
199/
200!-----------------------------------------------------------------------
201&nameos        !   ocean physical parameters
202!-----------------------------------------------------------------------
203   nn_eos      =  0       !  type of equation of state and Brunt-Vaisala frequency
204                                 !  =-1, TEOS-10
205                                 !  = 0, EOS-80
206                                 !  = 1, S-EOS   (simplified eos)
207   ln_useCT    = .false.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm
208   !                             !
209   !                      ! S-EOS coefficients :
210   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
211   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1)
212   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1)
213   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
214   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
215   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
216   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
217   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
218!!org GYRE   rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2)
219!!org GYRE   rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2)
220!!org  caution  now a0 = alpha / rau0   with rau0 = 1026
221/
222!-----------------------------------------------------------------------
223&namtra_adv    !   advection scheme for tracer
224!-----------------------------------------------------------------------
225   ln_traadv_cen2   =  .false.   !  2nd order centered scheme
226   ln_traadv_tvd    =  .true.    !  TVD scheme
227   ln_traadv_muscl  =  .false.   !  MUSCL scheme
228   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries
229   ln_traadv_ubs    =  .false.   !  UBS scheme
230   ln_traadv_qck    =  .false.   !  QUICKEST scheme
231   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl
232/
233!-----------------------------------------------------------------------
234&namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param)
235!-----------------------------------------------------------------------
236/
237!----------------------------------------------------------------------------------
238&namtra_ldf    !   lateral diffusion scheme for tracers
239!----------------------------------------------------------------------------------
240   !                       !  Operator type:
241   ln_traldf_lap   =  .true.   !    laplacian operator
242   ln_traldf_blp   =  .false.  !  bilaplacian operator
243   !                       !  Direction of action:
244   ln_traldf_lev   =  .false.  !  iso-level
245   ln_traldf_hor   =  .false.  !  horizontal (geopotential)
246   ln_traldf_iso   =  .true.   !  iso-neutral
247   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads
248   !
249   !                       !  iso-neutral options:       
250   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators)
251   rn_slpmax       =   0.01    !  slope limit                      (both operators)
252   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only)
253   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only)
254   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only)
255   !
256   !                       !  Coefficients:
257   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef
258   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
259   !                                !   =  0           constant
260   !                                !   = 10 F(k)      =ldf_c1d
261   !                                !   = 20 F(i,j)    =ldf_c2d
262   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
263   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d
264   !                                !   = 31 F(i,j,k,t)=F(local velocity)
265   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s]
266   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s]
267/
268!----------------------------------------------------------------------------------
269&namtra_ldfeiv !   eddy induced velocity param.
270!----------------------------------------------------------------------------------
271   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization
272/
273!-----------------------------------------------------------------------
274&namtra_dmp    !   tracer: T & S newtonian damping
275!-----------------------------------------------------------------------
276   ln_tradmp   =  .false.  !  add a damping termn (T) or not (F)
277/
278!-----------------------------------------------------------------------
279&namdyn_adv    !   formulation of the momentum advection
280!-----------------------------------------------------------------------
281/
282!-----------------------------------------------------------------------
283&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
284!-----------------------------------------------------------------------
285   ln_dynvor_ene = .true.  !  energy    conserving scheme 
286   ln_dynvor_ens = .false. !  enstrophy conserving scheme   
287   ln_dynvor_een = .false. !  energy & enstrophy scheme
288/
289!-----------------------------------------------------------------------
290&namdyn_hpg    !   Hydrostatic pressure gradient option
291!-----------------------------------------------------------------------
292   ln_hpg_zco  = .true.    !  z-coordinate - full steps
293   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
294/
295!-----------------------------------------------------------------------
296!namdyn_spg    !   surface pressure gradient   (CPP key only)
297!-----------------------------------------------------------------------
298
299!-----------------------------------------------------------------------
300&namdyn_ldf    !   lateral diffusion on momentum
301!-----------------------------------------------------------------------
302   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s]
303/
304!-----------------------------------------------------------------------
305&namzdf        !   vertical physics
306!-----------------------------------------------------------------------
307   nn_evdm     =    1      !  evd apply on tracer (=0) or on tracer and momentum (=1)
308/
309!-----------------------------------------------------------------------
310&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
311!-----------------------------------------------------------------------
312/
313!-----------------------------------------------------------------------
314&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
315!-----------------------------------------------------------------------
316   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves
317/
318!-----------------------------------------------------------------------
319&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls")
320!-----------------------------------------------------------------------
321/
322!-----------------------------------------------------------------------
323&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
324!-----------------------------------------------------------------------
325/
326!-----------------------------------------------------------------------
327&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
328!-----------------------------------------------------------------------
329   ln_tmx_itf  = .false.   !  ITF specific parameterisation
330/
331!-----------------------------------------------------------------------
332&namsol        !   elliptic solver / island / free surface
333!-----------------------------------------------------------------------
334   nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
335   nn_nmin     =    210    !  minimum of iterations for the SOR solver
336   rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain)
337/
338!-----------------------------------------------------------------------
339&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
340!-----------------------------------------------------------------------
341/
342!-----------------------------------------------------------------------
343&namctl        !   Control prints & Benchmark
344!-----------------------------------------------------------------------
345/
346!-----------------------------------------------------------------------
347&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
348!-----------------------------------------------------------------------
349/
350!-----------------------------------------------------------------------
351&namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra")
352!              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor")
353!-----------------------------------------------------------------------
354   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
355   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
356   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
357   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet)
358   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
359   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
360   ln_tra_trd  = .false.    ! (T) 3D tracer trend output
361   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
362   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
363/
364!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
365!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
366!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
367!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
368!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
369!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
370!!gm
371!-----------------------------------------------------------------------
372&namflo       !   float parameters                                      ("key_float")
373!-----------------------------------------------------------------------
374/
375!-----------------------------------------------------------------------
376&namptr       !   Poleward Transport Diagnostic
377!-----------------------------------------------------------------------
378/
379!-----------------------------------------------------------------------
380&namhsb       !  Heat and salt budgets
381!-----------------------------------------------------------------------
382/
383!-----------------------------------------------------------------------
384&namdct        ! transports through sections
385!-----------------------------------------------------------------------
386    nn_dct      = 60       !  time step frequency for transports computing
387    nn_dctwri   = 60       !  time step frequency for transports writing
388    nn_secdebug = 0        !      0 : no section to debug
389/
390!-----------------------------------------------------------------------
391&namobs       !  observation usage switch                               ('key_diaobs')
392!-----------------------------------------------------------------------
393/
394!-----------------------------------------------------------------------
395&nam_asminc   !   assimilation increments                               ('key_asminc')
396!-----------------------------------------------------------------------
397/
398!-----------------------------------------------------------------------
399&namsbc_wave   ! External fields from wave model
400!-----------------------------------------------------------------------
401/
402!-----------------------------------------------------------------------
403&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
404!-----------------------------------------------------------------------
405   ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water
406/
Note: See TracBrowser for help on using the repository browser.