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/ORCA2_LIM/EXP00 – NEMO

source: branches/nemo_v3_3_beta/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist @ 2340

Last change on this file since 2340 was 2340, checked in by rblod, 13 years ago

Add immediate send by default and key_nproci=1 key_nprocj=1 to avoid to add it each time, remove obsolete keys

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