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_UKMO_NOC_reproducibililty_env_bathy/NEMOGCM/CONFIG/AMM12/EXP00 – NEMO

source: branches/2013/dev_UKMO_NOC_reproducibililty_env_bathy/NEMOGCM/CONFIG/AMM12/EXP00/namelist @ 4179

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

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