New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
namelist in branches/nemo_v3_3_beta/NEMOGCM/CONFIG/GYRE/EXP00 – NEMO

source: branches/nemo_v3_3_beta/NEMOGCM/CONFIG/GYRE/EXP00/namelist @ 2375

Last change on this file since 2375 was 2375, checked in by gm, 13 years ago

v3.3beta: zdftke namelist simplification associated with an update of DOC

  • Property svn:keywords set to Author Date Id Revision
File size: 55.3 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  :  1 - run manager      (namrun)
3!! namelists    2 - Domain           (namzgr, namzgr_sco, namdom)
4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core
5!!                                    namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namsbc_alb)
6!!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide)
7!!              5 - bottom  boundary (nambfr, nambbc, nambbl)
8!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp)
9!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
10!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx)
11!!              9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr)
12!!              9 - miscellaneous    (namsol, nammpp, nammpp_dyndist, namctl)
13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
14!  CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE.
15
16!!======================================================================
17!!                   ***  Run management namelists  ***
18!!======================================================================
19!!   namrun            parameters of the run
20!!======================================================================
21
22!-----------------------------------------------------------------------
23&namrun        !   parameters of the run
24!-----------------------------------------------------------------------
25   nn_no       =       0   !  job number
26   cn_exp      =  "GYRE"   !  experience name
27   nn_it000    =       1   !  first time step
28   nn_itend    =    4320   !  last  time step (std 5475)
29   nn_date0    =  010101   !  initial calendar date yymmdd (used if nrstdt=1)
30   nn_leapy    =      30   !  Leap year calendar (1) or not (0)
31   nn_istate   =       0   !  output the initial state (1) or not (0)
32   nn_stock    =    4320   !  frequency of creation of a restart file (modulo referenced to 1)
33   nn_write    =      60   !  frequency of write in the output file   (modulo referenced to nit000)
34   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T)
35   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%)
36   ln_clobber  = .false.   !  clobber (overwrite) an existing file
37   nn_chunksz  = 0         !  chunksize (bytes) for NetCDF file (working only with iom_nf90 routines)
38   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
39   nn_rstctl   =       0       !  restart control = 0 nit000 is not compared to the restart file value
40                               !                  = 1 use ndate0 in namelist (not the value in the restart file)
41                               !                  = 2 calendar parameters read in the restart file
42   cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input)
43   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output)
44/
45!!======================================================================
46!!   namnc4            netcdf4 chunking and compression settings
47!!======================================================================
48!-----------------------------------------------------------------------
49&namnc4         !  netcdf4 chunking and compression settings     ("key_netcdf4")
50                !  (ignored if "key_netcdf4" is not used)
51!-----------------------------------------------------------------------
52   nn_nchunks_i   =   4       !  number of chunks in i-dimension
53   nn_nchunks_j   =   4       !  number of chunks in j-dimension
54   nn_nchunks_k   =   31      !  number of chunks in k-dimension
55                              !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
56                              !  is optimal for postprocessing which works exclusively with horizontal slabs
57   ln_nc4zip      =   .TRUE.  !  (T) use netcdf4 chunking and compression
58                              !  (F) ignore chunking information and produce netcdf3-compatible files
59/
60!!======================================================================
61!!                      ***  Domain namelists  ***
62!!======================================================================
63!!   namzgr       vertical coordinate
64!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate
65!!   namdom       space and time domain (bathymetry, mesh, timestep)
66!!======================================================================
67
68!-----------------------------------------------------------------------
69&namzgr        !   vertical coordinate
70!-----------------------------------------------------------------------
71   ln_zco      = .true.    !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined)
72   ln_zps      = .false.   !  z-coordinate - partial steps   (T/F)
73   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F)
74/
75!-----------------------------------------------------------------------
76&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate
77!-----------------------------------------------------------------------
78   rn_sbot_min =  300.     !  minimum depth of s-bottom surface (>0) (m)
79   rn_sbot_max = 5250.     !  maximum depth of s-bottom surface (= ocean depth) (>0) (m)
80   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20)
81   rn_thetb    =    0.75   !  bottom control parameter  (0<=thetb<= 1)
82   rn_rmax     =    0.15   !  maximum cut-off r-value allowed (0<r_max<1)
83   ln_s_sigma  = .false.   !  hybrid s-sigma coordinates
84   rn_bb       =    0.8    !  stretching with s-sigma
85   rn_hc       =  150.0    !  critical depth with s-sigma
86/
87!-----------------------------------------------------------------------
88&namdom        !   space and time domain (bathymetry, mesh, timestep)
89!-----------------------------------------------------------------------
90   nn_bathy    =    0      !  compute (=0) or read(=1) the bathymetry file
91   nn_closea    =   0      !  closed seas and lakes are removed (=0) or kept (=1) from the ORCA domain
92   nn_msh      =    0      !  create (=1) a mesh file (coordinates, scale factors, masks) or not (=0)
93   rn_e3zps_min=    5.     !  the thickness of the partial step is set larger than the minimum
94   rn_e3zps_rat=    0.1    !  of e3zps_min and e3zps_rat * e3t   (N.B. 0<e3zps_rat<1)
95                           !
96   rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nacc=0)   ==> 5760
97   nn_baro     =   60      !  number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts")
98   rn_atfp     =    0.1    !  asselin time filter parameter
99   nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k)
100                                 !                          =0, not used, rdt = rdttra
101   rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nacc=1)
102   rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nacc=1)
103   rn_rdth     =  800.           !  depth variation of tracer time step  (used if nacc=1)
104/
105!!======================================================================
106!!            ***  Surface Boundary Condition namelists  ***
107!!======================================================================
108!!   namsbc        surface boundary condition
109!!   namsbc_ana    analytical         formulation
110!!   namsbc_flx    flux               formulation
111!!   namsbc_clio   CLIO bulk formulea formulation
112!!   namsbc_core   CORE bulk formulea formulation
113!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled")
114!!   namtra_qsr    penetrative solar radiation
115!!   namsbc_rnf    river runoffs
116!!   namsbc_ssr    sea surface restoring term (for T and/or S)
117!!   namsbc_alb    albedo parameters
118!!======================================================================
119
120!-----------------------------------------------------------------------
121&namsbc        !   Surface Boundary Condition (surface module)
122!-----------------------------------------------------------------------
123   nn_fsbc     = 1         !  frequency of surface boundary condition computation
124                           !               (= the frequency of sea-ice model call)
125   ln_ana      = .true.    !  analytical formulation (T => fill namsbc_ana )
126   ln_flx      = .false.   !  flux formulation       (T => fill namsbc_flx )
127   ln_blk_clio = .false.   !  CLIO bulk formulation  (T => fill namsbc_clio)
128   ln_blk_core = .false.   !  CORE bulk formulation  (T => fill namsbc_core)
129   ln_cpl      = .false.   !  Coupled formulation    (T => fill namsbc_cpl )
130   nn_ice      = 0         !  =0 no ice boundary condition   ,
131                           !  =1 use observed ice-cover      ,
132                           !  =2 ice-model used                             ("key_lim3" or "key_lim2)
133   ln_dm2dc    = .false.   !  daily mean to diurnal cycle short wave (qsr)
134   ln_rnf      = .false.   !  runoffs (T => fill namsbc_rnf)
135   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
136   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked                              ,
137                           !                     =1 global mean of e-p-r set to zero at each nn_fsbc time step   ,
138                           !                     =2 annual global mean of e-p-r set to zero
139                           !                     =3 global emp set to zero and spread out over erp area
140/
141!-----------------------------------------------------------------------
142&namsbc_ana    !   analytical surface boundary condition
143!-----------------------------------------------------------------------
144   nn_tau000   =   100     !  gently increase the stress over the first ntau_rst time-steps
145   rn_utau0    =   0.1e0   !  uniform value for the i-stress
146   rn_vtau0    =   0.e0    !  uniform value for the j-stress
147   rn_qns0     =   0.e0    !  uniform value for the total heat flux
148   rn_qsr0     =   0.e0    !  uniform value for the solar radiation
149   rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P)
150/
151!-----------------------------------------------------------------------
152&namsbc_flx    !   surface boundary condition : flux formulation
153!-----------------------------------------------------------------------
154!              !   file name  ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
155!              !              !  (if <0  months)  !    name    !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
156   sn_utau     = 'utau'       ,        24         ,    'utau'  ,     .false.    , .false. , 'yearly'  , ''       , ''
157   sn_vtau     = 'vtau'       ,        24         ,    'vtau'  ,     .false.    , .false. , 'yearly'  , ''       , ''
158   sn_qtot     = 'qtot'       ,        24         ,    'qtot'  ,     .false.    , .false. , 'yearly'  , ''       , ''
159   sn_qsr      = 'qsr'        ,        24         ,    'qsr'   ,     .false.    , .false. , 'yearly'  , ''       , ''
160   sn_emp      = 'emp'        ,        24         ,    'emp'   ,     .false.    , .false. , 'yearly'  , ''       , ''
161!
162   cn_dir      = './'      !  root directory for the location of the flux files
163/     
164!-----------------------------------------------------------------------
165&namsbc_clio   !   namsbc_clio  CLIO bulk formulea
166!-----------------------------------------------------------------------
167!              !   file name  ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
168!              !              !  (if <0  months)  !    name    !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
169   sn_utau     = 'taux_1m'    ,       -1          , 'sozotaux' ,    .true.      , .true.  , 'yearly'  , ''       , ''
170   sn_vtau     = 'tauy_1m'    ,       -1          , 'sometauy' ,    .true.      , .true.  , 'yearly'  , ''       , ''
171   sn_wndm     = 'flx'        ,       -1          , 'socliowi' ,    .true.      , .true.  , 'yearly'  , ''       , ''
172   sn_tair     = 'flx'        ,       -1          , 'socliot2' ,    .true.      , .true.  , 'yearly'  , ''       , ''
173   sn_humi     = 'flx'        ,       -1          , 'socliohu' ,    .true.      , .true.  , 'yearly'  , ''       , ''
174   sn_ccov     = 'flx'        ,       -1          , 'socliocl' ,    .false.     , .true.  , 'yearly'  , ''       , ''
175   sn_prec     = 'flx'        ,       -1          , 'socliopl' ,    .false.     , .true.  , 'yearly'  , ''       , ''
176!
177   cn_dir      = './'      !  root directory for the location of the bulk files are
178/
179!-----------------------------------------------------------------------
180&namsbc_core   !   namsbc_core  CORE bulk formulea
181!-----------------------------------------------------------------------
182!              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights              ! rotation !
183!              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F)  ! 'monthly' ! filename             ! pairing  !
184   sn_wndi     = 'u10_core'   ,       -1          , 'u10'      ,    .true.      , .true.  , 'yearly'  ,'bicubic_weights_orca2.nc' , 'U1'
185   sn_wndj     = 'v10_core'   ,       -1          , 'v10'      ,    .true.      , .true.  , 'yearly'  ,'bicubic_weights_orca2.nc' , 'V1'
186   sn_qsr      = 'qsw_core'   ,       -1          , 'swdn'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', ''
187   sn_qlw      = 'qlw_core'   ,       -1          , 'lwdn'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', ''
188   sn_tair     = 't2_core'    ,       -1          , 't2'       ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', ''
189   sn_humi     = 'q2_core'    ,       -1          , 'q2'       ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', ''
190   sn_prec     = 'precip_core',       -1          , 'precip'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', ''
191   sn_snow     = 'snow_core'  ,       -1          , 'snow'     ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', ''
192   sn_tdif     = 'taudif_core',       24          , 'taudif'   ,    .true.      , .true.  , 'yearly'  ,'bilinear_weights_orca2.nc', ''
193!
194   cn_dir      = './'      !  root directory for the location of the bulk files
195   ln_2m       = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F)
196   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data ?
197   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow)
198/
199!-----------------------------------------------------------------------
200&namsbc_cpl    !   coupled ocean/atmosphere model                        ("key_coupled")
201!-----------------------------------------------------------------------
202/
203!-----------------------------------------------------------------------
204&namtra_qsr    !   penetrative solar radiation
205!-----------------------------------------------------------------------
206!              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
207!              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
208   sn_chl      = 'chlorophyll',        -1         , 'CHLA'     ,    .true.      , .true.  , 'yearly'  , ''       , ''
209 
210   cn_dir      = './'      !  root directory for the location of the runoff files
211   ln_traqsr   = .true.    !  Light penetration (T) or not (F)
212   ln_qsr_rgb  = .false.   !  RGB (Red-Green-Blue) light penetration
213   ln_qsr_2bd  = .true.    !  2 bands              light penetration
214   ln_qsr_bio  = .false.   !  bio-model light penetration
215   nn_chldta   =      0    !  RGB : Chl data (=1) or cst value (=0)
216   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
217   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
218   rn_si1      =   23.0    !  2 bands: longest depth of extinction
219/
220!-----------------------------------------------------------------------
221&namsbc_rnf    !   runoffs namelist surface boundary condition
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   sn_rnf      = 'runoff_1m_nomask' ,  -1         , 'sorunoff' ,    .true.      , .true.  , 'yearly'  , ''       , ''
226   sn_cnf      = 'runoff_1m_nomask' ,   0         , 'socoefr'  ,    .false.     , .true.  , 'yearly'  , ''       , ''
227   sn_s_rnf    = 'runoffs'          ,  24         , 'rosaline' ,    .true.      , .true.  ,   'yearly', ''       , ''
228   sn_t_rnf    = 'runoffs'          ,  24         , 'rotemper' ,    .true.      , .true.  ,   'yearly', ''       , ''
229   sn_dep_rnf  = 'runoffs'          ,   0         , 'rodepth'  ,    .false.     , .true.  ,   'yearly', ''       , ''
230   cn_dir       = './'      !  root directory for the location of the runoff files
231   ln_rnf_emp   =   .false. !  runoffs included into precipitation field (T) or into a file (F)
232   ln_rnf_mouth =   .false. !  specific treatment at rivers mouths
233   rn_hrnf      =   0.e0    !  depth over which enhanced vertical mixing is used
234   rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s]
235   rn_rfact     =   1.e0    !  multiplicative factor for runoff
236   ln_rnf_depth =  .false.  !  read in depth information for runoff
237   ln_rnf_tem   =  .false.  !  read in temperature information for runoff
238   ln_rnf_sal   =  .false.  !  read in salinity information for runoff
239/
240!-----------------------------------------------------------------------
241&namsbc_ssr    !   surface boundary condition : sea surface restoring
242!-----------------------------------------------------------------------
243!              !   file name  ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation !
244!              !              !  (if <0  months)  !   name     !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  !
245   sn_sst      = 'sst_data'   ,        24         ,  'sst'     ,     .false.    , .false. , 'yearly'  , ''       , ''
246   sn_sss      = 'sss_data'   ,        -1         ,  'sss'     ,     .true.     , .true.  , 'yearly'  , ''       , ''
247   
248   cn_dir      = './'      !  root directory for the location of the runoff files
249   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
250   nn_sssr     =     0     !  add a damping     term in the surface freshwater flux (=2)
251                           !  or to SSS only (=1) or no damping term (=0)
252   rn_dqdt     =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K]
253   rn_deds     =   -27.7   !  magnitude of the damping on salinity   [mm/day]
254   ln_sssr_bnd =   .false. !  flag to bound erp term (associated with nn_sssr=2)
255   rn_sssr_bnd =   0.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
256/     
257!-----------------------------------------------------------------------
258&namsbc_alb    !   albedo parameters
259!-----------------------------------------------------------------------
260   rn_cloud    =    0.06   !  cloud correction to snow and ice albedo
261   rn_albice   =    0.53   !  albedo of melting ice in the arctic and antarctic
262   rn_alphd    =    0.80   !  coefficients for linear interpolation used to
263   rn_alphc    =    0.65   !  compute albedo between two extremes values
264   rn_alphdi   =    0.72   !  (Pyane, 1972)
265/
266
267!!======================================================================
268!!               ***  Lateral boundary condition  ***
269!!======================================================================
270!!   namlbc        lateral momentum boundary condition
271!!   namcla        cross land advection
272!!   namobc        open boundaries parameters                           ("key_obc")
273!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")
274!!   nambdy        Unstructured open boundaries                         ("key_bdy")
275!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides")
276!!======================================================================
277
278!-----------------------------------------------------------------------
279&namlbc        !   lateral momentum boundary condition
280!-----------------------------------------------------------------------
281   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
282                           !  free slip  !   partial slip  !   no slip   ! strong slip
283/
284!-----------------------------------------------------------------------
285&namcla        !   cross land advection
286!-----------------------------------------------------------------------
287   nn_cla      =    0      !  advection between 2 ocean pts separates by land
288/
289!-----------------------------------------------------------------------
290&namobc        !   open boundaries parameters                           ("key_obc")
291!-----------------------------------------------------------------------
292    ln_obc_clim= .false.   !  climatological obc data files (T) or not (F)
293    ln_vol_cst = .true.    !  impose the total volume conservation (T) or not (F)
294    ln_obc_fla = .false.   !  Flather open boundary condition
295    nn_obcdta  =    1      !  = 0 the obc data are equal to the initial state
296                           !  = 1 the obc data are read in 'obc.dta' files
297    cn_obcdta  = 'annual'  !  set to annual if obc datafile hold 1 year of data
298                           !  set to monthly if obc datafile hold 1 month of data
299    rn_dpein   =    1.     !  damping time scale for inflow at east  open boundary
300    rn_dpwin   =    1.     !     -           -         -       west    -      -
301    rn_dpnin   =    1.     !     -           -         -       north   -      -
302    rn_dpsin   =    1.     !     -           -         -       south   -      -
303    rn_dpeob   = 3000.     !  time relaxation (days) for the east  open boundary
304    rn_dpwob   =   15.     !     -           -         -     west    -      -
305    rn_dpnob   = 3000.     !     -           -         -     north   -      -
306    rn_dpsob   =   15.     !     -           -         -     south   -      -
307    rn_volemp  =  1.       !  = 0 the total volume change with the surface flux (E-P-R)
308                           !  = 1 the total volume remains constant
309/
310!-----------------------------------------------------------------------
311&namagrif      !  AGRIF zoom                                            ("key_agrif")
312!-----------------------------------------------------------------------
313    nn_cln_update = 3      !  baroclinic update frequency
314    ln_spc_dyn    = .true. !  use 0 as special value for dynamics
315    rn_sponge_tra = 2880.  !  coefficient for tracer   sponge layer [s]
316    rn_sponge_dyn = 2880.  !  coefficient for dynamics sponge layer [s]
317/
318!-----------------------------------------------------------------------
319&nambdy        !  unstructured open boundaries                          ("key_bdy")
320!-----------------------------------------------------------------------
321    cn_mask    =  ''                        !  name of mask file (if ln_bdy_mask=.TRUE.)
322    cn_dta_frs_T  = 'bdydata_grid_T.nc'     !  name of data file (T-points)
323    cn_dta_frs_U  = 'bdydata_grid_U.nc'     !  name of data file (U-points)
324    cn_dta_frs_V  = 'bdydata_grid_V.nc'     !  name of data file (V-points)
325    cn_dta_fla_T  = 'bdydata_bt_grid_T.nc'  !  name of data file for Flather condition (T-points)
326    cn_dta_fla_U  = 'bdydata_bt_grid_U.nc'  !  name of data file for Flather condition (U-points)
327    cn_dta_fla_V  = 'bdydata_bt_grid_V.nc'  !  name of data file for Flather condition (V-points)
328    ln_clim    = .false.                    !  contain 1 (T) or 12 (F) time dumps and be cyclic
329    ln_vol     = .false.                    !  total volume correction (see volbdy parameter)
330    ln_mask    = .false.                    !  boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F)
331    ln_tides   = .false.                    !  Apply tidal harmonic forcing with Flather condition
332    ln_dyn_fla = .false.                    !  Apply Flather condition to velocities
333    ln_tra_frs = .false.                    !  Apply FRS condition to temperature and salinity
334    ln_dyn_frs = .false.                    !  Apply FRS condition to velocities
335    nn_rimwidth    =  9                     !  width of the relaxation zone
336    nn_dtactl      =  1                     !  = 0, bdy data are equal to the initial state
337                                            !  = 1, bdy data are read in 'bdydata   .nc' files
338    nn_volctl      =  0                     !  = 0, the total water flux across open boundaries is zero
339                                            !  = 1, the total volume of the system is conserved
340/
341!-----------------------------------------------------------------------
342&nambdy_tide     ! tidal forcing at unstructured boundaries             
343!-----------------------------------------------------------------------
344    filtide      = 'bdytide_'           !  file name root of tidal forcing files
345    tide_cpt     = 'M2','S1'            !  names of tidal components used
346    tide_speed   = 28.984106, 15.000001 !  phase speeds of tidal components (deg/hour)
347    ln_tide_date = .false.              !  adjust tidal harmonics for start date of run
348/
349
350!!======================================================================
351!!                 ***  Bottom boundary condition  ***
352!!======================================================================
353!!   nambfr        bottom friction
354!!   nambbc        bottom temperature boundary condition               
355!!   nambbl        bottom boundary layer scheme                         ("key_trabbl")
356!!======================================================================
357
358!-----------------------------------------------------------------------
359&nambfr        !   bottom friction
360!-----------------------------------------------------------------------
361   nn_bfr      =    2      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction
362                           !                              = 2 : nonlinear friction
363   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case)
364   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case)
365   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m^2/s^2)
366   ln_bfr2d    =   .false. !  horizontal variation of the bottom friction coef (read a 2D mask file )
367   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d = .true.)
368/
369!-----------------------------------------------------------------------
370&nambbc        !   bottom temperature boundary condition
371!-----------------------------------------------------------------------
372   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
373   nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux
374                           !     = 1 constant flux
375                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
376   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
377/
378!-----------------------------------------------------------------------
379&nambbl        !   bottom boundary layer scheme
380!-----------------------------------------------------------------------
381   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0)
382   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0)
383   rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s]
384   rn_gambbl   =  10.    !  advective bbl coefficient                 [s]
385/
386!!======================================================================
387!!                        Tracer (T & S ) namelists
388!!======================================================================
389!!   nameos        equation of state
390!!   namtra_adv    advection scheme
391!!   namtra_ldf    lateral diffusion scheme
392!!   namtra_dmp    T & S newtonian damping                              ("key_tradmp")
393!!======================================================================
394
395!-----------------------------------------------------------------------
396&nameos        !   ocean physical parameters
397!-----------------------------------------------------------------------
398   nn_eos      =    2      !  type of equation of state and Brunt-Vaisala frequency
399                           !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) )
400                           !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T )
401                           !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T )
402   rn_alpha    =    2.e-4  !  thermal expension coefficient (neos= 1 or 2)
403   rn_beta     =    7.7e-4 !  saline  expension coefficient (neos= 2)
404/
405!-----------------------------------------------------------------------
406&namtra_adv    !   advection scheme for tracer
407!-----------------------------------------------------------------------
408   ln_traadv_cen2   =  .false.  !  2nd order centered scheme   
409   ln_traadv_tvd    =  .true.   !  TVD scheme               
410   ln_traadv_muscl  =  .false.  !  MUSCL scheme             
411   ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
412   ln_traadv_ubs    =  .false.  !  UBS scheme                 
413/
414!-----------------------------------------------------------------------
415&namtra_ldf    !   lateral diffusion scheme for tracer
416!-----------------------------------------------------------------------
417                           !  Type of the operator :
418   ln_traldf_lap    =  .true.   !     laplacian operator       
419   ln_traldf_bilap  =  .false.  !     bilaplacian operator     
420                           !  Direction of action  :
421   ln_traldf_level  =  .false.  !     iso-level               
422   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T)
423   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp")
424   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp")   ! UNDER TEST, DO NOT USE
425   ln_traldf_gdia   =  .false.  !     griffies operator strfn diagnostics (require "key_ldfslp") ! UNDER TEST, DO NOT USE
426                           !  Coefficient
427   rn_aht_0         =  1000.    !     horizontal eddy diffusivity for tracers [m2/s]
428   rn_ahtb_0        =     0.    !     background eddy diffusivity for ldf_iso [m2/s]
429   rn_aeiv_0        =     0.    !     eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv")
430/
431!-----------------------------------------------------------------------
432&namtra_dmp    !   tracer: T & S newtonian damping                      ('key_tradmp')
433!-----------------------------------------------------------------------
434   nn_hdmp     =   -1      !  horizontal shape =-1, damping in Med and Red Seas only
435                           !                   =XX, damping poleward of XX degrees (XX>0)
436                           !                      + F(distance-to-coast) + Red and Med Seas
437   nn_zdmp     =    1      !  vertical   shape =0    damping throughout the water column
438                           !                   =1 no damping in the mixing layer (kz  criteria)
439                           !                   =2 no damping in the mixed  layer (rho crieria)
440   rn_surf     =   50.     !  surface time scale of damping   [days]
441   rn_bot      =  360.     !  bottom  time scale of damping   [days]
442   rn_dep      =  800.     !  depth of transition between rn_surf and rn_bot [meters]
443   nn_file     =    1      !  create a damping.coeff NetCDF file (=1) or not (=0)
444/
445!!======================================================================
446!!                      ***  Dynamics namelists  ***
447!!======================================================================
448!!   namdyn_adv    formulation of the momentum advection
449!!   namdyn_vor    advection scheme
450!!   namdyn_hpg    hydrostatic pressure gradient
451!!   namdyn_spg    surface pressure gradient                            (CPP key only)
452!!   namdyn_ldf    lateral diffusion scheme
453!!======================================================================
454
455!-----------------------------------------------------------------------
456&namdyn_adv    !   formulation of the momentum advection
457!-----------------------------------------------------------------------
458   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
459   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
460   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
461
462!-----------------------------------------------------------------------
463&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
464!-----------------------------------------------------------------------
465   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
466   ln_dynvor_ens = .false. !  energy conserving scheme   
467   ln_dynvor_mix = .false. !  mixed scheme               
468   ln_dynvor_een = .false. !  energy & enstrophy scheme 
469/
470!-----------------------------------------------------------------------
471&namdyn_hpg    !   Hydrostatic pressure gradient option
472!-----------------------------------------------------------------------
473   ln_hpg_zco  = .true.    !  z-coordinate - full steps                   
474   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
475   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
476   ln_hpg_hel  = .false.   !  s-coordinate (helsinki modification)
477   ln_hpg_wdj  = .false.   !  s-coordinate (weighted density jacobian)
478   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
479   ln_hpg_rot  = .false.   !  s-coordinate (ROTated axes scheme)
480   rn_gamma    = 0.e0      !  weighting coefficient (wdj scheme)
481   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T)
482                                 !           centered      time scheme  (F)
483   nn_dynhpg_rst =  0            ! =1 dynhpg restartable restart or not (=0)
484/
485!-----------------------------------------------------------------------
486!namdyn_spg    !   surface pressure gradient   (CPP key only)
487!-----------------------------------------------------------------------
488!                          !  explicit free surface                     ("key_dynspg_exp")
489!                          !  filtered free surface                     ("key_dynspg_flt")
490!                          !  split-explicit free surface               ("key_dynspg_ts")
491
492!-----------------------------------------------------------------------
493&namdyn_ldf    !   lateral diffusion on momentum
494!-----------------------------------------------------------------------
495                           !  Type of the operator :
496   ln_dynldf_lap    =  .true.   !  laplacian operator         
497   ln_dynldf_bilap  =  .false.  !  bilaplacian operator   
498                           !  Direction of action  :
499   ln_dynldf_level  =  .false.  !  iso-level               
500   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.)
501   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp")
502                           !  Coefficient
503   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s]
504   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s]
505   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s]
506/
507!!======================================================================
508!!             Tracers & Dynamics vertical physics namelists
509!!======================================================================
510!!       namzdf        vertical physics
511!!       namzdf_ric    richardson number dependent vertical mixing      ("key_zdfric"      )
512!!       namzdf_tke    TKE dependent vertical mixing                    ("key_zdftke"      )
513!!       namzdf_kpp    KPP dependent vertical mixing                    ("key_zdfkpp"      )
514!!       namzdf_ddm    double diffusive mixing parameterization         ("key_zdfddm"      )
515!!       namzdf_tmx    tidal mixing parameterization                    ("key_zdftmx"      )
516!!======================================================================
517
518!-----------------------------------------------------------------------
519&namzdf        !   vertical physics
520!-----------------------------------------------------------------------
521   rn_avm0     =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
522   rn_avt0     =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
523   nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0)
524   nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0)
525   ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F)
526   nn_evdm     =    1      !  evd apply on tracer (=0) or on tracer and momentum (=1)
527   rn_avevd    =  100.     !  evd mixing coefficient [m2/s]
528   ln_zdfnpc   = .false.   !  Non-Penetrative algorithm (T) or not (F)
529   nn_npc      =    1            !  frequency of application of npc
530   nn_npcp     =  365            !  npc control print frequency
531   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
532   nn_zdfexp   =    3            !  number of sub-timestep for ln_zdfexp=T
533/
534!-----------------------------------------------------------------------
535&namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" )
536!-----------------------------------------------------------------------
537   rn_avmri    = 100.e-4   !  maximum value of the vertical viscosity
538   rn_alp      =   5.      !  coefficient of the parameterization
539   nn_ric      =   2       !  coefficient of the parameterization
540/
541!-----------------------------------------------------------------------
542&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
543!-----------------------------------------------------------------------
544   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
545   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
546   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
547   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
548   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
549   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
550                           !                 = 1 bounded by the local vertical scale factor
551                           !                 = 2 first vertical derivative of mixing length bounded by 1
552                           !                 = 3 as =2 with distinct disspipative an mixing length scale
553   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
554   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
555   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
556   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
557   rn_lc       =   0.15    !  coef. associated to Langmuir cells
558   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves
559                           !        = 0 no penetration
560                           !        = 1 add a tke source below the ML
561                           !        = 2 add a tke source just at the base of the ML
562                           !        = 3 as = 1 applied on HF part of the stress    ("key_coupled")
563   rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
564   nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
565                           !        = 0  constant 10 m length scale
566                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees
567/
568!------------------------------------------------------------------------
569&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionnally:
570!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb")
571   ln_kpprimix = .true.    !  shear instability mixing
572   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s]
573   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s]
574   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability
575   rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s]
576   rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]
577   rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s]
578   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv
579   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt
580/
581!-----------------------------------------------------------------------
582&namzdf_gls                !   GLS vertical diffusion  ("key_zdfgls")
583!-----------------------------------------------------------------------
584   rn_emin        = 1.e-6  !  minimum value of e   [m2/s2]
585   rn_epsmin      = 1.e-12 !  minimum value of eps [m2/s3]
586   ln_length_lim  = .true. !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
587   rn_clim_galp   = 0.53   !  galperin limit
588   ln_crban       = .TRUE. !  Use Craig & Banner (1994) surface wave mixing parametrisation
589   ln_sigpsi      = .TRUE. !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
590   rn_crban       = 100.   !  Craig and Banner 1994 constant for wb tke flux
591   rn_charn       = 70000. !  Charnock constant for wb induced roughness length
592   nn_tkebc_surf  =   1    !  surface tke condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg)
593   nn_tkebc_bot   =   1    !  bottom tke condition (0/1=Dir/Neum)
594   nn_psibc_surf  =   1    !  surface psi condition (0/1/2=Dir/Neum/Dir Mellor-Blumberg)
595   nn_psibc_bot   =   1    !  bottom psi condition (0/1=Dir/Neum)
596   nn_stab_func   =   2    !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
597   nn_clos        =   1    !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
598/
599!-----------------------------------------------------------------------
600&namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm")
601!-----------------------------------------------------------------------
602   rn_avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
603   rn_hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
604/
605!-----------------------------------------------------------------------
606&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx")
607!-----------------------------------------------------------------------
608   rn_htmx     = 500.      !  vertical decay scale for turbulence (meters)
609   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1)
610   rn_tfe      = 0.333     !  tidal dissipation efficiency
611   rn_me       = 0.2       !  mixing efficiency
612   ln_tmx_itf  = .FALSE.   !  ITF specific parameterisation
613   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency
614/
615!!======================================================================
616!!                  ***  Miscelaneous namelists  ***
617!!======================================================================
618!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi)
619!!   nammpp_dyndist    Massively Parallel domain decomposition          ("key_agrif" && "key_mpp_dyndist")
620!!   namctl            Control prints & Benchmark
621!!   namsol            elliptic solver / island / free surface
622!!======================================================================
623
624!-----------------------------------------------------------------------
625&namsol        !   elliptic solver / island / free surface
626!-----------------------------------------------------------------------
627   nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
628                           !                   =2 successive-over-relaxation (sor)
629   nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test
630   rn_eps      =  1.e-6    !  absolute precision of the solver
631   nn_nmin     =    210    !  minimum of iterations for the SOR solver
632   nn_nmax     =    800    !  maximum of iterations for the SOR solver
633   nn_nmod     =     10    !  frequency of test for the SOR solver
634   rn_resmax   =  1.e-10   !  absolute precision for the SOR solver
635   rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain)
636/
637!-----------------------------------------------------------------------
638&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
639!-----------------------------------------------------------------------
640   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
641                           !  buffer blocking send or immediate non-blocking sends, resp.
642   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
643/
644!-----------------------------------------------------------------------
645&nammpp_dyndist !   Massively Parallel Distribution for AGRIF zoom      ("key_agrif" && "key_mpp_dyndist")
646!-----------------------------------------------------------------------
647   jpni        =    1      !  jpni   number of processors following i
648   jpnj        =    1      !  jpnj   number of processors following j
649   jpnij       =    1      !  jpnij  number of local domains
650/
651!-----------------------------------------------------------------------
652&namctl        !   Control prints & Benchmark
653!-----------------------------------------------------------------------
654   ln_ctl      = .false.   !  trends control print (expensive!)
655   nn_print    =    0      !  level of print (0 no extra print)
656   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
657   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
658   nn_jctls    =    0      !  start j indice of control               over a subdomain)
659   nn_jctle    =    0      !  end   j indice of control
660   nn_isplt    =    1      !  number of processors in i-direction
661   nn_jsplt    =    1      !  number of processors in j-direction
662   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench
663                           !     (no physical validity of the results)
664/
665
666!!======================================================================
667!!                       ***  Diagnostics namelists  ***
668!!======================================================================
669!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld")
670!!   namgap       level mean model-data gap                             ("key_diagap")
671!!   namflo       float parameters                                      ("key_float")
672!!   namptr       Poleward Transport Diagnostics
673!!======================================================================
674
675!-----------------------------------------------------------------------
676&namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra")
677!              !       or mixed-layer trends or barotropic vorticity    ('key_trdmld' or  "key_trdvor")
678!-----------------------------------------------------------------------
679   nn_trd      = 365       !  time step frequency dynamics and tracers trends
680   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
681   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
682   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
683   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
684   ln_trdmld_restart = .false.         !  restart for ML diagnostics
685   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
686/
687!-----------------------------------------------------------------------
688&namgap       !   level mean model-data gap                             ('key_diagap')
689!-----------------------------------------------------------------------
690   nn_gap     =  15        !  time-step frequency of model-data gap computation
691   nn_prg     =  10        !  time-step frequency of gap print in model output
692/
693!-----------------------------------------------------------------------
694&namflo       !   float parameters                                      ("key_float")
695!-----------------------------------------------------------------------
696    ln_rstflo = .false.    !  float restart (T) or not (F)
697    nn_writefl=      75    !  frequency of writing in float output file
698    nn_stockfl=    5475    !  frequency of creation of the float restart file
699    ln_argo   = .false.    !  Argo type floats (stay at the surface each 10 days)
700    ln_flork4 = .false.    !  trajectories computed with a 4th order Runge-Kutta (T)
701                           !  or computed with Blanke' scheme (F)
702/
703!-----------------------------------------------------------------------
704&namptr       !   Poleward Transport Diagnostic
705!-----------------------------------------------------------------------
706   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F)
707   ln_diaznl  = .false.    !  Add zonal means and meridional stream functions
708   ln_subbas  = .false.    !  Atlantic/Pacific/Indian basins computation (T) or not
709                           !  (orca configuration only, need input basins mask file named "subbasins.nc"
710   nf_ptr     =  1         !  Frequency of ptr computation [time step]
711   nf_ptr_wri =  15        !  Frequency of ptr outputs
712/
713!-----------------------------------------------------------------------
714&namhsb       !  Heat and salt budgets
715!-----------------------------------------------------------------------
716   ln_diahsb  = .false.     !  check the heat and salt budgets (T) or not (F)
717/
718!-----------------------------------------------------------------------
719!       namobs    observation usage switch
720!-----------------------------------------------------------------------
721!
722!  ln_t3d                  Logical switch for T profile observations         
723!  ln_s3d                  Logical switch for S profile observations         
724!  ln_ena                  Logical switch for ENACT insitu data set           
725!  ln_cor                  Logical switch for Coriolis insitu data set       
726!  ln_profb                Logical switch for feedback insitu data set     
727!  ln_sla                  Logical switch for SLA observations               
728!  ln_sladt                Logical switch for AVISO SLA data             
729!  ln_slafb                Logical switch for feedback SLA data           
730!  ln_ssh                  Logical switch for SSH observations             
731!  ln_sst                  Logical switch for SST observations             
732!  ln_reysst               Logical switch for Reynolds observations       
733!  ln_ghrsst               Logical switch for GHRSST observations         
734!  ln_sstfb                Logical switch for feedback SST data         
735!  ln_sss                  Logical switch for SSS observations             
736!  ln_seaice               Logical switch for Sea Ice observations       
737!  ln_vel3d                Logical switch for velocity observations         
738!  ln_velavcur             Logical switch for velocity daily av. cur.   
739!  ln_velhrcur             Logical switch for velocity high freq. cur.   
740!  ln_velavadcp            Logical switch for velocity daily av. ADCP 
741!  ln_velhradcp            Logical switch for velocity high freq. ADCP
742!  ln_velfb                Logical switch for feedback velocity data       
743!  ln_grid_global          Global distribtion of observations         
744!  ln_grid_search_lookup   Logical switch for obs grid search w/lookup table 
745!  grid_search_file        Grid search lookup file header
746!  enactfiles              ENACT input observation file names
747!  coriofiles              Coriolis input observation file name 
748!  profbfiles              Profile feedback input observation file name
749!  ln_profb_enatim         Enact feedback input time setting switch   
750!  slafilesact             Active SLA input observation file name
751!  slafilespas             Passive SLA input observation file name
752!  slafbfiles              Feedback SLA input observation file name
753!  sstfiles                GHRSST input observation file name       
754!  sstfbfiles              Feedback SST input observation file name
755!  seaicefiles             Sea Ice input observation file name
756!  velavcurfiles           Vel. cur. daily av. input file name 
757!  velhvcurfiles           Vel. cur. high freq. input file name 
758!  velavadcpfiles          Vel. ADCP daily av. input file name   
759!  velhvadcpfiles          Vel. ADCP high freq. input file name
760!  velfbfiles              Vel. feedback input observation file name
761!  dobsini                 Initial date in window YYYYMMDD.HHMMSS       
762!  dobsend                 Final date in window YYYYMMDD.HHMMSS         
763!  n1dint                  Type of vertical interpolation method       
764!  n2dint                  Type of horizontal interpolation method       
765!  ln_nea                  Rejection of observations near land switch   
766!  nmsshc                  MSSH correction scheme                         
767!  mdtcorr                 MDT  correction                               
768!  mdtcutoff               MDT cutoff for computed correction         
769!  ln_altbias              Logical switch for alt bias               
770!  ln_ignmis               Logical switch for ignoring missing files   
771!  endailyavtypes          ENACT daily average types                   
772 &namobs
773   ln_t3d = .false.
774   ln_s3d = .false.
775   ln_ena = .false.
776   ln_profb = .false.
777   ln_sla = .false.
778   ln_sladt = .false.
779   ln_slafb = .false.
780   ln_sst = .false.
781   ln_sstfb = .false.
782   nmsshc = 0
783   profbfiles = 'profiles_01.nc'
784   slafbfiles = 'sla_01.nc'
785   sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc'
786   ln_altbias = .false.
787   ln_grid_global = .true.
788   ln_grid_search_lookup = .false.
789   ln_ignmis = .true. 
790/
791!-----------------------------------------------------------------------
792!       nam_asminc    assimilation increments namelist
793!-----------------------------------------------------------------------
794! ln_bkgwri   Logical switch for writing out background state
795! ln_trjwri   Logical switch for writing out state trajectory
796! ln_trainc   Logical switch for applying tracer increments
797! ln_dyninc   Logical switch for applying velocity increments
798! ln_sshinc   Logical switch for applying SSH increments
799! ln_asmdin   Logical switch for Direct Initialization (DI)
800! ln_asmiau   Logical switch for Incremental Analysis Updating (IAU)
801! nitbkg      Timestep of background in [0,nitend-nit000-1]
802! nitdin      Timestep of background for DI in [0,nitend-nit000-1]
803! nitiaustr   Timestep of start of IAU interval in [0,nitend-nit000-1]
804! nitiaufin   Timestep of end of IAU interval in [0,nitend-nit000-1]
805! niaufn      Type of IAU weighting function
806! nittrjfrq   Frequency of trajectory output for 4D-VAR
807! ln_salfix   Logical switch for ensuring that the sa > salfixmin
808! salfixmin   Minimum salinity after applying the increments
809&nam_asminc
810    ln_bkgwri = .false.
811    ln_trjwri = .false.
812    ln_trainc = .false.
813    ln_dyninc = .false.
814    ln_sshinc = .false.
815    ln_asmdin = .false.
816    ln_asmiau = .false.
817    nitbkg = 0
818    nitdin = 0
819    nitiaustr = 1
820    nitiaufin = 15
821    niaufn = 0
822    nittrjfrq = 0
823    ln_salfix = .false.
824    salfixmin = -9999
825/
Note: See TracBrowser for help on using the repository browser.