source: NEMO/branches/2020/dev_r12973_AGRIF_CMEMS/cfgs/SHARED/namelist_ref @ 13141

Last change on this file since 13141 was 13141, checked in by jchanut, 3 months ago

#2129, corrections/add ons to initial state interpolation with AGRIF
1) add namelist flag for child grid initial state interpolation - ice not considered yet
2) provide depths and not thicknesses as inputs to vertical linear interpolation
3) extend initial state interpolation to a restart scenario for parent grid (warning should be added in that case in order to prevent users doing this at each model restart…)
The online interpolation seems to work fine in the VORTEX case (provided 0. is not considered as a special value in the initial velocity field, i.e. ln_spc_dyn=F)

  • Property svn:keywords set to Id
  • Property svn:mime-type set to text/x-fortran
File size: 110.6 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OCE :   Reference namelist_ref                                !!
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namtsd, namcrs, namc1d, namc1d_uvd)
5!! namelists    2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl,
6!!                                    namsbc_sas, namtra_qsr, namsbc_rnf,
7!!                                    namisf, namsbc_apr,
8!!                                    namsbc_ssr, namsbc_wave, namberg)
9!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide)
10!!              4 - top/bot boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl)
11!!              5 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_eiv, namtra_dmp)
12!!              6 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
13!!              7 - Vertical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm)
14!!              8 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb)
15!!              9 - Obs & Assim      (namobs, nam_asminc)
16!!             10 - miscellaneous    (nammpp, namctl, namsto)
17!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
18
19!!======================================================================
20!!              ***  Domain & Run management namelists  ***           !!
21!!                                                                    !!
22!!   namrun       parameters of the run
23!!   namdom       space and time domain
24!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
25!!   namwad       Wetting and drying                                    (default: OFF)
26!!   namtsd       data: temperature & salinity                          (default: OFF)
27!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
28!!   namc1d       1D configuration options                              ("key_c1d")
29!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d")
30!!   namc1d_uvd   1D data (currents)                                    ("key_c1d")
31!!======================================================================
32!
33!-----------------------------------------------------------------------
34&namrun        !   parameters of the run
35!-----------------------------------------------------------------------
36   nn_no       =       0   !  Assimilation cycle index
37   cn_exp      =  "ORCA2"  !  experience name
38   nn_it000    =       1   !  first time step
39   nn_itend    =    5840   !  last  time step (std 5840)
40   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
41   nn_time0    =       0   !  initial time of day in hhmm
42   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
43   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
44      ln_1st_euler = .false.  !  =T force a start with forward time step (ln_rstart=T)
45      nn_rstctl    =    0     !  restart control ==> activated only if ln_rstart=T
46      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
47      !                          !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
48      !                          !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart
49      cn_ocerst_in    = "restart"   !  suffix of ocean restart name (input)
50      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts
51      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output)
52      cn_ocerst_outdir = "."        !  directory in which to write output ocean restarts
53   nn_istate   =       0   !  output the initial state (1) or not (0)
54   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
55   nn_stock    =       0   !  used only if ln_rst_list = F: output restart freqeuncy (modulo referenced to 1)
56      !                          !    =  0 force to write restart files only at the end of the run
57      !                          !    = -1 do not do any restart
58   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
59   nn_write    =       0   !  used only if key_iomput is not defined: output frequency (modulo referenced to nn_it000)
60      !                          !    =  0 force to write output files only at the end of the run
61      !                          !    = -1 do not do any output file
62   ln_mskland  = .false.   !  mask land points in NetCDF outputs
63   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard
64   ln_clobber  = .true.    !  clobber (overwrite) an existing file
65   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
66   ln_xios_read = .false.  !  use XIOS to read restart file (only for a single file restart)
67   nn_wxios = 0      !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output
68/
69!-----------------------------------------------------------------------
70&namdom        !   time and space domain
71!-----------------------------------------------------------------------
72   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time
73   !
74   rn_Dt       = 5400.     !  time step for the dynamics and tracer
75   rn_atfp     =    0.1    !  asselin time filter parameter
76   !
77   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs)
78   !
79   ln_meshmask = .true.   !  =T create a mesh file
80/
81!-----------------------------------------------------------------------
82&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
83!-----------------------------------------------------------------------
84   ln_read_cfg = .false.     !  (=T) read the domain configuration file
85      !                      !  (=F) user defined configuration           (F => create/check namusr_def)
86      cn_domcfg = "domain_cfg"  ! domain configuration filename
87      !
88      ln_closea    = .false. !  (=T => fill namclo)
89      !                      !  (=F) no control of net precip/evap over closed sea
90      !
91   ln_write_cfg = .false.    !  (=T) create the domain configuration file
92      cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename
93      !
94   ln_use_jattr = .false.    !  use (T) the file attribute: open_ocean_jstart, if present
95   !                         !  in netcdf input files, as the start j-row for reading
96/
97!-----------------------------------------------------------------------
98&namclo        !   parameters of the closed sea (cs) behavior                (default: OFF)
99!-----------------------------------------------------------------------
100   ln_maskcs = .false.        ! (=T) cs are masked ; So, in this case ln_mask_csundef and ln_clo_rnf have no effect.
101      !                       ! (=F => set ln_mask_csundef and ln_clo_rnf)
102      !                       ! cs masks are read and net evap/precip over closed sea spread out depending on domain_cfg.nc masks.
103      !                       ! See ln_mask_csundef and ln_clo_rnf for specific option related to this case
104      !
105      ln_mask_csundef = .true.   ! (=T) undefined closed seas are masked ;
106      !                          ! (=F) undefined closed seas are kept and no specific treatment is done for these closed seas
107      !
108      ln_clo_rnf = .true.        ! (=T) river mouth specified in domain_cfg.nc masks (rnf and emp case) are added to the runoff mask.
109      !                          !      allow the treatment of closed sea outflow grid-points to be the same as river mouth grid-points
110/
111!-----------------------------------------------------------------------
112&namtsd        !    Temperature & Salinity Data  (init/dmp)             (default: OFF)
113!-----------------------------------------------------------------------
114   !                       ! =T  read T-S fields for:
115   ln_tsd_init = .false.         !  ocean initialisation
116   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp)
117
118   cn_dir      = './'      !  root directory for the T-S data location
119   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
120   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
121   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
122   sn_tem = 'data_1m_potential_temperature_nomask',  -1.     , 'votemper',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
123   sn_sal = 'data_1m_salinity_nomask'             ,  -1.     , 'vosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
124/
125!-----------------------------------------------------------------------
126&namwad        !   Wetting and Drying (WaD)                             (default: OFF)
127!-----------------------------------------------------------------------
128   ln_wd_il    = .false.   !  T/F activation of iterative   limiter
129   ln_wd_dl    = .false.   !  T/F activation of directional limiter
130   ln_wd_dl_bc = .false.   !  T/F Directional limiteer Baroclinic option
131   ln_wd_dl_rmp = .false.  !  T/F Turn on directional limiter ramp
132   rn_wdmin0   =  0.30     !  depth at which WaD starts
133   rn_wdmin1   =  0.2      !  Minimum wet depth on dried cells
134   rn_wdmin2   =  0.0001   !  Tolerance of min wet depth on dried cells
135   rn_wdld     =  2.5      !  Land elevation below which WaD is allowed
136   nn_wdit     =   20      !  Max iterations for WaD limiter
137   rn_wd_sbcdep =  5.0     !  Depth at which to taper sbc fluxes
138   rn_wd_sbcfra =  0.999   !  Fraction of SBC fluxes at taper depth (Must be <1)
139/
140!-----------------------------------------------------------------------
141&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T)
142!-----------------------------------------------------------------------
143   nn_factx    = 3         !  Reduction factor of x-direction
144   nn_facty    = 3         !  Reduction factor of y-direction
145   nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT
146      !                    !  0, coarse grid is binned with preferential treatment of the north fold
147      !                    !  1, coarse grid is binned with centering at the equator
148      !                    !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.
149   ln_msh_crs  = .false.   ! =T create a mesh & mask file
150   nn_crs_kz   = 0         ! 0, MEAN of volume boxes
151      !                    ! 1, MAX of boxes
152      !                    ! 2, MIN of boxes
153   ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F )
154/
155!-----------------------------------------------------------------------
156&namc1d        !   1D configuration options                             ("key_c1d" default: PAPA station)
157!-----------------------------------------------------------------------
158   rn_lat1d    =      50   !  Column latitude
159   rn_lon1d    =    -145   !  Column longitude
160   ln_c1d_locpt =  .true.   ! Localization of 1D config in a grid (T) or independant point (F)
161/
162!-----------------------------------------------------------------------
163&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
164!-----------------------------------------------------------------------
165   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F)
166/
167!-----------------------------------------------------------------------
168&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
169!-----------------------------------------------------------------------
170   !                       !  =T read U-V fields for:
171   ln_uvd_init   = .false.       !  ocean initialisation
172   ln_uvd_dyndmp = .false.       !  U-V restoring
173
174   cn_dir      = './'      !  root directory for the U-V data location
175   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
176   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
177   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
178   sn_ucur     = 'ucurrent'              ,         -1.       ,'u_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Ume'   , ''
179   sn_vcur     = 'vcurrent'              ,         -1.       ,'v_current',   .false.   , .true. , 'monthly' ,  ''              ,  'Vme'   , ''
180/
181
182!!======================================================================
183!!            ***  Surface Boundary Condition namelists  ***          !!
184!!                                                                    !!
185!!   namsbc          surface boundary condition manager                 (default: NO selection)
186!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
187!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
188!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
189!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
190!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
191!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
192!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
193!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
194!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
195!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
196!!   namberg         iceberg floats                                     (ln_icebergs=T)
197!!======================================================================
198!
199!-----------------------------------------------------------------------
200&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
201!-----------------------------------------------------------------------
202   nn_fsbc     = 2         !  frequency of SBC module call
203      !                    !  (control sea-ice & iceberg model call)
204                     ! Type of air-sea fluxes
205   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc)
206   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx )
207   ln_blk      = .false.   !  Bulk formulation                          (T => fill namsbc_blk )
208   ln_abl      = .false.   !  ABL  formulation                          (T => fill namsbc_abl )
209      !              ! Type of coupling (Ocean/Ice/Atmosphere) :
210   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 )
211   ln_mixcpl   = .false.   !  forced-coupled mixed formulation          ( requires key_oasis3 )
212   nn_components = 0       !  configuration of the opa-sas OASIS coupling
213      !                    !  =0 no opa-sas OASIS coupling: default single executable config.
214      !                    !  =1 opa-sas OASIS coupling: multi executable config., OPA component
215      !                    !  =2 opa-sas OASIS coupling: multi executable config., SAS component
216                     ! Sea-ice :
217   nn_ice      = 0         !  =0 no ice boundary condition
218      !                    !  =1 use observed ice-cover                 (  => fill namsbc_iif )
219      !                    !  =2 or 3 for SI3 and CICE, respectively
220   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges)
221      !                    !  =F levitating ice (no pressure, mass and salt exchanges)
222                     ! Misc. options of sbc :
223   ln_traqsr   = .false.   !  Light penetration in the ocean            (T => fill namtra_qsr)
224   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave
225   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
226   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked
227      !                    !     =1 global mean of e-p-r set to zero at each time step
228      !                    !     =2 annual global mean of e-p-r set to zero
229   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf)
230   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr )
231   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave)
232   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)
233   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)
234   nn_sdrift   =  0        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
235      !                    !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
236      !                    !   = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
237      !                    !   = 2 Phillips as (1) but using the wave frequency from a wave model
238   ln_tauwoc   = .false.   !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)
239   ln_tauw     = .false.   !  Activate ocean stress components from wave model
240   ln_stcor    = .false.   !  Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
241   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
242                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
243/
244!-----------------------------------------------------------------------
245&namsbc_flx    !   surface boundary condition : flux formulation        (ln_flx =T)
246!-----------------------------------------------------------------------
247   cn_dir      = './'      !  root directory for the fluxes data location
248   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
249   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
250   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
251   sn_utau     = 'utau'                  ,        24.        , 'utau'    , .false.     , .false., 'yearly'  , ''               , ''       , ''
252   sn_vtau     = 'vtau'                  ,        24.        , 'vtau'    , .false.     , .false., 'yearly'  , ''               , ''       , ''
253   sn_qtot     = 'qtot'                  ,        24.        , 'qtot'    , .false.     , .false., 'yearly'  , ''               , ''       , ''
254   sn_qsr      = 'qsr'                   ,        24.        , 'qsr'     , .false.     , .false., 'yearly'  , ''               , ''       , ''
255   sn_emp      = 'emp'                   ,        24.        , 'emp'     , .false.     , .false., 'yearly'  , ''               , ''       , ''
256/
257!-----------------------------------------------------------------------
258&namsbc_blk    !   namsbc_blk  generic Bulk formula          (ln_blk =T)
259!-----------------------------------------------------------------------
260   !                    !  bulk algorithm :
261   ln_NCAR      = .true.    ! "NCAR"      algorithm   (Large and Yeager 2008)
262   ln_COARE_3p0 = .false.   ! "COARE 3.0" algorithm   (Fairall et al. 2003)
263   ln_COARE_3p6 = .false.   ! "COARE 3.6" algorithm   (Edson et al. 2013)
264   ln_ECMWF     = .false.   ! "ECMWF"     algorithm   (IFS cycle 45r1)
265      !
266      rn_zqt     = 10.      !  Air temperature & humidity reference height (m)
267      rn_zu      = 10.      !  Wind vector reference height (m)
268      ln_Cd_L12  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2012)
269      ln_Cd_L15  = .false.  !  air-ice drags = F(ice conc.) (Lupkes et al. 2015)
270      !                     !  - module of the mean stress" data
271      rn_pfac    = 1.       !  multipl. factor for precipitation (total & snow)
272      rn_efac    = 1.       !  multipl. factor for evaporation (0. or 1.)
273      rn_vfac    = 0.       !  multipl. factor for ocean & ice velocity
274      !                     !  used to calculate the wind stress
275      !                     ! (0. => absolute or 1. => relative winds)
276      ln_skin_cs = .false.  !  use the cool-skin parameterization
277      ln_skin_wl = .false.  !  use the warm-layer parameterization
278      !                     !   ==> only available in ECMWF and COARE algorithms
279      ln_humi_sph = .true.  !  humidity "sn_humi" is specific humidity  [kg/kg]
280      ln_humi_dpt = .false. !  humidity "sn_humi" is dew-point temperature [K]
281      ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%]
282   !
283   cn_dir      = './'      !  root directory for the bulk data location
284   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
285   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename               ! rotation ! land/sea mask !
286   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                                      ! 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_hpgi     = 'NOT USED'                   ,   24.        , 'uhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG'     , ''
294   sn_hpgj     = 'NOT USED'                   ,   24.        , 'vhpg'    ,   .false.   , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG'     , ''
295   sn_prec     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'PRC_MOD1',   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
296   sn_snow     = 'ncar_precip.15JUNE2009_fill',   -1.        , 'SNOW'    ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
297   sn_slp      = 'slp.15JUNE2009_fill'        ,    6.        , 'SLP'     ,   .false.   , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc' , ''       , ''
298/
299!-----------------------------------------------------------------------
300&namsbc_abl    !   Atmospheric Boundary Layer formulation           (ln_abl = T)
301!-----------------------------------------------------------------------
302   cn_dir           = './'      !  root directory for the location of the ABL grid file
303   cn_dom           = 'dom_cfg_abl.nc'
304
305   cn_ablrst_in     = "restart_abl"   !  suffix of abl restart name (input)
306   cn_ablrst_out    = "restart_abl"   !  suffix of abl restart name (output)
307   cn_ablrst_indir  = "."             !  directory to read   input abl restarts
308   cn_ablrst_outdir = "."             !  directory to write output abl restarts
309
310   ln_hpgls_frc   = .false.
311   ln_geos_winds  = .false.
312   nn_dyn_restore = 2         ! restoring option for dynamical ABL variables: = 0 no restoring
313                              !                                               = 1 equatorial restoring
314                              !                                               = 2 global restoring
315   rn_ldyn_min   =  4.5       !  magnitude of the nudging on ABL dynamics at the bottom of the ABL   [hour]
316   rn_ldyn_max   =  1.5       !  magnitude of the nudging on ABL dynamics at the top of the ABL   [hour]
317   rn_ltra_min   =  4.5       !  magnitude of the nudging on ABL tracers  at the bottom of the ABL   [hour]
318   rn_ltra_max   =  1.5       !  magnitude of the nudging on ABL tracers  at the top of the ABL   [hour]
319   nn_amxl       =  0         ! mixing length: = 0 Deardorff 80 length-scale
320                              !                = 1 length-scale based on the distance to the PBL height
321                              !                = 2 Bougeault & Lacarrere 89 length-scale
322   rn_Cm         = 0.0667     ! 0.126 in MesoNH
323   rn_Ct         = 0.1667     ! 0.143 in MesoNH
324   rn_Ce         = 0.4        ! 0.4   in MesoNH
325   rn_Ceps       = 0.7        ! 0.85  in MesoNH
326   rn_Rod        = 0.15       ! c0 in RMCA17 mixing length formulation (not yet implemented)
327   rn_Ric        = 0.139      !  Critical Richardson number (to compute PBL height and diffusivities)
328/
329!-----------------------------------------------------------------------
330&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
331!-----------------------------------------------------------------------
332   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentially sending/receiving data
333   ln_usecplmask = .false. !  use a coupling mask file to merge data received from several models
334   !                       !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
335   nn_cats_cpl   =     5   !  Number of sea ice categories over which coupling is to be carried out (if not 1)
336   !_____________!__________________________!____________!_____________!______________________!________!
337   !             !        description       !  multiple  !    vector   !       vector         ! vector !
338   !             !                          ! categories !  reference  !     orientation      ! grids  !
339!***   send    ***
340   sn_snd_temp   =   'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   ''
341   sn_snd_alb    =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
342   sn_snd_thick  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
343   sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T'
344   sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
345   sn_snd_crtw   =   'none'                 ,    'no'    ,     ''      ,         ''           , 'U,V'
346   sn_snd_ifrac  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
347   sn_snd_wlev   =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
348   sn_snd_cond   =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
349   sn_snd_thick1 =   'ice and snow'         ,    'no'    ,     ''      ,         ''           ,   ''
350   sn_snd_mpnd   =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
351   sn_snd_sstfrz =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
352   sn_snd_ttilyr =   'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''
353!***  receive  ***
354   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
355   sn_rcv_taumod =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
356   sn_rcv_tau    =   'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward' ,  'U,V'
357   sn_rcv_dqnsdt =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
358   sn_rcv_qsr    =   'oce and ice'          ,    'no'    ,     ''      ,         ''           ,   ''
359   sn_rcv_qns    =   'oce and ice'          ,    'no'    ,     ''      ,         ''           ,   ''
360   sn_rcv_emp    =   'conservative'         ,    'no'    ,     ''      ,         ''           ,   ''
361   sn_rcv_rnf    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
362   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
363   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''
364   sn_rcv_hsig   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
365   sn_rcv_iceflx =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
366   sn_rcv_mslp   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
367   sn_rcv_phioc  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
368   sn_rcv_sdrfx  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
369   sn_rcv_sdrfy  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
370   sn_rcv_wper   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
371   sn_rcv_wnum   =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
372   sn_rcv_wfreq  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
373   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
374   sn_rcv_ts_ice =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
375   sn_rcv_isf    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
376   sn_rcv_icb    =   'none'                 ,    'no'    ,     ''      ,         ''           ,   ''
377   sn_rcv_tauwoc =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
378   sn_rcv_tauw   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
379   sn_rcv_wdrag  =   'none'                 ,    'no'    ,     ''      ,         ''          ,   ''
380/
381!-----------------------------------------------------------------------
382&namsbc_sas    !   Stand-Alone Surface module: ocean data               (SAS_SRC  only)
383!-----------------------------------------------------------------------
384   l_sasread   = .true.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
385      ln_3d_uve   = .false.   !  specify whether we are supplying a 3D u,v and e3 field
386      ln_read_frq = .false.   !  specify whether we must read frq or not
387
388   cn_dir      = './'      !  root directory for the ocean data location
389   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
390   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
391   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
392   sn_usp      = 'sas_grid_U'            ,       120.        , 'uos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
393   sn_vsp      = 'sas_grid_V'            ,       120.        , 'vos'     ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
394   sn_tem      = 'sas_grid_T'            ,       120.        , 'sosstsst',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
395   sn_sal      = 'sas_grid_T'            ,       120.        , 'sosaline',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
396   sn_ssh      = 'sas_grid_T'            ,       120.        , 'sossheig',   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
397   sn_e3t      = 'sas_grid_T'            ,       120.        , 'e3t_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
398   sn_frq      = 'sas_grid_T'            ,       120.        , 'frq_m'   ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,    ''
399/
400!-----------------------------------------------------------------------
401&namsbc_iif    !   Ice-IF : use observed ice cover                      (nn_ice = 1)
402!-----------------------------------------------------------------------
403   cn_dir      = './'      !  root directory for the ice cover data location
404   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
405   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
406   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
407   sn_ice      ='ice_cover_clim.nc'      ,        -12.       ,'ice_cover',   .true.    , .true.  , 'yearly'  ,   ''            ,   ''     ,   ''
408/
409!-----------------------------------------------------------------------
410&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T)
411!-----------------------------------------------------------------------
412   !                       !  type of penetration                        (default: NO selection)
413   ln_qsr_rgb  = .false.      !  RGB light penetration (Red-Green-Blue)
414   ln_qsr_2bd  = .false.      !  2BD light penetration (two bands)
415   ln_qsr_bio  = .false.      !  bio-model light penetration
416   !                       !  RGB & 2BD choices:
417   rn_abs      =   0.58       !  RGB & 2BD: fraction absorbed in the very near surface
418   rn_si0      =   0.35       !  RGB & 2BD: shortess depth of extinction
419   nn_chldta   =      0       !  RGB : Chl data (=1) or cst value (=0)
420   rn_si1      =   23.0       !  2BD : longest depth of extinction
421
422   cn_dir      = './'      !  root directory for the chlorophyl data location
423   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
424   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
425   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
426   sn_chl      ='chlorophyll'            ,        -1.        , 'CHLA'    ,   .true.    , .true. , 'yearly'  , ''               , ''       , ''
427/
428!-----------------------------------------------------------------------
429&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T)
430!-----------------------------------------------------------------------
431   nn_sstr     =     0     !  add a retroaction term to the surface heat flux (=1) or not (=0)
432      rn_dqdt     = -40.      !  magnitude of the retroaction on temperature   [W/m2/K]
433   nn_sssr     =     0     !  add a damping term to the surface freshwater flux (=2)
434      !                    !  or to SSS only (=1) or no damping term (=0)
435      rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
436      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2)
437      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
438      nn_sssr_ice =   1       ! control of sea surface restoring under sea-ice
439                              ! 0 = no restoration under ice : * (1-icefrac)
440                              ! 1 = restoration everywhere
441                              ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac
442
443   cn_dir      = './'      !  root directory for the SST/SSS data location
444   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
445   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
446   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
447   sn_sst      = 'sst_data'              ,        24.        ,  'sst'    ,    .false.  , .false., 'yearly'  ,    ''            ,    ''    ,     ''
448   sn_sss      = 'sss_data'              ,        -1.        ,  'sss'    ,    .true.   , .true. , 'yearly'  ,    ''            ,    ''    ,     ''
449/
450!-----------------------------------------------------------------------
451&namsbc_rnf    !   runoffs                                              (ln_rnf =T)
452!-----------------------------------------------------------------------
453   ln_rnf_mouth = .false.   !  specific treatment at rivers mouths
454      rn_hrnf     =  15.e0    !  depth over which enhanced vertical mixing is used    (ln_rnf_mouth=T)
455      rn_avt_rnf  =   1.e-3   !  value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
456   rn_rfact    =   1.e0    !  multiplicative factor for runoff
457   ln_rnf_depth = .false.   !  read in depth information for runoff
458   ln_rnf_tem  = .false.   !  read in temperature information for runoff
459   ln_rnf_sal  = .false.   !  read in salinity information for runoff
460   ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file
461      rn_rnf_max  = 5.735e-4  !  max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
462      rn_dep_max  = 150.      !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
463      nn_rnf_depth_file = 0   !  create (=1) a runoff depth file or not (=0)
464
465   cn_dir      = './'      !  root directory for the runoff data location
466   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
467   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
468   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
469   sn_rnf      = 'runoff_core_monthly'   ,        -1.        , 'sorunoff',   .true.    , .true. , 'yearly'  , ''               , ''       , ''
470   sn_cnf      = 'runoff_core_monthly'   ,         0.        , 'socoefr0',   .false.   , .true. , 'yearly'  , ''               , ''       , ''
471   sn_s_rnf    = 'runoffs'               ,        24.        , 'rosaline',   .true.    , .true. , 'yearly'  , ''               , ''       , ''
472   sn_t_rnf    = 'runoffs'               ,        24.        , 'rotemper',   .true.    , .true. , 'yearly'  , ''               , ''       , ''
473   sn_dep_rnf  = 'runoffs'               ,         0.        , 'rodepth' ,   .false.   , .true. , 'yearly'  , ''               , ''       , ''
474/
475!-----------------------------------------------------------------------
476&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T)
477!-----------------------------------------------------------------------
478   rn_pref     = 101000.   !  reference atmospheric pressure   [N/m2]/
479   ln_ref_apr  = .false.   !  ref. pressure: global mean Patm (T) or a constant (F)
480   ln_apr_obc  = .false.   !  inverse barometer added to OBC ssh data
481
482   cn_dir = './'        !  root directory for the Patm data location
483   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
484   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
485   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
486   sn_apr      = 'patm'                  ,         -1.       ,'somslpre' ,   .true.    , .true. , 'yearly'  ,    ''            ,    ''    ,      ''
487/
488!-----------------------------------------------------------------------
489&namisf       !  Top boundary layer (ISF)                               (default: OFF)
490!-----------------------------------------------------------------------
491   !
492   ! ---------------- ice shelf load -------------------------------
493   !
494   cn_isfload = 'uniform'      ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc)
495      rn_isfload_T = -1.9
496      rn_isfload_S = 34.4
497   !
498   ! ---------------- ice shelf melt formulation -------------------------------
499   !
500   ln_isf = .false.           ! activate ice shelf module
501      ln_isfdebug = .false.      ! add debug print in ISF code (global min/max/sum of specific variable)
502      cn_isfdir   = './'         ! directory for all ice shelf input file
503      !
504      ! ---------------- cavities opened -------------------------------
505      !
506      ln_isfcav_mlt = .false.    ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc)
507         cn_isfcav_mlt = '3eq'   ! ice shelf melting formulation (spe/2eq/3eq/oasis)
508         !                       ! spe = fwfisf is read from a forcing field
509         !                       ! 2eq = ISOMIP  like: 2 equations formulation (Hunter et al., 2006 for a short description)
510         !                       ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description)
511         !                       ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf
512         !              !  cn_isfcav_mlt = 2eq or 3eq cases:
513         cn_gammablk = 'vel'     ! scheme to compute gammat/s (spe,ad15,hj99)
514         !                       ! spe      = constant transfert velocity (rn_gammat0, rn_gammas0)
515         !                       ! vel      = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description)
516         !                       ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description)
517         rn_gammat0  = 1.4e-2    ! gammat coefficient used in spe, vel and vel_stab gamma computation method
518         rn_gammas0  = 4.0e-4    ! gammas coefficient used in spe, vel and vel_stab gamma computation method
519         !
520         rn_htbl     =  30.      ! thickness of the top boundary layer    (Losh et al. 2008)
521         !                       ! 0 => thickness of the tbl = thickness of the first wet cell
522         !
523         !* 'spe' and 'oasis' case
524         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________!
525         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
526         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
527         sn_isfcav_fwf = 'isfmlt_cav',      -12.      , 'fwflisf'  ,  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
528      !
529      ! ---------------- cavities parametrised -------------------------------
530      !
531      ln_isfpar_mlt = .false.   ! ice shelf melting parametrised
532         cn_isfpar_mlt = 'spe'  ! ice shelf melting parametrisation (spe/bg03/oasis)
533         !                      ! spe   = fwfisf is read from a forcing field
534         !                      ! bg03  = melt computed using Beckmann and Goosse parametrisation
535         !                      ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf
536         !
537         !* all cases
538         !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________!
539         !           !  file name  ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
540         !           !             !  (if <0  months)  !   name    !  (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
541         sn_isfpar_zmax = 'isfmlt_par',       0        ,'sozisfmax',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
542         sn_isfpar_zmin = 'isfmlt_par',       0        ,'sozisfmin',  .false.    , .true.  , 'yearly'  ,    ''    ,   ''     ,    ''
543         !* 'spe' and 'oasis' case
544         sn_isfpar_fwf = 'isfmlt_par' ,      -12.      ,'sofwfisf' ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    ''
545         !* 'bg03' case
546         sn_isfpar_Leff = 'isfmlt_par',       0.       ,'Leff'     ,  .false.    , .true.  , 'yearly'   ,    ''    ,   ''     ,    ''
547      !
548      ! ---------------- ice sheet coupling -------------------------------
549      !
550      ln_isfcpl = .false.
551         nn_drown       = 10        ! number of iteration of the extrapolation loop (fill the new wet cells)
552         ln_isfcpl_cons = .false.
553/
554!-----------------------------------------------------------------------
555&namsbc_wave   ! External fields from wave model                        (ln_wave=T)
556!-----------------------------------------------------------------------
557   cn_dir      = './'      !  root directory for the waves data location
558   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
559   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
560   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
561   sn_cdg      =  'sdw_ecwaves_orca2'    ,        6.         , 'drag_coeff' ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
562   sn_usd      =  'sdw_ecwaves_orca2'    ,        6.         , 'u_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
563   sn_vsd      =  'sdw_ecwaves_orca2'    ,        6.         , 'v_sd2d'     ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
564   sn_hsw      =  'sdw_ecwaves_orca2'    ,        6.         , 'hs'         ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
565   sn_wmp      =  'sdw_ecwaves_orca2'    ,        6.         , 'wmp'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
566   sn_wfr      =  'sdw_ecwaves_orca2'    ,        6.         , 'wfr'        ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
567   sn_wnum     =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_num'   ,  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
568   sn_tauwoc   =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
569   sn_tauwx    =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
570   sn_tauwy    =  'sdw_ecwaves_orca2'    ,        6.         , 'wave_stress',  .true.  , .true. , 'yearly'  ,  ''              , ''       , ''
571/
572!-----------------------------------------------------------------------
573&namberg       !   iceberg parameters                                   (default: OFF)
574!-----------------------------------------------------------------------
575   ln_icebergs = .false.      ! activate iceberg floats (force =F with "key_agrif")
576   !
577   !                          ! restart
578   cn_icbrst_in     = "restart_icb" !  suffix of iceberg restart name (input)
579   cn_icbrst_indir  = "./"          !  directory from which to read input ocean restarts
580   cn_icbrst_out    = "restart_icb" !  suffix of ocean restart name (output)
581   cn_icbrst_outdir = "./"          !  directory from which to read output ocean restarts         
582   !
583   !                          ! diagnostics:
584   ln_bergdia        = .true.        ! Calculate budgets
585   nn_verbose_level  = 0             ! Turn on more verbose output if level > 0
586   nn_verbose_write  = 15            ! Timesteps between verbose messages
587   nn_sample_rate    = 1             ! Timesteps between sampling for trajectory storage
588   !
589   !                          ! iceberg setting:
590   !                                 ! Initial mass required for an iceberg of each class
591   rn_initial_mass   = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
592   !                                 ! Proportion of calving mass to apportion to each class
593   rn_distribution   = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
594   !                                 ! Ratio between effective and real iceberg mass (non-dim)
595   !                                 ! i.e. number of icebergs represented at a point
596   rn_mass_scaling   = 2000., 200., 50., 20., 10., 5., 2., 1., 1., 1.
597                                     ! thickness of newly calved bergs (m)
598   rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
599   !
600   rn_rho_bergs            = 850.    ! Density of icebergs
601   rn_LoW_ratio            = 1.5     ! Initial ratio L/W for newly calved icebergs
602   ln_operator_splitting   = .true.  ! Use first order operator splitting for thermodynamics
603   rn_bits_erosion_fraction = 0.     ! Fraction of erosion melt flux to divert to bergy bits
604   rn_sicn_shift           = 0.      ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
605   ln_passive_mode         = .false. ! iceberg - ocean decoupling
606   nn_test_icebergs        =  10     ! Create test icebergs of this class (-1 = no)
607   !                                 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
608   rn_test_box             = 108.0,  116.0, -66.0, -58.0
609   ln_use_calving          = .false. ! Use calving data even when nn_test_icebergs > 0
610   rn_speed_limit          = 0.      ! CFL speed limit for a berg
611
612   cn_dir      = './'      !  root directory for the calving data location
613   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
614   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
615   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
616   sn_icb     =  'calving'              ,         -1.        ,'calvingmask',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
617/
618
619!!======================================================================
620!!               ***  Lateral boundary condition  ***                 !!
621!!                                                                    !!
622!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
623!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
624!!   nam_tide      Tidal forcing                                        (default: OFF)
625!!   nambdy        Unstructured open boundaries                         (default: OFF)
626!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
627!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
628!!======================================================================
629!
630!-----------------------------------------------------------------------
631&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
632!-----------------------------------------------------------------------
633   !                       !  free slip  !   partial slip  !   no slip   ! strong slip
634   rn_shlat    =  -9999.   !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
635   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs.
636/
637!-----------------------------------------------------------------------
638&namagrif      !  AGRIF zoom                                            ("key_agrif")
639!-----------------------------------------------------------------------
640   ln_agrif_2way   = .true.  !  activate two way nesting
641   ln_init_chfrpar = .false. !  initialize child grids from parent
642   ln_spc_dyn      = .true.  !  use 0 as special value for dynamics
643   rn_sponge_tra   = 2880.   !  coefficient for tracer   sponge layer [m2/s]
644   rn_sponge_dyn   = 2880.   !  coefficient for dynamics sponge layer [m2/s]
645   rn_trelax_tra   = 0.01    !  inverse of relaxation time (in steps) for tracers []
646   rn_trelax_dyn   = 0.01    !  inverse of relaxation time (in steps) for dynamics []
647   ln_chk_bathy    = .false. !  =T  check the parent bathymetry
648   ln_bry_south    = .true.  !  =T  south boundary open
649/
650!-----------------------------------------------------------------------
651&nam_tide      !   tide parameters                                      (default: OFF)
652!-----------------------------------------------------------------------
653   ln_tide     = .false.      ! Activate tides
654      nn_tide_var   = 1          !  Variant of tidal parameter set and tide-potential computation
655      !                          !     (1: default; 0: compatibility with previous versions)
656      ln_tide_dia   = .false.    !  Enable tidal diagnostic output
657      ln_tide_pot   = .false.               !  use tidal potential forcing
658         rn_tide_gamma = 0.7                   ! Tidal tilt factor
659         ln_scal_load  = .false.               ! Use scalar approximation for
660            rn_scal_load = 0.094               !     load potential
661         ln_read_load  = .false.               ! Or read load potential from file
662            cn_tide_load = 'tide_LOAD_grid_T.nc'  ! filename for load potential
663            !
664      ln_tide_ramp  = .false.               !  Use linear ramp for tides at startup
665         rn_tide_ramp_dt = 0.               !  ramp duration in days
666      sn_tide_cnames(1) = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg
667/
668!-----------------------------------------------------------------------
669&nambdy        !  unstructured open boundaries                          (default: OFF)
670!-----------------------------------------------------------------------
671   ln_bdy         = .false.   !  Use unstructured open boundaries
672   nb_bdy         = 0         !  number of open boundary sets
673   ln_coords_file = .true.    !  =T : read bdy coordinates from file
674      cn_coords_file = 'coordinates.bdy.nc'  !  bdy coordinates files
675   ln_mask_file   = .false.   !  =T : read mask from file
676      cn_mask_file = ''        !  name of mask file (if ln_mask_file=.TRUE.)
677   cn_dyn2d    = 'none'       !
678   nn_dyn2d_dta   =  0        !  = 0, bdy data are equal to the initial state
679      !                       !  = 1, bdy data are read in 'bdydata   .nc' files
680      !                       !  = 2, use tidal harmonic forcing data from files
681      !                       !  = 3, use external data AND tidal harmonic forcing
682   cn_dyn3d      =  'none'    !
683   nn_dyn3d_dta  =  0         !  = 0, bdy data are equal to the initial state
684   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
685   cn_tra        =  'none'    !
686   nn_tra_dta    =  0         !  = 0, bdy data are equal to the initial state
687   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
688   cn_ice        =  'none'    !
689   nn_ice_dta    =  0         !  = 0, bdy data are equal to the initial state
690   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
691   !
692   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers
693   ln_dyn3d_dmp  =.false.     !  open boundary condition for baroclinic velocities
694   rn_time_dmp   =  1.        !  Damping time scale in days
695   rn_time_dmp_out = 1.       !  Outflow damping time scale
696   nn_rimwidth   = 10         !  width of the relaxation zone
697   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter)
698   nn_volctl     =  1         !  = 0, the total water flux across open boundaries is zero
699/
700!-----------------------------------------------------------------------
701&nambdy_dta    !  open boundaries - external data                       (see nam_bdy)
702!-----------------------------------------------------------------------
703   ln_zinterp  = .false.      !  T if a vertical interpolation is required. Variables gdep[tuv] and e3[tuv] must exist in the file
704   !                          !  automatically defined to T if the number of vertical levels in bdy dta /= jpk
705   ln_full_vel = .false.      !  T if [uv]3d are "full" velocities and not only its baroclinic components
706   !                          !  in this case, baroclinic and barotropic velocities will be recomputed -> [uv]2d not needed
707   !
708   cn_dir      = 'bdydta/'    !  root directory for the BDY data location
709   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
710   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
711   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
712   bn_ssh      = 'amm12_bdyT_u2d'        ,         24.       , 'sossheig',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
713   bn_u2d      = 'amm12_bdyU_u2d'        ,         24.       , 'vobtcrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
714   bn_v2d      = 'amm12_bdyV_u2d'        ,         24.       , 'vobtcrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
715   bn_u3d      = 'amm12_bdyU_u3d'        ,         24.       , 'vozocrtx',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
716   bn_v3d      = 'amm12_bdyV_u3d'        ,         24.       , 'vomecrty',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
717   bn_tem      = 'amm12_bdyT_tra'        ,         24.       , 'votemper',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
718   bn_sal      = 'amm12_bdyT_tra'        ,         24.       , 'vosaline',    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
719!* for si3
720   bn_a_i      = 'amm12_bdyT_ice'        ,         24.       , 'siconc'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
721   bn_h_i      = 'amm12_bdyT_ice'        ,         24.       , 'sithic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
722   bn_h_s      = 'amm12_bdyT_ice'        ,         24.       , 'snthic'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
723   bn_t_i      = 'NOT USED'              ,         24.       , 'sitemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
724   bn_t_s      = 'NOT USED'              ,         24.       , 'sntemp'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
725   bn_tsu      = 'NOT USED'              ,         24.       , 'sittop'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
726   bn_s_i      = 'NOT USED'              ,         24.       , 'sisalt'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
727   ! melt ponds (be careful, bn_aip is the pond concentration (not fraction), so it differs from rn_iceapnd)
728   bn_aip      = 'NOT USED'              ,         24.       , 'siapnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
729   bn_hip      = 'NOT USED'              ,         24.       , 'sihpnd'  ,    .true.   , .false.,  'daily'  ,    ''            ,   ''     ,     ''
730   ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds
731   rn_ice_tem  = 270.         !  arbitrary temperature               of incoming sea ice
732   rn_ice_sal  = 10.          !       --   salinity                            --
733   rn_ice_age  = 30.          !       --   age                                 --
734   rn_ice_apnd = 0.2          !       --   pond fraction = a_ip/a_i            --
735   rn_ice_hpnd = 0.05         !       --   pond depth                          --
736/
737!-----------------------------------------------------------------------
738&nambdy_tide   !  tidal forcing at open boundaries                      (default: OFF)
739!-----------------------------------------------------------------------
740   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files
741   ln_bdytide_2ddta = .false.                   !
742/
743
744!!======================================================================
745!!                ***  Top/Bottom boundary condition  ***             !!
746!!                                                                    !!
747!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
748!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
749!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
750!!   nambbc        bottom temperature boundary condition                (default: OFF)
751!!   nambbl        bottom boundary layer scheme                         (default: OFF)
752!!======================================================================
753!
754!-----------------------------------------------------------------------
755&namdrg        !   top/bottom drag coefficient                          (default: NO selection)
756!-----------------------------------------------------------------------
757   ln_OFF      = .false.   !  free-slip       : Cd = 0                  (F => fill namdrg_bot
758   ln_lin      = .false.   !      linear  drag: Cd = Cd0 Uc0                   &   namdrg_top)
759   ln_non_lin  = .false.   !  non-linear  drag: Cd = Cd0 |U|
760   ln_loglayer = .false.   !  logarithmic drag: Cd = vkarmn/log(z/z0) |U|
761   !
762   ln_drgimp   = .true.    !  implicit top/bottom friction flag
763/
764!-----------------------------------------------------------------------
765&namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T)
766!-----------------------------------------------------------------------
767   rn_Cd0      =  1.e-3    !  drag coefficient [-]
768   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)
769   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag)
770   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases)
771   rn_z0       =  3.0e-3   !  roughness [m] (ln_loglayer=T)
772   ln_boost    = .false.   !  =T regional boost of Cd0 ; =F constant
773      rn_boost =  50.         !  local boost factor  [-]
774/
775!-----------------------------------------------------------------------
776&namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F)
777!-----------------------------------------------------------------------
778   rn_Cd0      =  1.e-3    !  drag coefficient [-]
779   rn_Uc0      =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)
780   rn_Cdmax    =  0.1      !  drag value maximum [-] (logarithmic drag)
781   rn_ke0      =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases)
782   rn_z0       =  3.e-3    !  roughness [m] (ln_loglayer=T)
783   ln_boost    = .false.   !  =T regional boost of Cd0 ; =F constant
784      rn_boost =  50.         !  local boost factor  [-]
785/
786!-----------------------------------------------------------------------
787&nambbc        !   bottom temperature boundary condition                (default: OFF)
788!-----------------------------------------------------------------------
789   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
790      nn_geoflx     = 2       !  geothermal heat flux: = 1 constant flux
791      !                       !                        = 2 read variable flux [mW/m2]
792      rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux       [mW/m2]
793
794   cn_dir      = './'      !  root directory for the geothermal data location
795   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
796   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
797   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
798   sn_qgh      ='geothermal_heating.nc'  ,       -12.        , 'heatflow',   .false.   , .true. , 'yearly'  ,   ''             ,   ''     ,   ''
799/
800!-----------------------------------------------------------------------
801&nambbl        !   bottom boundary layer scheme                         (default: OFF)
802!-----------------------------------------------------------------------
803   ln_trabbl   = .false.   !  Bottom Boundary Layer parameterisation flag
804      nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0)
805      nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0)
806      rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s]
807      rn_gambbl   =  10.      !  advective bbl coefficient                 [s]
808/
809
810!!======================================================================
811!!                        Tracer (T-S) namelists                      !!
812!!                                                                    !!
813!!   nameos        equation of state                                    (default: NO selection)
814!!   namtra_adv    advection scheme                                     (default: NO selection)
815!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
816!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
817!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
818!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
819!!======================================================================
820!
821!-----------------------------------------------------------------------
822&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
823!-----------------------------------------------------------------------
824   ln_teos10   = .false.         !  = Use TEOS-10
825   ln_eos80    = .false.         !  = Use EOS80
826   ln_seos     = .false.         !  = Use S-EOS (simplified Eq.)
827                                 !
828   !                     ! S-EOS coefficients (ln_seos=T):
829   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS
830   rn_a0       =  1.6550e-1      !  thermal expension coefficient
831   rn_b0       =  7.6554e-1      !  saline  expension coefficient
832   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos)
833   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos)
834   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos)
835   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos)
836   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos)
837/
838!-----------------------------------------------------------------------
839&namtra_adv    !   advection scheme for tracer                          (default: NO selection)
840!-----------------------------------------------------------------------
841   ln_traadv_OFF = .false. !  No tracer advection
842   ln_traadv_cen = .false. !  2nd order centered scheme
843      nn_cen_h   =  4            !  =2/4, horizontal 2nd order CEN / 4th order CEN
844      nn_cen_v   =  4            !  =2/4, vertical   2nd order CEN / 4th order COMPACT
845   ln_traadv_fct = .false. !  FCT scheme
846      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order
847      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order
848   ln_traadv_mus = .false. !  MUSCL scheme
849      ln_mus_ups = .false.       !  use upstream scheme near river mouths
850   ln_traadv_ubs = .false. !  UBS scheme
851      nn_ubs_v   =  2            !  =2  , vertical 2nd order FCT / COMPACT 4th order
852   ln_traadv_qck = .false. !  QUICKEST scheme
853/
854!-----------------------------------------------------------------------
855&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection)
856!-----------------------------------------------------------------------
857   !                       !  Operator type:
858   ln_traldf_OFF   = .false.   !  No explicit diffusion
859   ln_traldf_lap   = .false.   !    laplacian operator
860   ln_traldf_blp   = .false.   !  bilaplacian operator
861   !
862   !                       !  Direction of action:
863   ln_traldf_lev   = .false.   !  iso-level
864   ln_traldf_hor   = .false.   !  horizontal  (geopotential)
865   ln_traldf_iso   = .false.   !  iso-neutral (standard operator)
866   ln_traldf_triad = .false.   !  iso-neutral (triad    operator)
867   !
868   !                       !  iso-neutral options:
869   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators)
870   rn_slpmax       =  0.01     !  slope limit                           (both operators)
871   ln_triad_iso    = .false.   !  pure horizontal mixing in ML              (triad only)
872   rn_sw_triad     = 1         !  =1 switching triad ; =0 all 4 triads used (triad only)
873   ln_botmix_triad = .false.   !  lateral mixing on bottom                  (triad only)
874   !
875   !                       !  Coefficients:
876   nn_aht_ijk_t    = 0         !  space/time variation of eddy coefficient:
877      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file
878      !                             !   =  0           constant
879      !                             !   = 10 F(k)      =ldf_c1d
880      !                             !   = 20 F(i,j)    =ldf_c2d
881      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
882      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
883      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
884      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)
885      !                             !                           or   = 1/12 Ud*Ld^3 (blp case)
886      rn_Ud        = 0.01           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
887      rn_Ld        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10)
888/
889!-----------------------------------------------------------------------
890&namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF)
891!-----------------------------------------------------------------------
892   ln_mle      = .false.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation
893   rn_ce       = 0.06      ! magnitude of the MLE (typical value: 0.06 to 0.08)
894   nn_mle      = 1         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
895   rn_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_mle=0)
896   rn_time     = 172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0)
897   rn_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
898   nn_mld_uv   = 0         ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
899   nn_conv     = 0         ! =1 no MLE in case of convection ; =0 always MLE
900   rn_rho_c_mle = 0.01      ! delta rho criterion used to calculate MLD for FK
901/
902!-----------------------------------------------------------------------
903&namtra_eiv    !   eddy induced velocity param.                         (default: OFF)
904!-----------------------------------------------------------------------
905   ln_ldfeiv   = .false.   ! use eddy induced velocity parameterization
906      !
907      !                        !  Coefficients:
908      nn_aei_ijk_t    = 0           !  space/time variation of eddy coefficient:
909      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
910      !                             !   =  0           constant
911      !                             !   = 10 F(k)      =ldf_c1d
912      !                             !   = 20 F(i,j)    =ldf_c2d
913      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
914      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
915      !                        !  time invariant coefficients:  aei0 = 1/2  Ue*Le
916      rn_Ue        = 0.02           !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
917      rn_Le        = 200.e+3        !  lateral diffusive length   [m]   (nn_aht_ijk_t= 0, 10)
918      !
919      ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities
920/
921!-----------------------------------------------------------------------
922&namtra_dmp    !   tracer: T & S newtonian damping                      (default: OFF)
923!-----------------------------------------------------------------------
924   ln_tradmp   =  .false.  !  add a damping term (using resto.nc coef.)
925      nn_zdmp  =    0         !  vertical shape =0    damping throughout the water column
926      !                       !                 =1 no damping in the mixing layer (kz  criteria)
927      !                       !                 =2 no damping in the mixed  layer (rho crieria)
928      cn_resto = 'resto.nc'   !  Name of file containing restoration coeff. field (use dmp_tools to create this)
929/
930
931!!======================================================================
932!!                      ***  Dynamics namelists  ***                  !!
933!!                                                                    !!
934!!   nam_vvl       vertical coordinate options                          (default: z-star)
935!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
936!!   namdyn_vor    advection scheme                                     (default: NO selection)
937!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
938!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
939!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
940!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
941!!======================================================================
942!
943!-----------------------------------------------------------------------
944&nam_vvl       !   vertical coordinate options                          (default: z-star)
945!-----------------------------------------------------------------------
946   ln_vvl_zstar  = .true.           !  z-star vertical coordinate
947   ln_vvl_ztilde = .false.          !  z-tilde vertical coordinate: only high frequency variations
948   ln_vvl_layer  = .false.          !  full layer vertical coordinate
949   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar
950   ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator
951   rn_ahe3       =  0.0             !  thickness diffusion coefficient
952   rn_rst_e3t    = 30.0             !  ztilde to zstar restoration timescale [days]
953   rn_lf_cutoff  =  5.0             !  cutoff frequency for low-pass filter  [days]
954   rn_zdef_max   =  0.9             !  maximum fractional e3t deformation
955   ln_vvl_dbg    = .false.          !  debug prints    (T/F)
956/
957!-----------------------------------------------------------------------
958&namdyn_adv    !   formulation of the momentum advection                (default: NO selection)
959!-----------------------------------------------------------------------
960   ln_dynadv_OFF = .false. !  linear dynamics (no momentum advection)
961   ln_dynadv_vec = .false. !  vector form - 2nd centered scheme
962     nn_dynkeg     = 0        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction
963   ln_dynadv_cen2 = .false. !  flux form - 2nd order centered scheme
964   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
965/
966!-----------------------------------------------------------------------
967&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection)
968!-----------------------------------------------------------------------
969   ln_dynvor_ene = .false. !  energy    conserving scheme
970   ln_dynvor_ens = .false. !  enstrophy conserving scheme
971   ln_dynvor_mix = .false. !  mixed scheme
972   ln_dynvor_enT = .false. !  energy conserving scheme (T-point)
973   ln_dynvor_eeT = .false. !  energy conserving scheme (een using e3t)
974   ln_dynvor_een = .false. !  energy & enstrophy scheme
975      nn_een_e3f = 0          ! =0  e3f = mi(mj(e3t))/4
976      !                       ! =1  e3f = mi(mj(e3t))/mi(mj( tmask))
977   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T)        ==>>> PLEASE DO NOT ACTIVATE
978      !                    !  (f-point vorticity schemes only)
979/
980!-----------------------------------------------------------------------
981&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: NO selection)
982!-----------------------------------------------------------------------
983   ln_hpg_zco  = .false.   !  z-coordinate - full steps
984   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
985   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
986   ln_hpg_isf  = .false.   !  s-coordinate (sco ) adapted to isf
987   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
988   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme)
989/
990!-----------------------------------------------------------------------
991&namdyn_spg    !   surface pressure gradient                            (default: NO selection)
992!-----------------------------------------------------------------------
993   ln_dynspg_exp  = .false.   ! explicit free surface
994   ln_dynspg_ts   = .false.   ! split-explicit free surface
995      ln_bt_fw      = .true.     ! Forward integration of barotropic Eqs.
996      ln_bt_av      = .true.     ! Time filtering of barotropic variables
997         nn_bt_flt     = 1          ! Time filter choice  = 0 None
998         !                          !                     = 1 Boxcar over   nn_e sub-steps
999         !                          !                     = 2 Boxcar over 2*nn_e  "    "
1000      ln_bt_auto    = .true.     ! Number of sub-step defined from:
1001         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed
1002         nn_e         = 30          ! =F : the number of sub-step in rn_Dt seconds
1003      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F)
1004/
1005!-----------------------------------------------------------------------
1006&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO selection)
1007!-----------------------------------------------------------------------
1008   !                       !  Type of the operator :
1009   ln_dynldf_OFF = .false.     !  No operator (i.e. no explicit diffusion)
1010   ln_dynldf_lap = .false.     !    laplacian operator
1011   ln_dynldf_blp = .false.     !  bilaplacian operator
1012   !                       !  Direction of action  :
1013   ln_dynldf_lev = .false.     !  iso-level
1014   ln_dynldf_hor = .false.     !  horizontal  (geopotential)
1015   ln_dynldf_iso = .false.     !  iso-neutral (lap only)
1016   !                       !  Coefficient
1017   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coefficient :
1018      !                             !  =-30  read in eddy_viscosity_3D.nc file
1019      !                             !  =-20  read in eddy_viscosity_2D.nc file
1020      !                             !  =  0  constant
1021      !                             !  = 10  F(k)=c1d
1022      !                             !  = 20  F(i,j)=F(grid spacing)=c2d
1023      !                             !  = 30  F(i,j,k)=c2d*c1d
1024      !                             !  = 31  F(i,j,k)=F(grid spacing and local velocity)
1025      !                             !  = 32  F(i,j,k)=F(local gridscale and deformation rate)
1026      !                        !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case)
1027      !                             !                            or  = 1/12 Uv*Lv^3 (blp case)
1028      rn_Uv      = 0.1              !  lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30)
1029      rn_Lv      = 10.e+3           !  lateral viscous length   [m]   (nn_ahm_ijk_t= 0, 10)
1030      !                       !  Smagorinsky settings  (nn_ahm_ijk_t= 32) :
1031      rn_csmc       = 3.5         !  Smagorinsky constant of proportionality
1032      rn_minfac     = 1.0         !  multiplier of theorectical lower limit
1033      rn_maxfac     = 1.0         !  multiplier of theorectical upper limit
1034      !                       !  iso-neutral laplacian operator (ln_dynldf_iso=T) :
1035      rn_ahm_b      = 0.0         !  background eddy viscosity  [m2/s]
1036/
1037!-----------------------------------------------------------------------
1038&namdta_dyn    !   offline ocean input files                            (OFF_SRC only)
1039!-----------------------------------------------------------------------
1040   ln_dynrnf       =  .false.    !  runoffs option enabled (T) or not (F)
1041   ln_dynrnf_depth =  .false.    !  runoffs is spread in vertical (T) or not (F)
1042!   fwbcorr        = 3.786e-06   !  annual global mean of empmr for ssh correction
1043
1044   cn_dir      = './'      !  root directory for the ocean data location
1045   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
1046   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
1047   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
1048   sn_tem      = 'dyna_grid_T'           ,       120.        , 'votemper'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1049   sn_sal      = 'dyna_grid_T'           ,       120.        , 'vosaline'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1050   sn_mld      = 'dyna_grid_T'           ,       120.        , 'somixhgt'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1051   sn_emp      = 'dyna_grid_T'           ,       120.        , 'sowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1052   sn_fmf      = 'dyna_grid_T'           ,       120.        , 'iowaflup'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1053   sn_ice      = 'dyna_grid_T'           ,       120.        , 'soicecov'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1054   sn_qsr      = 'dyna_grid_T'           ,       120.        , 'soshfldo'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1055   sn_wnd      = 'dyna_grid_T'           ,       120.        , 'sowindsp'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1056   sn_uwd      = 'dyna_grid_U'           ,       120.        , 'uocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1057   sn_vwd      = 'dyna_grid_V'           ,       120.        , 'vocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1058   sn_wwd      = 'dyna_grid_W'           ,       120.        , 'wocetr_eff',  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1059   sn_avt      = 'dyna_grid_W'           ,       120.        , 'voddmavs'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1060   sn_ubl      = 'dyna_grid_U'           ,       120.        , 'sobblcox'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1061   sn_vbl      = 'dyna_grid_V'           ,       120.        , 'sobblcoy'  ,  .true.   , .true. , 'yearly'  , ''               , ''       , ''
1062/
1063
1064!!======================================================================
1065!!                     vertical physics namelists                     !!
1066!!                                                                    !!
1067!!    namzdf        vertical physics manager                            (default: NO selection)
1068!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
1069!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
1070!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
1071!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
1072!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
1073!!======================================================================
1074!
1075!-----------------------------------------------------------------------
1076&namzdf        !   vertical physics manager                             (default: NO selection)
1077!-----------------------------------------------------------------------
1078   !                       ! adaptive-implicit vertical advection
1079   ln_zad_Aimp = .false.      !  Courant number dependent scheme (Shchepetkin 2015)
1080   !
1081   !                       ! type of vertical closure (required)
1082   ln_zdfcst   = .false.      !  constant mixing
1083   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric)
1084   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke)
1085   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls)
1086   ln_zdfosm   = .false.      !  OSMOSIS BL closure                     (T =>   fill namzdf_osm)
1087   !
1088   !                       ! convection
1089   ln_zdfevd   = .false.      !  enhanced vertical diffusion
1090      nn_evdm     =    0         ! apply on tracer (=0) or on tracer and momentum (=1)
1091      rn_evd      =  100.        ! mixing coefficient [m2/s]
1092   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm
1093      nn_npc      =    1         ! frequency of application of npc
1094      nn_npcp     =  365         ! npc control print frequency
1095   !
1096   ln_zdfddm   = .false.   ! double diffusive mixing
1097      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity)
1098      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio
1099   !
1100   !                       ! gravity wave-driven vertical mixing
1101   ln_zdfiwm   = .false.      ! internal wave-induced mixing            (T =>   fill namzdf_iwm)
1102   ln_zdfswm   = .false.      ! surface  wave-induced mixing            (T => ln_wave=ln_sdw=T )
1103   !
1104   !                       ! coefficients
1105   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F)
1106   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F)
1107   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0)
1108   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0)
1109/
1110!-----------------------------------------------------------------------
1111&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T)
1112!-----------------------------------------------------------------------
1113   rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity
1114   rn_alp      =    5.     !  coefficient of the parameterization
1115   nn_ric      =    2      !  coefficient of the parameterization
1116   ln_mldw     =  .false.  !  enhanced mixing in the Ekman layer
1117      rn_ekmfc    =    0.7    !  Factor in the Ekman depth Equation
1118      rn_mldmin   =    1.0    !  minimum allowable mixed-layer depth estimate (m)
1119      rn_mldmax   = 1000.0    !  maximum allowable mixed-layer depth estimate (m)
1120      rn_wtmix    =   10.0    !  vertical eddy viscosity coeff [m2/s] in the mixed-layer
1121      rn_wvmix    =   10.0    !  vertical eddy diffusion coeff [m2/s] in the mixed-layer
1122/
1123!-----------------------------------------------------------------------
1124&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T)
1125!-----------------------------------------------------------------------
1126   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
1127   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
1128   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
1129   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
1130   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
1131   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it)
1132   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
1133   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
1134   !                       !                 = 1 bounded by the local vertical scale factor
1135   !                       !                 = 2 first vertical derivative of mixing length bounded by 1
1136   !                       !                 = 3 as =2 with distinct dissipative an mixing length scale
1137   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F)
1138      nn_mxlice    = 0        ! type of scaling under sea-ice
1139                              !    = 0 no scaling under sea-ice
1140                              !    = 1 scaling with constant sea-ice thickness
1141                              !    = 2  scaling with mean sea-ice thickness ( only with SI3 sea-ice model )
1142                              !    = 3  scaling with maximum sea-ice thickness
1143      rn_mxlice   = 10.       ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1)
1144   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value
1145   ln_drg      = .false.   !  top/bottom friction added as boundary condition of TKE
1146   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002)
1147      rn_lc       =   0.15    !  coef. associated to Langmuir cells
1148   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to NIWs
1149                              !        = 0 none ; = 1 add a tke source below the ML
1150                              !        = 2 add a tke source just at the base of the ML
1151                              !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T)
1152      rn_efr      =   0.05    !  fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
1153      nn_htau     =   1       !  type of exponential decrease of tke penetration below the ML
1154                              !        = 0  constant 10 m length scale
1155                              !        = 1  0.5m at the equator to 30m poleward of 40 degrees
1156      rn_eice     =   4       !  below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4
1157/
1158!-----------------------------------------------------------------------
1159&namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls =T)
1160!-----------------------------------------------------------------------
1161   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2]
1162   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
1163   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
1164   rn_clim_galp  = 0.267   !  galperin limit
1165   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
1166   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
1167   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
1168   rn_hsro       =  0.02   !  Minimum surface roughness
1169   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met>1)
1170   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3)
1171   !                             ! =3 requires ln_wave=T
1172   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
1173   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
1174   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
1175   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
1176/
1177!-----------------------------------------------------------------------
1178&namzdf_osm    !   OSM vertical diffusion                               (ln_zdfosm =T)
1179!-----------------------------------------------------------------------
1180   ln_use_osm_la = .false.      !  Use namelist  rn_osm_la
1181   rn_osm_la     = 0.3         !  Turbulent Langmuir number
1182   rn_osm_dstokes     = 5.     !  Depth scale of Stokes drift (m)
1183   nn_ave = 0                  ! choice of horizontal averaging on avt, avmu, avmv
1184   ln_dia_osm = .true.         ! output OSMOSIS-OBL variables
1185   rn_osm_hbl0 = 10.           ! initial hbl value
1186   ln_kpprimix = .true.        ! Use KPP-style Ri# mixing below BL
1187   rn_riinfty  = 0.7           ! Highest local Ri_g permitting shear instability
1188   rn_difri  =  0.005          ! max Ri# diffusivity at Ri_g = 0 (m^2/s)
1189   ln_convmix  = .true.        ! Use convective instability mixing below BL
1190   rn_difconv = 1.             ! diffusivity when unstable below BL  (m2/s)
1191   nn_osm_wave = 0             ! Method used to calculate Stokes drift
1192      !                        !  = 2: Use ECMWF wave fields
1193      !                        !  = 1: Pierson Moskowitz wave spectrum
1194      !                        !  = 0: Constant La# = 0.3
1195/
1196!-----------------------------------------------------------------------
1197&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T)
1198!-----------------------------------------------------------------------
1199   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
1200   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency
1201   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F)
1202
1203   cn_dir      = './'      !  root directory for the iwm data location                                                                           
1204   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
1205   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
1206   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
1207   sn_mpb      = 'NOT USED'              , -12               , 'mixing_power_bot' , .false.  , .true. , 'yearly' , '' , ''  , ''
1208   sn_mpp      = 'NOT USED'              , -12               , 'mixing_power_pyc' , .false.  , .true. , 'yearly' , '' , ''  , ''
1209   sn_mpc      = 'NOT USED'              , -12               , 'mixing_power_cri' , .false.  , .true. , 'yearly' , '' , ''  , ''
1210   sn_dsb      = 'NOT USED'              , -12               , 'decay_scale_bot'  , .false.  , .true. , 'yearly' , '' , ''  , ''
1211   sn_dsc      = 'NOT USED'              , -12               , 'decay_scale_cri'  , .false.  , .true. , 'yearly' , '' , ''  , ''
1212/
1213!!======================================================================
1214!!                  ***  Diagnostics namelists  ***                   !!
1215!!                                                                    !!
1216!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
1217!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
1218!!   namhsb       Heat and salt budgets                                 (default: OFF)
1219!!   namdiu       Cool skin and warm layer models                       (default: OFF)
1220!!   namdiu       Cool skin and warm layer models                       (default: OFF)
1221!!   namflo       float parameters                                      (default: OFF)
1222!!   nam_diadct   transports through some sections                      (default: OFF)
1223!!   nam_dia25h   25h Mean Output                                       (default: OFF)
1224!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
1225!!======================================================================
1226!
1227!-----------------------------------------------------------------------
1228&namtrd        !   trend diagnostics                                    (default: OFF)
1229!-----------------------------------------------------------------------
1230   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE
1231   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output
1232   ln_dyn_mxl  = .false.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1233   ln_vor_trd  = .false.   ! (T) 2D barotropic vorticity trends (not coded yet)
1234   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends
1235   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends
1236   ln_tra_trd  = .false.   ! (T) 3D tracer trend output
1237   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1238   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step)
1239/
1240!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
1241!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1242!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input)
1243!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output)
1244!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics
1245!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S
1246!!gm
1247/
1248!-----------------------------------------------------------------------
1249&namhsb        !  Heat and salt budgets                                 (default: OFF)
1250!-----------------------------------------------------------------------
1251   ln_diahsb   = .false.   !  output the heat and salt budgets (T) or not (F)
1252/
1253!-----------------------------------------------------------------------
1254&namdiu        !   Cool skin and warm layer models                      (default: OFF)
1255!-----------------------------------------------------------------------
1256   ln_diurnal      = .false.   !
1257   ln_diurnal_only = .false.   !
1258/
1259!-----------------------------------------------------------------------
1260&namflo        !   float parameters                                     (default: OFF)
1261!-----------------------------------------------------------------------
1262   ln_floats   = .false.      ! activate floats or not
1263      jpnfl       = 1         !    total number of floats during the run
1264      jpnnewflo   = 0         !    number of floats for the restart
1265      ln_rstflo   = .false.   !    float restart (T) or not (F)
1266      nn_writefl  =      75   !    frequency of writing in float output file
1267      nn_stockfl  =    5475   !    frequency of creation of the float restart file
1268      ln_argo     = .false.   !    Argo type floats (stay at the surface each 10 days)
1269      ln_flork4   = .false.   !    trajectories computed with a 4th order Runge-Kutta (T)
1270      !                       !    or computed with Blanke' scheme (F)
1271      ln_ariane   = .true.    !    Input with Ariane tool convention(T)
1272      ln_flo_ascii= .true.    !    Output with Ariane tool netcdf convention(F) or ascii file (T)
1273/
1274!-----------------------------------------------------------------------
1275&nam_diadct    !   transports through some sections                     (default: OFF)
1276!-----------------------------------------------------------------------
1277    ln_diadct  = .false.   ! Calculate transport thru sections or not
1278       nn_dct     = 15     !  time step frequency for transports computing
1279       nn_dctwri  = 15     !  time step frequency for transports writing
1280       nn_secdebug = 112   !      0 : no section to debug
1281       !                   !     -1 : debug all section
1282       !                   !  0 < n : debug section number n
1283/
1284!-----------------------------------------------------------------------
1285&nam_dia25h    !  25h Mean Output                                       (default: OFF)
1286!-----------------------------------------------------------------------
1287   ln_dia25h   = .false.   ! Choose 25h mean output or not
1288/
1289!-----------------------------------------------------------------------
1290&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
1291!-----------------------------------------------------------------------
1292   nn_nchunks_i =   4       !  number of chunks in i-dimension
1293   nn_nchunks_j =   4       !  number of chunks in j-dimension
1294   nn_nchunks_k =   31      !  number of chunks in k-dimension
1295   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1296   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs
1297   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
1298   !                       !  (F) ignore chunking information and produce netcdf3-compatible files
1299/
1300
1301!!======================================================================
1302!!               ***  Observation & Assimilation  ***                 !!
1303!!                                                                    !!
1304!!   namobs       observation and model comparison                      (default: OFF)
1305!!   nam_asminc   assimilation increments                               ('key_asminc')
1306!!======================================================================
1307!
1308!-----------------------------------------------------------------------
1309&namobs        !  observation usage switch                              (default: OFF)
1310!-----------------------------------------------------------------------
1311   ln_diaobs   = .false.             ! Logical switch for the observation operator
1312   !
1313   ln_t3d      = .false.             ! Logical switch for T profile observations
1314   ln_s3d      = .false.             ! Logical switch for S profile observations
1315   ln_sla      = .false.             ! Logical switch for SLA observations
1316   ln_sst      = .false.             ! Logical switch for SST observations
1317   ln_sss      = .false.             ! Logical swithc for SSS observations
1318   ln_sic      = .false.             ! Logical switch for Sea Ice observations
1319   ln_vel3d    = .false.             ! Logical switch for velocity observations
1320   ln_altbias  = .false.             ! Logical switch for altimeter bias correction
1321   ln_sstbias  = .false.             ! Logical switch for SST bias correction
1322   ln_nea      = .false.             ! Logical switch for rejection of observations near land
1323   ln_grid_global = .true.           ! Logical switch for global distribution of observations
1324   ln_grid_search_lookup = .false.   ! Logical switch for obs grid search w/lookup table
1325   ln_ignmis   = .true.              ! Logical switch for ignoring missing files
1326   ln_s_at_t   = .false.             ! Logical switch for computing model S at T obs if not there
1327   ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs
1328   ln_bound_reject  = .false.        ! Logical to remove obs near boundaries in LAMs.
1329   ln_sla_fp_indegs = .true.         ! Logical for SLA: T=> averaging footprint is in degrees, F=> in metres
1330   ln_sst_fp_indegs = .true.         ! Logical for SST: T=> averaging footprint is in degrees, F=> in metres
1331   ln_sss_fp_indegs = .true.         ! Logical for SSS: T=> averaging footprint is in degrees, F=> in metres
1332   ln_sic_fp_indegs = .true.         ! Logical for SIC: T=> averaging footprint is in degrees, F=> in metres
1333! All of the *files* variables below are arrays. Use namelist_cfg to add more files
1334   cn_profbfiles = 'profiles_01.nc'  ! Profile feedback input observation file names
1335   cn_slafbfiles = 'sla_01.nc'       ! SLA feedback input observation file names
1336   cn_sstfbfiles = 'sst_01.nc'       ! SST feedback input observation file names
1337   cn_sssfbfiles = 'sss_01.nc'       ! SSS feedback input observation file names
1338   cn_sicfbfiles = 'sic_01.nc'       ! SIC feedback input observation file names
1339   cn_velfbfiles = 'vel_01.nc'       ! Velocity feedback input observation file names
1340   cn_altbiasfile = 'altbias.nc'     ! Altimeter bias input file name
1341   cn_sstbiasfiles = 'sstbias.nc'    ! SST bias input file name
1342   cn_gridsearchfile ='gridsearch.nc' ! Grid search file name
1343   rn_gridsearchres = 0.5            ! Grid search resolution
1344   rn_mdtcorr  = 1.61                ! MDT  correction
1345   rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction
1346   rn_dobsini  = 00010101.000000     ! Initial date in window YYYYMMDD.HHMMSS
1347   rn_dobsend  = 00010102.000000     ! Final date in window YYYYMMDD.HHMMSS
1348   rn_sla_avglamscl = 0.             ! E/W diameter of SLA observation footprint (metres/degrees)
1349   rn_sla_avgphiscl = 0.             ! N/S diameter of SLA observation footprint (metres/degrees)
1350   rn_sst_avglamscl = 0.             ! E/W diameter of SST observation footprint (metres/degrees)
1351   rn_sst_avgphiscl = 0.             ! N/S diameter of SST observation footprint (metres/degrees)
1352   rn_sss_avglamscl = 0.             ! E/W diameter of SSS observation footprint (metres/degrees)
1353   rn_sss_avgphiscl = 0.             ! N/S diameter of SSS observation footprint (metres/degrees)
1354   rn_sic_avglamscl = 0.             ! E/W diameter of SIC observation footprint (metres/degrees)
1355   rn_sic_avgphiscl = 0.             ! N/S diameter of SIC observation footprint (metres/degrees)
1356   nn_1dint = 0                      ! Type of vertical interpolation method
1357   nn_2dint = 0                      ! Default horizontal interpolation method
1358   nn_2dint_sla = 0                  ! Horizontal interpolation method for SLA
1359   nn_2dint_sst = 0                  ! Horizontal interpolation method for SST
1360   nn_2dint_sss = 0                  ! Horizontal interpolation method for SSS
1361   nn_2dint_sic = 0                  ! Horizontal interpolation method for SIC
1362   nn_msshc     = 0                  ! MSSH correction scheme
1363   nn_profdavtypes = -1              ! Profile daily average types - array
1364/
1365!-----------------------------------------------------------------------
1366&nam_asminc    !   assimilation increments                              ('key_asminc')
1367!-----------------------------------------------------------------------
1368    ln_bkgwri  = .false.   !  Logical switch for writing out background state
1369    ln_trainc  = .false.   !  Logical switch for applying tracer increments
1370    ln_dyninc  = .false.   !  Logical switch for applying velocity increments
1371    ln_sshinc  = .false.   !  Logical switch for applying SSH increments
1372    ln_asmdin  = .false.   !  Logical switch for Direct Initialization (DI)
1373    ln_asmiau  = .false.   !  Logical switch for Incremental Analysis Updating (IAU)
1374    nitbkg     = 0         !  Timestep of background in [0,nitend-nit000-1]
1375    nitdin     = 0         !  Timestep of background for DI in [0,nitend-nit000-1]
1376    nitiaustr  = 1         !  Timestep of start of IAU interval in [0,nitend-nit000-1]
1377    nitiaufin  = 15        !  Timestep of end of IAU interval in [0,nitend-nit000-1]
1378    niaufn     = 0         !  Type of IAU weighting function
1379    ln_salfix  = .false.   !  Logical switch for ensuring that the sa > salfixmin
1380    salfixmin  = -9999     !  Minimum salinity after applying the increments
1381    nn_divdmp  = 0         !  Number of iterations of divergence damping operator
1382/
1383
1384!!======================================================================
1385!!                  ***  Miscellaneous namelists  ***                 !!
1386!!                                                                    !!
1387!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
1388!!   namctl            Control prints                                   (default: OFF)
1389!!   namsto            Stochastic parametrization of EOS                (default: OFF)
1390!!======================================================================
1391!
1392!-----------------------------------------------------------------------
1393&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
1394!-----------------------------------------------------------------------
1395   ln_listonly =  .false.  !  do nothing else than listing the best domain decompositions (with land domains suppression)
1396   !                       !  if T: the largest number of cores tested is defined by max(mppsize, jpni*jpnj)
1397   ln_nnogather =  .true.  !  activate code to avoid mpi_allgather use at the northfold
1398   jpni        =   0       !  number of processors following i (set automatically if < 1), see also ln_listonly = T
1399   jpnj        =   0       !  number of processors following j (set automatically if < 1), see also ln_listonly = T
1400/
1401!-----------------------------------------------------------------------
1402&namctl        !   Control prints                                       (default: OFF)
1403!-----------------------------------------------------------------------
1404   sn_cfctl%l_runstat = .TRUE.    ! switches and which areas produce reports with the proc integer settings.
1405   sn_cfctl%l_trcstat = .FALSE.   ! The default settings for the proc integers should ensure
1406   sn_cfctl%l_oceout  = .FALSE.   ! that  all areas report.
1407   sn_cfctl%l_layout  = .FALSE.   !
1408   sn_cfctl%l_prtctl  = .FALSE.   !
1409   sn_cfctl%l_prttrc  = .FALSE.   !
1410   sn_cfctl%l_oasout  = .FALSE.   !
1411   sn_cfctl%procmin   = 0         ! Minimum area number for reporting [default:0]
1412   sn_cfctl%procmax   = 1000000   ! Maximum area number for reporting [default:1000000]
1413   sn_cfctl%procincr  = 1         ! Increment for optional subsetting of areas [default:1]
1414   sn_cfctl%ptimincr  = 1         ! Timestep increment for writing time step progress info
1415   nn_print    =    0             !  level of print (0 no extra print)
1416   nn_ictls    =    0             !  start i indice of control sum (use to compare mono versus
1417   nn_ictle    =    0             !  end   i indice of control sum        multi processor runs
1418   nn_jctls    =    0             !  start j indice of control               over a subdomain)
1419   nn_jctle    =    0             !  end   j indice of control
1420   nn_isplt    =    1             !  number of processors in i-direction
1421   nn_jsplt    =    1             !  number of processors in j-direction
1422   ln_timing   = .false.          !  timing by routine write out in timing.output file
1423   ln_diacfl   = .false.          !  CFL diagnostics write out in cfl_diagnostics.ascii
1424/
1425!-----------------------------------------------------------------------
1426&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
1427!-----------------------------------------------------------------------
1428   ln_sto_eos  = .false.   ! stochastic equation of state
1429   nn_sto_eos  = 1         ! number of independent random walks
1430   rn_eos_stdxy = 1.4       ! random walk horz. standard deviation (in grid points)
1431   rn_eos_stdz = 0.7       ! random walk vert. standard deviation (in grid points)
1432   rn_eos_tcor = 1440.     ! random walk time correlation (in timesteps)
1433   nn_eos_ord  = 1         ! order of autoregressive processes
1434   nn_eos_flt  = 0         ! passes of Laplacian filter
1435   rn_eos_lim  = 2.0       ! limitation factor (default = 3.0)
1436   ln_rststo   = .false.   ! start from mean parameter (F) or from restart file (T)
1437   ln_rstseed  = .true.    ! read seed of RNG from restart file
1438   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
1439   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)
1440/
Note: See TracBrowser for help on using the repository browser.