source: branches/2017/dev_r7663_ISOMIP_TEST_CASE/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/EXP00/namelist_cfg @ 7668

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

Change location of ISOMIP, add plot routine and expected melt/moc/psi figure and README

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