New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
namelist_ref in NEMO/branches/2019/ENHANCE-02_ISF_domcfg – NEMO

source: NEMO/branches/2019/ENHANCE-02_ISF_domcfg/namelist_ref @ 11609

Last change on this file since 11609 was 11308, checked in by mathiot, 5 years ago

ENHANCE-02_ISF_domcfg: create domisf to isolate isf related work and add new isf geometry computation used by UKESM (ticket #2142)

File size: 19.6 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OCE :   Reference namelist_ref                                !!
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namzgr, namzgr_sco )
5!!              8 - diagnostics      (namnc4)
6!!             10 - miscellaneous    (nammpp, namctl)
7!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
8!-----------------------------------------------------------------------
9&namrun        !   parameters of the run
10!-----------------------------------------------------------------------
11   nn_no       =       0   !  Assimilation cycle index
12   cn_exp      =  "ORCA2"  !  experience name
13   nn_it000    =       1   !  first time step
14   nn_itend    =    5840   !  last  time step (std 5840)
15   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
16   nn_time0    =       0   !  initial time of day in hhmm
17   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
18   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
19      nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T
20      nn_rstctl   =    0      !  restart control ==> activated only if ln_rstart=T
21      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
22      !                          !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
23      !                          !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart
24      cn_ocerst_in    = "restart"   !  suffix of ocean restart name (input)
25      cn_ocerst_indir = "."         !  directory from which to read input ocean restarts
26      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output)
27      cn_ocerst_outdir = "."         !  directory in which to write output ocean restarts
28   ln_iscpl    = .false.   !  cavity evolution forcing or coupling to ice sheet model
29   nn_istate   =       0   !  output the initial state (1) or not (0)
30   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
31   nn_stock    =    5840   !  frequency of creation of a restart file (modulo referenced to 1)
32   nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
33   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000)
34   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%)
35   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard
36   ln_clobber  = .true.    !  clobber (overwrite) an existing file
37   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
38/
39!-----------------------------------------------------------------------
40&namdom        !   space and time domain (bathymetry, mesh, timestep)
41!-----------------------------------------------------------------------
42   nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file
43   rn_bathy    =    0.     !  value of the bathymetry. if (=0) bottom flat at jpkm1
44   nn_msh      =    0      !  create (=1) a mesh file or not (=0)
45   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0)
46   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of
47   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1
48                           !
49   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0)
50   rn_atfp     =    0.1    !  asselin time filter parameter
51   ln_crs      = .false.      !  Logical switch for coarsening module
52   jphgr_msh   =       0               !  type of horizontal mesh
53                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc
54                                       !  = 1 geographical mesh on the sphere with regular grid-spacing
55                                       !  = 2 f-plane with regular grid-spacing
56                                       !  = 3 beta-plane with regular grid-spacing
57                                       !  = 4 Mercator grid with T/U point at the equator
58   ppglam0     =       0.0             !  longitude of first raw and column T-point (jphgr_msh = 1)
59   ppgphi0     =     -35.0             ! latitude  of first raw and column T-point (jphgr_msh = 1)
60   ppe1_deg    =       1.0             !  zonal      grid-spacing (degrees)
61   ppe2_deg    =       0.5             !  meridional grid-spacing (degrees)
62   ppe1_m      =    5000.0             !  zonal      grid-spacing (degrees)
63   ppe2_m      =    5000.0             !  meridional grid-spacing (degrees)
64   ppsur       =    -4762.96143546300  !  ORCA r4, r2 and r05 coefficients
65   ppa0        =      255.58049070440  ! (default coefficients)
66   ppa1        =      245.58132232490  !
67   ppkth       =       21.43336197938  !
68   ppacr       =        3.0            !
69   ppdzmin     =       10.             !  Minimum vertical spacing
70   pphmax      =     5000.             !  Maximum depth
71   ldbletanh   =    .TRUE.             !  Use/do not use double tanf function for vertical coordinates
72   ppa2        =      100.760928500000 !  Double tanh function parameters
73   ppkth2      =       48.029893720000 !
74   ppacr2      =       13.000000000000 !
75/
76!-----------------------------------------------------------------------
77&namcfg        !   parameters of the configuration
78!-----------------------------------------------------------------------
79   !
80   ln_e3_dep   = .true.    ! =T : e3=dk[depth] in discret sens.
81   !                       !      ===>>> will become the only possibility in v4.0
82   !                       ! =F : e3 analytical derivative of depth function
83   !                       !      only there for backward compatibility test with v3.6
84   !
85   cp_cfg      = "default" !  name of the configuration
86   cp_cfz      = "no zoom" !  name of the zoom of configuration
87   jp_cfg      =      0    !  resolution of the configuration
88   jpkdta      =     31    !  number of levels      ( >= jpk )
89   jpiglo      =     10    !  1st dimension of global domain --> i =jpidta
90   jpjglo      =     12    !  2nd    -                  -    --> j =jpjdta
91   jperio      =      0    !  lateral cond. type (between 0 and 6)
92                                 !  = 0 closed                 ;   = 1 cyclic East-West
93                                 !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot
94                                 !  = 4 cyclic East-West AND North fold T-point pivot
95                                 !  = 5 North fold F-point pivot
96                                 !  = 6 cyclic East-West AND North fold F-point pivot
97   ln_use_jattr = .false.  !  use (T) the file attribute: open_ocean_jstart, if present
98                           !  in netcdf input files, as the start j-row for reading
99   ln_domclo = .false.     ! computation of closed sea masks (see namclo)
100/
101!-----------------------------------------------------------------------
102&namzgr        !   vertical coordinate                                  (default: NO selection)
103!-----------------------------------------------------------------------
104   ln_zco      = .false.   !  z-coordinate - full    steps
105   ln_zps      = .false.   !  z-coordinate - partial steps
106   ln_sco      = .false.   !  s- or hybrid z-s-coordinate
107   ln_isfcav   = .false.   !  ice shelf cavity             (T: see namzgr_isf)
108   ln_linssh   = .false.   !  linear free surface
109/
110!-----------------------------------------------------------------------
111&namzgr_isf    !   isf cavity geometry definition                       (default: OFF)
112!-----------------------------------------------------------------------
113   rn_isfdep_min    = 10.         ! minimum isf draft tickness (if lower, isf draft set to this value)
114   rn_glhw_min      = 1.e-3       ! minimum water column thickness to define the grounding line
115   rn_isfhw_min     = 10          ! minimum water column thickness in the cavity once the grounding line defined.
116   ln_isfchannel    = .false.     ! remove channel (based on 2d mask build from isfdraft-bathy)
117   ln_isfconnect    = .false.     ! force connection under the ice shelf (based on 2d mask build from isfdraft-bathy)
118      nn_kisfmax       = 999         ! limiter in level on the previous condition. (if change larger than this number, get back to value before we enforce the connection)
119      rn_zisfmax       = 7000.       ! limiter in m     on the previous condition. (if change larger than this number, get back to value before we enforce the connection)
120   ln_isfcheminey   = .false.     ! close cheminey
121   ln_isfsubgl      = .false.     ! remove subglacial lake created by the remapping process
122      rn_isfsubgllon   =    0.0      !  longitude of the seed to determine the open ocean
123      rn_isfsubgllat   =    0.0      !  latitude  of the seed to determine the open ocean
124/
125!-----------------------------------------------------------------------
126&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate                (default: OFF)
127!-----------------------------------------------------------------------
128   ln_s_sh94   = .false.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)|
129   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied
130   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch
131                           !  stretching coefficients for all functions
132   rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m)
133   rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m)
134   rn_hc       =  150.0    !  critical depth for transition to stretched coordinates
135                        !!!!!!!  Envelop bathymetry
136   rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1)
137                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.)
138   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20)
139   rn_bb       =    0.8    !  stretching with SH94 s-sigma
140                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.)
141   rn_alpha    =    4.4    !  stretching with SF12 s-sigma
142   rn_efold    =    0.0    !  efold length scale for transition to stretched coord
143   rn_zs       =    1.0    !  depth of surface grid box
144                           !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b
145   rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb
146   rn_zb_b     =   -0.2    !  offset for calculating Zb
147                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above]
148   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)
149/
150!-----------------------------------------------------------------------
151&namclo ! (closed sea : need ln_domclo = .true. in namcfg)              (default: OFF)
152!-----------------------------------------------------------------------
153   rn_lon_opnsea = -2.0     ! longitude seed of open ocean
154   rn_lat_opnsea = -2.0     ! latitude  seed of open ocean
155   nn_closea = 8           ! number of closed seas ( = 0; only the open_sea mask will be computed)
156   !                name   ! lon_src ! lat_src ! lon_trg ! lat_trg ! river mouth area   ! net evap/precip correction scheme ! radius tgt   ! id trg
157   !                       ! (degree)! (degree)! (degree)! (degree)! local/coast/global ! (glo/rnf/emp)                     !     (m)      !
158   ! North American lakes
159   sn_lake(1) = 'superior' ,  -86.57 ,  47.30  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2   
160   sn_lake(2) = 'michigan' ,  -87.06 ,  42.74  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2   
161   sn_lake(3) = 'huron'    ,  -82.51 ,  44.74  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2   
162   sn_lake(4) = 'erie'     ,  -81.13 ,  42.25  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2   
163   sn_lake(5) = 'ontario'  ,  -77.72 ,  43.62  , -66.49  , 50.45   , 'local'            , 'rnf'                             ,   550000.0 , 2   
164   ! African Lake
165   sn_lake(6) = 'victoria' ,   32.93 ,  -1.08  ,  30.44  , 31.37   , 'coast'            , 'emp'                             ,   100000.0 , 3   
166   ! Asian Lakes
167   sn_lake(7) = 'caspian'  ,   50.0  ,  44.0   ,   0.0   ,  0.0    , 'global'           , 'glo'                             ,        0.0 , 1     
168   sn_lake(8) = 'aral'     ,   60.0  ,  45.0   ,   0.0   ,  0.0    , 'global'           , 'glo'                             ,        0.0 , 1   
169/
170!-----------------------------------------------------------------------
171&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
172!-----------------------------------------------------------------------
173   !                       !  free slip  !   partial slip  !   no slip   ! strong slip
174   rn_shlat    =  0        !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
175   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs.
176/
177!-----------------------------------------------------------------------
178&namagrif      !  AGRIF zoom                                            ("key_agrif")
179!-----------------------------------------------------------------------
180   ln_spc_dyn    = .true.  !  use 0 as special value for dynamics
181   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s]
182   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s]
183   ln_chk_bathy  = .false. !  =T  check the parent bathymetry
184/
185!-----------------------------------------------------------------------
186&nambdy        !  unstructured open boundaries                          (default: OFF)
187!-----------------------------------------------------------------------
188   ln_bdy         = .false.   !  Use unstructured open boundaries
189   nb_bdy         = 0         !  number of open boundary sets
190   ln_coords_file = .true.    !  =T : read bdy coordinates from file
191      cn_coords_file = 'coordinates.bdy.nc'  !  bdy coordinates files
192   ln_mask_file   = .false.   !  =T : read mask from file
193      cn_mask_file = ''        !  name of mask file (if ln_mask_file=.TRUE.)
194   cn_dyn2d    = 'none'       !
195   nn_dyn2d_dta   =  0        !  = 0, bdy data are equal to the initial state
196      !                       !  = 1, bdy data are read in 'bdydata   .nc' files
197      !                       !  = 2, use tidal harmonic forcing data from files
198      !                       !  = 3, use external data AND tidal harmonic forcing
199   cn_dyn3d      =  'none'    !
200   nn_dyn3d_dta  =  0         !  = 0, bdy data are equal to the initial state
201   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
202   cn_tra        =  'none'    !
203   nn_tra_dta    =  0         !  = 0, bdy data are equal to the initial state
204   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
205   cn_ice        =  'none'    !
206   nn_ice_dta    =  0         !  = 0, bdy data are equal to the initial state
207   !                          !  = 1, bdy data are read in 'bdydata   .nc' files
208   rn_ice_tem    = 270.       !  si3 only: arbitrary temperature of incoming sea ice
209   rn_ice_sal    = 10.        !  si3 only:      --   salinity           --
210   rn_ice_age    = 30.        !  si3 only:      --   age                --
211   !
212   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers
213   ln_dyn3d_dmp  =.false.     !  open boundary condition for baroclinic velocities
214   rn_time_dmp   =  1.        !  Damping time scale in days
215   rn_time_dmp_out = 1.        !  Outflow damping time scale
216   nn_rimwidth   = 10         !  width of the relaxation zone
217   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter)
218   nn_volctl     =  1         !  = 0, the total water flux across open boundaries is zero
219   nb_jpk_bdy    = -1         ! number of levels in the bdy data (set < 0 if consistent with planned run)
220/
221!-----------------------------------------------------------------------
222&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
223!-----------------------------------------------------------------------
224   nn_nchunks_i =   4       !  number of chunks in i-dimension
225   nn_nchunks_j =   4       !  number of chunks in j-dimension
226   nn_nchunks_k =   31      !  number of chunks in k-dimension
227   !                       !  setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
228   !                       !  is optimal for postprocessing which works exclusively with horizontal slabs
229   ln_nc4zip   = .true.    !  (T) use netcdf4 chunking and compression
230   !                       !  (F) ignore chunking information and produce netcdf3-compatible files
231/
232!-----------------------------------------------------------------------
233&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
234!-----------------------------------------------------------------------
235   cn_mpi_send =  'I'      !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
236   !                       !  buffer blocking send or immediate non-blocking sends, resp.
237   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation
238   ln_nnogather =  .true.  !  activate code to avoid mpi_allgather use at the northfold
239   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1)
240   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1)
241/
242!-----------------------------------------------------------------------
243&namctl        !   Control prints                                       (default: OFF)
244!-----------------------------------------------------------------------
245   ln_ctl = .TRUE.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T
246     sn_cfctl%l_config = .TRUE.     ! IF .true. then control which reports are written with the following
247       sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings.
248       sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure
249       sn_cfctl%l_oceout  = .FALSE. ! that  all areas report.
250       sn_cfctl%l_layout  = .FALSE. !
251       sn_cfctl%l_mppout  = .FALSE. !
252       sn_cfctl%l_mpptop  = .FALSE. !
253       sn_cfctl%procmin   = 0       ! Minimum area number for reporting [default:0]
254       sn_cfctl%procmax   = 1000000 ! Maximum area number for reporting [default:1000000]
255       sn_cfctl%procincr  = 1       ! Increment for optional subsetting of areas [default:1]
256       sn_cfctl%ptimincr  = 1       ! Timestep increment for writing time step progress info
257   nn_print    =    0      !  level of print (0 no extra print)
258   nn_ictls    =    0      !  start i indice of control sum (use to compare mono versus
259   nn_ictle    =    0      !  end   i indice of control sum        multi processor runs
260   nn_jctls    =    0      !  start j indice of control               over a subdomain)
261   nn_jctle    =    0      !  end   j indice of control
262   nn_isplt    =    1      !  number of processors in i-direction
263   nn_jsplt    =    1      !  number of processors in j-direction
264   ln_timing   = .false.   !  timing by routine write out in timing.output file
265   ln_diacfl   = .false.   !  CFL diagnostics write out in cfl_diagnostics.ascii
266/
Note: See TracBrowser for help on using the repository browser.