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_r5836_NOC3_vvl_by_default/NEMOGCM/CONFIG/GYRE/EXP00 – NEMO

source: branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg @ 5866

Last change on this file since 5866 was 5866, checked in by gm, 8 years ago

#1613: vvl by default: add ln_linssh and remove key_vvl

File size: 24.7 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
35   ln_linssh   = .true.    !  linear free surface
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&namagrif      !  AGRIF zoom                                            ("key_agrif")
159!-----------------------------------------------------------------------
160/
161!-----------------------------------------------------------------------
162&nam_tide      !    tide parameters (#ifdef key_tide)
163!-----------------------------------------------------------------------
164/
165!-----------------------------------------------------------------------
166&nambdy        !  unstructured open boundaries                          ("key_bdy")
167!-----------------------------------------------------------------------
168/
169!-----------------------------------------------------------------------
170&nambdy_dta      !  open boundaries - external data           ("key_bdy")
171!-----------------------------------------------------------------------
172/
173!-----------------------------------------------------------------------
174&nambdy_tide     ! tidal forcing at open boundaries
175!-----------------------------------------------------------------------
176/
177!-----------------------------------------------------------------------
178&nambfr        !   bottom friction
179!-----------------------------------------------------------------------
180   nn_bfr      =    2      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
181/
182!-----------------------------------------------------------------------
183&nambbc        !   bottom temperature boundary condition
184!-----------------------------------------------------------------------
185   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
186   nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux
187/
188!-----------------------------------------------------------------------
189&nambbl        !   bottom boundary layer scheme
190!-----------------------------------------------------------------------
191/
192!-----------------------------------------------------------------------
193&nameos        !   ocean physical parameters
194!-----------------------------------------------------------------------
195   nn_eos      =  0       !  type of equation of state and Brunt-Vaisala frequency
196                                 !  =-1, TEOS-10
197                                 !  = 0, EOS-80
198                                 !  = 1, S-EOS   (simplified eos)
199   ln_useCT    = .false.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm
200   !                             !
201   !                      ! S-EOS coefficients :
202   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
203   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1)
204   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1)
205   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
206   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
207   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
208   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
209   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
210!!org GYRE   rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2)
211!!org GYRE   rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2)
212!!org  caution  now a0 = alpha / rau0   with rau0 = 1026
213/
214!-----------------------------------------------------------------------
215&namtra_adv    !   advection scheme for tracer
216!-----------------------------------------------------------------------
217   ln_traadv_fct =  .true.   !  FCT scheme
218      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order
219      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order
220      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping
221      !                             !        (number of sub-timestep = nn_fct_zts)
222/
223!-----------------------------------------------------------------------
224&namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param)
225!-----------------------------------------------------------------------
226/
227!----------------------------------------------------------------------------------
228&namtra_ldf    !   lateral diffusion scheme for tracers
229!----------------------------------------------------------------------------------
230   !                       !  Operator type:
231   ln_traldf_lap   =  .true.   !    laplacian operator
232   ln_traldf_blp   =  .false.  !  bilaplacian operator
233   !                       !  Direction of action:
234   ln_traldf_lev   =  .false.  !  iso-level
235   ln_traldf_hor   =  .false.  !  horizontal (geopotential)
236   ln_traldf_iso   =  .true.   !  iso-neutral
237   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads
238   !
239   !                       !  iso-neutral options:       
240   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators)
241   rn_slpmax       =   0.01    !  slope limit                      (both operators)
242   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only)
243   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only)
244   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only)
245   !
246   !                       !  Coefficients:
247   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef
248   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
249   !                                !   =  0           constant
250   !                                !   = 10 F(k)      =ldf_c1d
251   !                                !   = 20 F(i,j)    =ldf_c2d
252   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
253   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d
254   !                                !   = 31 F(i,j,k,t)=F(local velocity)
255   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s]
256   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s]
257/
258!----------------------------------------------------------------------------------
259&namtra_ldfeiv !   eddy induced velocity param.
260!----------------------------------------------------------------------------------
261   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization
262/
263!-----------------------------------------------------------------------
264&namtra_dmp    !   tracer: T & S newtonian damping
265!-----------------------------------------------------------------------
266   ln_tradmp   =  .false.  !  add a damping termn (T) or not (F)
267/
268!-----------------------------------------------------------------------
269&namdyn_adv    !   formulation of the momentum advection
270!-----------------------------------------------------------------------
271/
272!-----------------------------------------------------------------------
273&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
274!-----------------------------------------------------------------------
275   ln_dynvor_ene = .true.  !  enstrophy conserving scheme
276   ln_dynvor_ens = .false. !  energy conserving scheme
277   ln_dynvor_mix = .false. !  mixed scheme
278   ln_dynvor_een = .false. !  energy & enstrophy scheme
279      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)
280/
281!-----------------------------------------------------------------------
282&namdyn_hpg    !   Hydrostatic pressure gradient option
283!-----------------------------------------------------------------------
284   ln_hpg_zco  = .true.    !  z-coordinate - full steps
285   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
286/
287!-----------------------------------------------------------------------
288!namdyn_spg    !   surface pressure gradient   (CPP key only)
289!-----------------------------------------------------------------------
290
291!-----------------------------------------------------------------------
292&namdyn_ldf    !   lateral diffusion on momentum
293!-----------------------------------------------------------------------
294   !                       !  Type of the operator :
295   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F
296   ln_dynldf_lap =  .true.     !    laplacian operator
297   ln_dynldf_blp =  .false.    !  bilaplacian operator
298   !                       !  Direction of action  :
299   ln_dynldf_lev =  .true.     !  iso-level
300   ln_dynldf_hor =  .false.    !  horizontal (geopotential)
301   ln_dynldf_iso =  .false.    !  iso-neutral
302   !                       !  Coefficient
303   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef
304   !                                !  =-30  read in eddy_viscosity_3D.nc file
305   !                                !  =-20  read in eddy_viscosity_2D.nc file
306   !                                !  =  0  constant
307   !                                !  = 10  F(k)=c1d
308   !                                !  = 20  F(i,j)=F(grid spacing)=c2d
309   !                                !  = 30  F(i,j,k)=c2d*c1d
310   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity)
311   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s]
312   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s]
313   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s]
314/
315!-----------------------------------------------------------------------
316&namzdf        !   vertical physics
317!-----------------------------------------------------------------------
318   nn_evdm     =    1      !  evd apply on tracer (=0) or on tracer and momentum (=1)
319/
320!-----------------------------------------------------------------------
321&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
322!-----------------------------------------------------------------------
323/
324!-----------------------------------------------------------------------
325&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
326!-----------------------------------------------------------------------
327   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves
328/
329!-----------------------------------------------------------------------
330&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls")
331!-----------------------------------------------------------------------
332/
333!-----------------------------------------------------------------------
334&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
335!-----------------------------------------------------------------------
336/
337!-----------------------------------------------------------------------
338&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
339!-----------------------------------------------------------------------
340   ln_tmx_itf  = .false.   !  ITF specific parameterisation
341/
342!-----------------------------------------------------------------------
343&namsol        !   elliptic solver / island / free surface
344!-----------------------------------------------------------------------
345   nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
346   nn_nmin     =    210    !  minimum of iterations for the SOR solver
347   rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain)
348/
349!-----------------------------------------------------------------------
350&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
351!-----------------------------------------------------------------------
352/
353!-----------------------------------------------------------------------
354&namctl        !   Control prints & Benchmark
355!-----------------------------------------------------------------------
356/
357!-----------------------------------------------------------------------
358&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
359!-----------------------------------------------------------------------
360/
361!-----------------------------------------------------------------------
362&namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra")
363!              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor")
364!-----------------------------------------------------------------------
365   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
366   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
367   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
368   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet)
369   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
370   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
371   ln_tra_trd  = .false.    ! (T) 3D tracer trend output
372   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
373   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
374/
375!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
376!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
377!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
378!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
379!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
380!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
381!!gm
382!-----------------------------------------------------------------------
383&namflo       !   float parameters                                      ("key_float")
384!-----------------------------------------------------------------------
385/
386!-----------------------------------------------------------------------
387&namptr       !   Poleward Transport Diagnostic
388!-----------------------------------------------------------------------
389/
390!-----------------------------------------------------------------------
391&namhsb       !  Heat and salt budgets
392!-----------------------------------------------------------------------
393/
394!-----------------------------------------------------------------------
395&namdct        ! transports through sections
396!-----------------------------------------------------------------------
397    nn_dct      = 60       !  time step frequency for transports computing
398    nn_dctwri   = 60       !  time step frequency for transports writing
399    nn_secdebug = 0        !      0 : no section to debug
400/
401!-----------------------------------------------------------------------
402&namobs       !  observation usage switch                               ('key_diaobs')
403!-----------------------------------------------------------------------
404/
405!-----------------------------------------------------------------------
406&nam_asminc   !   assimilation increments                               ('key_asminc')
407!-----------------------------------------------------------------------
408/
409!-----------------------------------------------------------------------
410&namsbc_wave   ! External fields from wave model
411!-----------------------------------------------------------------------
412/
413!-----------------------------------------------------------------------
414&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
415!-----------------------------------------------------------------------
416   ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water
417/
Note: See TracBrowser for help on using the repository browser.