source: CONFIG/IPSLCM/IPSLCM5/trunk/EXP00/PARAM/namelist_ORCA2 @ 555

Last change on this file since 555 was 555, checked in by aclsce, 15 years ago

Modified orca2 namelist parameter concerning wind stress coupling
variable (no impact on the results)

File size: 46.6 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/OPA  :  1 - run manager      (namrun)
3!! namelists    2 - Domain           (nam_zgr, nam_zgr_sco, namdom)
4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core
5!!                                    namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb)
6!!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, namtide)
7!!              5 - bottom  boundary (nambfr, nambbc, nambbl)
8!!              6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp)
9!!              7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf)
10!!              8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm)
11!!              9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr)
12!!              9 - miscellaneous    (namsol, nam_mpp, nam_mpp_dyndist, namctl)
13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
14!  CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE.
15
16!!======================================================================
17!!                   ***  Run management namelists  ***
18!!======================================================================
19!!   namrun            parameters of the run
20!!======================================================================
21
22!-----------------------------------------------------------------------
23&namrun        !   parameters of the run
24!-----------------------------------------------------------------------
25   no          =       0       !  job number
26   cexper      =  "ORCA2"      !  experience name
27   cn_ocerst_in  = "restartopa"!  suffix of ocean restart name (input)
28   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output)
29   ln_rstart   = .false.       !  start from rest (F) or from a restart file (T)
30   nrstdt      =       1       !  restart control = 0 nit000 is not compared to the restart file value
31                               !                  = 1 use ndate0 in namelist (not the value in the restart file)
32                               !                  = 2 calendar parameters read in the restart file
33   nit000      =      76       !  first time step
34   nitend      =     150       !  last  time step
35   ndate0      =18600106       !  initial calendar date yymmdd (used if nrstdt=1)
36   nleapy      =      30       !  Leap year calendar (1) or not (0)
37   ninist      =       0       !  output the initial state (1) or not (0)
38   nstock      =    5400       !  frequency of creation of a restart file
39   nwrite      =       1       !  frequency of write in the output file
40   ln_dimgnnn  = .false.       !  DIMG file format: 1 file for all processors (F) or by processor (T)
41/
42!!======================================================================
43!!                      ***  Domain namelists  ***
44!!======================================================================
45!!   nam_zgr       vertical coordinate
46!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
47!!   namdom        space and time domain (bathymetry, mesh, timestep)
48!!======================================================================
49
50!-----------------------------------------------------------------------
51&nam_zgr       !   vertical coordinate
52!-----------------------------------------------------------------------
53   ln_zco      = .true.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined)
54   ln_zps      = .false.    !  z-coordinate - partial steps   (T/F)
55   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F)
56/
57!-----------------------------------------------------------------------
58&nam_zgr_sco   !   s-coordinate or hybrid z-s-coordinate
59!-----------------------------------------------------------------------
60   sbot_min    =  300.     !  minimum depth of s-bottom surface (>0) (m)
61   sbot_max    = 5250.     !  maximum depth of s-bottom surface (= ocean depth) (>0) (m)
62   theta       =    6.0    !  surface control parameter (0<=theta<=20)
63   thetb       =    0.75   !  bottom control parameter  (0<=thetb<= 1)
64   r_max       =    0.15   !  maximum cut-off r-value allowed (0<r_max<1)
65/
66!-----------------------------------------------------------------------
67&namdom        !   space and time domain (bathymetry, mesh, timestep)
68!-----------------------------------------------------------------------
69   ntopo       =    1      !  compute (=0) or read(=1) the bathymetry file
70   e3zps_min   =    5.     !  the thickness of the partial step is set larger than the minimum
71   e3zps_rat   =    0.1    !  of e3zps_min and e3zps_rat * e3t   (N.B. 0<e3zps_rat<1)
72   nmsh        =    0      !  create (=1) a mesh file (coordinates, scale factors, masks) or not (=0)
73   nacc        =    0      !  =1 acceleration of convergence method used, rdt < rdttra(k)
74                           !  =0, no acceleration, rdt = rdttra
75   atfp        =    0.1    !  asselin time filter parameter
76   rdt         = 5760.     !  time step for the dynamics (and tracer if nacc=0)
77   rdtmin      = 5760.     !  minimum time step on tracers (used if nacc=1)
78   rdtmax      = 5760.     !  maximum time step on tracers (used if nacc=1)
79   rdth        =  800.     !  depth variation of tracer time step  (used if nacc=1)
80   nn_baro     =   64      !  number of barotropic time step (for the split explicit algorithm) ("key_dynspg_ts")
81   nclosea     =    1      !  = 0 no closed sea in the model domain
82                           !  = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...)
83/
84!!======================================================================
85!!            ***  Surface Boundary Condition namelists  ***
86!!======================================================================
87!!   namsbc        surface boundary condition
88!!   namsbc_ana    analytical         formulation
89!!   namsbc_flx    flux               formulation
90!!   namsbc_clio   CLIO bulk formulea formulation
91!!   namsbc_core   CORE bulk formulea formulation
92!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled")
93!!   namqsr        penetrative solar radiation
94!!   namsbc_rnf    river runoffs
95!!   namsbc_ssr    sea surface restoring term (for T and/or S)
96!!   namalb        albedo parameters
97!!======================================================================
98
99!-----------------------------------------------------------------------
100&namsbc        !   Surface Boundary Condition (surface module)
101!-----------------------------------------------------------------------
102   nn_fsbc     = 5         !  frequency of surface boundary condition computation
103                           !               (= the frequency of sea-ice model call)
104   ln_ana      = .false.   !  analytical formulation (T => fill namsbc_ana )
105   ln_flx      = .false.   !  flux formulation       (T => fill namsbc_flx )
106   ln_blk_clio = .false.    !  CLIO bulk formulation  (T => fill namsbc_clio)
107   ln_blk_core = .false.   !  CORE bulk formulation  (T => fill namsbc_core)
108   ln_cpl      = .true.   !  Coupled formulation    (T => fill namsbc_cpl )
109   nn_ice      = 2         !  =0 no ice boundary condition   ,
110                           !  =1 use observed ice-cover      ,
111                           !  =2 ice-model used                             ("key_lim3" or "key_lim2)
112   nn_ico_cpl  = 0         !  ice-ocean coupling : =0 each nn_fsbc
113                           !                       =1 stresses recomputed each ocean time step ("key_lim3" only)
114                           !                       =2 combination of 0 and 1 cases             ("key_lim3" only)
115   ln_dm2dc    = .false.   !  daily mean to diurnal cycle short wave (qsr)
116   ln_rnf      = .true.    !  runoffs (T => fill namsbc_rnf)
117   ln_ssr      = .false.    !  Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
118   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked                              ,
119                           !                     =1 global mean of e-p-r set to zero at each nn_fsbc time step   ,
120                           !                     =2 annual global mean of e-p-r set to zero
121/
122!-----------------------------------------------------------------------
123&namsbc_ana    !   analytical surface boundary condition
124!-----------------------------------------------------------------------
125   nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps
126   rn_utau0    =   0.5     !  uniform value for the i-stress
127   rn_vtau0    =   0.e0    !  uniform value for the j-stress
128   rn_q0       =   0.e0    !  uniform value for the total heat flux
129   rn_qsr0     =   0.e0    !  uniform value for the solar radiation
130   rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P)
131/
132!-----------------------------------------------------------------------
133&namsbc_flx    !   surface boundary condition : flux formulation
134!-----------------------------------------------------------------------
135!              !      file name     ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly' or !
136!              !                    !  (if <0  months)  !    name    !    (logical)   !  (T/F)  !  'monthly'  !
137   sn_utau     =      'utau'        ,        24.        ,    'utau'  ,     .false.    , .false. ,   'yearly'
138   sn_vtau     =      'vtau'        ,        24.        ,    'vtau'  ,     .false.    , .false. ,   'yearly'
139   sn_qtot     =      'qtot'        ,        24.        ,    'qtot'  ,     .false.    , .false. ,   'yearly'
140   sn_qsr      =      'qsr'         ,        24.        ,    'qsr'   ,     .false.    , .false. ,   'yearly'
141   sn_emp      =      'emp'         ,        24.        ,    'emp'   ,     .false.    , .false. ,   'yearly'
142!
143   cn_dir      = './'      !  root directory for the location of the flux files
144/     
145!-----------------------------------------------------------------------
146&namsbc_clio   !   namsbc_clio  CLIO bulk formulea
147!-----------------------------------------------------------------------
148!              !      file name     ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly' or !
149!              !                    !  (if <0  months)  !    name    !    (logical)   !  (T/F)  !  'monthly'  !
150   sn_utau     =    'taux_1m'       ,       -1.         , 'sozotaux' ,    .true.      , .true.  ,   'yearly'   
151   sn_vtau     =    'tauy_1m'       ,       -1.         , 'sometauy' ,    .true.      , .true.  ,   'yearly'   
152   sn_wndm     =    'flx'           ,       -1.         , 'socliowi' ,    .true.      , .true.  ,   'yearly'     
153   sn_tair     =    'flx'           ,       -1.         , 'socliot2' ,    .true.      , .true.  ,   'yearly'     
154   sn_humi     =    'flx'           ,       -1.         , 'socliohu' ,    .true.      , .true.  ,   'yearly'   
155   sn_ccov     =    'flx'           ,       -1.         , 'socliocl' ,    .false.     , .true.  ,   'yearly'     
156   sn_prec     =    'flx'           ,       -1.         , 'socliopl' ,    .false.     , .true.  ,   'yearly'       
157!
158   cn_dir      = './'      !  root directory for the location of the bulk files are
159/
160!-----------------------------------------------------------------------
161&namsbc_core   !   namsbc_core  CORE bulk formulea
162!-----------------------------------------------------------------------
163!              !   file name        ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly' or !
164!              !                    !  (if <0  months)  !   name     !    (logical)   !  (T/F)  !  'monthly'  !
165   sn_wndi     =    'u10'           ,       24.         , 'U_10_MOD' ,    .false.     , .true.  ,   'yearly'   
166   sn_wndj     =    'v10'           ,       24.         , 'V_10_MOD' ,    .false.     , .true.  ,   'yearly'   
167   sn_qsr      =    'rad'           ,       24.         , 'SWDN_MOD' ,    .false.     , .true.  ,   'yearly'     
168   sn_qlw      =    'rad'           ,       24.         , 'LWDN_MOD' ,    .false.     , .true.  ,   'yearly'     
169   sn_tair     =    't10'           ,       24.         , 'T_10_MOD' ,    .false.     , .true.  ,   'yearly'     
170   sn_humi     =    'q10'           ,       24.         , 'Q_10_MOD' ,    .false.     , .true.  ,   'yearly'   
171   sn_prec     =    'precip'        ,       -1.         , 'RAIN'     ,    .false.     , .true.  ,   'yearly'     
172   sn_snow     =    'precip'        ,       -1.         , 'SNOW'     ,    .false.     , .true.  ,   'yearly'       
173!
174   cn_dir      = './'      !  root directory for the location of the bulk files
175   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F)
176   alpha_precip= 1.        !  multiplicative factor for precipitation (total & snow)
177/
178!-----------------------------------------------------------------------
179&namsbc_cpl    !   coupled ocean/atmosphere model                        ("key_coupled")
180!-----------------------------------------------------------------------
181! SEND
182cn_snd_temperature= 'weighted oce and ice'  ! 'oce only' 'weighted oce and ice' 'mixed oce-ice'
183cn_snd_albedo     = 'weighted ice'          ! 'none' 'weighted ice' 'mixed oce-ice'
184cn_snd_thickness  = 'none'                  ! 'none' 'weighted ice and snow'
185cn_snd_crt_nature = 'none'                  ! 'none' 'oce only' 'weighted oce and ice' 'mixed oce-ice'
186cn_snd_crt_refere = 'spherical'             ! 'spherical' 'cartesian'
187cn_snd_crt_orient = 'eastward-northward'    ! 'eastward-northward' or 'local grid'
188cn_snd_crt_grid   = 'T'                     ! 'T'
189! RECEIVE
190cn_rcv_w10m       = 'none'                  ! 'none' 'coupled'
191cn_rcv_tau_nature = 'mixed oce-ice'         ! 'oce only' 'oce and ice' 'mixed oce-ice'
192cn_rcv_tau_refere = 'cartesian'             ! 'spherical' 'cartesian'
193cn_rcv_tau_orient = 'eastward-northward'    ! 'eastward-northward' or 'local grid'
194cn_rcv_tau_grid   = 'U,V'                   ! 'T' 'U,V' 'U,V,F' 'U,V,I' 'T,F' 'T,I' 'T,U,V'
195cn_rcv_dqnsdt     = 'coupled'               ! 'none' 'coupled'
196cn_rcv_qsr        = 'oce and ice'           ! 'conservative' 'oce and ice' 'mixed oce-ice'
197cn_rcv_qns        = 'oce and ice'           ! 'conservative' 'oce and ice' 'mixed oce-ice'
198cn_rcv_emp        = 'oce and ice'           ! 'conservative' 'oce and ice' 'mixed oce-ice'
199cn_rcv_rnf        = 'coupled'               ! 'coupled' 'climato' 'mixed'
200cn_rcv_cal        = 'coupled'               ! 'none' 'coupled'
201/
202!-----------------------------------------------------------------------
203&namqsr        !   penetrative solar radiation
204!-----------------------------------------------------------------------
205   ln_traqsr   = .true.    !  penetrative solar radiation (T) or not (F)
206   rabs        =   0.58    !  fraction of qsr associated with xsi1
207   xsi1        =   0.35    !  first depth of extinction
208   xsi2        =   23.0    !  second depth of extinction
209/
210!-----------------------------------------------------------------------
211&namsbc_rnf    !   runoffs namelist surface boundary condition
212!-----------------------------------------------------------------------
213!              !     file name      ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly' or !
214!              !                    !  (if <0  months)  !   name     !    (logical)   !  (T/F)  !  'monthly'  !
215   sn_rnf    =   'runoff_1m_nomask' ,        -1.        , 'sorunoff' ,    .true.      , .true.  ,   'yearly'     
216   sn_cnf    =   'runoff_1m_nomask' ,         0.        , 'socoefr'  ,    .false.     , .true.  ,   'yearly'             
217!
218   cn_dir       = './'      !  root directory for the location of the runoff files
219   ln_rnf_emp   =   .true. !  runoffs included into precipitation field (T) or into a file (F)
220   ln_rnf_mouth =   .true. !  specific treatment at rivers mouths
221   rn_hrnf      =   0.e0    !  depth over which enhanced vertical mixing is used
222   rn_avt_rnf   =   1.e-3   !  value of the additional vertical mixing coef. [m2/s]
223/
224!-----------------------------------------------------------------------
225&namsbc_ssr    !   surface boundary condition : sea surface restoring
226!-----------------------------------------------------------------------
227!              !     file name      ! frequency (hours) ! variable   ! time interpol. !  clim   ! 'yearly' or !
228!              !                    !  (if <0  months)  !   name     !    (logical)   !  (T/F)  !  'monthly'  !
229   sn_sst      =    'sst_data'     ,        24.        ,  'sst'     ,     .false.    , .false. ,   'yearly'
230   sn_sss      =    'sss_data'     ,        -1.        ,  'sss'     ,     .true.     , .true.  ,   'yearly'
231!   
232   cn_dir      = './'      !  root directory for the location of the runoff files
233   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0)
234   nn_sssr     =     1     !  add a damping     term in the surface freshwater flux (=1) or not (=0)
235   dqdt        =   -40.    !  magnitude of the retroaction on temperature   [W/m2/K]
236   deds        =   -27.7   !  magnitude of the damping on salinity   [mm/day/psu]
237/     
238!-----------------------------------------------------------------------
239&namalb        !   albedo parameters
240!-----------------------------------------------------------------------
241   cgren       =    0.0    !  correction of the snow or ice albedo to take into account the
242   albice      =    0.5    !  albedo of melting ice in the arctic and antarctic
243   alphd       =    0.80   !  coefficients for linear interpolation used to
244   alphc       =    0.65   !  compute albedo between two extremes values
245   alphdi      =    0.72   !  (Pyane, 1972)
246/
247!!======================================================================
248!!               ***  Lateral boundary condition  ***
249!!======================================================================
250!!   namlbc        lateral momentum boundary condition
251!!   namcla        cross land advection
252!!   namobc        open boundaries parameters                           ("key_obc")
253!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")
254!!   nambdy        Unstructured open boundaries                         ("key_bdy")
255!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides")
256!!======================================================================
257
258!-----------------------------------------------------------------------
259&namlbc        !   lateral momentum boundary condition
260!-----------------------------------------------------------------------
261   shlat       =    2.     !      shlat = 0 : free slip
262                           !  0 < shlat < 2 : partial slip
263                           !      shlat = 2 : no slip
264                           !  2 < shlat     : strong slip
265/
266!-----------------------------------------------------------------------
267&namcla        !   cross land advection
268!-----------------------------------------------------------------------
269   n_cla       =    0      !  advection between 2 ocean pts separates by land
270/
271!-----------------------------------------------------------------------
272&namobc        !   open boundaries parameters                           ("key_obc")
273!-----------------------------------------------------------------------
274    nobc_dta   =    1      !  = 0 the obc data are equal to the initial state
275                           !  = 1 the obc data are read in 'obc.dta' files
276    cffile     = 'annual'  !  set to annual if obc datafile hold 1 year of data
277                           !  set to monthly if obc datafile hold 1 month of data
278    rdpein     =    1.     !  ???
279    rdpwin     =    1.     !  ???
280    rdpnin     =    1.     !  ???
281    rdpsin     =    1.     !  ???
282    rdpeob     = 3000.     !  time relaxation (days) for the east  open boundary
283    rdpwob     =   15.     !    "        "             "     west         "
284    rdpnob     = 3000.     !    "        "             "     north        "
285    rdpsob     =   15.     !    "        "             "     south        "
286    zbsic1     =  140.e+6  !   barotropic stream function on first  isolated coastline
287    zbsic2     =    1.e+6  !    "                   "        second       "
288    zbsic3     =    0.     !    "                   "        thrid        "
289    ln_obc_clim= .false.   !  climatological obc data files (T) or not (F)
290    ln_vol_cst = .true.    !  impose the total volume conservation (T) or not (F)
291/
292!-----------------------------------------------------------------------
293&namagrif      !                                                        ("key_agrif")
294!-----------------------------------------------------------------------
295    nbclineupdate = 3      !  baroclinic update frequency
296    ln_spc_dyn    = .true. !  use 0 as special value for dynamics
297    visc_tra      = 2880.  !  viscosity coeeficient for tracers sponge layer
298    visc_dyn      = 2880.  !  viscosity coeeficient for dynamics sponge layer
299/
300!-----------------------------------------------------------------------
301&nambdy        !  unstructured open boundaries parameters               ("key_bdy")
302!-----------------------------------------------------------------------
303    filbdy_mask    =  ''                  !  name of mask file (if ln_bdy_mask=.TRUE.)
304    filbdy_data_T  = 'bdydata_grid_T.nc'  !  name of data file (T-points)
305    filbdy_data_U  = 'bdydata_grid_U.nc'  !  name of data file (U-points)
306    filbdy_data_V  = 'bdydata_grid_V.nc'  !  name of data file (V-points)
307    ln_bdy_clim    = .false.              !  contain 1 (T) or 12 (F) time dumps and be cyclic
308    ln_bdy_vol     = .true.               !  total volume correction (see volbdy parameter)
309    ln_bdy_mask    = .false.              !  boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F)
310    ln_bdy_tides   = .true.               !  Apply tidal harmonic forcing with Flather condition
311    ln_bdy_dyn_fla = .true.               !  Apply Flather condition to velocities
312    ln_bdy_tra_frs = .false.              !  Apply FRS condition to temperature and salinity
313    ln_bdy_dyn_frs = .false.              !  Apply FRS condition to velocities
314    nbdy_dta       =  1                   !  = 0, bdy data are equal to the initial state
315                                          !  = 1, bdy data are read in 'bdydata   .nc' files
316    nb_rimwidth    = 9                    !  width of the relaxation zone
317    volbdy         = 0                    !  = 0, the total water flux across open boundaries is zero
318                                          !  = 1, the total volume of the system is conserved
319/
320!-----------------------------------------------------------------------
321&namtide        ! tidal forcing at unstructured boundaries             
322!-----------------------------------------------------------------------
323    filtide      = 'bdytide_'           !  file name root of tidal forcing files
324    tide_cpt     = 'M2','S1'            !  names of tidal components used
325    tide_speed   = 28.984106, 15.000001 !  phase speeds of tidal components (deg/hour)
326    ln_tide_date = .false.              !  adjust tidal harmonics for start date of run
327/
328!!======================================================================
329!!                 ***  Bottom boundary condition  ***
330!!======================================================================
331!!   nambfr        bottom friction
332!!   nambbc        bottom temperature boundary condition                ("key_trabbc")
333!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv")
334!!======================================================================
335
336!-----------------------------------------------------------------------
337&nambfr        !   bottom friction
338!-----------------------------------------------------------------------
339   nbotfr      =    1      !  type of bottom friction :   = 0 : no   slip,  = 2 : nonlinear friction
340                           !                              = 3 : free slip,  = 1 :    linear friction
341   bfri1       =    4.e-4  !  bottom drag coefficient (linear case)
342   bfri2       =    1.e-3  !  bottom drag coefficient (non linear case)
343   bfeb2       =    2.5e-3 !  bottom turbulent kinetic energy background  (m^2/s^2)
344/
345!-----------------------------------------------------------------------
346&nambbc        !   bottom temperature boundary condition
347!-----------------------------------------------------------------------
348   ngeo_flux   =    2      !  geothermal heat flux = 0 no flux considered
349                           !                       = 1 constant flux
350                           !                       = 2 variable flux (read in geothermal_heating.nc in mW/m2)
351   ngeo_flux_const = 86.4e-3  !  Constant value of geothermal heat flux [W/m2]
352/
353!-----------------------------------------------------------------------
354&nambbl        !   bottom boundary layer scheme
355!-----------------------------------------------------------------------
356!                          !  diffusive bbl                             ("key_trabbl")
357!                          !  advective bbl                             ("key_trabbl_adv")
358   atrbbl      =  10000.   !  lateral mixing coefficient in the bbl  [m2/s]
359/
360!!======================================================================
361!!                        Tracer (T & S ) namelists
362!!======================================================================
363!!   nameos        equation of state
364!!   nam_traadv    advection scheme
365!!   nam_traldf    lateral diffusion scheme
366!!   namtdp        tracer newtonian damping                             ("key_tradmp")
367!!======================================================================
368
369!-----------------------------------------------------------------------
370&nameos        !   ocean physical parameters
371!-----------------------------------------------------------------------
372   neos        =    0      !  type of equation of state and Brunt-Vaisala frequency
373                           !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) )
374                           !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T )
375                           !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T )
376   ralpha      =    2.e-4  !  thermal expension coefficient (neos= 1 or 2)
377   rbeta       =    0.001  !  saline  expension coefficient (neos= 2)
378/
379!-----------------------------------------------------------------------
380&nam_traadv    !   advection scheme for tracer
381!-----------------------------------------------------------------------
382   ln_traadv_cen2   =  .true.   !  2nd order centered scheme   
383   ln_traadv_tvd    =  .false.  !  TVD scheme               
384   ln_traadv_muscl  =  .false.  !  MUSCL scheme             
385   ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
386   ln_traadv_ubs    =  .false.  !  UBS scheme                 
387/
388!-----------------------------------------------------------------------
389&nam_traldf    !   lateral diffusion scheme for tracer
390!-----------------------------------------------------------------------
391!                               !  Type of the operator :
392   ln_traldf_lap    =  .true.   !     laplacian operator       
393   ln_traldf_bilap  =  .false.  !     bilaplacian operator     
394                                !  Direction of action  :
395   ln_traldf_level  =  .false.  !     iso-level               
396   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T)
397   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp")
398!                               !  Coefficient
399   aht0        =  2000.         !     horizontal eddy diffusivity for tracers [m2/s]
400   ahtb0       =     0.         !     background eddy diffusivity for ldf_iso [m2/s]
401   aeiv0       =  2000.         !     eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv")
402/
403!-----------------------------------------------------------------------
404&namtdp        !   tracer newtonian damping                             ('key_tradmp')
405!-----------------------------------------------------------------------
406   ndmp        =   -1      !  type of damping in temperature and salinity
407                           !     ='latitude', damping poleward of 'ndmp' degrees and function
408                           !                  of the distance-to-coast. Red and Med Seas as ndmp=-1
409                           !     =-1 damping only in Med and Red Seas
410   ndmpf       =    1      !  create a damping.coeff NetCDF file (=1) or not (=0)
411   nmldmp      =    1      !  type of damping: =0 damping throughout the water column
412                           !                   =1 no damping in the mixed layer defined by avt >5cm2/s )
413                           !                   =2 no damping in the mixed layer defined rho<rho(surf)+.01 )
414   sdmp        =   50.     !  surface time scale for internal damping (days)
415   bdmp        =  360.     !  bottom  time scale for internal damping (days)
416   hdmp        =  800.     !  depth of transition between sdmp and bdmp (meters)
417/
418!!======================================================================
419!!                      ***  Dynamics namelists  ***
420!!======================================================================
421!!   nam_dynadv    formulation of the momentum advection
422!!   nam_dynvor    advection scheme
423!!   nam_dynhpg    hydrostatic pressure gradient
424!!   namflg        hydrostatic pressure gradient time stepping
425!!   nam_dynspg    surface pressure gradient                            (CPP key only)
426!!   nam_dynldf    lateral diffusion scheme
427!!======================================================================
428
429!-----------------------------------------------------------------------
430&nam_dynadv    !   formulation of the momentum advection
431!-----------------------------------------------------------------------
432   ln_dynadv_vec = .false.  !  vector form (T) or flux form (F) 
433   ln_dynadv_cen2= .true. !  flux form - 2nd order centered scheme
434   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme
435
436!-----------------------------------------------------------------------
437&nam_dynvor    !   option of physics/algorithm (not control by CPP keys)
438!-----------------------------------------------------------------------
439   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
440   ln_dynvor_ens = .true.  !  energy conserving scheme   
441   ln_dynvor_mix = .false. !  mixed scheme               
442   ln_dynvor_een = .false. !  energy & enstrophy scheme 
443/
444!-----------------------------------------------------------------------
445&nam_dynhpg    !   Hydrostatic pressure gradient option
446!-----------------------------------------------------------------------
447   ln_hpg_zco  = .true.   !  z-coordinate - full steps                   
448   ln_hpg_zps  = .false.    !  z-coordinate - partial steps (interpolation)
449   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation)
450   ln_hpg_hel  = .false.   !  s-coordinate (helsinki modification)
451   ln_hpg_wdj  = .false.   !  s-coordinate (weighted density jacobian)
452   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial)
453   ln_hpg_rot  = .false.   !  s-coordinate (ROTated axes scheme)
454   gamm        = 0.e0      !  weighting coefficient (wdj scheme)
455/
456!-----------------------------------------------------------------------
457&namflg        !   algorithm flags (algorithm not control by CPP keys)
458!-----------------------------------------------------------------------
459   ln_dynhpg_imp = .false. !  hydrostatic pressure gradient: semi-implicit time scheme  (T)
460                           !                                 centered      time scheme  (F)
461   nn_dynhpg_rst =  0      !  add dynhpg implicit variables in restart ot not (1/0)
462/
463!-----------------------------------------------------------------------
464!nam_dynspg    !   surface pressure gradient   (CPP key only)
465!-----------------------------------------------------------------------
466!                          !  explicit free surface                     ("key_dynspg_exp")
467!                          !  filtered free surface                     ("key_dynspg_flt")
468!                          !  split-explicit free surface               ("key_dynspg_ts")
469!                          !  rigid-lid                                 ("key_dynspg_rl")
470
471!-----------------------------------------------------------------------
472&nam_dynldf    !   lateral diffusion on momentum
473!-----------------------------------------------------------------------
474!                               !  Type of the operator :
475   ln_dynldf_lap    =  .true.   !     laplacian operator         
476   ln_dynldf_bilap  =  .false.  !     bilaplacian operator   
477!                               !  Direction of action  :
478   ln_dynldf_level  =  .false.  !     iso-level               
479   ln_dynldf_hor    =  .true.   !     horizontal (geopotential)        (require "key_ldfslp" in s-coord.)
480   ln_dynldf_iso    =  .false.  !     iso-neutral                      (require "key_ldfslp")
481                                !  Coefficient
482   ahm0    = 40000.             !     horizontal eddy viscosity   [m2/s]
483   ahmb0   =     0.             !     background eddy viscosity for ldf_iso [m2/s]
484/
485!!======================================================================
486!!             Tracers & Dynamics vertical physics namelists
487!!======================================================================
488!!       namzdf        vertical physics
489!!       namnpc        non penetrative convection                       
490!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      )
491!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      )
492!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      )
493!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      )
494!!======================================================================
495
496!-----------------------------------------------------------------------
497&namzdf        !   vertical physics
498!-----------------------------------------------------------------------
499   avm0        =   1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst")
500   avt0        =   1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst")
501   ln_zdfnpc   = .false.   !  convection: Non-Penetrative algorithm (T) or not (F)
502   ln_zdfevd   = .true.    !  convection: enhanced vertical diffusion (T) or not (F)   
503   avevd       = 100.      !  vertical coefficient for enhanced diffusion scheme [m2/s]
504   n_evdm      =   1       !  enhanced mixing apply on tracer (=0) or on tracer and momentum (=1)
505   ln_zdfexp   =  .false.  !  split explicit (T) or implicit (F) time stepping
506   n_zdfexp    =   3       !  number of sub-timestep for ln_zdfexp=T
507/
508!-----------------------------------------------------------------------
509&namnpc        !   non penetrative convection
510!-----------------------------------------------------------------------
511   nnpc1       =    1      !  non penetrative convective scheme computation frequency
512   nnpc2       =  365      !  non penetrative convective scheme print frequency
513/
514!-----------------------------------------------------------------------
515&namric        !   richardson number dependent vertical diffusion       ("key_zdfric" )
516!-----------------------------------------------------------------------
517   avmri       = 100.e-4   !  maximum value of the vertical viscosity
518   alp         =   5.      !  coefficient of the parameterization
519   nric        =   2       !  coefficient of the parameterization
520/
521!-----------------------------------------------------------------------
522&namtke        !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke")
523!-----------------------------------------------------------------------
524   ln_rstke    = .false.   !  restart with tke from a run without tke (T) or not (F)
525   nn_itke     =  50       !  number of iterative loops if ln_rstke=T
526   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
527   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation
528   rn_ebb      =   3.75    !  coef. of the surface input of tke
529   rn_efave    =   1.      !  boost of the tke diffusion ( avtke=rn_efave*avm )
530   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2]
531   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2]
532   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom
533                           !                 = 1 bounded by the local vertical scale factor
534                           !                 = 2 first vertical derivative of mixing length bounded by 1
535                           !                 = 3 same criteria as case 2 but applied in a different way
536   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
537   nn_avb      =   0       !  profile for constant background used on avt & avm (=1) or not (=0)
538   nn_ave      =   1       !  horizontal averaged on avt (=1) or not (=0)
539   ln_mxl0     = .false.   !  mixing length scale surface value as function of wind stress (T) or not (F)
540   rn_lmin     =   0.4     !  interior buoyancy lenght scale minimum value
541   rn_lmin0    =   0.4     !  surface  buoyancy lenght scale minimum value
542   nn_etau     =   0       !  exponentially deceasing penetration of tke due to internal & intertial waves
543                           !        = 0 no penetration ( O(2 km) resolution)
544                           !        = 1 additional tke source
545                           !        = 2 additional tke source applied only at the base of the mixed layer
546   nn_htau     =   2       !  type of exponential decrease of tke penetration
547                           !        = 0  constant 10 m length scale
548                           !        = 1  ???
549                           !        = 2  ???
550   rn_efr      =   0.05    !  fraction of surface tke value which penetrates inside the ocean
551   ln_lc       = .false.   !  Langmuir cell effect
552   rn_lc       =   0.15    !  coef. associated to Langmuir cells
553   nn_havtb    =   0       !  horizontal shape for avtb (=1) or not (=0)
554/
555!------------------------------------------------------------------------
556&namkpp        !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionnally:
557!                                                                         "key_kppcustom" or "key_kpplktb")
558!------------------------------------------------------------------------
559   ln_kpprimix = .true.    !  shear instability mixing
560   difmiw      =  1.0e-04  !  constant internal wave viscosity [m2/s]
561   difsiw      =  0.1e-04  !  constant internal wave diffusivity [m2/s]
562   Riinfty     =  0.8      !  local Richardson Number limit for shear instability
563   difri       =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s]
564   bvsqcon     = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]
565   difcon      =  1.       !  maximum mixing in interior convection [m2/s]
566   navb        =  0        !  horizontal averaged (=1) or not (=0) on avt and amv
567   nave        =  1        !  constant (=0) or profile (=1) background on avt
568/
569!-----------------------------------------------------------------------
570&namddm        !   double diffusive mixing parameterization                  ("key_zdfddm")
571!-----------------------------------------------------------------------
572      avts     = 1.e-4     !  maximum avs (vertical mixing on salinity)
573      hsbfr    = 1.6       !  heat/salt buoyancy flux ratio
574/
575!!======================================================================
576!!                  ***  Miscelaneous namelists  ***
577!!======================================================================
578!!   nam_mpp           Massively Parallel Processing                    ("key_mpp_mpi)
579!!   nam_mpp_dyndist   Massively Parallel domain decomposition          ("key_agrif" && "key_mpp_dyndist")
580!!   namctl            Control prints & Benchmark
581!!   namsol            elliptic solver / island / free surface
582!!======================================================================
583
584!-----------------------------------------------------------------------
585&namsol        !   elliptic solver / island / free surface
586!-----------------------------------------------------------------------
587   nsolv       =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg)
588                           !                   =2 successive-over-relaxation (sor)
589                           !                   =3 FETI (fet)                               ("key_feti")
590                           !                   =4 sor with extra outer halo
591   nsol_arp    =      0    !  absolute/relative (0/1) precision convergence test
592   nmin        =    300    !  minimum of iterations for the SOR solver
593   nmax        =    800    !  maximum of iterations for the SOR solver
594   nmod        =     10    !  frequency of test for the SOR solver
595   eps         =  1.e-6    !  absolute precision of the solver
596   resmax      =  1.e-10   !  absolute precision for the SOR solver
597   sor         =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain)
598   epsisl      =  1.e-10   !  absolute precision on stream function solver
599   nmisl       =   4000    !  maximum pcg iterations for island                            ("key_islands")
600   rnu         =      1.   !  strength of the additional force used in filtered free surface
601/
602!-----------------------------------------------------------------------
603&nam_mpp      !   Massively Parallel Processing                         ("key_mpp_mpi)
604!-----------------------------------------------------------------------
605   c_mpi_send =  'S'       !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
606                           !  buffer blocking send or immediate non-blocking sends, resp.
607   nn_buffer  =   0        !  size in bytes of exported buffer ('B' case), 0 no exportation
608/
609!-----------------------------------------------------------------------
610&nam_mpp_dyndist !   Massively Parallel Distribution                    ("key_agrif" && "key_mpp_dyndist")
611!-----------------------------------------------------------------------
612   jpni       =    1       !  jpni   number of processors following i
613   jpnj       =    1       !  jpnj   number of processors following j
614   jpnij      =    1       !  jpnij  number of local domains
615/
616!-----------------------------------------------------------------------
617&namctl       !   Control prints & Benchmark
618!-----------------------------------------------------------------------
619   ln_ctl     = .false.    !  trends control print (expensive!)
620   nprint     =    0       !  level of print (0 no extra print)
621   nictls     =    0       !  start i indice of control sum (use to compare mono versus
622   nictle     =    0       !  end   i indice of control sum        multi processor runs
623   njctls     =    0       !  start j indice of control               over a subdomain)
624   njctle     =    0       !  end   j indice of control
625   isplt      =    1       !  number of processors in i-direction
626   jsplt      =    1       !  number of processors in j-direction
627   nbench     =    0       !  Bench mode (1/0): CAUTION use zero except for bench
628                           !     (no physical validity of the results)
629   nbit_cmp   =    0       !  bit comparison mode (1/0): CAUTION use zero except for test
630                           !     of comparison between single and multiple processor runs
631/
632!!======================================================================
633!!                       ***  Diagnostics namelists  ***
634!!======================================================================
635!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld")
636!!   namgap       level mean model-data gap                             ("key_diagap")
637!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr")
638!!   namflo       float parameters                                      ("key_float")
639!!   namptr       Poleward Transport Diagnostics
640!!======================================================================
641
642!-----------------------------------------------------------------------
643&namtrd       !   diagnostics on dynamics and/or tracer trends             ("key_trddyn" and/or "key_trdtra")
644!                          or mixed-layer trends                        ('key_trdmld')
645!                          or barotropic vorticity                      ("key_trdvor")
646!-----------------------------------------------------------------------
647   ntrd       = 365        !  time step frequency dynamics and tracers trends
648   nctls      =   0        !  control surface type in mixed-layer trends (0,1 or n<jpk)
649   ucf        =   1.       !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
650   cn_trdrst_in  = "restart_mld"   ! suffix of ocean restart name (input)
651   cn_trdrst_out = "restart_mld"   ! suffix of ocean restart name (output)
652   ln_trdmld_restart = .false.     !  restart for ML diagnostics
653   ln_trdmld_instant = .false.     !  flag to diagnose trends of instantantaneous or mean ML T/S
654/
655!-----------------------------------------------------------------------
656&namgap       !   level mean model-data gap                             ('key_diagap')
657!-----------------------------------------------------------------------
658   ngap       =  15        !  time-step frequency of model-data gap computation
659   nprg       =  10        !  time-step frequency of gap print in model output
660/
661!-----------------------------------------------------------------------
662&namspr       !   surface pressure diagnostic
663!-----------------------------------------------------------------------
664   nmaxp      =   1000     !  maximum of iterations for the solver
665   epsp       =  1.e-3     !  absolute precision of the solver
666   niterp     =    400     !  number of iteration done by the solver
667/
668!-----------------------------------------------------------------------
669&namflo       !   float parameters                                      ("key_float")
670!-----------------------------------------------------------------------
671    ln_rstflo = .false.    !  float restart (T) or not (F)
672    nwritefl  =      75    !  frequency of writing in float output file
673    nstockfl  =    5475    !  frequency of creation of the float restart file
674    ln_argo   = .false.    !  Argo type floats (stay at the surface each 10 days)
675    ln_flork4 = .false.    !  trajectories computed with a 4th order Runge-Kutta (T)
676                           !  or computed with Blanke' scheme (F)
677/
678!-----------------------------------------------------------------------
679&namptr       !   Poleward Transport Diagnostic
680!-----------------------------------------------------------------------
681   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F)
682   ln_subbas  = .false.    !  Atlantic/Pacific/Indian basins computation (T) or not
683                           !  (orca configuration only, need input basins mask file named "subbasins.nc"
684   nf_ptr     =  15        !  Frequency of ptr computation [time step]
685/
Note: See TracBrowser for help on using the repository browser.