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.
1_namelist in branches/2012/dev_r3452_UKMO9_RESTART/NEMOGCM/CONFIG/ORCA2_LIM/EXP00 – NEMO

source: branches/2012/dev_r3452_UKMO9_RESTART/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist @ 3594

Last change on this file since 3594 was 3594, checked in by rfurner, 11 years ago

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