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/2013/dev_r3987_METO1_MERCATOR6_OBC_BDY_merge/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00 – NEMO

source: branches/2013/dev_r3987_METO1_MERCATOR6_OBC_BDY_merge/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/namelist @ 9260

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