source: CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA1_ref @ 2265

Last change on this file since 2265 was 2265, checked in by omamce, 10 years ago

O.M.

Add namelists for ORCA1 LIM2. No diurnal cycle presently.
Add namelists for ORCA1 LIM3 : still not working.

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