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/2012/dev_r3389_INGV4_stokes/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00 – NEMO

source: branches/2012/dev_r3389_INGV4_stokes/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist @ 3529

Last change on this file since 3529 was 3529, checked in by adani, 11 years ago

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