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_cfg in branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/EXP00 – NEMO

source: branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/EXP00/namelist_cfg @ 9089

Last change on this file since 9089 was 9089, checked in by gm, 6 years ago

dev_merge_2017: bug correction in zdfdrg + ISOMIP cfg

  • Property svn:mime-type set to text/x-fortran
File size: 23.6 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  : ISOMIP Configuration namelist to overwrite reference dynamical namelist
3!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4!-----------------------------------------------------------------------
5&namusr_def    !   ISOMIP user defined namelist 
6!-----------------------------------------------------------------------
7   ln_zps      = .true.    ! z-partial-step coordinate
8   rn_e1deg    =   0.3     !  zonal      grid-spacing (degrees)
9   rn_e2deg    =   0.1     !  meridional grid-spacing (degrees)
10   rn_e3       =   30.     !  vertical resolution
11/
12!-----------------------------------------------------------------------
13&namrun        !   parameters of the run
14!-----------------------------------------------------------------------
15   cn_exp      =  "ISOMIP" !  experience name
16   nn_it000    =       1   !  first time step
17   nn_itend    =  525600   !  last  time step
18   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
19   ln_clobber  = .true.    !  clobber (overwrite) an existing file
20   nn_istate   =       0   !  output the initial state (1) or not (0)
21   nn_stock    =  99999999 !  frequency of creation of a restart file (modulo referenced to 1)
22   nn_write    =      48   !  frequency of write in the output file   (modulo referenced to nn_it000)
23   nn_istate   =       0   !  output the initial state (1) or not (0)
24/
25!-----------------------------------------------------------------------
26&namcfg     !   parameters of the configuration   
27!-----------------------------------------------------------------------
28   ln_read_cfg = .false.   !  (=T) read the domain configuration file
29      !                    !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules
30      cn_domcfg = "ISOMIP_zps_domcfg"         ! domain configuration filename
31      !
32   ln_write_cfg= .false.   !  (=T) create the domain configuration file
33      cn_domcfg_out = "ISOMIP_cfg_out"    ! filename of the created file
34/
35!-----------------------------------------------------------------------
36&namdom        !   space and time domain (bathymetry, mesh, timestep)
37!-----------------------------------------------------------------------
38   ln_linssh   = .false.    !  =T  linear free surface  ==>>  model level are fixed in time
39   !
40   nn_msh      =    0      !  create (>0) a mesh file or not (=0)
41   !
42   rn_rdt      = 1800.     !  time step for the dynamics (and tracer if nn_acc=0)
43/
44!-----------------------------------------------------------------------
45&namcrs        !   Grid coarsening for dynamics output and/or
46!              !   passive tracer coarsened online simulations
47!-----------------------------------------------------------------------
48/
49!-----------------------------------------------------------------------
50&namtsd    !   data : Temperature  & Salinity
51!-----------------------------------------------------------------------
52   cn_dir        = './'      !  root directory for the location of the runoff files
53   ln_tsd_init   = .false.   !  Initialisation of ocean T & S with T &S input data (T) or not (F)
54   ln_tsd_tradmp = .false.   !  damping of ocean T & S toward T &S input data (T) or not (F)
55/
56!-----------------------------------------------------------------------
57&namsbc        !   Surface Boundary Condition (surface module)
58!-----------------------------------------------------------------------
59   nn_fsbc     = 1         !  frequency of surface boundary condition computation
60   !                       !     (also = the frequency of sea-ice model call)
61   ln_blk      = .false.    !  Bulk formulation                          (T => fill namsbc_blk )
62   ln_usr      = .true.    !  user defined formulation                  (T => check usrdef_sbc)
63   nn_ice      = 0         !  =0 no ice boundary condition   ,
64                           !  =1 use observed ice-cover      ,
65                           !  =2 ice-model used               
66   ln_ice_embd = .false.   !  =F levitating ice with mass and salt exchange but no presure effect
67                           !  =T embedded sea-ice (full salt and mass exchanges and pressure)
68                     ! Misc. options of sbc :
69   ln_traqsr   = .false.   !  Light penetration (T) or not (F)
70   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf)
71   ln_isf      = .true.    !  ice shelf melting/freezing                (T => fill namsbc_isf)
72   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
73   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked
74                           !     =1 global mean of e-p-r set to zero at each time step
75                           !     =2 annual global mean of e-p-r set to zero
76/
77!-----------------------------------------------------------------------
78&namsbc_flx    !   surface boundary condition : flux formulation
79!-----------------------------------------------------------------------
80/
81!-----------------------------------------------------------------------
82&namsbc_clio   !   namsbc_clio  CLIO bulk formulae
83!-----------------------------------------------------------------------
84/
85!-----------------------------------------------------------------------
86&namsbc_core   !   namsbc_core  CORE bulk formulae
87!-----------------------------------------------------------------------
88/
89!-----------------------------------------------------------------------
90&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae
91!-----------------------------------------------------------------------
92/
93!-----------------------------------------------------------------------
94&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
95!-----------------------------------------------------------------------
96/
97!-----------------------------------------------------------------------
98&namtra_qsr    !   penetrative solar radiation
99!-----------------------------------------------------------------------
100/
101!-----------------------------------------------------------------------
102&namsbc_rnf    !   runoffs namelist surface boundary condition
103!-----------------------------------------------------------------------
104/
105!-----------------------------------------------------------------------
106&namsbc_isf    !  Top boundary layer (ISF)
107!-----------------------------------------------------------------------
108!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
109!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
110! nn_isf == 4
111   sn_fwfisf   = 'rnfisf'  ,         -12       ,'sowflisf',     .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , ''
112! nn_isf == 3
113   sn_rnfisf   = 'rnfisf'  ,         -12       ,'sofwfisf',     .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , ''
114! nn_isf == 2 and 3
115   sn_depmax_isf = 'rnfisf' ,        -12       ,'sozisfmax' ,   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , ''
116   sn_depmin_isf = 'rnfisf' ,        -12       ,'sozisfmin' ,   .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , ''
117! nn_isf == 2
118   sn_Leff_isf = 'rnfisf'  ,         -12       ,'Leff'    ,     .false.    , .true.  , 'yearly'  ,  ''      ,   ''     , ''
119! for all case
120   nn_isf      = 1         !  ice shelf melting/freezing
121                           !  1 = presence of ISF    2 = bg03 parametrisation
122                           !  3 = rnf file for isf   4 = ISF fwf specified
123                           !  option 1 and 4 need ln_isfcav = .true. (domzgr)
124! only for nn_isf = 1 or 2
125   rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula
126   rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula
127! only for nn_isf = 1 or 4
128   rn_hisf_tbl =  30.      ! thickness of the top boundary layer    (Losh et al. 2008)
129                          ! 0 => thickness of the tbl = thickness of the first wet cell
130! only for nn_isf = 1
131   nn_isfblk   = 1        ! 1 ISOMIP  like: 2 equations formulation (Hunter et al., 2006)
132                          ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015)
133   nn_gammablk = 0        ! 0 = cst Gammat (= gammat/s)
134                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010)
135                          ! 2 = velocity and stability dependent Gamma    (Holland et al. 1999)
136/
137!-----------------------------------------------------------------------
138&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk
139!-----------------------------------------------------------------------
140/
141!-----------------------------------------------------------------------
142&namsbc_ssr    !   surface boundary condition : sea surface restoring
143!-----------------------------------------------------------------------
144/
145!-----------------------------------------------------------------------
146&namsbc_alb    !   albedo parameters
147!-----------------------------------------------------------------------
148/
149!-----------------------------------------------------------------------
150&namberg       !   iceberg parameters
151!-----------------------------------------------------------------------
152/
153!-----------------------------------------------------------------------
154&namlbc        !   lateral momentum boundary condition
155!-----------------------------------------------------------------------
156                           !  free slip  !   partial slip  !   no slip   ! strong slip
157   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
158   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs.
159/
160!-----------------------------------------------------------------------
161&namagrif      !  AGRIF zoom                                            ("key_agrif")
162!-----------------------------------------------------------------------
163/
164!-----------------------------------------------------------------------
165&nam_tide      !    tide parameters (#ifdef key_tide)
166!-----------------------------------------------------------------------
167/
168!-----------------------------------------------------------------------
169&nambdy        !  unstructured open boundaries                          ("key_bdy")
170!-----------------------------------------------------------------------
171/
172!-----------------------------------------------------------------------
173&nambdy_dta      !  open boundaries - external data           ("key_bdy")
174!-----------------------------------------------------------------------
175/
176!-----------------------------------------------------------------------
177&nambdy_tide     ! tidal forcing at open boundaries
178!-----------------------------------------------------------------------
179/
180!-----------------------------------------------------------------------
181&namdrg            !   top/bottom drag coefficient                      (default: NO selection)
182!-----------------------------------------------------------------------
183   ln_NONE    = .false.    !  free-slip       : Cd = 0                  (F => fill namdrg_bot
184   ln_lin     = .false.    !      linear  drag: Cd = Cd0 Uc0                   &   namdrg_top)
185   ln_non_lin = .true.     !  non-linear  drag: Cd = Cd0 |U|
186   ln_loglayer= .false.    !  logarithmic drag: Cd = vkarmn/log(z/z0) |U|
187   !
188   ln_drgimp  = .false.    !  implicit top/bottom friction flag
189/
190!-----------------------------------------------------------------------
191&namdrg_top        !   TOP friction                                     (ln_isfcav=T)
192!-----------------------------------------------------------------------
193   rn_Cd0     =  2.5e-3    !  drag coefficient [-]
194   rn_Uc0     =  0.16      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)
195   rn_Cdmax   =  0.1       !  drag value maximum [-] (logarithmic drag)
196   rn_ke0     =  0.0e-0    !  background kinetic energy  [m2/s2] (non-linear cases)
197   rn_z0      =  3.0e-3    !  roughness [m] (ln_loglayer=T)
198   ln_boost   = .false.    !  =T regional boost of Cd0 ; =F constant
199      rn_boost=  50.          !  local boost factor  [-]
200/
201!-----------------------------------------------------------------------
202&namdrg_bot        !   BOTTOM friction                                 
203!-----------------------------------------------------------------------
204   rn_Cd0     =  1.e-3    !  drag coefficient [-]
205   rn_Uc0     =  0.4      !  ref. velocity [m/s] (linear drag=Cd0*Uc0)
206   rn_Cdmax   =  0.1      !  drag value maximum [-] (logarithmic drag)
207   rn_ke0     =  2.5e-3   !  background kinetic energy  [m2/s2] (non-linear cases)
208   rn_z0      =  3.e-3    !  roughness [m] (ln_loglayer=T)
209   ln_boost   = .false.   !  =T regional boost of Cd0 ; =F constant
210      rn_boost=  50.         !  local boost factor  [-]
211/
212!-----------------------------------------------------------------------
213&nambbc        !   bottom temperature boundary condition
214!-----------------------------------------------------------------------
215   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom
216/
217!-----------------------------------------------------------------------
218&nambbl        !   bottom boundary layer scheme
219!-----------------------------------------------------------------------
220   nn_bbl_ldf  =  0      !  diffusive bbl (=1)   or not (=0)
221   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0)
222/
223!-----------------------------------------------------------------------
224&nameos        !   ocean physical parameters
225!-----------------------------------------------------------------------
226   ln_teos10   = .false.         !  = Use TEOS-10 equation of state
227   ln_eos80    = .true.          !  = Use EOS80 equation of state
228/
229!-----------------------------------------------------------------------
230&namtra_adv    !   advection scheme for tracer
231!-----------------------------------------------------------------------
232   ln_traadv_fct =  .true.   !  FCT scheme
233      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order
234      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order
235/
236!-----------------------------------------------------------------------
237&namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param)
238!-----------------------------------------------------------------------
239/
240!----------------------------------------------------------------------------------
241&namtra_ldf    !   lateral diffusion scheme for tracers
242!----------------------------------------------------------------------------------
243   !                       !  Operator type:
244   ln_traldf_lap   =  .true.   !    laplacian operator
245   ln_traldf_blp   =  .false.  !  bilaplacian operator
246   !                       !  Direction of action:
247   ln_traldf_lev   =  .false.  !  iso-level
248   ln_traldf_hor   =  .true.   !  horizontal (geopotential)
249   ln_traldf_iso   =  .false.  !  iso-neutral (standard operator)
250   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator)
251   !
252   !                       !  iso-neutral options:       
253   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators)
254   rn_slpmax       =   0.01    !  slope limit                      (both operators)
255   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only)
256   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only)
257   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only)
258   !
259   !                       !  Coefficients:
260   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef
261   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
262   !                                !   =  0           constant
263   !                                !   = 10 F(k)      =ldf_c1d
264   !                                !   = 20 F(i,j)    =ldf_c2d
265   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
266   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
267   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
268   rn_aht_0        =  100.     !  lateral eddy diffusivity   (lap. operator) [m2/s]
269   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s]
270/
271!----------------------------------------------------------------------------------
272&namtra_ldfeiv !   eddy induced velocity param.
273!----------------------------------------------------------------------------------
274   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization
275/
276!-----------------------------------------------------------------------
277&namtra_dmp    !   tracer: T & S newtonian damping
278!-----------------------------------------------------------------------
279   ln_tradmp   =  .false.  !  add a damping termn (T) or not (F)
280/
281!-----------------------------------------------------------------------
282&namdyn_adv    !   formulation of the momentum advection
283!-----------------------------------------------------------------------
284   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)
285   nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction
286   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme
287   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
288/
289!-----------------------------------------------------------------------
290/
291!-----------------------------------------------------------------------
292&namdyn_vor    !   option of physics/algorithm (not control by CPP keys)
293!-----------------------------------------------------------------------
294   ln_dynvor_ene = .true.  !  enstrophy conserving scheme
295   ln_dynvor_ens = .false. !  energy conserving scheme
296   ln_dynvor_mix = .false. !  mixed scheme
297   ln_dynvor_een = .false. !  energy & enstrophy scheme
298      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)
299/
300!-----------------------------------------------------------------------
301&namdyn_hpg    !   Hydrostatic pressure gradient option
302!-----------------------------------------------------------------------
303   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation)
304   ln_hpg_isf  = .true.    !  s-coordinate adapted for isf (standard jacobian formulation)
305/
306!-----------------------------------------------------------------------
307&namdyn_spg    !   surface pressure gradient
308!-----------------------------------------------------------------------
309   ln_dynspg_ts  = .true.  !  split-explicit free surface
310/
311!-----------------------------------------------------------------------
312&namdyn_ldf    !   lateral diffusion on momentum
313!-----------------------------------------------------------------------
314   !                       !  Type of the operator :
315   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F
316   ln_dynldf_lap =  .true.     !    laplacian operator
317   ln_dynldf_blp =  .false.    !  bilaplacian operator
318   !                       !  Direction of action  :
319   ln_dynldf_lev =  .false.    !  iso-level
320   ln_dynldf_hor =  .true.     !  horizontal (geopotential)
321   ln_dynldf_iso =  .false.    !  iso-neutral
322   !                       !  Coefficient
323   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef
324   !                                !  =-30  read in eddy_viscosity_3D.nc file
325   !                                !  =-20  read in eddy_viscosity_2D.nc file
326   !                                !  =  0  constant
327   !                                !  = 10  F(k)=c1d
328   !                                !  = 20  F(i,j)=F(grid spacing)=c2d
329   !                                !  = 30  F(i,j,k)=c2d*c1d
330   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity)
331   rn_ahm_0      =    600.     !  horizontal laplacian eddy viscosity   [m2/s]
332   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s]
333   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s]
334   !
335   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km)
336/
337!-----------------------------------------------------------------------
338&namzdf        !   vertical physics                                     (default: NO selection)
339!-----------------------------------------------------------------------
340   !                       ! type of vertical closure
341   ln_zdfcst   = .true.       !  constant mixing
342   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric)
343   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke)
344   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls)
345   ln_zdfosm   = .false.      !  OSMOSIS BL closure                     (T =>   fill namzdf_osm)
346   !
347   !                       ! convection
348   ln_zdfevd   = .true.       !  enhanced vertical diffusion
349      nn_evdm     =    1         ! apply on tracer (=0) or on tracer and momentum (=1)
350      rn_evd      =   0.1        ! mixing coefficient [m2/s]
351   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm
352      nn_npc      =    1         ! frequency of application of npc
353      nn_npcp     =  365         ! npc control print frequency
354   !
355   ln_zdfddm   = .false.   ! double diffusive mixing
356   !
357   !                       ! gravity wave-driven vertical mixing
358   ln_zdfiwm   = .false.      ! internal wave-induced mixing            (T =>   fill namzdf_iwm)
359   ln_zdfswm   = .false.      ! surface  wave-induced mixing            (T => ln_wave=ln_sdw=T )
360   !
361   !                       ! coefficients
362   rn_avm0     =   1.e-3     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F)
363   rn_avt0     =   5.e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F)
364   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0)
365   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0)
366/
367!-----------------------------------------------------------------------
368&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi)
369!-----------------------------------------------------------------------
370/
371!-----------------------------------------------------------------------
372&namctl        !   Control prints & Benchmark
373!-----------------------------------------------------------------------
374/
375!-----------------------------------------------------------------------
376&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
377!-----------------------------------------------------------------------
378/
379!-----------------------------------------------------------------------
380&namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra")
381!              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor")
382!-----------------------------------------------------------------------
383/
384!-----------------------------------------------------------------------
385&namhsb       !  Heat and salt budgets
386!-----------------------------------------------------------------------
387/
388!-----------------------------------------------------------------------
389&namobs       !  observation usage switch                               (ln_diaobs =T)
390!-----------------------------------------------------------------------
391/
392!-----------------------------------------------------------------------
393&nam_asminc   !   assimilation increments                               ('key_asminc')
394!-----------------------------------------------------------------------
395/
396!-----------------------------------------------------------------------
397&namsbc_wave   ! External fields from wave model
398!-----------------------------------------------------------------------
399/
Note: See TracBrowser for help on using the repository browser.