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/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/AMM12/EXP00 – NEMO

source: branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg @ 3993

Last change on this file since 3993 was 3993, checked in by clevy, 11 years ago

Configuration setting/Step3 bugfixes,doc, and redistribute variables in namcfg and namdom see ticket:#1074

File size: 25.5 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  : AMM12 configuration namelist used to overwrite defaults values defined in SHARED/namelist_ref
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!-----------------------------------------------------------------------
5&namrun        !   parameters of the run
6!-----------------------------------------------------------------------
7   cn_exp      =  "AMM12"  !  experience name
8   nn_itend    =    2880   !  last  time step (std 1 day = 288)
9   nn_date0    =  20070101 !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
10   nn_leapy    =       1   !  Leap year calendar (1) or not (0)
11   ln_rstart   =  .true.  !  start from rest (F) or from a restart file (T)
12   cn_ocerst_in   = "amm12_restart_oce"   !  suffix of ocean restart name (input)
13   cn_ocerst_out  = "AMM12.restart.out"   !  suffix of ocean restart name (input)
14   nn_istate   =       1   !  output the initial state (1) or not (0)
15   nn_stock    =     2880  !  frequency of creation of a restart file (modulo referenced to 1)
16   nn_write    =     144   !  frequency of write in the output file   (modulo referenced to nit000)
17/
18!-----------------------------------------------------------------------
19&namcfg        !   parameters of the configuration
20!-----------------------------------------------------------------------
21   cp_cfg      =  "amm"                !  name of the configuration
22   jp_cfg      =     011               !  resolution of the configuration
23   jpidta      =     198               !  1st lateral dimension ( >= jpi )
24   jpjdta      =     224               !  2nd    "         "    ( >= jpj )
25   jpkdta      =      51               !  number of levels      ( >= jpk )
26   jpiglo      =     198               !  1st dimension of global domain --> i =jpidta
27   jpjglo      =     224               !  2nd    -                  -    --> j  =jpjdta
28   jpizoom     =       1               !  left bottom (i,j) indices of the zoom
29   jpjzoom     =       1               !  in data domain indices
30   jperio      =       0               !  lateral cond. type (between 0 and 6)
31/
32!-----------------------------------------------------------------------
33&namzgr        !   vertical coordinate
34!-----------------------------------------------------------------------
35   ln_zps      = .false.   !  z-coordinate - partial steps   (T/F)
36   ln_sco      = .true.    !  s- or hybrid z-s-coordinate    (T/F)
37/
38!-----------------------------------------------------------------------
39&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate
40!-----------------------------------------------------------------------
41   ln_s_sh94   = .false.   !  Song & Haidvogel 1994 hybrid S-sigma   (T)|
42   ln_s_sf12   = .true.    !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied
43   ln_sigcrit  = .true.    !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch
44                           !  stretching coefficients for all functions
45   rn_hc       =   50.0    !  critical depth for transition to stretched coordinates
46/
47!-----------------------------------------------------------------------
48&namdom        !   space and time domain (bathymetry, mesh, timestep)
49!-----------------------------------------------------------------------
50   rn_rdt      =  300.     !  time step for the dynamics (and tracer if nn_acc=0)
51   nn_baro     =   30      !  number of barotropic time step            ("key_dynspg_ts")
52   rn_rdtmin   =   300.          !  minimum time step on tracers (used if nn_acc=1)
53   rn_rdtmax   =   300.          !  maximum time step on tracers (used if nn_acc=1)
54   rn_rdth     =  300.           !  depth variation of tracer time step  (used if nn_acc=1)
55   jphgr_msh   =       0               !  type of horizontal mesh
56   ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1)
57   ppgphi0     =  999999.0             ! latitude  of first raw and column T-point (jphgr_msh = 1)
58   ppe1_deg    =  999999.0             !  zonal      grid-spacing (degrees)
59   ppe2_deg    =  999999.0             !  meridional grid-spacing (degrees)
60   ppe1_m      =  999999.0             !  zonal      grid-spacing (degrees)
61   ppe2_m      =  999999.0             !  meridional grid-spacing (degrees)
62   ppsur       =  999998.0             !  ORCA r4, r2 and r05 coefficients
63   ppa0        =  999998.0             ! (default coefficients)
64   ppa1        =  999998.0             !
65   ppkth       =      23.563           !
66   ppacr       =       9.0             !
67   ppdzmin     =       6.0             !  Minimum vertical spacing
68   pphmax      =    5720.              !  Maximum depth
69   ldbletanh   =  .FALSE.              !  Use/do not use double tanf function for vertical coordinates
70   ppa2        =  999999.              !  Double tanh function parameters
71   ppkth2      =  999999.              !
72   ppacr2      =  999999.
73/
74!-----------------------------------------------------------------------
75&namtsd    !   data : Temperature  & Salinity
76!-----------------------------------------------------------------------
77   ln_tsd_init   = .false.   !  Initialisation of ocean T & S with T &S input data (T) or not (F)
78   ln_tsd_tradmp = .false.   !  damping of ocean T & S toward T &S input data (T) or not (F)
79/
80!-----------------------------------------------------------------------
81&namsbc        !   Surface Boundary Condition (surface module)
82!-----------------------------------------------------------------------
83   nn_fsbc     = 1         !  frequency of surface boundary condition computation
84                           !     (also = the frequency of sea-ice model call)
85   ln_flx      = .true.    !  flux formulation       (T => fill namsbc_flx )
86   ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core)
87   nn_ice      = 0         !  =0 no ice boundary condition   ,
88   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
89   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked
90/
91!-----------------------------------------------------------------------
92&namsbc_ana    !   analytical surface boundary condition
93!-----------------------------------------------------------------------
94/
95!-----------------------------------------------------------------------
96&namsbc_flx    !   surface boundary condition : flux formulation
97!-----------------------------------------------------------------------
98!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation !
99!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
100   sn_utau     = 'amm12_utau'     ,          1        ,  'utau'      , .false.      , .false. , 'daily'   ,  ''      ,  ''
101   sn_vtau     = 'amm12_vtau'     ,          1        ,  'vtau'      , .false.      , .false. , 'daily'   ,  ''      ,  ''
102   sn_qtot     = 'amm12_flx'      ,          3        ,  'sonsfldo'  ,  .true.      , .false. , 'daily'   ,  ''      ,  ''
103   sn_qsr      = 'amm12_flx'      ,          3        ,  'soshfldo'  ,  .true.      , .false. , 'daily'   ,  ''      ,  ''
104   sn_emp      = 'amm12_flx'      ,          3        ,  'sowafldo'  ,  .true.      , .false. , 'daily'   ,  ''      ,  ''
105   cn_dir      = './fluxes/'        !  root directory for the location of the flux files
106/
107!-----------------------------------------------------------------------
108&namsbc_clio   !   namsbc_clio  CLIO bulk formulae
109!-----------------------------------------------------------------------
110/
111!-----------------------------------------------------------------------
112&namsbc_core   !   namsbc_core  CORE bulk formulae
113!-----------------------------------------------------------------------
114/
115!-----------------------------------------------------------------------
116&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae
117!-----------------------------------------------------------------------
118/
119!-----------------------------------------------------------------------
120&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_coupled")
121!-----------------------------------------------------------------------
122/
123!-----------------------------------------------------------------------
124&namtra_qsr    !   penetrative solar radiation
125!-----------------------------------------------------------------------
126   ln_traqsr   = .false.   !  Light penetration (T) or not (F)
127   nn_chldta   =      0    !  RGB : Chl data (=1) or cst value (=0)
128/
129!-----------------------------------------------------------------------
130&namsbc_rnf    !   runoffs namelist surface boundary condition
131!-----------------------------------------------------------------------
132!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation !
133!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
134   sn_rnf      = 'amm12_rivers'       ,        24         , 'rorunoff',   .false.    , .true. , 'yearly'  , ''       , ''
135   sn_cnf      = 'runoff_1m_nomask'   ,         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''
136   sn_s_rnf    = 'amm12_rivers'       ,        24         , 'rosaline',   .false.    , .true. , 'yearly'  , ''       , ''
137   sn_t_rnf    = 'amm12_rivers'       ,        24         , 'rotemper',   .false.    , .true. , 'yearly'  , ''       , ''
138   sn_dep_rnf  = 'amm12_rivers'       ,        24         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''
139   ln_rnf_mouth = .false.   !  specific treatment at rivers mouths
140   ln_rnf_depth = .true.    !  read in depth information for runoff
141   ln_rnf_tem   = .true.    !  read in temperature information for runoff
142   ln_rnf_sal   = .true.    !  read in salinity information for runoff
143/
144!-----------------------------------------------------------------------
145&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk
146!-----------------------------------------------------------------------
147!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
148!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
149   sn_apr      = 'patm'    ,         -1        ,'somslpre',    .true.      , .true.  , 'yearly'  ,  ''      ,   ''
150
151   cn_dir      = './'       !  root directory for the location of the bulk files
152   rn_pref     = 101000._wp !  reference atmospheric pressure   [N/m2]/
153   ln_ref_apr  = .false.    !  ref. pressure: global mean Patm (T) or a constant (F)
154   ln_apr_obc  = .false.    !  inverse barometer added to OBC ssh data
155/
156!-----------------------------------------------------------------------
157&namsbc_ssr    !   surface boundary condition : sea surface restoring
158!-----------------------------------------------------------------------
159!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation !
160!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !
161   sn_sst      = 'amm12_sstref'    ,        24         ,  'sst'     ,     .true.     , .false. , 'daily'   , .false. , ''       , ''
162   sn_sss      = 'sss_data'        ,        -1         ,  'sss'     ,     .true.     , .true.  , 'yearly'  , .false. , ''       , ''
163
164   cn_dir      = 'fluxes/' !  root directory for the location of the runoff files
165   nn_sstr     =     1     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
166   nn_sssr     =     0     !  add a damping     term in the surface freshwater flux (=2)
167                           !  or to SSS only (=1) or no damping term (=0)
168   rn_deds     =  -27.7    !  magnitude of the damping on salinity   [mm/day]
169/
170!-----------------------------------------------------------------------
171&namsbc_alb    !   albedo parameters
172!-----------------------------------------------------------------------
173/
174!-----------------------------------------------------------------------
175&namberg       !   iceberg parameters
176!-----------------------------------------------------------------------
177/
178!-----------------------------------------------------------------------
179&namlbc        !   lateral momentum boundary condition
180!-----------------------------------------------------------------------
181   rn_shlat    =     0     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
182                           !  free slip  !   partial slip  !   no slip   ! strong slip
183/
184!-----------------------------------------------------------------------
185&namcla        !   cross land advection
186!-----------------------------------------------------------------------
187/
188!-----------------------------------------------------------------------
189&namobc        !   open boundaries parameters                           ("key_obc")
190!-----------------------------------------------------------------------
191/
192!-----------------------------------------------------------------------
193&namagrif      !  AGRIF zoom                                            ("key_agrif")
194!-----------------------------------------------------------------------
195/
196!-----------------------------------------------------------------------
197&nam_tide      !   tide parameters (#ifdef key_tide)
198!-----------------------------------------------------------------------
199   ln_tide_pot   = .false.   !  use tidal potential forcing
200   clname(1)     =   'Q1'   !  name of constituent
201   clname(2)     =   'O1'
202   clname(3)     =   'P1'
203   clname(4)     =   'S1'
204   clname(5)     =   'K1'
205   clname(6)     =   '2N2'
206   clname(7)     =   'MU2'
207   clname(8)     =   'N2'
208   clname(9)     =   'NU2'
209   clname(10)    =   'M2'
210   clname(11)    =   'L2'
211   clname(12)    =   'T2'
212   clname(13)    =   'S2'
213   clname(14)    =   'K2'
214   clname(15)    =   'M4'
215/
216!-----------------------------------------------------------------------
217&nambdy        !  unstructured open boundaries                          ("key_bdy")
218!-----------------------------------------------------------------------
219/
220!-----------------------------------------------------------------------
221&nambdy_dta      !  open boundaries - external data           ("key_bdy")
222!-----------------------------------------------------------------------
223!              !   file name    ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
224!              !                !  (if <0  months)  !    name    !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
225   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''
226   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''
227   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''
228   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''
229   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''
230   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''
231   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''
232   cn_dir  =    './bdydta/'
233   ln_full_vel = .false.
234/
235!-----------------------------------------------------------------------
236&nambdy_tide     ! tidal forcing at open boundaries
237!-----------------------------------------------------------------------
238   filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files
239/
240!-----------------------------------------------------------------------
241&nambfr        !   bottom friction
242!-----------------------------------------------------------------------
243   nn_bfr      =    2      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
244                           !                              = 2 : nonlinear friction
245   rn_bfri2    =    2.5e-3 !  bottom drag coefficient (non linear case)
246   rn_bfeb2    =    0.0e0  !  bottom turbulent kinetic energy background  (m2/s2)
247   ln_loglayer = .false.   !  loglayer bottom friction (only effect when nn_bfr = 2)
248   rn_bfrz0    =    0.003  !  bottom roughness (only effect when ln_loglayer = .true.)
249/
250!-----------------------------------------------------------------------
251&nambbc        !   bottom temperature boundary condition
252!-----------------------------------------------------------------------
253   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
254/
255!-----------------------------------------------------------------------
256&nambbl        !   bottom boundary layer scheme
257!-----------------------------------------------------------------------
258   nn_bbl_ldf  =  0      !  diffusive bbl (=1)   or not (=0)
259/
260!-----------------------------------------------------------------------
261&nameos        !   ocean physical parameters
262!-----------------------------------------------------------------------
263/
264!-----------------------------------------------------------------------
265&namtra_adv    !   advection scheme for tracer
266!-----------------------------------------------------------------------
267/
268!----------------------------------------------------------------------------------
269&namtra_ldf    !   lateral diffusion scheme for tracers
270!----------------------------------------------------------------------------------
271   ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T)
272   ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp")
273   rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]
274   rn_aht_0         =    50.    !  horizontal eddy diffusivity for tracers [m2/s]
275/
276!-----------------------------------------------------------------------
277&namtra_dmp    !   tracer: T & S newtonian damping
278!-----------------------------------------------------------------------
279   ln_tradmp   =  .false.  !  add a damping termn (T) or not (F)
280/
281!-----------------------------------------------------------------------
282&namdyn_adv    !   formulation of the momentum advection
283!-----------------------------------------------------------------------
284/
285!-----------------------------------------------------------------------
286&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
287!-----------------------------------------------------------------------
288/
289!-----------------------------------------------------------------------
290&namdyn_hpg    !   Hydrostatic pressure gradient option
291!-----------------------------------------------------------------------
292   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
293   ln_hpg_prj  = .true.   !  s-coordinate (Pressure Jacobian scheme)
294/
295!-----------------------------------------------------------------------
296!namdyn_spg    !   surface pressure gradient   (CPP key only)
297!-----------------------------------------------------------------------
298!                          !  explicit free surface                     ("key_dynspg_exp")
299!                          !  filtered free surface                     ("key_dynspg_flt")
300!                          !  split-explicit free surface               ("key_dynspg_ts")
301
302!-----------------------------------------------------------------------
303&namdyn_ldf    !   lateral diffusion on momentum
304!-----------------------------------------------------------------------
305   !                       !  Type of the operator :
306   ln_dynldf_bilap  =  .true.   !  bilaplacian operator
307                           !  Coefficient
308   rn_ahm_0_lap     = 60.0      !  horizontal laplacian eddy viscosity   [m2/s]
309   rn_ahmb_0        =  0.0      !  background eddy viscosity for ldf_iso [m2/s]
310   rn_ahm_0_blp     = -1.0e+10  !  horizontal bilaplacian eddy viscosity [m4/s]
311/
312!-----------------------------------------------------------------------
313&namzdf        !   vertical physics
314!-----------------------------------------------------------------------
315   rn_avm0     =   0.1e-6  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
316   rn_avt0     =   0.1e-6  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
317   ln_zdfevd   = .false.   !  enhanced vertical diffusion (evd) (T) or not (F)
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/
328!------------------------------------------------------------------------
329&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally:
330!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
331/
332!-----------------------------------------------------------------------
333&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls")
334!-----------------------------------------------------------------------
335   rn_charn =  100000.     !  Charnock constant for wb induced roughness length
336/
337!-----------------------------------------------------------------------
338&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
339!-----------------------------------------------------------------------
340/
341!-----------------------------------------------------------------------
342&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
343!-----------------------------------------------------------------------
344   ln_tmx_itf  = .FALSE.   !  ITF specific parameterisation
345/
346!-----------------------------------------------------------------------
347&namsol        !   elliptic solver / island / free surface
348!-----------------------------------------------------------------------
349/
350!-----------------------------------------------------------------------
351&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
352!-----------------------------------------------------------------------
353/
354!-----------------------------------------------------------------------
355&namctl        !   Control prints & Benchmark
356!-----------------------------------------------------------------------
357/
358!-----------------------------------------------------------------------
359&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
360!-----------------------------------------------------------------------
361/
362!-----------------------------------------------------------------------
363&namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra")
364!              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor")
365!-----------------------------------------------------------------------
366/
367!-----------------------------------------------------------------------
368&namflo       !   float parameters                                      ("key_float")
369!-----------------------------------------------------------------------
370/
371!-----------------------------------------------------------------------
372&namptr       !   Poleward Transport Diagnostic
373!-----------------------------------------------------------------------
374   ln_diaznl  = .false.    !  Add zonal means and meridional stream functions
375   ln_subbas  = .false.    !  Atlantic/Pacific/Indian basins computation (T) or not
376                           !  (orca configuration only, need input basins mask file named "subbasins.nc"
377   ln_ptrcomp = .false.    !  Add decomposition : overturning
378/
379!-----------------------------------------------------------------------
380&namhsb       !  Heat and salt budgets
381!-----------------------------------------------------------------------
382/
383!-----------------------------------------------------------------------
384&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm')
385!-----------------------------------------------------------------------
386/
387!-----------------------------------------------------------------------
388&namdct        ! transports through sections
389!-----------------------------------------------------------------------
390/
391!-----------------------------------------------------------------------
392&namobs       !  observation usage switch                               ('key_diaobs')
393!-----------------------------------------------------------------------
394/
395!-----------------------------------------------------------------------
396&nam_asminc   !   assimilation increments                               ('key_asminc')
397!-----------------------------------------------------------------------
398/
399!-----------------------------------------------------------------------
400&namsbc_wave   ! External fields from wave model
401!-----------------------------------------------------------------------
402/
403!-----------------------------------------------------------------------
404&namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed)
405!-----------------------------------------------------------------------
406/
Note: See TracBrowser for help on using the repository browser.