Index: /trunk/LMDZ.VENUS/deftank/context_lmdz_physics-78niv.xml
===================================================================
--- /trunk/LMDZ.VENUS/deftank/context_lmdz_physics-78niv.xml	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/context_lmdz_physics-78niv.xml	(revision 2279)
@@ -0,0 +1,619 @@
+<context id="LMDZ">
+
+    <!-- Definition of a calendar is done in code via wxios_set_cal -->
+    
+    <!-- Definition of domains and groups of domains -->
+    <domain_definition>
+        <domain id="dom_glo" data_dim="2" />
+    </domain_definition>
+        
+    <!-- Definition of vertical axes 
+         NB: these must be defined and set in code via
+         wxios_add_vaxis -->
+    <axis_definition>
+        <axis id="presnivs" 
+              standard_name="Pseudo-pressure of model vertical levels" 
+              unit="Pa">
+        </axis>
+        <axis id="altitude" 
+              standard_name="Pseudo-altitude of model vertical levels" 
+              unit="km">
+        </axis>
+    </axis_definition>
+    
+    <grid_definition>
+        <!-- toggle axis id below to change output vertical axis -->
+        <grid id="grid_3d">
+            <domain id="dom_glo" />
+                <axis id="presnivs" />
+                <!-- <axis id="altitude" /> -->
+        </grid>
+        <!-- to output scalars -->
+        <grid id="scalarpoint" >
+            <scalar />
+        </grid>
+    </grid_definition>
+    
+    <!-- Definition of model variables -->
+    <field_definition prec="4" 
+                      freq_op="1ts"
+                      enabled=".TRUE.">
+        
+        <!-- 0D variables 
+        <field_group id="fields_0D" grid_ref="scalarpoint">
+            <field id="ls"
+                   long_name="Solar longitude"
+                   unit="degrees" />
+        </field_group>    -->
+        
+        <!-- 2D variables -->
+        <field_group id="fields_2D" domain_ref="dom_glo">
+            <field id="phis"
+                   long_name="Surface geopotential (gz)" 
+                   unit="m2/s2" />
+            <field id="aire"
+                   long_name="Mesh area"
+                   unit="m2" />
+            <field id="tsol"
+                   long_name="Surface Temperature" 
+                   unit="K" />
+            <field id="psol"
+                   long_name="Surface Pressure" 
+                   unit="Pa" />
+            <field id="cdragh"
+                   long_name="Drag coef on T" 
+                   unit="-" />
+            <field id="cdragm"
+                   long_name="Drag coef on U" 
+                   unit="-" />
+            <field id="tops"
+                   long_name="Solar rad. at TOA" 
+                   unit="W/m2" />
+            <field id="topl"
+                   long_name="IR rad. at TOA" 
+                   unit="W/m2" />
+            <field id="sols"
+                   long_name="Solar rad. at surface" 
+                   unit="W/m2" />
+            <field id="soll"
+                   long_name="IR rad. at surface" 
+                   unit="W/m2" />
+        </field_group>
+
+        <!-- 3D variables -->
+        <field_group id="fields_3D" grid_ref="grid_3d">
+            <field id="temp"
+                   long_name="Atmospheric temperature"
+                   unit="K" />
+            <field id="pres"
+                   long_name="Atmospheric pressure" 
+                   unit="Pa" />
+            <field id="geop"
+                   long_name="Geopotential (gz)" 
+                   unit="m2/s2" />
+            <field id="vitu"
+                   long_name="Westward Zonal Wind"
+                   unit="m/s" />
+            <field id="vitv" 
+                   long_name="Northward Meridional Wind" 
+                   unit="m/s" />
+            <field id="vitw"
+                   long_name="Vertical wind" 
+                   unit="Pa/s" />
+            <field id="Kz"
+                   long_name="vertical diffusion coef" 
+                   unit="m2/s" />
+            <field id="mmean"
+                   long_name="Mean molecular mass" 
+                   unit="g/mol" />
+            <field id="rho"
+                   long_name="Atmospheric density" 
+                   unit="kg/m3" />
+            <field id="dudyn"
+                   long_name="Dynamics du" 
+                   unit="m/s2" />
+            <field id="duvdf"
+                   long_name="Boundary layer du" 
+                   unit="m/s2" />
+            <field id="dvvdf"
+                   long_name="Boundary layer dv" 
+                   unit="m/s2" />
+            <field id="duajs"
+                   long_name="Dry adj. du" 
+                   unit="m/s2" />
+            <field id="dugwo"
+                   long_name="GW oro du" 
+                   unit="m/s2" />
+            <field id="dugwno"
+                   long_name="GW non-oro du" 
+                   unit="m/s2" />
+            <field id="dumolvis"
+                   long_name="molec viscosity (u)" 
+                   unit="m/s2" />
+            <field id="dvmolvis"
+                   long_name="molec viscosity (v)" 
+                   unit="m/s2" />
+            <field id="dtdyn"
+                   long_name="Dynamics dT" 
+                   unit="K/s" />
+            <field id="dtphy"
+                   long_name="Physics dT" 
+                   unit="K/s" />
+            <field id="dtvdf"
+                   long_name="Boundary layer dT" 
+                   unit="K/s" />
+            <field id="dtajs"
+                   long_name="Dry adj. dT" 
+                   unit="K/s" />
+            <field id="dtswr"
+                   long_name="SW radiation dT" 
+                   unit="K/s" />
+            <field id="dtswrNLTE"
+                   long_name="SW radiation (NLTE) dT" 
+                   unit="K/s" />
+            <field id="dtswrLTE"
+                   long_name="SW radiation (LTE) dT" 
+                   unit="K/s" />
+            <field id="dtlwr"
+                   long_name="LW radiation dT" 
+                   unit="K/s" />
+            <field id="dtlwrNLTE"
+                   long_name="LW radiation (NLTE) dT" 
+                   unit="K/s" />
+            <field id="dtlwrLTE"
+                   long_name="LW radiation (LTE) dT" 
+                   unit="K/s" />
+            <field id="dteuv"
+                   long_name="EUV radiation dT" 
+                   unit="K/s" />
+            <field id="dtcond"
+                   long_name="Thermal conduction dT" 
+                   unit="K/s" />
+            <field id="dtec"
+                   long_name="Kinetic dissip dT" 
+                   unit="K/s" />
+            <field id="SWnet"
+                   long_name="Net SW flux" 
+                   unit="W/m2" />
+            <field id="LWnet"
+                   long_name="Net LW flux" 
+                   unit="W/m2" />
+            <field id="fluxvdf"
+                   long_name="PBL net flux" 
+                   unit="W/m2" />
+            <field id="fluxdyn"
+                   long_name="Dyn net flux" 
+                   unit="W/m2" />
+            <field id="fluxajs"
+                   long_name="Dry adj. net flux" 
+                   unit="W/m2" />
+            <field id="fluxec"
+                   long_name="Kinetic Energy net flux" 
+                   unit="W/m2" />
+            <field id="co2"
+                   long_name="CO2" 
+                   unit="mol/mol" />
+            <field id="co"
+                   long_name="CO" 
+                   unit="mol/mol" />
+            <field id="h2"
+                   long_name="H2" 
+                   unit="mol/mol" />
+            <field id="h2o"
+                   long_name="H2O" 
+                   unit="mol/mol" />
+            <field id="o1d"
+                   long_name="O(1D)" 
+                   unit="mol/mol" />
+            <field id="o"
+                   long_name="O" 
+                   unit="mol/mol" />
+            <field id="o2"
+                   long_name="O2" 
+                   unit="mol/mol" />
+            <field id="o2dg"
+                   long_name="O2(DG)" 
+                   unit="mol/mol" />
+            <field id="o3"
+                   long_name="O3" 
+                   unit="mol/mol" />
+            <field id="h"
+                   long_name="H" 
+                   unit="mol/mol" />
+            <field id="oh"
+                   long_name="OH" 
+                   unit="mol/mol" />
+            <field id="ho2"
+                   long_name="HO2" 
+                   unit="mol/mol" />
+            <field id="h2o2"
+                   long_name="H2O2" 
+                   unit="mol/mol" />
+            <field id="cl"
+                   long_name="Cl" 
+                   unit="mol/mol" />
+            <field id="clo"
+                   long_name="ClO" 
+                   unit="mol/mol" />
+            <field id="cl2"
+                   long_name="Cl2" 
+                   unit="mol/mol" />
+            <field id="hcl"
+                   long_name="HCl" 
+                   unit="mol/mol" />
+            <field id="hocl"
+                   long_name="HOCl" 
+                   unit="mol/mol" />
+            <field id="clco"
+                   long_name="ClCO" 
+                   unit="mol/mol" />
+            <field id="clco3"
+                   long_name="ClCO3" 
+                   unit="mol/mol" />
+            <field id="cocl2"
+                   long_name="COCl2" 
+                   unit="mol/mol" />
+            <field id="s"
+                   long_name="S" 
+                   unit="mol/mol" />
+            <field id="so"
+                   long_name="SO" 
+                   unit="mol/mol" />
+            <field id="so2"
+                   long_name="SO2" 
+                   unit="mol/mol" />
+            <field id="so3"
+                   long_name="SO3" 
+                   unit="mol/mol" />
+            <field id="s2o2"
+                   long_name="S2O2" 
+                   unit="mol/mol" />
+            <field id="ocs"
+                   long_name="OCS" 
+                   unit="mol/mol" />
+            <field id="hso3"
+                   long_name="HSO3" 
+                   unit="mol/mol" />
+            <field id="h2so4"
+                   long_name="H2SO4" 
+                   unit="mol/mol" />
+            <field id="s2"
+                   long_name="S2" 
+                   unit="mol/mol" />
+            <field id="clso2"
+                   long_name="ClSO2" 
+                   unit="mol/mol" />
+            <field id="oscl"
+                   long_name="OSCl" 
+                   unit="mol/mol" />
+            <field id="n2"
+                   long_name="N2" 
+                   unit="mol/mol" />
+            <field id="h2oliq"
+                   long_name="H2Oliq" 
+                   unit="mol/mol" />
+            <field id="h2so4liq"
+                   long_name="H2SO4liq" 
+                   unit="mol/mol" />
+            <field id="Fsedim"
+                   long_name="tendency from sedim" 
+                   unit="kg.m-2.s-1" />
+            <field id="d_qmoldifCO2"
+                   long_name="Diff molec CO2" 
+                   unit="kg/kg" />
+            <field id="d_qmoldifO3p"
+                   long_name="Diff molec 0" 
+                   unit="kg/kg" />
+            <field id="d_qmoldifN2"
+                   long_name="Diff molec N2" 
+                   unit="kg/kg" />
+            <field id="BV2"
+		           long_name="Freq Brunt Väisälä"
+		           unit="Hz^2" />
+       </field_group>
+
+    </field_definition>
+    
+    <!-- Definition of output files
+   output_freq: ts is physics time steps, d is day
+   enabled=".true." to output a file
+      -->
+
+    <file_definition>
+        <!-- Instantaneous outputs; every physics time steps -->
+        <file id="ins"
+              name="Xins"
+              output_freq="1ts" 
+              type="one_file"
+              enabled=".false.">
+                    
+            <!-- VARS 0D
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="ls" />
+            </field_group> -->
+
+            <!-- VARS 2D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="phis" operation="once" />
+                <field field_ref="aire" operation="once" />
+                <field field_ref="tops" />
+                <field field_ref="tsol" />
+                <field field_ref="psol" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="temp" />
+                <field field_ref="pres" />
+                <field field_ref="geop" />
+                <field field_ref="vitu" />
+                <field field_ref="vitv" />
+                <field field_ref="vitw" />
+                <field field_ref="dudyn" />
+                <field field_ref="duvdf" />
+                <field field_ref="duajs" />
+                <field field_ref="dtdyn" />
+                <field field_ref="dtvdf" />
+                <field field_ref="dtajs" />
+                <field field_ref="dteuv" />
+                <field field_ref="dtcond" />
+                <field field_ref="dtswrNLTE" />
+                <field field_ref="dtswrLTE" />
+                <field field_ref="dtlwrNLTE" />
+                <field field_ref="dtlwrLTE" />
+                <field field_ref="rho"   />
+                <field field_ref="mmean" />
+                <field field_ref="h"    />
+                <field field_ref="o"    />
+                <field field_ref="o2"   />
+                <field field_ref="n2"   />
+                <field field_ref="co2"  />
+                <field field_ref="dugwno" />
+                <field field_ref="d_qmoldifCO2"  />
+                <field field_ref="d_qmoldifO3p"  />
+                <field field_ref="d_qmoldifN2"  />
+            </field_group>
+        </file>
+
+        <!-- Instantaneous outputs; 500/Vd (if day_step=480000) -->
+        <file id="ins-500perVd"
+              name="Xins-500perVd"
+              output_freq="192ts" 
+              type="one_file"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="phis" operation="once" />
+                <field field_ref="aire" operation="once" />
+                <field field_ref="tops" />
+                <field field_ref="tsol" />
+                <field field_ref="psol" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="temp" />
+                <field field_ref="pres" />
+                <field field_ref="geop" />
+                <field field_ref="vitu" />
+                <field field_ref="vitv" />
+                <field field_ref="vitw" />
+                <field field_ref="dudyn" />
+                <field field_ref="duvdf" />
+                <field field_ref="duajs" />
+                <field field_ref="dtdyn" />
+                <field field_ref="dtvdf" />
+                <field field_ref="dtajs" />
+                <field field_ref="dteuv" />
+                <field field_ref="dtcond" />
+                <field field_ref="dtswrNLTE" />
+                <field field_ref="dtswrLTE" />
+                <field field_ref="dtlwrNLTE" />
+                <field field_ref="dtlwrLTE" />
+                <field field_ref="rho"  />
+                <field field_ref="mmean" />
+                <field field_ref="h"    />
+                <field field_ref="o"    />
+                <field field_ref="o2"   />
+                <field field_ref="n2"   />
+                <field field_ref="co2"  />
+                <field field_ref="dugwno" />
+                <field field_ref="d_qmoldifCO2"  />
+                <field field_ref="d_qmoldifO3p"  />
+                <field field_ref="d_qmoldifN2"  />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; every Vday -->
+        <file id="diurnalave"
+              name="Xdiurnalave"
+              output_freq="1d" 
+              type="one_file"
+              output_level="1"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp" level="1" />
+                <field field_ref="pres" level="1" />
+                <field field_ref="geop" level="1" />
+                <field field_ref="vitu" level="1" />
+                <field field_ref="vitv" level="1" />
+                <field field_ref="vitw" level="2" />
+                <field field_ref="rho"  level="3" />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; 1000/Vd (if day_step=480000) -->
+        <file id="ave-1kperVd"
+              name="Xave-1kperVd"
+              output_freq="96ts" 
+              type="one_file"
+              output_level="3"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp" level="1" />
+                <field field_ref="pres" level="1" />
+                <field field_ref="geop" level="1" />
+                <field field_ref="vitu" level="1" />
+                <field field_ref="vitv" level="1" />
+                <field field_ref="vitw" level="2" />
+                <field field_ref="rho"  level="3" />
+                <field field_ref="mmean" level="3" />
+                <field field_ref="h"     level="3" />
+                <field field_ref="o"     level="3" />
+                <field field_ref="o2"    level="3" />
+                <field field_ref="n2"    level="3" />
+                <field field_ref="co2"   level="3" />
+                <field field_ref="dugwno" level="3" />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; 100/Vd (if day_step=480000) -->
+        <file id="ave-100perVd"
+              name="Xave-100perVd"
+              output_freq="960ts" 
+              type="one_file"
+              output_level="5"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp"  level="1" />
+                <field field_ref="pres"  level="1" />
+                <field field_ref="geop"  level="1" />
+                <field field_ref="vitu"  level="1" />
+                <field field_ref="vitv"  level="1" />
+                <field field_ref="vitw"  level="2" />
+                <field field_ref="rho"   level="2" />
+                <field field_ref="mmean" level="2" />
+                <field field_ref="dugwno" level="3" />
+                <field field_ref="h"     level="3" />
+                <field field_ref="o"     level="3" />
+                <field field_ref="o2"    level="3" />
+                <field field_ref="n2"    level="3" />
+                <field field_ref="co2"   level="3" />
+                <field field_ref="co"    level="4" />
+                <field field_ref="h2o"   level="4" />
+                <field field_ref="o3"    level="4" />
+                <field field_ref="hcl"   level="4" />
+                <field field_ref="so"    level="4" />
+                <field field_ref="so2"   level="4" />
+                <field field_ref="ocs"   level="4" />
+                <field field_ref="o2dg"  level="4" />
+                <field field_ref="h2so4" level="5" />
+                <field field_ref="h2oliq" level="5" />
+                <field field_ref="h2so4liq" level="5" />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; 24/Vd (if day_step=480000) -->
+        <file id="ave-24perVd"
+              name="Xave-24perVd"
+              output_freq="4000ts" 
+              type="one_file"
+              output_level="5"
+              enabled=".true.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp"  level="1" />
+                <field field_ref="pres"  level="1" />
+                <field field_ref="geop"  level="1" />
+                <field field_ref="vitu"  level="1" />
+                <field field_ref="vitv"  level="1" />
+                <field field_ref="vitw"  level="2" />
+                <field field_ref="dudyn" level="2" />
+                <field field_ref="duvdf" level="2" />
+                <field field_ref="duajs" level="2" />
+                <field field_ref="dtdyn" level="2" />
+                <field field_ref="dtvdf" level="2" />
+                <field field_ref="dtajs" level="2" />
+                <field field_ref="dtswr" level="2" />
+                <field field_ref="dtlwr" level="2" />
+                <field field_ref="rho"   level="3" />
+                <field field_ref="mmean" level="3" />
+                <field field_ref="dugwno" level="3" />
+                <field field_ref="h"     level="3" />
+                <field field_ref="o"     level="3" />
+                <field field_ref="o2"    level="3" />
+                <field field_ref="n2"    level="3" />
+                <field field_ref="co2"   level="3" />
+                <field field_ref="co"    level="4" />
+                <field field_ref="h2o"   level="4" />
+                <field field_ref="o3"    level="4" />
+                <field field_ref="hcl"   level="4" />
+                <field field_ref="so"    level="4" />
+                <field field_ref="so2"   level="4" />
+                <field field_ref="ocs"   level="4" />
+                <field field_ref="h2so4" level="4" />
+                <field field_ref="o2dg"  level="4" />
+                <field field_ref="so3"   level="5" />
+                <field field_ref="s2"    level="5" />
+                <field field_ref="h2oliq" level="5" />
+                <field field_ref="h2so4liq" level="5" />
+                <field field_ref="dteuv" level="5" />
+                <field field_ref="dtcond" level="5" />
+                <field field_ref="LWnet" level="6" />
+                <field field_ref="SWnet" level="6" />
+                <field field_ref="fluxvdf" level="6" />
+                <field field_ref="fluxajs" level="6" />
+                <field field_ref="Kz"    level="6" />
+            </field_group>
+        </file>
+
+    </file_definition>
+
+</context>
Index: /trunk/LMDZ.VENUS/deftank/context_lmdz_physics-notrac.xml
===================================================================
--- /trunk/LMDZ.VENUS/deftank/context_lmdz_physics-notrac.xml	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/context_lmdz_physics-notrac.xml	(revision 2279)
@@ -0,0 +1,535 @@
+<context id="LMDZ">
+
+    <!-- Definition of a calendar is done in code via wxios_set_cal -->
+    
+    <!-- Definition of domains and groups of domains -->
+    <domain_definition>
+        <domain id="dom_glo" data_dim="2" />
+    </domain_definition>
+        
+    <!-- Definition of vertical axes 
+         NB: these must be defined and set in code via
+         wxios_add_vaxis -->
+    <axis_definition>
+        <axis id="presnivs" 
+              standard_name="Pseudo-pressure of model vertical levels" 
+              unit="Pa">
+        </axis>
+        <axis id="altitude" 
+              standard_name="Pseudo-altitude of model vertical levels" 
+              unit="km">
+        </axis>
+    </axis_definition>
+    
+    <grid_definition>
+        <!-- toggle axis id below to change output vertical axis -->
+        <grid id="grid_3d">
+            <domain id="dom_glo" />
+                <axis id="presnivs" />
+                <!-- <axis id="altitude" /> -->
+        </grid>
+        <!-- to output scalars -->
+        <grid id="scalarpoint" >
+            <scalar />
+        </grid>
+    </grid_definition>
+    
+    <!-- Definition of model variables -->
+    <field_definition prec="4" 
+                      freq_op="1ts"
+                      enabled=".TRUE.">
+        
+        <!-- 0D variables 
+        <field_group id="fields_0D" grid_ref="scalarpoint">
+            <field id="ls"
+                   long_name="Solar longitude"
+                   unit="degrees" />
+        </field_group>    -->
+        
+        <!-- 2D variables -->
+        <field_group id="fields_2D" domain_ref="dom_glo">
+            <field id="phis"
+                   long_name="Surface geopotential (gz)" 
+                   unit="m2/s2" />
+            <field id="aire"
+                   long_name="Mesh area"
+                   unit="m2" />
+            <field id="tsol"
+                   long_name="Surface Temperature" 
+                   unit="K" />
+            <field id="psol"
+                   long_name="Surface Pressure" 
+                   unit="Pa" />
+            <field id="cdragh"
+                   long_name="Drag coef on T" 
+                   unit="-" />
+            <field id="cdragm"
+                   long_name="Drag coef on U" 
+                   unit="-" />
+            <field id="tops"
+                   long_name="Solar rad. at TOA" 
+                   unit="W/m2" />
+            <field id="topl"
+                   long_name="IR rad. at TOA" 
+                   unit="W/m2" />
+            <field id="sols"
+                   long_name="Solar rad. at surface" 
+                   unit="W/m2" />
+            <field id="soll"
+                   long_name="IR rad. at surface" 
+                   unit="W/m2" />
+        </field_group>
+
+        <!-- 3D variables -->
+        <field_group id="fields_3D" grid_ref="grid_3d">
+            <field id="temp"
+                   long_name="Atmospheric temperature"
+                   unit="K" />
+            <field id="pres"
+                   long_name="Atmospheric pressure" 
+                   unit="Pa" />
+            <field id="geop"
+                   long_name="Geopotential (gz)" 
+                   unit="m2/s2" />
+            <field id="vitu"
+                   long_name="Westward Zonal Wind"
+                   unit="m/s" />
+            <field id="vitv" 
+                   long_name="Northward Meridional Wind" 
+                   unit="m/s" />
+            <field id="vitw"
+                   long_name="Vertical wind" 
+                   unit="Pa/s" />
+            <field id="Kz"
+                   long_name="vertical diffusion coef" 
+                   unit="m2/s" />
+            <field id="mmean"
+                   long_name="Mean molecular mass" 
+                   unit="g/mol" />
+            <field id="rho"
+                   long_name="Atmospheric density" 
+                   unit="kg/m3" />
+            <field id="dudyn"
+                   long_name="Dynamics du" 
+                   unit="m/s2" />
+            <field id="duvdf"
+                   long_name="Boundary layer du" 
+                   unit="m/s2" />
+            <field id="dvvdf"
+                   long_name="Boundary layer dv" 
+                   unit="m/s2" />
+            <field id="duajs"
+                   long_name="Dry adj. du" 
+                   unit="m/s2" />
+            <field id="dugwo"
+                   long_name="GW oro du" 
+                   unit="m/s2" />
+            <field id="dugwno"
+                   long_name="GW non-oro du" 
+                   unit="m/s2" />
+            <field id="dumolvis"
+                   long_name="molec viscosity (u)" 
+                   unit="m/s2" />
+            <field id="dvmolvis"
+                   long_name="molec viscosity (v)" 
+                   unit="m/s2" />
+            <field id="dtdyn"
+                   long_name="Dynamics dT" 
+                   unit="K/s" />
+            <field id="dtphy"
+                   long_name="Physics dT" 
+                   unit="K/s" />
+            <field id="dtvdf"
+                   long_name="Boundary layer dT" 
+                   unit="K/s" />
+            <field id="dtajs"
+                   long_name="Dry adj. dT" 
+                   unit="K/s" />
+            <field id="dtswr"
+                   long_name="SW radiation dT" 
+                   unit="K/s" />
+            <field id="dtswrNLTE"
+                   long_name="SW radiation (NLTE) dT" 
+                   unit="K/s" />
+            <field id="dtswrLTE"
+                   long_name="SW radiation (LTE) dT" 
+                   unit="K/s" />
+            <field id="dtlwr"
+                   long_name="LW radiation dT" 
+                   unit="K/s" />
+            <field id="dtlwrNLTE"
+                   long_name="LW radiation (NLTE) dT" 
+                   unit="K/s" />
+            <field id="dtlwrLTE"
+                   long_name="LW radiation (LTE) dT" 
+                   unit="K/s" />
+            <field id="dteuv"
+                   long_name="EUV radiation dT" 
+                   unit="K/s" />
+            <field id="dtcond"
+                   long_name="Thermal conduction dT" 
+                   unit="K/s" />
+            <field id="dtec"
+                   long_name="Kinetic dissip dT" 
+                   unit="K/s" />
+            <field id="SWnet"
+                   long_name="Net SW flux" 
+                   unit="W/m2" />
+            <field id="LWnet"
+                   long_name="Net LW flux" 
+                   unit="W/m2" />
+            <field id="fluxvdf"
+                   long_name="PBL net flux" 
+                   unit="W/m2" />
+            <field id="fluxdyn"
+                   long_name="Dyn net flux" 
+                   unit="W/m2" />
+            <field id="fluxajs"
+                   long_name="Dry adj. net flux" 
+                   unit="W/m2" />
+            <field id="fluxec"
+                   long_name="Kinetic Energy net flux" 
+                   unit="W/m2" />
+            <field id="co2"
+                   long_name="CO2" 
+                   unit="mol/mol" />
+            <field id="co"
+                   long_name="CO" 
+                   unit="mol/mol" />
+            <field id="h2"
+                   long_name="H2" 
+                   unit="mol/mol" />
+            <field id="h2o"
+                   long_name="H2O" 
+                   unit="mol/mol" />
+            <field id="o1d"
+                   long_name="O(1D)" 
+                   unit="mol/mol" />
+            <field id="o"
+                   long_name="O" 
+                   unit="mol/mol" />
+            <field id="o2"
+                   long_name="O2" 
+                   unit="mol/mol" />
+            <field id="o2dg"
+                   long_name="O2(DG)" 
+                   unit="mol/mol" />
+            <field id="o3"
+                   long_name="O3" 
+                   unit="mol/mol" />
+            <field id="h"
+                   long_name="H" 
+                   unit="mol/mol" />
+            <field id="oh"
+                   long_name="OH" 
+                   unit="mol/mol" />
+            <field id="ho2"
+                   long_name="HO2" 
+                   unit="mol/mol" />
+            <field id="h2o2"
+                   long_name="H2O2" 
+                   unit="mol/mol" />
+            <field id="cl"
+                   long_name="Cl" 
+                   unit="mol/mol" />
+            <field id="clo"
+                   long_name="ClO" 
+                   unit="mol/mol" />
+            <field id="cl2"
+                   long_name="Cl2" 
+                   unit="mol/mol" />
+            <field id="hcl"
+                   long_name="HCl" 
+                   unit="mol/mol" />
+            <field id="hocl"
+                   long_name="HOCl" 
+                   unit="mol/mol" />
+            <field id="clco"
+                   long_name="ClCO" 
+                   unit="mol/mol" />
+            <field id="clco3"
+                   long_name="ClCO3" 
+                   unit="mol/mol" />
+            <field id="cocl2"
+                   long_name="COCl2" 
+                   unit="mol/mol" />
+            <field id="s"
+                   long_name="S" 
+                   unit="mol/mol" />
+            <field id="so"
+                   long_name="SO" 
+                   unit="mol/mol" />
+            <field id="so2"
+                   long_name="SO2" 
+                   unit="mol/mol" />
+            <field id="so3"
+                   long_name="SO3" 
+                   unit="mol/mol" />
+            <field id="s2o2"
+                   long_name="S2O2" 
+                   unit="mol/mol" />
+            <field id="ocs"
+                   long_name="OCS" 
+                   unit="mol/mol" />
+            <field id="hso3"
+                   long_name="HSO3" 
+                   unit="mol/mol" />
+            <field id="h2so4"
+                   long_name="H2SO4" 
+                   unit="mol/mol" />
+            <field id="s2"
+                   long_name="S2" 
+                   unit="mol/mol" />
+            <field id="clso2"
+                   long_name="ClSO2" 
+                   unit="mol/mol" />
+            <field id="oscl"
+                   long_name="OSCl" 
+                   unit="mol/mol" />
+            <field id="n2"
+                   long_name="N2" 
+                   unit="mol/mol" />
+            <field id="h2oliq"
+                   long_name="H2Oliq" 
+                   unit="mol/mol" />
+            <field id="h2so4liq"
+                   long_name="H2SO4liq" 
+                   unit="mol/mol" />
+            <field id="d_qmoldifCO2"
+                   long_name="Diff molec CO2" 
+                   unit="kg/kg" />
+            <field id="d_qmoldifO3p"
+                   long_name="Diff molec 0" 
+                   unit="kg/kg" />
+            <field id="d_qmoldifN2"
+                   long_name="Diff molec N2" 
+                   unit="kg/kg" />
+            <field id="BV2"
+		           long_name="Freq Brunt Väisälä"
+		           unit="Hz^2" />
+       </field_group>
+
+    </field_definition>
+    
+    <!-- Definition of output files
+   output_freq: ts is physics time steps, d is day
+   enabled=".true." to output a file
+      -->
+
+    <file_definition>
+        <!-- Instantaneous outputs; every physics time steps -->
+        <file id="ins"
+              name="Xins"
+              output_freq="1ts" 
+              type="one_file"
+              enabled=".false.">
+                    
+            <!-- VARS 0D
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="ls" />
+            </field_group> -->
+
+            <!-- VARS 2D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="phis" operation="once" />
+                <field field_ref="aire" operation="once" />
+                <field field_ref="tops" />
+                <field field_ref="tsol" />
+                <field field_ref="psol" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="temp" />
+                <field field_ref="pres" />
+                <field field_ref="geop" />
+                <field field_ref="vitu" />
+                <field field_ref="vitv" />
+                <field field_ref="vitw" />
+            </field_group>
+        </file>
+
+        <!-- Instantaneous outputs; 2000/Vd (if day_step=240000) -->
+        <file id="ins-2kperVd"
+              name="Xins-2kperVd"
+              output_freq="24ts" 
+              type="one_file"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="phis" operation="once" />
+                <field field_ref="aire" operation="once" />
+                <field field_ref="tops" />
+                <field field_ref="tsol" />
+                <field field_ref="psol" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="instant"
+                         freq_op="1ts">
+                <field field_ref="temp" />
+                <field field_ref="pres" />
+                <field field_ref="geop" />
+                <field field_ref="vitu" />
+                <field field_ref="vitv" />
+                <field field_ref="vitw" />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; every Vday -->
+        <file id="diurnalave"
+              name="Xdiurnalave"
+              output_freq="1d" 
+              type="one_file"
+              output_level="1"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp" level="1" />
+                <field field_ref="pres" level="1" />
+                <field field_ref="geop" level="1" />
+                <field field_ref="vitu" level="1" />
+                <field field_ref="vitv" level="1" />
+                <field field_ref="vitw" level="2" />
+                <field field_ref="rho"  level="3" />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; 1000/Vd (if day_step=240000) -->
+        <file id="ave-1kperVd"
+              name="Xave-1kperVd"
+              output_freq="48ts" 
+              type="one_file"
+              output_level="1"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp" level="1" />
+                <field field_ref="pres" level="1" />
+                <field field_ref="geop" level="1" />
+                <field field_ref="vitu" level="1" />
+                <field field_ref="vitv" level="1" />
+                <field field_ref="vitw" level="2" />
+                <field field_ref="rho"  level="3" />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; 100/Vd (if day_step=240000) -->
+        <file id="ave-100perVd"
+              name="Xave-100perVd"
+              output_freq="480ts" 
+              type="one_file"
+              output_level="1"
+              enabled=".false.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp"  level="1" />
+                <field field_ref="pres"  level="1" />
+                <field field_ref="geop"  level="1" />
+                <field field_ref="vitu"  level="1" />
+                <field field_ref="vitv"  level="1" />
+                <field field_ref="BV2"   level="2" />
+                <field field_ref="vitw"  level="2" />
+                <field field_ref="dudyn" level="2" />
+                <field field_ref="duvdf" level="2" />
+                <field field_ref="duajs" level="2" />
+                <field field_ref="dtdyn" level="2" />
+                <field field_ref="dtvdf" level="2" />
+                <field field_ref="dtajs" level="2" />
+                <field field_ref="dtswr" level="2" />
+                <field field_ref="dtlwr" level="2" />
+                <field field_ref="rho"   level="3" />
+                <field field_ref="mmean" level="3" />
+            </field_group>
+        </file>
+
+        <!-- averages outputs; 24/Vd (if day_step=240000) -->
+        <file id="ave-24perVd"
+              name="Xave-24perVd"
+              output_freq="2000ts" 
+              type="one_file"
+              output_level="2"
+              enabled=".true.">
+                    
+            <!-- VARS 2D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="phis" level="1" operation="once" />
+                <field field_ref="aire" level="1" operation="once" />
+                <field field_ref="tops" level="1" />
+                <field field_ref="tsol" level="1" />
+                <field field_ref="psol" level="1" />
+            </field_group>
+
+            <!-- VARS 3D -->
+            <field_group operation="average"
+                         freq_op="1ts">
+                <field field_ref="temp"  level="1" />
+                <field field_ref="pres"  level="1" />
+                <field field_ref="geop"  level="1" />
+                <field field_ref="vitu"  level="1" />
+                <field field_ref="vitv"  level="1" />
+                <field field_ref="BV2"   level="1" />
+                <field field_ref="vitw"  level="2" />
+                <field field_ref="dudyn" level="2" />
+                <field field_ref="duvdf" level="2" />
+                <field field_ref="duajs" level="2" />
+                <field field_ref="dtdyn" level="2" />
+                <field field_ref="dtvdf" level="2" />
+                <field field_ref="dtajs" level="2" />
+                <field field_ref="dtswr" level="2" />
+                <field field_ref="dtlwr" level="2" />
+                <field field_ref="rho"   level="3" />
+                <field field_ref="mmean" level="3" />
+                <field field_ref="LWnet" level="4" />
+                <field field_ref="SWnet" level="4" />
+                <field field_ref="fluxvdf" level="4" />
+                <field field_ref="fluxajs" level="4" />
+                <field field_ref="Kz"    level="4" />
+            </field_group>
+        </file>
+
+    </file_definition>
+
+</context>
Index: /trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml
===================================================================
--- /trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml	(revision 2278)
+++ /trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml	(revision 2279)
@@ -302,10 +302,10 @@
                    long_name="Diff molec 0" 
                    unit="kg/kg" />
-             <field id="d_qmoldifN2"
+            <field id="d_qmoldifN2"
                    long_name="Diff molec N2" 
                    unit="kg/kg" />
-             <field id="BV2"
-		   long_name="Freq Brunt Väisälä"
-		   unit="Hz^2" />
+            <field id="BV2"
+		           long_name="Freq Brunt Väisälä"
+		           unit="Hz^2" />
        </field_group>
 
@@ -470,4 +470,5 @@
                 <field field_ref="vitu"  level="1" />
                 <field field_ref="vitv"  level="1" />
+                <field field_ref="BV2"   level="2" />
                 <field field_ref="vitw"  level="2" />
                 <field field_ref="dudyn" level="2" />
@@ -479,5 +480,4 @@
                 <field field_ref="dtswr" level="2" />
                 <field field_ref="dtlwr" level="2" />
-                <field field_ref="BV2"   level="2" />
                 <field field_ref="rho"   level="3" />
                 <field field_ref="mmean" level="3" />
@@ -511,4 +511,5 @@
                 <field field_ref="vitu"  level="1" />
                 <field field_ref="vitv"  level="1" />
+                <field field_ref="BV2"   level="1" />
                 <field field_ref="vitw"  level="2" />
                 <field field_ref="dudyn" level="2" />
@@ -520,5 +521,4 @@
                 <field field_ref="dtswr" level="2" />
                 <field field_ref="dtlwr" level="2" />
-                <field field_ref="BV2"   level="2" />
                 <field field_ref="rho"   level="3" />
                 <field field_ref="mmean" level="3" />
@@ -535,4 +535,9 @@
                 <field field_ref="h2oliq" level="5" />
                 <field field_ref="h2so4liq" level="5" />
+                <field field_ref="LWnet" level="6" />
+                <field field_ref="SWnet" level="6" />
+                <field field_ref="fluxvdf" level="6" />
+                <field field_ref="fluxajs" level="6" />
+                <field field_ref="Kz"    level="6" />
             </field_group>
         </file>
Index: /trunk/LMDZ.VENUS/deftank/gcm-96x96x50-chemistry.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/gcm-96x96x50-chemistry.def	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/gcm-96x96x50-chemistry.def	(revision 2279)
@@ -0,0 +1,93 @@
+## $Header$
+#
+## Planet:
+planet_type=venus
+# 
+## Number of dynamical steps per day (must be a multiple of iperiod)      
+day_step=240000
+## 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=1.e4
+## dissipation time scale (s) for shortest wavelengths for u,v (nxgradrot)
+tetagrot=1.e4
+## dissipation time scale (s) for shortest wavelengths for u,v (divgrad) 
+tetatemp=1.e4
+## 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=5
+## runwith or without tracers                                                 
+iflag_trac=1
+## 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=50.
+# deltaz et hdelta in km
+dissip_deltaz=30.
+dissip_hdelta=5.
+# pupstart in Pa
+dissip_pupstart=1.e4
+## Sponge layer 
+# 1: over the topmost 4 levels
+# 2: over layers at pressures down to 100 times below the topmost layer pressure
+iflag_top_bound=1
+## Sponge layer 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=1
+#  Coefficient for the sponge layer (value in topmost layer)
+au_top_bound=1.e-4
+#
+###############################################
+### 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/gcm-96x96x50-notrac.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/gcm-96x96x50-notrac.def	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/gcm-96x96x50-notrac.def	(revision 2279)
@@ -0,0 +1,93 @@
+## $Header$
+#
+## Planet:
+planet_type=venus
+# 
+## Number of dynamical steps per day (must be a multiple of iperiod)      
+day_step=240000
+## 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=1.e4
+## dissipation time scale (s) for shortest wavelengths for u,v (nxgradrot)
+tetagrot=1.e4
+## dissipation time scale (s) for shortest wavelengths for u,v (divgrad) 
+tetatemp=1.e4
+## 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=5
+## 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=50.
+# deltaz et hdelta in km
+dissip_deltaz=30.
+dissip_hdelta=5.
+# pupstart in Pa
+dissip_pupstart=1.e4
+## Sponge layer 
+# 1: over the topmost 4 levels
+# 2: over layers at pressures down to 100 times below the topmost layer pressure
+iflag_top_bound=1
+## Sponge layer 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=1
+#  Coefficient for the sponge layer (value in topmost layer)
+au_top_bound=1.e-4
+#
+###############################################
+### 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/gcm-96x96x78-chemistry.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/gcm-96x96x78-chemistry.def	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/gcm-96x96x78-chemistry.def	(revision 2279)
@@ -0,0 +1,93 @@
+## $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=1
+## 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=5
+## runwith or without tracers                                                 
+iflag_trac=1
+## 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=10.
+# deltaz et hdelta in km
+dissip_deltaz=10.
+dissip_hdelta=5.
+# pupstart in Pa
+dissip_pupstart=10.
+## Sponge layer 
+# 1: over the topmost 4 levels
+# 2: over layers at pressures down to 100 times below the topmost layer pressure
+iflag_top_bound=1
+## Sponge layer 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)
+au_top_bound=1.e-4
+#
+###############################################
+### 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/gcm.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/gcm.def	(revision 2278)
+++ /trunk/LMDZ.VENUS/deftank/gcm.def	(revision 2279)
@@ -1,85 +1,93 @@
 ## $Header$
 #
-## Planete:
+## Planet:
 planet_type=venus
 # 
-## nombre de pas par jour (multiple de iperiod) ( ici pour  dt = 1 min )      
-day_step=120000
-## periode pour le pas Matsuno (en pas)
+## Number of dynamical steps per day (must be a multiple of iperiod)      
+day_step=240000
+## Apply a Matsuno step every iperiod dynamical step
 iperiod=5
-## periode de la dissipation (en pas)
-## dissip_period=5  PAR DEFAUT: 0 ie calcul autom. (=> 25 pour runs recents)
-## choix de l'operateur de dissipation (star ou  non star )
+## 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
-## nombre d'iterations de l'operateur de dissipation   gradiv
+## iterate lateral dissipation operator gradiv nitergdiv times 
 nitergdiv=1
-## nombre d'iterations de l'operateur de dissipation  nxgradrot
+## iterate lateral dissipation operator nxgradrot nitergrot times
 nitergrot=2
-## nombre d'iterations de l'operateur de dissipation  divgrad            
+## iterate lateral dissipation operator divgrad niterh times             
 niterh=2
-## temps de dissipation des plus petites long.d ondes pour u,v (gradiv)  
-tetagdiv=2.e4
-## temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
-tetagrot=2.e4
-## temps de dissipation des plus petites long.d ondes pour  h ( divgrad) 
-tetatemp=2.e4
-## coefficient pour gamdissip                                            
+## dissipation time scale (s) for shortest wavelengths for u,v (gradiv)  
+tetagdiv=1.e4
+## dissipation time scale (s) for shortest wavelengths for u,v (nxgradrot)
+tetagrot=1.e4
+## dissipation time scale (s) for shortest wavelengths for u,v (divgrad) 
+tetatemp=1.e4
+## coefficient for gamdissip                                            
 coefdis=0.
-## choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog) 
+## time marching scheme (Matsuno if purmats is y, else Matsuno-Leapfrog)
 purmats=n
-## avec ou sans physique                                                 
-## 0: pas de physique (e.g. en mode Shallow Water)
-## 1: avec physique (e.g. physique phylmd)
-## 2: avec rappel newtonien dans la dynamique                                         
+# 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
-## avec ou sans fichiers de demarrage (start.nc, startphy.nc) ?
-## (sans fichiers de demarrage, initialisation des champs par iniacademic
-##  dans la dynamique) PAS AU POINT POUR VENUS
+## 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
-## periode de la physique (en pas)                                       
-iphysiq=20
-## avec ou sans traceurs                                                 
+## call physics every iphysiq dynamical steps                                   
+iphysiq=5
+## runwith or without tracers                                                 
 iflag_trac=0
-##  Avec ou sans strato // i.e. Couche eponge et second palier pour dissip horiz.
+## run with or without stratosphere // i.e. a sponge layer and secondary
+## higher altitude level of horizontal dissipation
 ok_strato=y
-## Dissipation horizontale
+## Horizontal dissipation multipliers along the vertical
 dissip_fac_mid=2.
-dissip_fac_up=2.
-# deltaz et hdelta en km
-dissip_deltaz=10.
+dissip_fac_up=50.
+# deltaz et hdelta in km
+dissip_deltaz=30.
 dissip_hdelta=5.
-# pupstart en Pa
-dissip_pupstart=1.e3
-## Couche eponge 
-#   1: dans les 4 derniers niveaux
-#   2: dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+# pupstart in Pa
+dissip_pupstart=1.e4
+## Sponge layer 
+# 1: over the topmost 4 levels
+# 2: over layers at pressures down to 100 times below the topmost layer pressure
 iflag_top_bound=1
-## Mode Couche eponge 
-#   mode = 0 : pas de sponge
+## Sponge layer mode 
+#   mode = 0 : no sponge
 #   mode = 1 : u et v -> 0
-#   mode = 2 : u et v -> moyenne zonale
-#   mode = 3 : u, v et h -> moyenne zonale
+#   mode = 2 : u et v -> zonal average
+#   mode = 3 : u, v et h -> zonale average
 mode_top_bound=1
-#  Coefficient pour la couche eponge (valeur derniere couche)
-tau_top_bound=1.e-5
+#  Coefficient for the sponge layer (value in topmost layer)
+au_top_bound=1.e-4
 #
-## longitude en degres du centre du zoom                                 
+###############################################
+### Zoom parameters
+###############################################
+## longitude (degrees) of zoom center                              
 clon=0.
-## latitude en degres du centre du zoom                                  
+## latitude (degrees) of zoom center                              
 clat=0.
-## facteur de grossissement du zoom,selon longitude                      
+## enhancement factor of zoom, along longitudes                     
 grossismx=1.0
-## facteur de grossissement du zoom ,selon latitude                      
+## enhancement factor of zoom, along latitudes                      
 grossismy=1.0
-##  Fonction  f(y)  hyperbolique  si = .true.  , sinon  sinusoidale         
+##  Use an hyperbolic function f(y) if .true., else use a sine      
 fxyhypb=y
-## extension en longitude  de la zone du zoom  ( fraction de la zone totale)
+## extention along longitudes of zoom region (fraction of global domain)
 dzoomx=0.0
-## extension en latitude de la zone  du zoom  ( fraction de la zone totale)
+## extention along latitudes of zoom region (fraction of global domain)
 dzoomy=0.0
-##raideur du zoom en  X
+## zoom stiffness along longitudes
 taux=3.
-##raideur du zoom en  Y
+## zoom stiffness along latitudes
 tauy=3.
-##  Fonction  f(y) avec y = Sin(latit.) si = .true. , sinon y = latit.         
+##  Fonction  f(y) as y = Sin(latitude) if = .true. ,  else  y = latitude
 ysinus=y
Index: /trunk/LMDZ.VENUS/deftank/physiq-96x96x50-chemistry.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/physiq-96x96x50-chemistry.def	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/physiq-96x96x50-chemistry.def	(revision 2279)
@@ -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=n
+## 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=n
+## choice of species profile to use
+nltemodel=2
+## Compute CO2 IR absorption or not
+callnirco2=n
+##CO2 IR Absorption model
+nircorr=0
+##include thermosphere or not
+callthermos=n
+
+## Thermospheric options 
+##~~~~~~~~~~~~~~~~~~~~~~
+#Method to include solar variability?
+#0-> Old method   1-> Variability with E10.7 as observed
+solvarmod=0
+
+## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
+solarcondate = 1993.4
+
+# 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/physiq-96x96x50-notrac.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/physiq-96x96x50-notrac.def	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/physiq-96x96x50-notrac.def	(revision 2279)
@@ -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=n
+## 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 = 0
+#
+# Reinitialization of tracer abundances 
+reinit_trac=n
+#
+# parameters for chemistry and microphysics
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+#
+# use chemistry ?
+ok_chem=n
+# use clouds ? (needed for chemistry)
+ok_cloud=n
+# 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=n
+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=n
+## choice of species profile to use
+nltemodel=2
+## Compute CO2 IR absorption or not
+callnirco2=n
+##CO2 IR Absorption model
+nircorr=0
+##include thermosphere or not
+callthermos=n
+
+## Thermospheric options 
+##~~~~~~~~~~~~~~~~~~~~~~
+#Method to include solar variability?
+#0-> Old method   1-> Variability with E10.7 as observed
+solvarmod=0
+
+## (Solar min=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
+solarcondate = 1993.4
+
+# 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/physiq-96x96x78-chemistry.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/physiq-96x96x78-chemistry.def	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/physiq-96x96x78-chemistry.def	(revision 2279)
@@ -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=1996.4 ave=1993.4 max=1990.6) ; Only used if solvarmod=0
+solarcondate = 1993.4
+
+# 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       
+#
+#
+#
+#
+# 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/physiq.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/physiq.def	(revision 2278)
+++ /trunk/LMDZ.VENUS/deftank/physiq.def	(revision 2279)
@@ -1,55 +1,23 @@
-## $Header: /home/cvsroot/LMDZ4/physiq.def,v 1.2 2004/06/22 11:45:18 lmdzadmin Exp $
+## $Header$
 #
-# PARAMETRES ANCIENNEMENT DANS gcm.def
-##~~~~~~~~~~~~~~~~~~~~~~~~~
-##  Cycle diurne  ou non                 
+##  Include diurnal cycle or not
 cycle_diurne=y
-##  Soil Model  ou non               
+##  Include soil model or not
 soil_model=y
-##  Orodr  ou  non   pour l orographie              
+## Use orodr or not for orography
 ok_orodr=n
-##  Orolf  ou  non   pour l orographie              
+## Use orolf or not for orography
 ok_orolf=n
-##  Gravity Waves non-orographiques
+## Use non-orographic Gravity Waves of not
 ok_gw_nonoro=n
-## Nombre  d'appels des routines de rayonnements ( par jour)                 
-nbapp_rad=2000
-## Nombre  d'appels des routines de chimie ( par jour)                 
+## 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  pour la convection : 1 pour LMD, 2 pour Tiedtke, 3 KE(nvlle version JYG), 30 KE(version IPCC AR4), 4 KE vect
+##  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
 #
-#
-# Parametres fichiers de sortie
-##~~~~~~~~~~~~~~~~~~~~~~~~~
-#
-### OK_journe= y sortir fichier journalier histday.nc, =n pas de fichier histday.nc
-### Pour Venus, ca n'a pas de sens
-OK_journe=n
-### OK_mensuel= y sortir fichier mensuel histmth.nc, =n pas de fichier histmth.nc
-### Pour Venus, seul fichier de sorties moyennees
-OK_mensuel=y
-## frequence (en  jours ) de l'ecriture du fichier histmth pour Venus               
-# regle la frequence d'ecriture dans histmth (mensuel) et/ou histins 
-ecritphy=0.1
-### OK_instan=y, ecrire sorties "instantanees" (meme freq que histmth)
-OK_instan=n
-# 
-# parametres niveau de sorties differents fichiers
-#
-#niveau de sortie "day" lev_histday
-# - lev_hist*=1 => basiques, 2D
-# - lev_hist*=2 => basiques, 3D (defaut)
-# - lev_hist*=3 => transfert radiatif
-# - lev_hist*=4 => champs tendances 3d
-# - lev_hist*=5 => traceurs, autres
-lev_histday=2
-#niveau de sortie "mth" lev_histmth 
-lev_histmth=4
-#niveau de sortie "ins" lev_histins avec
-lev_histins=4
-#
-# parametres climatiques
-##~~~~~~~~~~~~~~~~~~~~~~~~~
+# orbital parameters
+##~~~~~~~~~~~~~~~~~~
 #
 # VENUS ##
@@ -59,5 +27,5 @@
 # R_incl = 177.4
 R_incl = 0.
-# solaire: effective, rapportee a 1 UA
+# solar: effective flux, given at 1 UA
 solaire = 328.
 #
@@ -71,14 +39,24 @@
 tr_scheme = 0
 #
+# Reinitialization of tracer abundances 
+reinit_trac=n
+#
 # parameters for chemistry and microphysics
-##~~~~~~~~~~~~~~~~~~~~~~~~~
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 #
-reinit_trac=n
+# use chemistry ?
+ok_chem=n
+# use clouds ? (needed for chemistry)
 ok_cloud=n
-ok_chem=n
+# 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=n
+nb_mode=3
 #
-# parametres boundary layer
-##~~~~~~~~~~~~~~~~~~~~~~~~~
+# parameters for the boundary layer
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 #
 iflag_pbl=8
@@ -88,29 +66,28 @@
 ok_kzmin=n
 #
+# Surface thermal inertia
 inertie=2000
 #
-# parametre convection seche
+# dry convection parameters
 ##~~~~~~~~~~~~~~~~~~~~~~~~~
 #
 iflag_ajs = 1
 #
-# parametres nlte
-##~~~~~~~~~~~~~~~~~~~~~~~~~
+# Thermosphere and nlte parameters
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 #
-##Pour avoir acces a tmoy
-moyzon_ch=y
-##Physique non-LTE ou non
+## Compute non-LTE or not
 callnlte=n
-## choix pour profil des especes
+## choice of species profile to use
 nltemodel=2
-##Absorption IR du CO2 ou non
+## Compute CO2 IR absorption or not
 callnirco2=n
-##Choix modele absorption IR du CO2
+##CO2 IR Absorption model
 nircorr=0
-##Thermosphere ou non
+##include thermosphere or not
 callthermos=n
 
 ## Thermospheric options 
-##~~~~~~~~~~~~~~~~~~~~~~~~~
+##~~~~~~~~~~~~~~~~~~~~~~
 #Method to include solar variability?
 #0-> Old method   1-> Variability with E10.7 as observed
@@ -124,3 +101,36 @@
 ## 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/run.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/run.def	(revision 2278)
+++ /trunk/LMDZ.VENUS/deftank/run.def	(revision 2279)
@@ -1,40 +1,68 @@
-#
-# $Header: /users/lmdz/cvsroot/LMDZ.3.3/run.def,v 1.2.2.3 2002/07/12 14:11:18 lmdzadmin Exp $
+## $Header$
 #
 INCLUDEDEF=physiq.def
 INCLUDEDEF=gcm.def
-## Calendrier
+## Calendar
 #calend=venus
-## Jour de l etat initial ( = 350  si 20 Decembre ,par expl. ,comme ici )
+## reference day of initial state
 dayref=1
-##  Annee de l etat  initial (   avec  4  chiffres   )
+##  year of the initial state (   4 figures number   )
 anneeref=1111
-## Remise a zero de la date initiale
+## Flag to reset initial date (0: no / 1: yes)
 raz_date=0
-## Reinit des variables de controle
+## Reinitialize control variables (in dynamics)
 resetvarc=n
-## Nombre de jours d integration
+#
+## Length of the simulation run
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+## Number of days of time integration
 nday=2
-## Si on veut moins d un jour :
+## If one wants to run less than 1 full day:
 less1day=n
-## si less1day T, fraction de jour du run :
-fractday=0.01
-## si less1day T et que les runs sont enchaines, 
-## indiquer ici la fraction de jour de depart de l etat initial precedent
-## Il faut mettre dans ce cas raz_date=1
-## Defaut=0.
-#starttime=0. 
-## periode de sortie des variables de controle (en pas)
+## if less1day is y, fraction of day to run :
+fractday=0.1
+## si less1day is y and if successive runs are made,
+## then indicate here the fraction of the initial day for the previous run
+## !! and also set raz_date=1 !!
+## Default=0.
+starttime=0. 
+#
+#
+## Control output information in the dynamics every iconser dynamical steps
+## Only used in LMDZ, in serial run (not in parallel)... 
 iconser=500
 #
-## sorties instantanees dans la dynamique (fichiers dyn_hist.nc and co.)
+## Outputs in the dynamical core (only LMDZ, and not used usually)
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+## time series outputs in the dynamics (files dyn_hist.nc and co.)
 ok_dyn_ins=n
-## periode d ecriture des sorties instantanees dans la dynamique
-## (en pas dynamiques)
+## frequency of time series outputs in the dynamics (in dynamical steps)
 iecri=500 
-## sorties de valeurs moyennes dans la dynamique (fichiers dyn_hist_ave.nc and co.)
+## time averaged outputs in the dynamics (files dyn_hist_ave.nc and co.)
 ok_dyn_ave=n
-## periode de stockage des moyennes dans la dynamique et dans dynzon (en jour)
+## rate of storage of averages in the dynamics and in dynzon (in days)
 periodav=0.1
-## flag de sortie dynzon
+## flag to trigger dynzon outputs
 ok_dynzon=n
+#
+## Parameters that can be used for newstart, when modifying initial state
+##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# modify topography (e.g. when changing resolution), needs a RELIEF.nc file
+# topoflag=y
+# Remove topography
+# notopo=y
+# Initialize winds to zero
+# razvitu=y
+# razvitv=y
+# Initialize zonal wind to superrotation (uini)
+# uini=y
+# Initialize temperature to smooth VIRA profile
+# raztemp=y
+# Initialize surface temperature to 735K
+# razTS=y
+# Reinitialize the albedo value (we may consider a map later)
+# albedoflag=y
+# if yes, value used for albedo (default: 0.1)
+# albedo=0.1
+
+
Index: /trunk/LMDZ.VENUS/deftank/z2sig-venus95niv.def
===================================================================
--- /trunk/LMDZ.VENUS/deftank/z2sig-venus95niv.def	(revision 2279)
+++ /trunk/LMDZ.VENUS/deftank/z2sig-venus95niv.def	(revision 2279)
@@ -0,0 +1,106 @@
+ 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
+35.22671758
+35.72503773
+36.22515687
+36.72527601
+37.22539515
+
+
+
+
+
+
+
+
+
+
