Index: trunk/LMDZ.VENUS/deftank/context_lmdz_physics-78niv.xml
===================================================================
--- trunk/LMDZ.VENUS/deftank/context_lmdz_physics-78niv.xml (revision 2464)
+++ trunk/LMDZ.VENUS/deftank/context_lmdz_physics-78niv.xml (revision 2465)
@@ -290,4 +290,7 @@
long_name="N2"
unit="mol/mol" />
+
+ long_name="Freq Brunt Väisälä"
+ unit="Hz^2" />
+
+
Index: trunk/LMDZ.VENUS/deftank/context_lmdz_physics-notrac.xml
===================================================================
--- trunk/LMDZ.VENUS/deftank/context_lmdz_physics-notrac.xml (revision 2464)
+++ trunk/LMDZ.VENUS/deftank/context_lmdz_physics-notrac.xml (revision 2465)
@@ -290,4 +290,7 @@
long_name="N2"
unit="mol/mol" />
+
+ long_name="Freq Brunt Väisälä"
+ unit="Hz^2" />
+
+
Index: trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml
===================================================================
--- trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml (revision 2464)
+++ trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml (revision 2465)
@@ -290,4 +290,7 @@
long_name="N2"
unit="mol/mol" />
+
+ long_name="Freq Brunt Väisälä"
+ unit="Hz^2" />
+
+
Index: trunk/LMDZ.VENUS/deftank/gcm-96x96x90-chemistry.def
===================================================================
--- trunk/LMDZ.VENUS/deftank/gcm-96x96x90-chemistry.def (revision 2465)
+++ trunk/LMDZ.VENUS/deftank/gcm-96x96x90-chemistry.def (revision 2465)
@@ -0,0 +1,110 @@
+## $Header$
+#
+## Planet:
+planet_type=venus
+#
+## Number of dynamical steps per day (must be a multiple of iperiod)
+day_step=480000
+## Apply a Matsuno step every iperiod dynamical step
+iperiod=5
+## dissipation is applied every dissip_period dynamical steps
+## DEFAULT: dissip_period=0 , meaning dissip_period is automatically computed
+## (in practice it is =>25 in recent runs)
+dissip_period=5
+## dissipation operator to use (star or non-star)
+lstardis=y
+## iterate lateral dissipation operator gradiv nitergdiv times
+nitergdiv=1
+## iterate lateral dissipation operator nxgradrot nitergrot times
+nitergrot=2
+## iterate lateral dissipation operator divgrad niterh times
+niterh=2
+## dissipation time scale (s) for shortest wavelengths for u,v (gradiv)
+tetagdiv=3000
+## dissipation time scale (s) for shortest wavelengths for u,v (nxgradrot)
+tetagrot=4000
+## dissipation time scale (s) for shortest wavelengths for u,v (divgrad)
+tetatemp=4000
+## coefficient for gamdissip
+coefdis=0.
+## time marching scheme (Matsuno if purmats is y, else Matsuno-Leapfrog)
+purmats=n
+# run with (true) or without (false) physics
+physic=y
+## Physics package type
+## 0: no physics (e.g. Shallow Water mode)
+## 1: with physics (e.g. phyvenus physics package)
+## 2: with a netwonian relaxation scheme in the dynamics
+iflag_phys=1
+## run with or without initial condition files (start.nc, startphy.nc) ?
+## (in the without case, initialization of fields is done via the iniacademic
+## routine in the dynamics => not available for Venus
+read_start=y
+## call physics every iphysiq dynamical steps
+iphysiq=10
+## runwith or without tracers
+iflag_trac=0
+## run with or without stratosphere // i.e. a sponge layer and secondary
+## higher altitude level of horizontal dissipation
+ok_strato=y
+## Horizontal dissipation multipliers along the vertical
+dissip_fac_mid=2.
+dissip_fac_up=30.
+# deltaz et hdelta in km
+dissip_deltaz=20.
+dissip_hdelta=10.
+# pupstart in Pa
+dissip_pupstart=1.
+
+## Sponge layer
+# 0: LMDZ.GENERIC style
+# 1 and 2: LMDZ.EARTH style
+# 1: in last 4 levels
+# 2: in levels with pressure less than 100 times the last layer pressure
+iflag_top_bound=0
+
+## sponge layer parameters LMDZ.EARTH style
+## Mode
+# mode = 0 : no sponge
+# mode = 1 : u et v -> 0
+# mode = 2 : u et v -> zonal average
+# mode = 3 : u, v et h -> zonale average
+mode_top_bound=3
+# Coefficient for the sponge layer (value in topmost layer)
+tau_top_bound=1.e-4
+
+# sponge layer parameters LMDZ.GENERIC style
+callsponge=y
+# mode_sponge ( 0: h -> h_mean , ucov -> 0 , vcov -> 0
+# 1: h -> h_mean , ucov -> ucov_mean , vcov -> 0
+# 2: h -> h_mean , ucov -> ucov_mean , vcov -> vcov_mean )
+mode_sponge=2
+# nsponge: number of topmost atmospheric layers over which extends the sponge
+nsponge=10
+# tetasponge characteristic time scale (seconds) at topmost layer
+# (time scale then doubles with decreasing layer index)
+tetasponge=1.e4
+
+###############################################
+### Zoom parameters
+###############################################
+## longitude (degrees) of zoom center
+clon=0.
+## latitude (degrees) of zoom center
+clat=0.
+## enhancement factor of zoom, along longitudes
+grossismx=1.0
+## enhancement factor of zoom, along latitudes
+grossismy=1.0
+## Use an hyperbolic function f(y) if .true., else use a sine
+fxyhypb=y
+## extention along longitudes of zoom region (fraction of global domain)
+dzoomx=0.0
+## extention along latitudes of zoom region (fraction of global domain)
+dzoomy=0.0
+## zoom stiffness along longitudes
+taux=3.
+## zoom stiffness along latitudes
+tauy=3.
+## Fonction f(y) as y = Sin(latitude) if = .true. , else y = latitude
+ysinus=y
Index: trunk/LMDZ.VENUS/deftank/physiq-96x96x78-chemistry.def
===================================================================
--- trunk/LMDZ.VENUS/deftank/physiq-96x96x78-chemistry.def (revision 2464)
+++ trunk/LMDZ.VENUS/deftank/physiq-96x96x78-chemistry.def (revision 2465)
@@ -94,11 +94,11 @@
solvarmod=0
-## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
-solarcondate = 1993.4
+## (Solar min=80 ave=140 max=300)
+fixed_euv_value = 140.
# value for the UV heating efficiency
##(experimental values between 0.19 and 0.23, lower values may
## be used to compensate for low 15 um cooling)
-euveff = 0.23
+euveff = 0.21
#
#
Index: trunk/LMDZ.VENUS/deftank/physiq-96x96x90-chemistry.def
===================================================================
--- trunk/LMDZ.VENUS/deftank/physiq-96x96x90-chemistry.def (revision 2465)
+++ trunk/LMDZ.VENUS/deftank/physiq-96x96x90-chemistry.def (revision 2465)
@@ -0,0 +1,136 @@
+## $Header$
+#
+## Include diurnal cycle or not
+cycle_diurne=y
+## Include soil model or not
+soil_model=y
+## Use orodr or not for orography
+ok_orodr=n
+## Use orolf or not for orography
+ok_orolf=n
+## Use non-orographic Gravity Waves of not
+ok_gw_nonoro=y
+## Number of calls to the radiative transfer (per day)
+nbapp_rad=24000
+## Number of calls to the chemistry routines (per day)
+nbapp_chim=1
+## Flag for convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(new version JYG), 30 KE(version IPCC AR4), 4 KE vect
+iflag_con=0
+#
+# orbital parameters
+##~~~~~~~~~~~~~~~~~~
+#
+# VENUS ##
+# R_ecc = 0.006787
+R_ecc = 0.
+R_peri = 0.
+# R_incl = 177.4
+R_incl = 0.
+# solar: effective flux, given at 1 UA
+solaire = 328.
+#
+# parameters for the tracers
+##~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+# 0: nothing (passive tracers)
+# 1: pseudo-chemistry relaxation (phytrac_relax)
+# 2: surface emission (phytrac_emiss)
+# 3: full chemistry (phytrac_chem)
+tr_scheme = 3
+#
+# Reinitialization of tracer abundances
+reinit_trac=n
+#
+# parameters for chemistry and microphysics
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+# use chemistry ?
+ok_chem=y
+# use clouds ? (needed for chemistry)
+ok_cloud=y
+# if yes, which scheme:
+# 1 => simple scheme (Aurelien Stolzenbach)
+# 2 => full microphysical scheme (Sabrina Guilbon)
+cl_scheme=1
+# use sedimentation (goes with clouds)
+ok_sedim=y
+nb_mode=3
+#
+# parameters for the boundary layer
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+iflag_pbl=8
+z0 = 0.01
+lmixmin = 35.
+ksta = 1.e-7
+ok_kzmin=n
+#
+# Surface thermal inertia
+inertie=2000
+#
+# dry convection parameters
+##~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+iflag_ajs = 1
+#
+# Thermosphere and nlte parameters
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+## Compute non-LTE or not
+callnlte=y
+## choice of species profile to use
+nltemodel=2
+## Compute CO2 IR absorption or not
+callnirco2=y
+##CO2 IR Absorption model
+nircorr=1
+##include thermosphere or not
+callthermos=y
+
+## Thermospheric options
+##~~~~~~~~~~~~~~~~~~~~~~
+#Method to include solar variability?
+#0-> Old method 1-> Variability with E10.7 as observed
+solvarmod=0
+
+## (Solar min=80 ave=140 max=300)
+fixed_euv_value = 140.
+
+# value for the UV heating efficiency
+##(experimental values between 0.19 and 0.23, lower values may
+## be used to compensate for low 15 um cooling)
+euveff = 0.21
+#
+#
+#
+#
+# Parameters for IOIPSL output files
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+## OLD. Now we use XIOS => see context_lmdz_physics.xml to taylor the output files
+#
+### OK_journe= y for daily output file histday.nc, =n no histday.nc output
+### Meaningless for Venus
+OK_journe=n
+### OK_mensuel= y for monthly output file histmth.nc, =n no histmth.nc
+### For Venus, only these averaged outputs
+OK_mensuel=n
+## rate (in days) at which the Venus histmth file is to be written
+# sets the output rate in histmth and/or histins
+ecritphy=0.1
+### OK_instan=y, make some "instantaneous" outputs (same rate as histmth)
+OK_instan=n
+#
+# Output levels for the various output files
+#
+# output level for "day" lev_histday
+# - lev_hist*=1 => baseline 2D fields
+# - lev_hist*=2 => baseline 3D fields (default)
+# - lev_hist*=3 => radiative transfert
+# - lev_hist*=4 => 3D tendencies
+# - lev_hist*=5 => tracers and others
+lev_histday=2
+#output level for "mth" lev_histmth
+lev_histmth=2
+#output level for "ins" lev_histins
+lev_histins=2
+
Index: trunk/LMDZ.VENUS/deftank/traceur-chemistry-withHe.def
===================================================================
--- trunk/LMDZ.VENUS/deftank/traceur-chemistry-withHe.def (revision 2465)
+++ trunk/LMDZ.VENUS/deftank/traceur-chemistry-withHe.def (revision 2465)
@@ -0,0 +1,37 @@
+36
+10 10 co2
+10 10 co
+10 10 h2
+10 10 h2o
+10 10 o1d
+10 10 o
+10 10 o2
+10 10 o2dg
+10 10 o3
+10 10 h
+10 10 oh
+10 10 ho2
+10 10 h2o2
+10 10 cl
+10 10 clo
+10 10 cl2
+10 10 hcl
+10 10 hocl
+10 10 clco
+10 10 clco3
+10 10 cocl2
+10 10 s
+10 10 so
+10 10 so2
+10 10 so3
+10 10 s2o2
+10 10 ocs
+10 10 hso3
+10 10 h2so4
+10 10 s2
+10 10 clso2
+10 10 oscl
+10 10 n2
+10 10 he
+10 10 h2oliq
+10 10 h2so4liq
Index: trunk/LMDZ.VENUS/deftank/z2sig-venus90niv.def
===================================================================
--- trunk/LMDZ.VENUS/deftank/z2sig-venus90niv.def (revision 2465)
+++ trunk/LMDZ.VENUS/deftank/z2sig-venus90niv.def (revision 2465)
@@ -0,0 +1,101 @@
+ 1.00000 H: atmospheric scale height (km) (used as a reference only)
+ 0.000500016264 Typical log(sigma) for 1st layer (z=H*log(sigma))
+ 0.00198207167 Typical log(sigma) for 2nd layer, etc...
+ 0.00790814403
+ 0.0212424006
+ 0.0449459739
+ 0.0819828957
+ 0.135316208
+ 0.207908899
+ 0.302724212
+ 0.422724098
+ 0.570872605
+ 0.750131905
+ 0.963465035
+ 1.213835
+ 1.50420499
+ 1.837538
+ 2.17087102
+ 2.50420403
+ 2.83753705
+ 3.17087007
+ 3.50420308
+ 3.8375361
+ 4.17086935
+ 4.50420237
+ 4.83753538
+ 5.1708684
+ 5.50420141
+ 5.83753443
+ 6.17086744
+ 6.50420046
+ 6.83753347
+ 7.17086649
+ 7.5041995
+ 7.83753252
+ 8.17086601
+ 8.50419903
+ 8.837533
+ 9.17086601
+ 9.50419998
+ 9.837533
+ 10.170867
+ 10.5042
+ 10.837534
+ 11.170867
+ 11.5042009
+ 11.837534
+ 12.4375343
+ 13.0375338
+ 13.8875341
+ 14.7375336
+15.22908758
+15.72908758
+16.22908758
+16.72908758
+17.22908758
+17.72908758
+18.22908758
+18.72908758
+19.22908758
+19.72908758
+20.22908758
+20.72908758
+21.22908758
+21.72908758
+22.22908758
+22.72908758
+23.22908758
+23.72908758
+24.22908758
+24.72908758
+25.22908758
+25.72908758
+26.22908758
+26.72908758
+27.22908758
+27.72908758
+28.22908758
+28.72908758
+29.22908758
+29.72880695
+30.22915096
+30.72949498
+31.22938577
+31.72927656
+32.22899593
+32.72871529
+33.22856881
+33.72842233
+34.22886155
+34.72886155
+
+
+
+
+
+
+
+
+
+