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 NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/eORCA1/EXPREF – NEMO

source: NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/eORCA1/EXPREF/namelist_cfg @ 12224

Last change on this file since 12224 was 12224, checked in by agn, 4 years ago

add eORCA1 configuration with EXPREF holding OSMOSIS xml & namelists
File size: 38.2 KB
Line 
1!-----------------------------------------------------------------------
2&namrun        !   parameters of the run
3!-----------------------------------------------------------------------
4   nn_no       =       0   !  job number (no more used...)
5   cn_exp      =  "eORCA1L3_OSM_N"  !  experience name
6   nn_it000    =       1 !8993 !33 !8992   !  first time step
7   nn_itend    =     35040 !72 !35040 !35040 !6 4 !11680 !32 !2464 !35040  !64 !9000 !116800  !  last  time step (std 5475)
8   nn_date0    =  19800101 !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
9   nn_leapy    =       0   !  Leap year calendar (1) or not (0)
10   ln_rstart   = .false. !.false.   !  start from rest (F) or from a restart file (T)
11   !
12   nn_rstctl   =   2 !0 ! 2            !  restart control ==> activated only if ln_rstart=T
13   cn_ocerst_in    =  "restart"  ! "eORCA1L3_OSM_no_w_00000032_restart" !  suffix of ocean restart name (input)
14   cn_ocerst_indir = "RESTARTS"         !  directory from which to read input ocean restarts
15   cn_ocerst_outdir = "RESTARTS"         !  directory from which to read input ocean restarts
16   nn_istate   =      0 !1   !  output the initial state (1) or not (0)
17   ln_rst_list = .false.   !  output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
18   nn_stock    =    11680 !8892 ! 32 !8892 !11680   !  frequency of creation of a restart file (modulo referenced to 1)
19   nn_stocklist = 8892,11680,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
20   nn_write    =    11680   !  frequency of write in the output file   (modulo referenced to nn_it000)
21   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%)
22   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard
23   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
24/
25!-----------------------------------------------------------------------
26&namcfg        !   parameters of the configuration                      (default: user defined GYRE)
27!-----------------------------------------------------------------------
28   ln_read_cfg = .true.    !  (=T) read the domain configuration file
29      !                    !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules
30      cn_domcfg = "eORCA_R1_zps_domcfg"    ! domain configuration filename
31/
32!-----------------------------------------------------------------------
33&namdom        !   time and space domain
34!-----------------------------------------------------------------------
35   rn_rdt      =  2700.     !  time step for the dynamics
36   ln_meshmask = .false. !.true.    ! create meshfile 
37/
38!-----------------------------------------------------------------------
39&namtsd        !   data : Temperature  & Salinity
40!-----------------------------------------------------------------------
41   !                       ! =T  read T-S fields for:
42   ln_tsd_init = .true.         !  ocean initialisation
43   ln_tsd_dmp  = .false.         !  T-S restoring   (see namtra_dmp)
44   
45   cn_dir      = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/'      !  root directory for the T-S data location
46   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
47   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
48   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
49   sn_tem  = 'conservative_temperature_WOA13_decav_Reg1L75_clim', -1 ,'votemper' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc'  ,   ''    ,    ''
50   sn_sal  = 'absolute_salinity_WOA13_decav_Reg1L75_clim'       , -1 ,'vosaline' , .true. , .true. , 'yearly' , 'weights_3D_WOA13d1_2_eorca1_bilinear.nc'  ,   ''    ,    ''
51/
52!-----------------------------------------------------------------------
53&namwad        !   Wetting and drying  default is no WAD
54!-----------------------------------------------------------------------
55/
56!-----------------------------------------------------------------------
57&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T)
58!-----------------------------------------------------------------------
59/
60!-----------------------------------------------------------------------
61&namc1d        !   1D configuration options                             ("key_c1d")
62!-----------------------------------------------------------------------
63/
64!-----------------------------------------------------------------------
65&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d")
66!-----------------------------------------------------------------------
67/
68!-----------------------------------------------------------------------
69&namc1d_uvd    !   data: U & V currents                                 ("key_c1d")
70!-----------------------------------------------------------------------
71/
72!-----------------------------------------------------------------------
73&namsbc        !   Surface Boundary Condition (surface module)
74!-----------------------------------------------------------------------
75   nn_fsbc     = 1         !  frequency of surface boundary condition computation
76                           !     (also = the frequency of sea-ice & iceberg model call)
77                     ! Type of air-sea fluxes
78   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
79   nn_ice      = 2         !  =0 no ice boundary condition   
80                           !  =1 use observed ice-cover                 (  => fill namsbc_iif )
81                           !  =2 or 3 automatically for LIM3 or CICE    ("key_lim3" or "key_cice")
82                           !          except in AGRIF zoom where it has to be specified
83   ln_ice_embd = .false.   !  =T embedded sea-ice (pressure + mass and salt exchanges)
84                           !  =F levitating ice (no pressure, mass and salt exchanges)
85                     ! Misc. options of sbc :
86   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr)
87   ln_dm2dc    = .true.   !  daily mean to diurnal cycle on short wave
88   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf)
89   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr)
90   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked
91                           !     =1 global mean of e-p-r set to zero at each time step
92                           !     =2 annual global mean of e-p-r set to zero
93   ln_isf      = .true.   !  ice shelf                                 (T   => fill namsbc_isf)
94   ln_wave     = .true.   !  Activate coupling with wave  (T => fill namsbc_wave)
95   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)
96   ln_sdw      = .true. !.false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)
97   nn_sdrift   =  1        !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
98                           !   = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
99                           !   = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
100                           !   = 2 Phillips as (1) but using the wave frequency from a wave model
101   ln_tauwoc    = .false.   !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)
102   ln_stcor    = .true.   !  Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
103   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
104                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
105/
106!-----------------------------------------------------------------------
107&namsbc_flx    !   surface boundary condition : flux formulation
108!-----------------------------------------------------------------------
109/
110!-----------------------------------------------------------------------
111&namsbc_blk   !   namsbc_blk  generic Bulk formula                      (ln_blk =T)
112!-----------------------------------------------------------------------
113!              !  file name                   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                              ! rotation ! land/sea mask !
114!              !                              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                             ! pairing  ! filename      !
115   sn_wndi     = 'u_10_core_oonly'       ,         6         , 'U_10_MOD',   .true.     ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bicubic.nc'  , 'U1'    , ''
116   sn_wndj     = 'v_10_core_oonly'       ,         6         , 'V_10_MOD',   .true.     ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bicubic.nc'  , 'V1'    , ''
117   sn_qsr      = 'rad_core_oonly'        ,        24         , 'SWDN_MOD',   .false.    ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bilinear.nc' , ''      , ''
118   sn_qlw      = 'rad_core_oonly'        ,        24         , 'LWDN_MOD',   .true.     ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bilinear.nc' , ''      , ''
119   sn_tair     = 't_10_core_oonly'       ,         6         , 'T_10_MOD',   .false.    ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bilinear.nc' , ''      , ''
120   sn_humi     = 'q_10_core_oonly'       ,         6         , 'Q_10_MOD',   .false.    ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bilinear.nc' , ''      , ''
121   sn_prec     = 'tprecip_core_oonly'    ,        -1         , 'TPRECIP' ,   .true.     ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bilinear.nc' , ''      , ''
122   sn_snow     = 'tprecip_core_oonly'    ,        -1         , 'SNOW'    ,   .true.     ,  .false. , 'yearly'  , 'weights_coreII_2_eorca1_bilinear.nc' , ''      , ''
123   sn_slp      = 'slp_core_oonly'        ,         6         , 'SLP'     ,   .true.     , .false.  , 'yearly'  , 'weights_coreII_2_eorca1_bilinear.nc' , ''      , ''
124   ln_NCAR     = .true.   ! "NCAR"      algorithm   (Large and Yeager 2008)
125   cn_dir      = '/hpcdata/scratch/agn/OSMOSIS/FORCING/IAF_CORE/'      !  root directory for the location of the flux files
126/
127!-----------------------------------------------------------------------
128&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
129!-----------------------------------------------------------------------
130/
131!-----------------------------------------------------------------------
132&namsbc_sas    !   Stand-Alone Surface boundary condition
133!-----------------------------------------------------------------------
134/
135!-----------------------------------------------------------------------
136&namsbc_iif    !   Ice-IF : use observed ice cover                      (nn_ice = 1)
137!-----------------------------------------------------------------------
138/
139!-----------------------------------------------------------------------
140&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T)
141!-----------------------------------------------------------------------
142!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
143!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
144   sn_chl      ='merged_ESACCI_BIOMER4V1R1_CHL_REG05',  -1  , 'CHLA' , .true. , .true. , 'yearly' , 'weights_reg05_2_eorca1_bilinear.nc' , '' , ''
145   cn_dir      = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/'      !  root directory for the location of the flux files
146   ln_qsr_rgb  = .false.   !  RGB (Red-Green-Blue) light penetration
147   ln_qsr_2bd  = .true.  !  2 bands              light penetration
148   ln_qsr_bio  = .false.  !  bio-model light penetration
149   nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0)
150   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
151   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
152   rn_si1      =   23.0    !  2 bands: longest depth of extinction
153/
154!-----------------------------------------------------------------------
155&namsbc_rnf    !   runoffs namelist surface boundary condition          (ln_rnf =T)
156!-----------------------------------------------------------------------
157!              !  file name                                       ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask!
158!              !                                                  !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename     !
159   sn_rnf      = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , ''
160   sn_cnf      = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc',         0         , 'socoeff' ,   .false.    , .true. , 'yearly'  , ''       , ''       , ''
161
162   cn_dir      = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/'      !  root directory for the location of the flux files
163   ln_rnf_mouth = .false.    !  specific treatment at rivers mouths
164   ln_rnf_depth = .false.    !  read in depth information for runoff
165   ln_rnf_tem   = .false.    !  read in temperature information for runoff
166   ln_rnf_sal   = .false.    !  read in salinity information for runoff
167   ln_rnf_depth_ini = .false. !.true. !  compute depth at initialisation from runoff file
168   rn_rnf_max   = 0.05       !  max value of the runoff climatology over global domain ( if ln_rnf_depth_ini = .true )
169   rn_dep_max = 150.         !  depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
170   nn_rnf_depth_file = 0     ! create (=1) a runoff depth file or not (=0)
171/
172!-----------------------------------------------------------------------
173&namsbc_isf    !  Top boundary layer (ISF)                              (nn_isf >0)
174!-----------------------------------------------------------------------
175!              ! file name                                          ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask!
176!              !                                                    !  (if <0  months)  !   name   !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename     !
177   sn_rnfisf     = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' ,   -12             ,'sornfisf'     , .false. , .true. , 'yearly'  ,  ''     ,   ''     ,  ''
178   sn_depmax_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' ,   -12             ,'sodepmax_isf' , .false. , .true. , 'yearly'  ,  ''     ,   ''     ,  ''
179   sn_depmin_isf = 'runoff-icb_DaiTrenberth_Depoorter_eORCA1_JD.nc' ,   -12             ,'sodepmin_isf' , .false. , .true. , 'yearly'  ,  ''     ,   ''     ,  ''
180   nn_isf      = 3         !  ice shelf melting/freezing
181                           !  1 = presence of ISF    2 = bg03 parametrisation
182                           !  3 = rnf file for isf   4 = ISF fwf specified
183                           !  option 1 and 4 need ln_isfcav = .true. (domzgr)
184/
185!-----------------------------------------------------------------------
186&namsbc_iscpl  !   land ice / ocean coupling option                     
187!-----------------------------------------------------------------------
188/
189!-----------------------------------------------------------------------
190&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T)
191!-----------------------------------------------------------------------
192/
193!-----------------------------------------------------------------------
194&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T)
195!-----------------------------------------------------------------------
196!              ! file name ! frequency (hours) ! variable ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
197!              !           !  (if <0  months)  !   name   !   (logical) !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
198   sn_sss      = 'sss_absolute_salinity_WOA13_decav_Reg1L75_clim', -1. , 'sosaline', .true. , .true. , 'yearly' , 'weights_WOA13d1_2_eorca1_bilinear.nc' , '' , ''
199   cn_dir      = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/'      !  root directory for the location of the flux files
200      rn_deds     =  -166.67  !  magnitude of the damping on salinity   [mm/day]
201      ln_sssr_bnd =  .true.   !  flag to bound erp term (associated with nn_sssr=2)
202      rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day]
203/
204!-----------------------------------------------------------------------
205&namsbc_wave   ! External fields from wave model                        (ln_wave=T)
206!-----------------------------------------------------------------------
207!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
208!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
209!   sn_cdg      =  'sdw_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'yearly'   ,  ''      , ''       , ''
210   sn_usd      =  'ust' ,        6          , 'ust'     ,     .true.   , .false. , 'monthly'   ,  'weights_grid05_eORCA1_bicubic.nc'    , 'U1'     , ''
211   sn_vsd      =  'vst' ,        6          , 'vst'     ,     .true.   , .false. , 'monthly'   ,  'weights_grid05_eORCA1_bicubic.nc'    , 'V1'     , ''
212   sn_hsw      =  'swh' ,        6          , 'swh'     ,     .true.   , .false. , 'monthly'   ,  'weights_grid05_eORCA1_bilin.nc'      , ''       , ''
213   sn_wmp      =  'mwp' ,        6          , 'mwp'     ,     .true.   , .false. , 'monthly'   ,  'weights_grid05_eORCA1_bilin.nc'      , ''       , ''
214!   sn_wnum     =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'yearly'   ,  ''      , ''       , ''
215!   sn_tauoc    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'yearly'   ,  ''      , ''       , ''
216!
217   cn_dir      = '/hpcdata/scratch/agn/OSMOSIS/FORCING/ERA_INTERIM/WAVES/'      !  root directory for the location of the flux files
218/
219!-----------------------------------------------------------------------
220&namberg       !   iceberg parameters                                   (default: No iceberg)
221!-----------------------------------------------------------------------
222/
223!-----------------------------------------------------------------------
224&namlbc        !   lateral momentum boundary condition
225!-----------------------------------------------------------------------
226   !                       !  free slip  !   partial slip  !   no slip   ! strong slip
227   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
228   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs.
229/
230!-----------------------------------------------------------------------
231&namagrif      !  AGRIF zoom                                            ("key_agrif")
232!-----------------------------------------------------------------------
233/
234!-----------------------------------------------------------------------
235&nam_tide      !   tide parameters
236!-----------------------------------------------------------------------
237/
238!-----------------------------------------------------------------------
239&nambdy        !  unstructured open boundaries
240!-----------------------------------------------------------------------
241/
242!-----------------------------------------------------------------------
243&nambdy_dta    !  open boundaries - external data
244!-----------------------------------------------------------------------
245/
246!-----------------------------------------------------------------------
247&nambdy_tide   !  tidal forcing at open boundaries
248!-----------------------------------------------------------------------
249/
250!-----------------------------------------------------------------------
251&namdrg        !   top/bottom drag coefficient                          (default: NO selection)
252!-----------------------------------------------------------------------
253   ln_lin = .true.    !   linear  drag: Cd = Cd0 Uc0
254/
255!-----------------------------------------------------------------------
256&namdrg_top    !   TOP friction                                         (ln_isfcav=T)
257!-----------------------------------------------------------------------
258/
259!-----------------------------------------------------------------------
260&namdrg_bot        !   BOTTOM friction
261!-----------------------------------------------------------------------
262/
263!-----------------------------------------------------------------------
264&nambbc        !   bottom temperature boundary condition                (default: NO)
265!-----------------------------------------------------------------------
266!              !  file name      ! frequency (hours) ! variable  ! time interp.!  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
267!              !                 !  (if <0  months)  !   name    !  (logical)  !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      !
268   sn_qgh      ='Goutorbe_ghflux.nc',  -12.  , 'gh_flux'    ,   .false.     , .true. , 'yearly'  , 'weights_Goutorbe1_2_eorca1_bilinear.nc'       , ''       , ''
269   !
270   ln_trabbc   = .true.   !  Apply a geothermal heating at the ocean bottom
271      nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux
272      !                       !     = 1 constant flux
273      !                       !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)
274      rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2]
275   cn_dir      = '/hpcdata/scratch/agn/OSMOSIS/FORCING/SHACONEMO_eORCA_LIM3/'      !  root directory for the location of the flux files
276/
277!-----------------------------------------------------------------------
278&nambbl        !   bottom boundary layer scheme                         (default: NO)
279!-----------------------------------------------------------------------
280/
281!-----------------------------------------------------------------------
282&nameos        !   ocean Equation Of Seawater                           (default: NO)
283!-----------------------------------------------------------------------
284   ln_teos10    = .true.         !  = Use TEOS-10 equation of state
285/
286!-----------------------------------------------------------------------
287&namtra_adv    !   advection scheme for tracer                          (default: NO selection)
288!-----------------------------------------------------------------------
289   ln_traadv_fct = .true. !  FCT scheme
290      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order
291      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order
292/
293!-----------------------------------------------------------------------
294&namtra_mle !   mixed layer eddy parametrisation (Fox-Kemper param) (default: NO)
295!-----------------------------------------------------------------------
296   ln_mle      = .true.   ! (T) use the Mixed Layer Eddy (MLE) parameterisation
297   rn_ce       = 0.04      ! magnitude of the MLE (typical value: 0.06 to 0.08)
298
299/
300!-----------------------------------------------------------------------
301&namtra_ldf    !   lateral diffusion scheme for tracers
302!----------------------------------------------------------------------------------
303   !                       !  Operator type:
304   ln_traldf_lap   =  .true.   !    laplacian operator
305   ln_traldf_blp   =  .false.  !  bilaplacian operator
306   !
307   !                       !  Direction of action:
308   ln_traldf_iso   =  .true.  !  iso-neutral (standard operator)
309   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator)
310   !
311   !                       !  iso-neutral options:       
312   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators)
313   rn_slpmax       =   0.01    !  slope limit                      (both operators)
314   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only)
315   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only)
316   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only)
317   !
318   !                       !  Coefficients:
319   nn_aht_ijk_t    = 20        !  space/time variation of eddy coefficient:
320      !                             !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file
321      !                             !   =  0           constant
322      !                             !   = 10 F(k)      =ldf_c1d
323      !                             !   = 20 F(i,j)    =ldf_c2d
324      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
325      !                             !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d
326      !                             !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
327      !                        !  time invariant coefficients:  aht0 = 1/2  Ud*Ld   (lap case)
328      !                             !                           or   = 1/12 Ud*Ld^3 (blp case)
329/
330!-----------------------------------------------------------------------
331&namtra_ldfeiv !   eddy induced velocity param.                         (default: NO)
332!-----------------------------------------------------------------------
333    ln_ldfeiv     = .true. ! use eddy induced velocity parameterization
334      rn_aeiv_0     = 1000.   ! eddy induced velocity coefficient   [m2/s]
335      nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient
336      !                             !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file
337      !                             !   =  0           constant
338      !                             !   = 10 F(k)      =ldf_c1d
339      !                             !   = 20 F(i,j)    =ldf_c2d
340      !                             !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation
341      !                             !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d
342      ln_ldfeiv_dia =.true.  ! diagnose eiv stream function and velocities
343/
344!-----------------------------------------------------------------------
345&namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO)
346!-----------------------------------------------------------------------
347   ln_tradmp   =  .false.
348/
349!-----------------------------------------------------------------------
350&nam_vvl       !   vertical coordinate options                          (default: z-star)
351!-----------------------------------------------------------------------
352   ln_vvl_zstar  = .true.           !  z-star vertical coordinate
353/
354!-----------------------------------------------------------------------
355&namdyn_adv    !   formulation of the momentum advection                (default: NO selection)
356!-----------------------------------------------------------------------
357   ln_dynadv_vec = .true. !  vector form - 2nd centered scheme
358     nn_dynkeg     = 0        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction
359   ln_dynadv_cen2 = .false. !  flux form - 2nd order centered scheme
360/
361!-----------------------------------------------------------------------
362&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO)
363!-----------------------------------------------------------------------
364   ln_dynvor_een = .true.  !  energy & enstrophy scheme
365      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)
366/
367!-----------------------------------------------------------------------
368&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: NO selection)
369!-----------------------------------------------------------------------
370   ln_hpg_sco  = .true.   !  s-coordinate (standard jacobian formulation)
371/
372!-----------------------------------------------------------------------
373&namdyn_spg    !   surface pressure gradient                            (default: NO)
374!-----------------------------------------------------------------------
375   ln_dynspg_ts   = .true.   ! split-explicit free surface
376      ln_bt_fw      = .false.     ! Forward integration of barotropic Eqs.
377      ln_bt_av      = .true.     ! Time filtering of barotropic variables
378         nn_bt_flt     = 1          ! Time filter choice  = 0 None
379         !                          !                     = 1 Boxcar over   nn_baro sub-steps
380         !                          !                     = 2 Boxcar over 2*nn_baro  "    "
381      ln_bt_auto    = .true.     ! Number of sub-step defined from:
382         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed
383         nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds
384      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F)
385/
386!-----------------------------------------------------------------------
387&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO selection)
388!-----------------------------------------------------------------------
389   !                       !  Type of the operator :
390   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F
391   ln_dynldf_lap =  .true.     !    laplacian operator
392   !                       !  Direction of action  :
393   ln_dynldf_lev =  .true.     !  iso-level
394   !                       !  Coefficient
395   nn_ahm_ijk_t  =  30         !  space/time variation of eddy coef
396   !                                !  =-30  read in eddy_viscosity_3D.nc file
397   !                                !  =-20  read in eddy_viscosity_2D.nc file
398   !                                !  =  0  constant
399   !                                !  = 10  F(k)=c1d
400   !                                !  = 20  F(i,j)=F(grid spacing)=c2d
401   !                                !  = 30  F(i,j,k)=c2d*c1d
402   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity)
403   !                                !  = 32  F(i,j,k)=F(local gridscale and deformation rate)
404   !                           !  time invariant coefficients :  ahm = 1/2  Uv*Lv   (lap case)
405   !                                !                            or  = 1/12 Uv*Lv^3 (blp case)
406/
407!-----------------------------------------------------------------------
408&namzdf        !   vertical physics                                     (default: NO selection)
409!-----------------------------------------------------------------------
410   ln_zdfosm   = .true.      !  OSMOSIS BL closure                     (T =>   fill namzdf_osm)
411   ln_zdftke = .false.
412   !
413   !                       ! convection
414   ln_zdfevd   = .false.    !  enhanced vertical diffusion
415   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F)
416   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F)
417/
418!-----------------------------------------------------------------------
419&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T)
420!-----------------------------------------------------------------------
421/
422!-----------------------------------------------------------------------
423&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T)
424!-----------------------------------------------------------------------
425   !
426   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to NIWs
427                              !        = 0 none ; = 1 add a tke source below the ML
428                              !        = 2 add a tke source just at the base of the ML
429                              !        = 3 as = 1 applied on HF part of the stress           (ln_cpl=T)
430/
431!-----------------------------------------------------------------------
432&namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls =T)
433!-----------------------------------------------------------------------
434/
435!-----------------------------------------------------------------------
436&namzdf_osm    !   OSM vertical diffusion                               (ln_zdfosm =T)
437!-----------------------------------------------------------------------
438   ln_use_osm_la = .false.      !  Use namelist  rn_osm_la
439   rn_osm_la     = 0.3         !  Turbulent Langmuir number
440   rn_osm_dstokes     = 5.     !  Depth scale of Stokes drift (m)
441   nn_ave = 0                  ! choice of horizontal averaging on avt, avmu, avmv
442   ln_dia_osm = .true.         ! output OSMOSIS-OBL variables
443   rn_osm_hbl0 = 10.           ! initial hbl value
444   ln_kpprimix = .true.        ! Use KPP-style Ri# mixing below BL
445   rn_riinfty  = 0.7           ! Highest local Ri_g permitting shear instability
446   rn_difri  =  0.005          ! max Ri# diffusivity at Ri_g = 0 (m^2/s)
447   ln_convmix  = .true.        ! Use convective instability mixing below BL
448   rn_difconv = 1. !0.01 !1.             ! diffusivity when unstable below BL  (m2/s)
449   nn_osm_wave = 2             ! Method used to calculate Stokes drift
450                               !  = 2: Use ECMWF wave fields
451                               !  = 1: Pierson Moskowitz wave spectrum
452                               !  = 0: Constant La# = 0.3
453!========================only if key_FK_osm is set =================================
454   ln_osm_mle = .true.        !  Use integrated FK-OSM model
455/
456!-----------------------------------------------------------------------
457&namosm_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF)
458!-----------------------------------------------------------------------
459   rn_osm_mle_ce       = 0.04      ! magnitude of the MLE (typical value: 0.06 to 0.08)
460   nn_osm_mle      = 0         ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
461   rn_osm_mle_lf       = 5.e+3     ! typical scale of mixed layer front (meters)                      (case rn_osm_mle=0)
462   rn_osm_mle_time     = 43200 !172800.   ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_osm_mle=0)
463   rn_osm_mle_lat      = 20.       ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1)
464   rn_osm_mle_rho_c = 0.03 ! 0.01      ! delta rho criterion used to calculate MLD for FK
465   rn_osm_mle_thresh = 0.0001      ! delta b criterion used for FK MLD criterion
466   rn_osm_mle_tau     = 172800.   ! time scale for FK-OSM (seconds)  (case rn_osm_mle=0)
467!========================only if key_FK_osm is set =================================
468/
469!-----------------------------------------------------------------------
470&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T)
471!-----------------------------------------------------------------------
472/
473!-----------------------------------------------------------------------
474&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
475!-----------------------------------------------------------------------
476  ln_nnogather = .true.  !  activate code to avoid mpi_allgather use at the northfold
477  jpni         = 16 !12 !16 !12 !16 !12 !16      ! 12 jpni   number of processors following i (set automatically if < 1)
478  jpnj         = 15 !6 !15 !6 !15 !6 ! 15      ! 8 jpnj   number of processors following j (set automatically if < 1)
479/
480!-----------------------------------------------------------------------
481&namctl        !   Control prints
482!-----------------------------------------------------------------------
483   ln_ctl = .false.                 ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T
484/
485!-----------------------------------------------------------------------
486&namsto        ! Stochastic parametrization of EOS                      (default: NO)
487!-----------------------------------------------------------------------
488/
489!-----------------------------------------------------------------------
490&namtrd        !   trend diagnostics                                    (default F)
491!-----------------------------------------------------------------------
492/
493!-----------------------------------------------------------------------
494&namptr        !   Poleward Transport Diagnostic                        (default F)
495!-----------------------------------------------------------------------
496/
497!-----------------------------------------------------------------------
498&namhsb        !  Heat and salt budgets                                 (default F)
499!-----------------------------------------------------------------------
500/
501!-----------------------------------------------------------------------
502&namhth        !  Extra ML depth and thermocline depth diagnostics      (default F)
503!-----------------------------------------------------------------------
504   ln_diahth   = .true.   !  produce extra diagnostics of BL and thermocline depth
505/
506!-----------------------------------------------------------------------
507&namdiu        !   Cool skin and warm layer models                      (default F)
508!-----------------------------------------------------------------------
509/
510!-----------------------------------------------------------------------
511&namflo        !   float parameters                                     ("key_float")
512!-----------------------------------------------------------------------
513/
514!-----------------------------------------------------------------------
515&nam_diaharm   !   Harmonic analysis of tidal constituents              ("key_diaharm")
516!-----------------------------------------------------------------------
517/
518!-----------------------------------------------------------------------
519&namdct        ! transports through some sections                       ("key_diadct")
520!-----------------------------------------------------------------------
521/
522!-----------------------------------------------------------------------
523&nam_diatmb    !  Top Middle Bottom Output                              (default F)
524!-----------------------------------------------------------------------
525/
526!-----------------------------------------------------------------------
527&nam_dia25h    !  25h Mean Output                                       (default F)
528!-----------------------------------------------------------------------
529   ln_dia25h   = .false.   ! Choose 25h mean output or not
530/
Note: See TracBrowser for help on using the repository browser.