Index: LMDZ6/branches/cirrus/DefLists/README
===================================================================
--- LMDZ6/branches/cirrus/DefLists/README	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/README	(revision 5202)
@@ -16,2 +16,4 @@
 physiq.def_NPv6.4          NPv6.3 without 2 CM6A bugs
 physiq.def_LUDOBest2       Ludovic' physic tuned by Frederic 
+physiq.def_NPv7.0.1b       Physique Ludovic debuguee + condensation grande echelle epuree E.Vignon
+physiq.def_NPv7.0.1c       physiq.def_NPv7.0.1b + ok_bug_ajs_cv=n
Index: LMDZ6/branches/cirrus/DefLists/context_input_lmdz_xios3.xml
===================================================================
--- LMDZ6/branches/cirrus/DefLists/context_input_lmdz_xios3.xml	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/context_input_lmdz_xios3.xml	(revision 5202)
@@ -13,6 +13,6 @@
        <field id="landice"      field_ref="landice_reg" domain_ref="from_landiceref"  />
        <field id="mask"         field_ref="mask_reg"    domain_ref="from_relief" />
-       <field id="sic_limit"    field_ref="sic_reg"     domain_ref="from_limit_amip"  />
-       <field id="sst_limit"    field_ref="sst_reg"     domain_ref="from_limit_amip" />
+       <field id="sic_limit"    field_ref="sic_reg"     axis_ref="time_sic"  domain_ref="from_limit_amip"  />
+       <field id="sst_limit"    field_ref="sst_reg"     axis_ref="time_sst" domain_ref="from_limit_amip" />
        <field id="rugos_limit"  field_ref="rugos_reg"   domain_ref="from_limit" > exp(this) </field>
        <field id="albedo_limit" field_ref="albedo_reg"  domain_ref="from_limit"  />
@@ -59,9 +59,9 @@
 <!-- Case with amip files --> 
       <file id="amipbc_sic_1x1" name="amipbc_sic_1x1" >
-        <field id="sic_reg" name="sicbcs" domain_ref="domain_limit_amip" axis_ref="time_sic" operation="instant" freq_offset="1ts"/>
+        <field id="sic_reg" name="sicbcs" domain_ref="domain_limit_amip" axis_ref="time_sic_read" operation="instant" freq_offset="1ts"/>
       </file>
  
       <file id="amipbc_sst_1x1" name="amipbc_sst_1x1" >
-         <field id="sst_reg" name="tosbcs" domain_ref="domain_limit_amip" axis_ref="time_sst" operation="instant" freq_offset="1ts">this+273.15</field> 
+         <field id="sst_reg" name="tosbcs" domain_ref="domain_limit_amip" axis_ref="time_sst_read" operation="instant" freq_offset="1ts">this+273.15</field> 
       </file>
 
@@ -69,9 +69,9 @@
 <!--  
       <file id="siconc" name="hadisst_siconc" >
-        <field id="sic_reg" name="sicbcs" domain_ref="domain_limit_amip" axis_ref="time_sic" operation="instant" freq_offset="1ts"/>
+        <field id="sic_reg" name="sicbcs" domain_ref="domain_limit_amip" axis_ref="time_sic_read" operation="instant" freq_offset="1ts"/>
       </file>
  
       <file id="tos" name="hadisst_tos" >
-         <field id="sst_reg" name="tosbcs" domain_ref="domain_limit_amip" axis_ref="time_sst" operation="instant" freq_offset="1ts">this+273.15</field> 
+         <field id="sst_reg" name="tosbcs" domain_ref="domain_limit_amip" axis_ref="time_sst_read" operation="instant" freq_offset="1ts">this+273.15</field> 
       </file>
 -->
@@ -144,6 +144,8 @@
 
   <axis_definition>
-    <axis id="time_sic"/>
-    <axis id="time_sst"/>
+    <axis id="time_sic_read"/>
+    <axis id="time_sst_read"/>
+    <axis id="time_sic" axis_ref="sic_reg::"/>
+    <axis id="time_sst" axis_ref="sst_reg::"/>
   </axis_definition>
 
Index: LMDZ6/branches/cirrus/DefLists/field_def_lmdz.xml
===================================================================
--- LMDZ6/branches/cirrus/DefLists/field_def_lmdz.xml	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/field_def_lmdz.xml	(revision 5202)
@@ -210,5 +210,5 @@
         <field id="flat_x"    long_name="flat off-wake"    unit="W/m2" />
         <field id="flat_w"    long_name="flat within-wake"    unit="W/m2" />
-        <field id="delta_tsurf"    long_name="Temperature difference (w-x)"    unit="K" enabled="false" /> <!-- ym --field disabled-- delta_tsurf is not conform to the associated grid (the is an arry containing all sub_surface) --> 
+        <field id="delta_tsurf"    long_name="Temperature difference (w-x)"    unit="K" enabled="false" /> <!-- ym -field disabled- delta_tsurf is not conform to the associated grid (the is an arry containing all sub_surface) --> 
         <field id="cdragh_x"    long_name="cdragh off-wake"    unit="-" />
         <field id="cdragh_w"    long_name="cdragh within-wake"    unit="-" />
@@ -223,4 +223,5 @@
         <field id="fqfonte"    long_name="Land ice melt"    unit="kg/m2/s" />
         <field id="runofflic"  long_name="Land ice melt to ocean" unit="kg/m2/s" />
+        <field id="runoffland" long_name="Runoff to ocean" unit="kg/m2/s" />
         <field id="mrroli"  long_name="Runoff flux over land ice" unit="kg/m2/s" />
         <field id="taux"    long_name="Zonal wind stress"    unit="Pa" />
@@ -587,9 +588,71 @@
         <field id="taur" long_name="momentum flux due to rain" 	       unit="Pa" 	       detect_missing_value=".true." />
         <field id="SSS" long_name="bulk sea-surface salinity" unit="ppt" detect_missing_value=".true." />
-    </field_group>
+            <!-- Begin Added SN isotopes 2D fields 07 2023 -->
+            <!-- water oxygen H216O H217O H218O -->
+     <field_group id="iso2D_H216O" enabled="False" > <!-- H216O -->
+        <field id="precipH216O" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
+        <field id="plulH216O"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
+        <field id="plucH216O"   long_name="Convective Precip."    unit="kg/(s*m2)" />
+        <field id="runofflandH216O" long_name="Runoff to ocean"   unit="kg/(s*m2)" /> 
+        <field id="evapH216O"   long_name="Evaporat."             unit="kg/(s*m2)" />
+        <field id="evap_terH216O"  long_name="evaporation at surface ter"    unit="kg/(s*m2)" />
+        <field id="evap_licH216O"  long_name="evaporation at surface lic"    unit="kg/(s*m2)" />
+        <field id="evap_oceH216O"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
+        <field id="evap_sicH216O"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
+     </field_group>
+     <field_group id="iso2D_H217O" enabled="False" > <!-- H217O -->
+        <field id="precipH217O" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
+        <field id="plulH217O"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
+        <field id="plucH217O"   long_name="Convective Precip."    unit="kg/(s*m2)" />
+        <field id="runofflandH217O" long_name="Runoff to ocean"   unit="kg/(s*m2)" /> 
+        <field id="evapH217O"   long_name="Evaporat."             unit="kg/(s*m2)" />
+        <field id="evap_terH217O"  long_name="evaporation at surface ter"    unit="kg/(s*m2)" />
+        <field id="evap_licH217O"  long_name="evaporation at surface lic"    unit="kg/(s*m2)" />
+        <field id="evap_oceH217O"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
+        <field id="evap_sicH217O"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
+     </field_group>
+     <field_group id="iso2D_H218O" enabled="False" > <!-- H218O -->
+        <field id="precipH218O" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
+        <field id="plulH218O"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
+        <field id="plucH218O"   long_name="Convective Precip."    unit="kg/(s*m2)" />
+        <field id="runofflandH218O" long_name="Runoff to ocean"   unit="kg/(s*m2)" />
+        <field id="evapH218O"   long_name="Evaporat."             unit="kg/(s*m2)" />
+        <field id="evap_terH218O"  long_name="evaporation at surface ter"    unit="kg/(s*m2)" />
+        <field id="evap_licH218O"  long_name="evaporation at surface lic"    unit="kg/(s*m2)" />
+        <field id="evap_oceH218O"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
+        <field id="evap_sicH218O"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
+    </field_group>
+            <!-- water hydrogen HDO HTO -->
+    <field_group id="iso2D_HDO" enabled="False" > <!-- HDO -->
+        <field id="precipHDO" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
+        <field id="plulHDO"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
+        <field id="plucHDO"   long_name="Convective Precip."    unit="kg/(s*m2)" />
+        <field id="runofflandHDO" long_name="Runoff to ocean"   unit="kg/(s*m2)" />
+        <field id="evapHDO"   long_name="Evaporat."             unit="kg/(s*m2)" />
+        <field id="evap_terHDO"  long_name="evaporation at surface ter"    unit="kg/(s*m2)" />
+        <field id="evap_licHDO"  long_name="evaporation at surface lic"    unit="kg/(s*m2)" />
+        <field id="evap_oceHDO"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
+        <field id="evap_sicHDO"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
+    </field_group>
+    <field_group id="iso2D_HTO" enabled="False" > <!-- HTO -->
+        <field id="precipHTO" long_name="Precip Totale liq+sol" unit="kg/(s*m2)" />
+        <field id="plulHTO"   long_name="Large-scale Precip."   unit="kg/(s*m2)" />
+        <field id="plucHTO"   long_name="Convective Precip."    unit="kg/(s*m2)" />
+        <field id="runofflandHTO" long_name="Runoff to ocean"   unit="kg/(s*m2)" />
+        <field id="evapHTO"   long_name="Evaporat."             unit="kg/(s*m2)" />
+        <field id="evap_terHTO"  long_name="evaporation at surface ter"    unit="kg/(s*m2)" />
+        <field id="evap_licHTO"  long_name="evaporation at surface lic"    unit="kg/(s*m2)" />
+        <field id="evap_oceHTO"  long_name="evaporation at surface oce"    unit="kg/(s*m2)" />
+        <field id="evap_sicHTO"  long_name="evaporation at surface sic"    unit="kg/(s*m2)" />
+   </field_group>
+          <!-- end Added SN isotopes 2D fields 07 2023 -->
+   </field_group>
 
     <field_group id="fields_3D" grid_ref="grid_glo_presnivs" >
         <field id="tke"    long_name="TKE"    unit="m2/s2" />
 	<field id="tke_dissip" long_name="TKE DISSIPATION" unit="m2/s3" />
+	<field id="tke_buoy" long_name="TKE Buoyancy term" unit="m2/s3" />
+	<field id="tke_shear" long_name="TKE Shear term" unit="m2/s3" />
+    <field id="tke_trans" long_name="TKE Transport term" unit="m2/s3" />
         <field id="tke_ter"    long_name="Max Turb. Kinetic Energy ter"    unit="m2/s2" />
         <field id="tke_lic"    long_name="Max Turb. Kinetic Energy lic"    unit="m2/s2" />
@@ -600,5 +663,20 @@
         <field id="tke_max_oce"    long_name="Max Turb. Kinetic Energy oce"    unit="m2/s2" operation="maximum"/>
         <field id="tke_max_sic"    long_name="Max Turb. Kinetic Energy sic"    unit="m2/s2" operation="maximum"/>
-        <field id="l_mix_ter"    long_name="PBL mixing length ter"    unit="m" />
+	
+	<field id="qrain_lsc"    long_name="LS specific rain content"    unit="kg/kg" />
+	<field id="qsnow_lsc"    long_name="LS specific snow content"    unit="kg/kg" />
+	<field id="dqreva"    long_name="LS rain tendency due to evaporation"    unit="kg/kg/s" />
+	<field id="dqssub"    long_name="LS snow tendency due to sublimation"    unit="kg/kg/s" />
+	<field id="dqrauto"    long_name="LS rain tendency due to autoconversion"    unit="kg/kg/s" />
+	<field id="dqrcol"    long_name="LS rain tendency due to collection"    unit="kg/kg/s" />
+	<field id="dqrmelt"    long_name="LS rain tendency due to melting"    unit="kg/kg/s" />
+	<field id="dqrfreez"    long_name="LS rain tendency due to freezing"    unit="kg/kg/s" />
+	<field id="dqsauto"    long_name="LS snow tendency due to autoconversion"    unit="kg/kg/s" />
+	<field id="dqsagg"    long_name="LS snow tendency due to aggregation"    unit="kg/kg/s" />
+	<field id="dqsrim"    long_name="LS snow tendency due to riming"    unit="kg/kg/s" />
+	<field id="dqsmelt"    long_name="LS snow tendency due to melting"    unit="kg/kg/s" />
+	<field id="dqsfreez"    long_name="LS snow tendency due to freezing"    unit="kg/kg/s" />
+
+	<field id="l_mix_ter"    long_name="PBL mixing length ter"    unit="m" />
         <field id="l_mix_lic"    long_name="PBL mixing length lic"    unit="m" />
         <field id="l_mix_oce"    long_name="PBL mixing length oce"    unit="m" />
@@ -645,5 +723,8 @@
         <field id="zfull"    long_name="Altitude of full pressure levels"    unit="m" />
         <field id="zhalf"    long_name="Altitude of half pressure levels"    unit="m" />
-        <field id="rneb"    long_name="Cloud fraction"    unit="-" />
+	<field id="rneb"    long_name="Cloud fraction"    unit="-" />
+	<field id="cldfraliq"    long_name="Fraction of liquid cloud"    unit="-" />
+        <field id="sigma2_icefracturb"    long_name="Variance of the diagnostic supersaturation distribution (icefrac_turb)"    unit="-" />
+        <field id="mean_icefracturb"     long_name="Mean of the diagnostic supersaturation distribution (icefrac_turb)"    unit="-" />
         <field id="rnebcon"    long_name="Convective Cloud Fraction"    unit="-" />
         <field id="rnebls"    long_name="LS Cloud fraction"    unit="-" />
@@ -835,5 +916,74 @@
 	<field id="fluxt"     long_name="flux h"     unit="W/m2" />
         <field id="fluxq"     long_name="flux q"     unit="-" />
-    </field_group>
+             <!-- Begin Added SN isotopes 3D fields 07 2023 -->
+             <!-- water oxygen H216O H217O H218O -->
+     <field_group id="iso3D_H216O" enabled="False" > <!-- H216O -->       
+        <field id="ovapH216O"   long_name="Specific humidity" unit="kg/kg" />
+        <field id="oliqH216O"   long_name="Liquid water"      unit="kg/kg" />
+        <field id="ocondH216O"  long_name="Condensed water"   unit="kg/kg" />
+        <field id="dqdynH216O"  long_name="Dynamics dQ"       unit="(kg/kg)/s" />
+        <field id="dqldynH216O" long_name="Dynamics dQL"      unit="(kg/kg)/s" />
+        <field id="dqconH216O"  long_name="Convection dQ"     unit="(kg/kg)/s" />
+        <field id="dqlscH216O"  long_name="Condensation dQ"   unit="(kg/kg)/s" />
+        <field id="dqevaH216O"  long_name="Reevaporation dQ"  unit="(kg/kg)/s" />
+        <field id="dqvdfH216O"  long_name="Boundary-layer dQ" unit="(kg/kg)/s" />
+        <field id="dqajsH216O"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
+        <field id="dqtheH216O"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
+     </field_group>
+     <field_group id="iso3D_H217O" enabled="False" > <!-- H217O -->
+        <field id="ovapH217O"   long_name="Specific humidity" unit="kg/kg" />
+        <field id="oliqH217O"   long_name="Liquid water"      unit="kg/kg" />
+        <field id="ocondH217O"  long_name="Condensed water"   unit="kg/kg" />
+        <field id="dqdynH217O"  long_name="Dynamics dQ"       unit="(kg/kg)/s" />
+        <field id="dqldynH217O" long_name="Dynamics dQL"      unit="(kg/kg)/s" />
+        <field id="dqconH217O"  long_name="Convection dQ"     unit="(kg/kg)/s" />
+        <field id="dqlscH217O"  long_name="Condensation dQ"   unit="(kg/kg)/s" />
+        <field id="dqevaH217O"  long_name="Reevaporation dQ"  unit="(kg/kg)/s" />
+        <field id="dqvdfH217O"  long_name="Boundary-layer dQ" unit="(kg/kg)/s" />
+        <field id="dqajsH217O"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
+        <field id="dqtheH217O"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
+     </field_group>
+     <field_group id="iso3D_H218O" enabled="False" > <!-- H218O -->
+        <field id="ovapH218O"   long_name="Specific humidity" unit="kg/kg" />
+        <field id="oliqH218O"   long_name="Liquid water"      unit="kg/kg" />
+        <field id="ocondH218O"  long_name="Condensed water"   unit="kg/kg" />
+        <field id="dqdynH218O"  long_name="Dynamics dQ"       unit="(kg/kg)/s" />
+        <field id="dqldynH218O" long_name="Dynamics dQL"      unit="(kg/kg)/s" />
+        <field id="dqconH218O"  long_name="Convection dQ"     unit="(kg/kg)/s" />
+        <field id="dqlscH218O"  long_name="Condensation dQ"   unit="(kg/kg)/s" />
+        <field id="dqevaH218O"  long_name="Reevaporation dQ"  unit="(kg/kg)/s" />
+        <field id="dqvdfH218O"  long_name="Boundary-layer dQ" unit="(kg/kg)/s" />
+        <field id="dqajsH218O"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
+        <field id="dqtheH218O"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
+     </field_group>
+            <!-- water hydrogen HDO HTO -->
+     <field_group id="iso3D_HDO" enabled="False" > <!-- HDO -->
+        <field id="ovapHDO"   long_name="Specific humidity" unit="kg/kg" />
+        <field id="oliqHDO"   long_name="Liquid water"      unit="kg/kg" />
+        <field id="ocondHDO"  long_name="Condensed water"   unit="kg/kg" />
+        <field id="dqdynHDO"  long_name="Dynamics dQ"       unit="(kg/kg)/s" />
+        <field id="dqldynHDO" long_name="Dynamics dQL"      unit="(kg/kg)/s" />
+        <field id="dqconHDO"  long_name="Convection dQ"     unit="(kg/kg)/s" />
+        <field id="dqlscHDO"  long_name="Condensation dQ"   unit="(kg/kg)/s" />
+        <field id="dqevaHDO"  long_name="Reevaporation dQ"  unit="(kg/kg)/s" />
+        <field id="dqvdfHDO"  long_name="Boundary-layer dQ" unit="(kg/kg)/s" />
+        <field id="dqajsHDO"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
+        <field id="dqtheHDO"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
+     </field_group>
+     <field_group id="iso3D_HTO" enabled="False" > <!-- HTO -->
+        <field id="ovapHTO"   long_name="Specific humidity" unit="kg/kg" />
+        <field id="oliqHTO"   long_name="Liquid water"      unit="kg/kg" />
+        <field id="ocondHTO"  long_name="Condensed water"   unit="kg/kg" />
+        <field id="dqdynHTO"  long_name="Dynamics dQ"       unit="(kg/kg)/s" />
+        <field id="dqldynHTO" long_name="Dynamics dQL"      unit="(kg/kg)/s" />
+        <field id="dqconHTO"  long_name="Convection dQ"     unit="(kg/kg)/s" />
+        <field id="dqlscHTO"  long_name="Condensation dQ"   unit="(kg/kg)/s" />
+        <field id="dqevaHTO"  long_name="Reevaporation dQ"  unit="(kg/kg)/s" />
+        <field id="dqvdfHTO"  long_name="Boundary-layer dQ" unit="(kg/kg)/s" />
+        <field id="dqajsHTO"  long_name="Dry adjust. dQ"    unit="(kg/kg)/s" />
+        <field id="dqtheHTO"  long_name="Thermal dQ"        unit="(kg/kg)/s" />
+     </field_group>
+             <!-- End Added SN isotopes 3D fields 07 2023 -->
+   </field_group>
 
     <field_group id="fields_3Dkp1" grid_ref="grid_glo_presinter" >
@@ -886,4 +1036,7 @@
         <field id="budg_3D_backgr_so2"   long_name="SO2 background tendency"                               unit="kg(S)/m2/layer/s" />
         <field id="R2SO4"                long_name="H2SO4 mass fraction in aerosol"                        unit="%" />
+        <field id="SAD_sulfate"                long_name="SAD WET sulfate aerosols"                        unit="cm2/cm3" />
+        <field id="reff_sulfate"                long_name="Effective	radius of WET sulfate aerosols"              unit="cm" />
+        <field id="sulfMMR"                long_name="Sulfate aerosol concentration (dry mass mixing ratio)"                        unit="kg(H2SO4)/kg(air)" />
         <field id="OCS_lifetime"         long_name="OCS lifetime"                                          unit="s" />
         <field id="SO2_lifetime"         long_name="SO2 lifetime"                                          unit="s" />
@@ -927,4 +1080,79 @@
         <field id="BIN35" long_name="Sulfur aersol bin" unit="kg kg-1" />
         <field id="BIN36" long_name="Sulfur aersol bin" unit="kg kg-1" />
+	
+	<field id="BIN01_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN02_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN03_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN04_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN05_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN06_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN07_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN08_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN09_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN10_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN11_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN12_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN13_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN14_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN15_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN16_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN17_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN18_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN19_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN20_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN21_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN22_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN23_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN24_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN25_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN26_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN27_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN28_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN29_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN30_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN31_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN32_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN33_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN34_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN35_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+        <field id="BIN36_nd_mode" long_name="Sulfur aersol bin" unit="kg kg-1" />
+	
+	<field id="BIN01_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN02_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN03_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN04_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN05_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN06_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN07_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN08_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN09_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN10_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN11_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN12_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN13_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN14_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN15_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN16_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN17_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN18_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN19_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN20_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN21_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN22_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN23_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN24_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN25_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN26_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN27_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN28_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN29_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN30_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN31_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN32_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN33_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN34_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN35_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+        <field id="BIN36_sulfmmr_mode" long_name="Dry sulfate MMR in bin" unit="kg(H2SO4)/kg(air)" />
+	
 	<field id="GASOCS"  long_name="Carbonyl sulfide gas" unit="kg kg-1" />
 	<field id="GASSO2"  long_name="Sulfur dioxide gas" unit="kg kg-1" />
Index: LMDZ6/branches/cirrus/DefLists/file_def_histday_lmdz.xml
===================================================================
--- LMDZ6/branches/cirrus/DefLists/file_def_histday_lmdz.xml	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/file_def_histday_lmdz.xml	(revision 5202)
@@ -2,5 +2,5 @@
 <file_definition type="one_file">
     <file_group id="defile">
-        <file id="histday" name="histday" output_freq="1d" output_level="_AUTO_" enabled="_AUTO_" compression_level="2" sync_freq="1d">
+        <file id="histday" name="histday" output_freq="1d" output_level="1" enabled=".TRUE." compression_level="0" sync_freq="1d">
             
             <!-- VARS 1D -->
@@ -131,4 +131,5 @@
                 <field field_ref="mrroli" level="10" />
                 <field field_ref="runofflic" level="10" />
+                <field field_ref="runoffland" level="1" /> 
                 <field field_ref="taux" level="10" />
                 <field field_ref="tauy" level="10" />
@@ -450,4 +451,53 @@
                 <field field_ref="dqthe2d" level="10" />
                 <field field_ref="dqajs2d" level="10" />
+            <!-- Begin Added SN isotopes 2D fields 07 2023 -->
+            <!-- water oxygen H216O H217O H218O -->
+                <field field_ref="precipH216O"     level="1" />
+                <field field_ref="plulH216O"       level="1" />
+                <field field_ref="plucH216O"       level="1" />
+                <field field_ref="runofflandH216O" level="1" />
+                <field field_ref="evapH216O"       level="1" />
+                <field field_ref="evap_terH216O"   level="1" />
+                <field field_ref="evap_licH216O"   level="1" />
+                <field field_ref="evap_oceH216O"   level="1" />
+                <field field_ref="evap_sicH216O"   level="1" />
+                <field field_ref="precipH217O"     level="1" />
+                <field field_ref="plulH217O"       level="1" />
+                <field field_ref="plucH217O"       level="1" />
+                <field field_ref="runofflandH217O" level="1" />
+                <field field_ref="evapH217O"       level="1" />
+                <field field_ref="evap_terH217O"   level="1" />
+                <field field_ref="evap_licH217O"   level="1" />
+                <field field_ref="evap_oceH217O"   level="1" />
+                <field field_ref="evap_sicH217O"   level="1" />
+                <field field_ref="precipH218O"     level="1" />
+                <field field_ref="plulH218O"       level="1" />
+                <field field_ref="plucH218O"       level="1" />
+                <field field_ref="runofflandH218O" level="1" />
+                <field field_ref="evapH218O"       level="1" />
+                <field field_ref="evap_terH218O"   level="1" />
+                <field field_ref="evap_licH218O"   level="1" />
+                <field field_ref="evap_oceH218O"   level="1" />
+                <field field_ref="evap_sicH218O"   level="1" />
+            <!-- water hydrogen HDO HTO -->
+                <field field_ref="precipHDO"    level="1" />
+                <field field_ref="plulHDO"      level="1" />
+                <field field_ref="plucHDO"      level="1" />
+                <field field_ref="runofflandHDO" level="1" />
+                <field field_ref="evapHDO"      level="1" />
+                <field field_ref="evap_terHDO"  level="1" />
+                <field field_ref="evap_licHDO"  level="1" />
+                <field field_ref="evap_oceHDO"  level="1" />
+                <field field_ref="evap_sicHDO"  level="1" />
+                <field field_ref="plulHTO"      level="1" />
+                <field field_ref="precipHTO"    level="1" />
+                <field field_ref="plucHTO"      level="1" />
+                <field field_ref="runofflandHTO" level="1" />
+                <field field_ref="evapHTO"      level="1" />
+                <field field_ref="evap_terHTO"  level="1" />
+                <field field_ref="evap_licHTO"  level="1" />
+                <field field_ref="evap_oceHTO"  level="1" />
+                <field field_ref="evap_sicHTO"  level="1" />
+             <!-- end Added SN isotopes 2D fields 07 2023 -->
             </field_group>
 
@@ -625,4 +675,63 @@
                 <field field_ref="l_mixmin_sic" level="10" />
                 <field field_ref="ozone_daylight" level="10" />
+             <!-- Begin Added SN isotopes 3D fields 07 2023 -->
+             <!-- water oxygen H216O H217O H218O -->
+                <field field_ref="ovapH216O"   level="1" />
+                <field field_ref="oliqH216O"   level="1" />
+                <field field_ref="ocondH216O"  level="1" />
+                <field field_ref="dqdynH216O"  level="1" />
+                <field field_ref="dqldynH216O" level="1" />
+                <field field_ref="dqconH216O"  level="1" />
+                <field field_ref="dqlscH216O"  level="1" />
+                <field field_ref="dqevaH216O"  level="1" />
+                <field field_ref="dqvdfH216O"  level="1" />
+                <field field_ref="dqajsH216O"  level="1" />
+                <field field_ref="dqtheH216O"  level="1" />
+                <field field_ref="ovapH217O"   level="1" />
+                <field field_ref="oliqH217O"   level="1" />
+                <field field_ref="ocondH217O"  level="1" />
+                <field field_ref="dqdynH217O"  level="1" />
+                <field field_ref="dqldynH217O" level="1" />
+                <field field_ref="dqconH217O"  level="1" />
+                <field field_ref="dqlscH217O"  level="1" />
+                <field field_ref="dqevaH217O"  level="1" />
+                <field field_ref="dqvdfH217O"  level="1" />
+                <field field_ref="dqajsH217O"  level="1" />
+                <field field_ref="dqtheH217O"  level="1" />
+                <field field_ref="ovapH218O"   level="1" />
+                <field field_ref="oliqH218O"   level="1" />
+                <field field_ref="ocondH218O"  level="1" />
+                <field field_ref="dqdynH218O"  level="1" />
+                <field field_ref="dqldynH218O" level="1" />
+                <field field_ref="dqconH218O"  level="1" />
+                <field field_ref="dqlscH218O"  level="1" />
+                <field field_ref="dqevaH218O"  level="1" />
+                <field field_ref="dqvdfH218O"  level="1" />
+                <field field_ref="dqajsH218O"  level="1" />
+                <field field_ref="dqtheH218O"  level="1" />
+            <!-- water hydrogen HDO HTO -->
+                <field field_ref="ovapHDO"   level="1" />
+                <field field_ref="oliqHDO"   level="1" />
+                <field field_ref="ocondHDO"  level="1" />
+                <field field_ref="dqdynHDO"  level="1" />
+                <field field_ref="dqldynHDO" level="1" />
+                <field field_ref="dqconHDO"  level="1" />
+                <field field_ref="dqlscHDO"  level="1" />
+                <field field_ref="dqevaHDO"  level="1" />
+                <field field_ref="dqvdfHDO"  level="1" />
+                <field field_ref="dqajsHDO"  level="1" />
+                <field field_ref="dqtheHDO"  level="1" />
+                <field field_ref="ovapHTO"   level="1" />
+                <field field_ref="oliqHTO"   level="1" />
+                <field field_ref="ocondHTO"  level="1" />
+                <field field_ref="dqdynHTO"  level="1" />
+                <field field_ref="dqldynHTO" level="1" />
+                <field field_ref="dqconHTO"  level="1" />
+                <field field_ref="dqlscHTO"  level="1" />
+                <field field_ref="dqevaHTO"  level="1" />
+                <field field_ref="dqvdfHTO"  level="1" />
+                <field field_ref="dqajsHTO"  level="1" />
+                <field field_ref="dqtheHTO"  level="1" />
+             <!-- End Added SN isotopes 3D fields 07 2023 -->
             </field_group>
             <field_group operation="average" grid_ref="grid_out_presinter">
Index: LMDZ6/branches/cirrus/DefLists/file_def_histdaystrataer_lmdz.xml
===================================================================
--- LMDZ6/branches/cirrus/DefLists/file_def_histdaystrataer_lmdz.xml	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/file_def_histdaystrataer_lmdz.xml	(revision 5202)
@@ -38,4 +38,7 @@
 	<field field_ref="budg_3D_backgr_so2" level="10" />
 	<field field_ref="R2SO4" level="1" />
+	<field field_ref="SAD_sulfate" level="1" />
+        <field field_ref="reff_sulfate" level="5" />
+	<field field_ref="sulfMMR" level="1" />
 	<field field_ref="OCS_lifetime" level="10" />
 	<field field_ref="SO2_lifetime" level="10" />
@@ -87,4 +90,78 @@
         <field field_ref="BIN35" level="2" />
         <field field_ref="BIN36" level="2" />
+	
+        <field field_ref="BIN01_nd_mode" level="3" />
+        <field field_ref="BIN02_nd_mode" level="3" />
+        <field field_ref="BIN03_nd_mode" level="3" />
+        <field field_ref="BIN04_nd_mode" level="3" />
+        <field field_ref="BIN05_nd_mode" level="3" />
+        <field field_ref="BIN06_nd_mode" level="3" />
+        <field field_ref="BIN07_nd_mode" level="3" />
+        <field field_ref="BIN08_nd_mode" level="3" />
+        <field field_ref="BIN09_nd_mode" level="3" />
+        <field field_ref="BIN10_nd_mode" level="3" />
+        <field field_ref="BIN11_nd_mode" level="3" />
+        <field field_ref="BIN12_nd_mode" level="3" />
+        <field field_ref="BIN13_nd_mode" level="3" />
+        <field field_ref="BIN14_nd_mode" level="3" />
+        <field field_ref="BIN15_nd_mode" level="3" />
+        <field field_ref="BIN16_nd_mode" level="3" />
+        <field field_ref="BIN17_nd_mode" level="3" />
+        <field field_ref="BIN18_nd_mode" level="3" />
+        <field field_ref="BIN19_nd_mode" level="3" />
+        <field field_ref="BIN20_nd_mode" level="3" />
+        <field field_ref="BIN21_nd_mode" level="3" />
+        <field field_ref="BIN22_nd_mode" level="3" />
+        <field field_ref="BIN23_nd_mode" level="3" />
+        <field field_ref="BIN24_nd_mode" level="3" />
+        <field field_ref="BIN25_nd_mode" level="3" />
+        <field field_ref="BIN26_nd_mode" level="3" />
+        <field field_ref="BIN27_nd_mode" level="3" />
+        <field field_ref="BIN28_nd_mode" level="3" />
+        <field field_ref="BIN29_nd_mode" level="3" />
+        <field field_ref="BIN30_nd_mode" level="3" />
+        <field field_ref="BIN31_nd_mode" level="3" />
+        <field field_ref="BIN32_nd_mode" level="3" />
+        <field field_ref="BIN33_nd_mode" level="3" />
+        <field field_ref="BIN34_nd_mode" level="3" />
+        <field field_ref="BIN35_nd_mode" level="3" />
+        <field field_ref="BIN36_nd_mode" level="3" />
+	
+        <field field_ref="BIN01_sulfmmr_mode" level="3" />
+        <field field_ref="BIN02_sulfmmr_mode" level="3" />
+        <field field_ref="BIN03_sulfmmr_mode" level="3" />
+        <field field_ref="BIN04_sulfmmr_mode" level="3" />
+        <field field_ref="BIN05_sulfmmr_mode" level="3" />
+        <field field_ref="BIN06_sulfmmr_mode" level="3" />
+        <field field_ref="BIN07_sulfmmr_mode" level="3" />
+        <field field_ref="BIN08_sulfmmr_mode" level="3" />
+        <field field_ref="BIN09_sulfmmr_mode" level="3" />
+        <field field_ref="BIN10_sulfmmr_mode" level="3" />
+        <field field_ref="BIN11_sulfmmr_mode" level="3" />
+        <field field_ref="BIN12_sulfmmr_mode" level="3" />
+        <field field_ref="BIN13_sulfmmr_mode" level="3" />
+        <field field_ref="BIN14_sulfmmr_mode" level="3" />
+        <field field_ref="BIN15_sulfmmr_mode" level="3" />
+        <field field_ref="BIN16_sulfmmr_mode" level="3" />
+        <field field_ref="BIN17_sulfmmr_mode" level="3" />
+        <field field_ref="BIN18_sulfmmr_mode" level="3" />
+        <field field_ref="BIN19_sulfmmr_mode" level="3" />
+        <field field_ref="BIN20_sulfmmr_mode" level="3" />
+        <field field_ref="BIN21_sulfmmr_mode" level="3" />
+        <field field_ref="BIN22_sulfmmr_mode" level="3" />
+        <field field_ref="BIN23_sulfmmr_mode" level="3" />
+        <field field_ref="BIN24_sulfmmr_mode" level="3" />
+        <field field_ref="BIN25_sulfmmr_mode" level="3" />
+        <field field_ref="BIN26_sulfmmr_mode" level="3" />
+        <field field_ref="BIN27_sulfmmr_mode" level="3" />
+        <field field_ref="BIN28_sulfmmr_mode" level="3" />
+        <field field_ref="BIN29_sulfmmr_mode" level="3" />
+        <field field_ref="BIN30_sulfmmr_mode" level="3" />
+        <field field_ref="BIN31_sulfmmr_mode" level="3" />
+        <field field_ref="BIN32_sulfmmr_mode" level="3" />
+        <field field_ref="BIN33_sulfmmr_mode" level="3" />
+        <field field_ref="BIN34_sulfmmr_mode" level="3" />
+        <field field_ref="BIN35_sulfmmr_mode" level="3" />
+        <field field_ref="BIN36_sulfmmr_mode" level="3" />
       </field_group>
       
Index: LMDZ6/branches/cirrus/DefLists/file_def_histmth_lmdz.xml
===================================================================
--- LMDZ6/branches/cirrus/DefLists/file_def_histmth_lmdz.xml	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/file_def_histmth_lmdz.xml	(revision 5202)
@@ -2,5 +2,5 @@
 <file_definition>
         <file_group id="defile">
-        <file id="histmth" name="histmth" output_freq="1mo" output_level="_AUTO_"  enabled="_AUTO_" compression_level="2">
+        <file id="histmth" name="histmth" output_freq="1mo" output_level="1"  enabled=".TRUE." compression_level="0">
 
             <!-- VARS 1D -->
@@ -165,4 +165,5 @@
                 <field field_ref="mrroli" level="1" />
                 <field field_ref="runofflic" level="5" />
+                <field field_ref="runoffland" level="1" />
                 <field field_ref="taux" level="1" />
                 <field field_ref="tauy" level="1" />
@@ -499,4 +500,53 @@
                 <field field_ref="dqthe2d" level="4" />
                 <field field_ref="dqajs2d" level="4" />
+            <!-- Begin Added SN isotopes 2D fields 07 2023 -->
+            <!-- water oxygen H216O H217O H218O -->
+                <field field_ref="precipH216O"     level="1" />
+                <field field_ref="plulH216O"       level="1" />
+                <field field_ref="plucH216O"       level="1" />
+                <field field_ref="runofflandH216O" level="1" />
+                <field field_ref="evapH216O"       level="1" />
+                <field field_ref="evap_terH216O"   level="1" />
+                <field field_ref="evap_licH216O"   level="1" />
+                <field field_ref="evap_oceH216O"   level="1" />
+                <field field_ref="evap_sicH216O"   level="1" />
+                <field field_ref="precipH217O"     level="1" />
+                <field field_ref="plulH217O"       level="1" />
+                <field field_ref="plucH217O"       level="1" />
+                <field field_ref="runofflandH217O" level="1" />
+                <field field_ref="evapH217O"       level="1" />
+                <field field_ref="evap_terH217O"   level="1" />
+                <field field_ref="evap_licH217O"   level="1" />
+                <field field_ref="evap_oceH217O"   level="1" />
+                <field field_ref="evap_sicH217O"   level="1" />
+                <field field_ref="precipH218O"     level="1" />
+                <field field_ref="plulH218O"       level="1" />
+                <field field_ref="plucH218O"       level="1" />
+                <field field_ref="runofflandH218O" level="1" />
+                <field field_ref="evapH218O"       level="1" />
+                <field field_ref="evap_terH218O"   level="1" />
+                <field field_ref="evap_licH218O"   level="1" />
+                <field field_ref="evap_oceH218O"   level="1" />
+                <field field_ref="evap_sicH218O"   level="1" />
+            <!-- water hydrogen HDO HTO -->
+                <field field_ref="precipHDO"     level="1" />
+                <field field_ref="plulHDO"       level="1" />
+                <field field_ref="plucHDO"       level="1" />
+                <field field_ref="runofflandHDO" level="1" />
+                <field field_ref="evapHDO"       level="1" />
+                <field field_ref="evap_terHDO"   level="1" />
+                <field field_ref="evap_licHDO"   level="1" />
+                <field field_ref="evap_oceHDO"   level="1" />
+                <field field_ref="evap_sicHDO"   level="1" />
+                <field field_ref="plulHTO"       level="1" />
+                <field field_ref="precipHTO"     level="1" />
+                <field field_ref="plucHTO"       level="1" />
+                <field field_ref="runofflandHTO" level="1" />
+                <field field_ref="evapHTO"       level="1" />
+                <field field_ref="evap_terHTO"   level="1" />
+                <field field_ref="evap_licHTO"   level="1" />
+                <field field_ref="evap_oceHTO"   level="1" />
+                <field field_ref="evap_sicHTO"   level="1" />
+         <!-- end Added SN isotopes 2D fields 07 2023 -->
             </field_group>
 
@@ -680,4 +730,63 @@
                 <field field_ref="heat_volc" level="1" />
 		<field field_ref="dqch4"     level="10" />
+             <!-- Begin Added SN isotopes 3D fields 07 2023 -->
+             <!-- water oxygen H216O H217O H218O -->
+                <field field_ref="ovapH216O"   level="1" />
+                <field field_ref="oliqH216O"   level="1" />
+                <field field_ref="ocondH216O"  level="1" />
+                <field field_ref="dqdynH216O"  level="1" />
+                <field field_ref="dqldynH216O" level="1" />
+                <field field_ref="dqconH216O"  level="1" />
+                <field field_ref="dqlscH216O"  level="1" />
+                <field field_ref="dqevaH216O"  level="1" />
+                <field field_ref="dqvdfH216O"  level="1" />
+                <field field_ref="dqajsH216O"  level="1" />
+                <field field_ref="dqtheH216O"  level="1" />
+                <field field_ref="ovapH217O"   level="1" />
+                <field field_ref="oliqH217O"   level="1" />
+                <field field_ref="ocondH217O"  level="1" />
+                <field field_ref="dqdynH217O"  level="1" />
+                <field field_ref="dqldynH217O" level="1" />
+                <field field_ref="dqconH217O"  level="1" />
+                <field field_ref="dqlscH217O"  level="1" />
+                <field field_ref="dqevaH217O"  level="1" />
+                <field field_ref="dqvdfH217O"  level="1" />
+                <field field_ref="dqajsH217O"  level="1" />
+                <field field_ref="dqtheH217O"  level="1" />
+                <field field_ref="ovapH218O"   level="1" />
+                <field field_ref="oliqH218O"   level="1" />
+                <field field_ref="ocondH218O"  level="1" />
+                <field field_ref="dqdynH218O"  level="1" />
+                <field field_ref="dqldynH218O" level="1" />
+                <field field_ref="dqconH218O"  level="1" />
+                <field field_ref="dqlscH218O"  level="1" />
+                <field field_ref="dqevaH218O"  level="1" />
+                <field field_ref="dqvdfH218O"  level="1" />
+                <field field_ref="dqajsH218O"  level="1" />
+                <field field_ref="dqtheH218O"  level="1" />
+            <!-- water hydrogen HDO HTO -->
+                <field field_ref="ovapHDO"   level="1" />
+                <field field_ref="oliqHDO"   level="1" />
+                <field field_ref="ocondHDO"  level="1" />
+                <field field_ref="dqdynHDO"  level="1" />
+                <field field_ref="dqldynHDO" level="1" />
+                <field field_ref="dqconHDO"  level="1" />
+                <field field_ref="dqlscHDO"  level="1" />
+                <field field_ref="dqevaHDO"  level="1" />
+                <field field_ref="dqvdfHDO"  level="1" />
+                <field field_ref="dqajsHDO"  level="1" />
+                <field field_ref="dqtheHDO"  level="1" />
+                <field field_ref="ovapHTO"   level="1" />
+                <field field_ref="oliqHTO"   level="1" />
+                <field field_ref="ocondHTO"  level="1" />
+                <field field_ref="dqdynHTO"  level="1" />
+                <field field_ref="dqldynHTO" level="1" />
+                <field field_ref="dqconHTO"  level="1" />
+                <field field_ref="dqlscHTO"  level="1" />
+                <field field_ref="dqevaHTO"  level="1" />
+                <field field_ref="dqvdfHTO"  level="1" />
+                <field field_ref="dqajsHTO"  level="1" />
+                <field field_ref="dqtheHTO"  level="1" />
+             <!-- End Added SN isotopes 3D fields 07 2023 -->
             </field_group>
 
Index: LMDZ6/branches/cirrus/DefLists/file_def_histstrataer_lmdz.xml
===================================================================
--- LMDZ6/branches/cirrus/DefLists/file_def_histstrataer_lmdz.xml	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/file_def_histstrataer_lmdz.xml	(revision 5202)
@@ -41,4 +41,7 @@
 	  <field field_ref="budg_3D_backgr_so2" level="1" />
 	  <field field_ref="R2SO4" level="1" />
+	  <field field_ref="SAD_sulfate" level="1" />
+	  <field field_ref="reff_sulfate" level="1" />
+	  <field field_ref="sulfMMR" level="1" />
 	  <field field_ref="OCS_lifetime" level="1" />
 	  <field field_ref="SO2_lifetime" level="1" />
@@ -90,4 +93,78 @@
 	  <field field_ref="GASSO2"  level="1" />
 	  <field field_ref="GASH2SO4"  level="1" />
+	  
+          <field field_ref="BIN01_nd_mode" level="2" />
+          <field field_ref="BIN02_nd_mode" level="2" />
+          <field field_ref="BIN03_nd_mode" level="2" />
+          <field field_ref="BIN04_nd_mode" level="2" />
+          <field field_ref="BIN05_nd_mode" level="2" />
+          <field field_ref="BIN06_nd_mode" level="2" />
+          <field field_ref="BIN07_nd_mode" level="2" />
+          <field field_ref="BIN08_nd_mode" level="2" />
+          <field field_ref="BIN09_nd_mode" level="2" />
+          <field field_ref="BIN10_nd_mode" level="2" />
+          <field field_ref="BIN11_nd_mode" level="2" />
+          <field field_ref="BIN12_nd_mode" level="2" />
+          <field field_ref="BIN13_nd_mode" level="2" />
+          <field field_ref="BIN14_nd_mode" level="2" />
+          <field field_ref="BIN15_nd_mode" level="2" />
+          <field field_ref="BIN16_nd_mode" level="2" />
+          <field field_ref="BIN17_nd_mode" level="2" />
+          <field field_ref="BIN18_nd_mode" level="2" />
+          <field field_ref="BIN19_nd_mode" level="2" />
+          <field field_ref="BIN20_nd_mode" level="2" />
+          <field field_ref="BIN21_nd_mode" level="2" />
+          <field field_ref="BIN22_nd_mode" level="2" />
+          <field field_ref="BIN23_nd_mode" level="2" />
+          <field field_ref="BIN24_nd_mode" level="2" />
+          <field field_ref="BIN25_nd_mode" level="2" />
+          <field field_ref="BIN26_nd_mode" level="2" />
+          <field field_ref="BIN27_nd_mode" level="2" />
+          <field field_ref="BIN28_nd_mode" level="2" />
+          <field field_ref="BIN29_nd_mode" level="2" />
+          <field field_ref="BIN30_nd_mode" level="2" />
+          <field field_ref="BIN31_nd_mode" level="2" />
+          <field field_ref="BIN32_nd_mode" level="2" />
+          <field field_ref="BIN33_nd_mode" level="2" />
+          <field field_ref="BIN34_nd_mode" level="2" />
+          <field field_ref="BIN35_nd_mode" level="2" />
+          <field field_ref="BIN36_nd_mode" level="2" />
+	  
+          <field field_ref="BIN01_sulfmmr_mode" level="2" />
+          <field field_ref="BIN02_sulfmmr_mode" level="2" />
+          <field field_ref="BIN03_sulfmmr_mode" level="2" />
+          <field field_ref="BIN04_sulfmmr_mode" level="2" />
+          <field field_ref="BIN05_sulfmmr_mode" level="2" />
+          <field field_ref="BIN06_sulfmmr_mode" level="2" />
+          <field field_ref="BIN07_sulfmmr_mode" level="2" />
+          <field field_ref="BIN08_sulfmmr_mode" level="2" />
+          <field field_ref="BIN09_sulfmmr_mode" level="2" />
+          <field field_ref="BIN10_sulfmmr_mode" level="2" />
+          <field field_ref="BIN11_sulfmmr_mode" level="2" />
+          <field field_ref="BIN12_sulfmmr_mode" level="2" />
+          <field field_ref="BIN13_sulfmmr_mode" level="2" />
+          <field field_ref="BIN14_sulfmmr_mode" level="2" />
+          <field field_ref="BIN15_sulfmmr_mode" level="2" />
+          <field field_ref="BIN16_sulfmmr_mode" level="2" />
+          <field field_ref="BIN17_sulfmmr_mode" level="2" />
+          <field field_ref="BIN18_sulfmmr_mode" level="2" />
+          <field field_ref="BIN19_sulfmmr_mode" level="2" />
+          <field field_ref="BIN20_sulfmmr_mode" level="2" />
+          <field field_ref="BIN21_sulfmmr_mode" level="2" />
+          <field field_ref="BIN22_sulfmmr_mode" level="2" />
+          <field field_ref="BIN23_sulfmmr_mode" level="2" />
+          <field field_ref="BIN24_sulfmmr_mode" level="2" />
+          <field field_ref="BIN25_sulfmmr_mode" level="2" />
+          <field field_ref="BIN26_sulfmmr_mode" level="2" />
+          <field field_ref="BIN27_sulfmmr_mode" level="2" />
+          <field field_ref="BIN28_sulfmmr_mode" level="2" />
+          <field field_ref="BIN29_sulfmmr_mode" level="2" />
+          <field field_ref="BIN30_sulfmmr_mode" level="2" />
+          <field field_ref="BIN31_sulfmmr_mode" level="2" />
+          <field field_ref="BIN32_sulfmmr_mode" level="2" />
+          <field field_ref="BIN33_sulfmmr_mode" level="2" />
+          <field field_ref="BIN34_sulfmmr_mode" level="2" />
+          <field field_ref="BIN35_sulfmmr_mode" level="2" />
+          <field field_ref="BIN36_sulfmmr_mode" level="2" />
 	</field_group>
 	
Index: LMDZ6/branches/cirrus/DefLists/gcm.def
===================================================================
--- LMDZ6/branches/cirrus/DefLists/gcm.def	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/gcm.def	(revision 5202)
@@ -33,12 +33,8 @@
 ##  dans la dynamique)
 read_start=y
+## choix du qsat pour le test sur l humidite pour l advection de q: q advecte< qsat aval
+adv_qsat_liq=n
 ## periode de la physique (en pas dynamiques, n'a de sens que si iflag_phys=1)                                       
 iphysiq=10
-##  Avec ou sans strato
-ok_strato=y 
-#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
-iflag_top_bound=2
-#  Coefficient pour la couche eponge (valeur derniere couche)
-tau_top_bound=5.e-5
 ## longitude en degres du centre du zoom                                 
 clon=0.
Index: LMDZ6/branches/cirrus/DefLists/gcm.def_144x142
===================================================================
--- LMDZ6/branches/cirrus/DefLists/gcm.def_144x142	(revision 4951)
+++ LMDZ6/branches/cirrus/DefLists/gcm.def_144x142	(revision 5202)
@@ -23,4 +23,6 @@
 ## choix du shema d'integration temporelle (Matsuno:y ou Matsuno-leapfrog:n) 
 purmats=n
+## choix du qsat pour le test sur l humidite pour l advection de q: q advecte< qsat aval
+adv_qsat_liq=n
 ## avec ou sans physique
 ## 0: pas de physique (e.g. en mode Shallow Water)
Index: LMDZ6/branches/cirrus/DefLists/lmdz.xml
===================================================================
--- LMDZ6/branches/cirrus/DefLists/lmdz.xml	(revision 5202)
+++ LMDZ6/branches/cirrus/DefLists/lmdz.xml	(revision 5202)
@@ -0,0 +1,4 @@
+<context_group>
+  <context id="LMDZ" src="./context_lmdz.xml"/>
+  <context id="LMDZ" src="./context_input_lmdz.xml"/>
+</context_group>
Index: LMDZ6/branches/cirrus/DefLists/physiq.def_NPv7.0.1c
===================================================================
--- LMDZ6/branches/cirrus/DefLists/physiq.def_NPv7.0.1c	(revision 5202)
+++ LMDZ6/branches/cirrus/DefLists/physiq.def_NPv7.0.1c	(revision 5202)
@@ -0,0 +1,632 @@
+#
+# $Id: physiq.def_NPv6.2_MR1 3867 2021-03-23 16:49:48Z fairhead $
+#
+
+##############################################################
+# PARAMETRES ET CLES DE LA PHYSIQUE LMDZ
+##############################################################
+# D: Valeur par default
+#
+
+###########################################################
+# divers. Devrait aller dans config.def ?
+###########################################################
+##   Si=.T. ,  lecture du fichier limit avec la bonne annee 
+ok_limitvrai=n
+
+###########################################################
+# Parametres couche limite de surface
+###########################################################
+
+          #  ...
+f_cdrag_stable=1
+
+          # Facteur sur les coefficients de trainee terrestres
+f_cdrag_ter=1
+
+          # Vitesse minimum dans les calculs de trainee Cd ( || V || + Vmin )
+min_wind_speed=0.1
+
+          #  Modelisation des bosquets
+ifl_pbltree=1
+Cd_frein=0.016
+
+###########################################################
+# parameter des flux air-mer
+###########################################################
+
+          # Facteur sur les flux sensibles et latent
+          # Introduit pour reduire l evaporation quand on avait z0h=z0m
+          # CONSEILLE : 1
+f_cdrag_oce=1
+
+          # Flag controlant le calcul des z0
+iflag_z0_oce=1
+
+          # Rapport entre les Z0 latent et sensible
+f_z0qh_oce=1.55
+
+          # Activation des rafales : Oui 1 ou non 0
+iflag_gusts=2
+
+          # facteur pour les rafales des thermiques (bounary layer)
+f_gust_bl=0.845
+
+          # facteur pour les poches (wake)
+f_gust_wk=0.2
+
+          # facteur sur qsat pour tenir compte de la salinite
+f_qsat_oce=0.98
+
+          #cdmmax : valeur max de Cd moment. Sur ocean seulement
+cdmmax=0.2
+
+          #cdhmax : valeur max de Cd evap/sensible. Sur ocean seulement
+cdhmax=0.2
+
+####################################################################
+# Couche limite turbulente
+####################################################################
+
+          # Nouvelle version de Mellor et Yamada
+new_yamada4=y
+
+          # Choix du schema numerique pour new_yamada4=y
+          #  1 schema de MAR. Bien pour CL stables mais detruit les stratocu
+          #  5 schema de MAR modifie. Precaunise.
+yamada4_num=5
+
+          # Flag de controle de la couche limite stable
+iflag_corr_sta=4
+
+          # min sur les fonction de stabilite a la surface
+f_ri_cd_min=0.01
+
+          # max du Ric pour Kz. Decouplage plus grand pour Ric plus grand.
+yamada4_ric=0.18
+
+          # Longueur de melange minimum pour Kz
+lmixmin=0
+
+          #shema de la couche de surface (D:1, 1:LMD, 8:Mellor-Yamada) 
+iflag_pbl=12
+
+          # Seuils pour la diffusion turbulente
+ksta_ter=1e-07
+ksta=1e-10
+
+          #ok_kzmin : calcul Kzmin dans la CL de surface (D: y)
+ok_kzmin=n
+
+          # Creation de TKE par orographie sous-maille
+addtkeoro=2
+alphatkeoro=1
+smallscales_tkeoro=y
+
+
+####################################################################
+# Ondes de gravite / Orographie sous maille
+####################################################################
+
+          #  Orodr  ou  non   pour l orographie
+ok_orodr=y
+
+          #  Orolf  ou  non   pour l orographie              
+ok_orolf=y
+
+          #Rugoro : inclusion de l orographie sous maille dans le Cd de surface
+          # Obsolete avec l activation des ondes de gravite
+f_rugoro=0
+
+          # Parametres controlant l effet de l orographie sous maille
+sso_gfrcri=0.7
+sso_gkdrag=0.6
+sso_gkwake=0.4
+sso_gklift=0.1
+
+          # Random gravity waves
+          #Activitate the stochastic parameterization of convective GWs
+          # (should be equal to ok_strato)
+ok_gwd_rando=y
+
+          # Scale the momentum fluxes due to convective GWs
+gwd_rando_ruwmax=0.5
+
+          #Scale the saturated momentum flux of  convective GWs 
+gwd_rando_sat=0.6
+
+          # Important pour controler la periode de la QBO
+gwd_front_ruwmax=3
+
+          # zstd_orodr_t : minimum subgrid orography standard deviation for activation of drag parameterisation (D:10, 0.:CMIP7)
+zstd_orodr_t = 0.
+
+          # zpmm_orodr_t : minimum zpic-zval value for activation of drag parameterisation (D:100., 0.:CMIP7)
+zpmm_orodr_t = 0.
+
+          # zpmm_orolf_t : minimum zpic-zval value for activation of lift parameterisation (D:100., 0.:CMIP7)
+zpmm_orolf_t =0.
+
+###############################################################
+# Rayonnement
+###############################################################
+          # activation du nouveau code de rayonnement RRTM
+          # 0 : Ancien code, 1 : RRTM et 2 : ECRAD (D=0)
+iflag_rrtm=1
+
+          # Nombre de bandes pour le SW. Mettre 2 si iflag_rrtm=0
+NSW=6
+
+    # Albedo et inertie des calottes
+alb_vis_sno_lic=0.96
+alb_nir_sno_lic=0.68
+    # Controle de l'inertie des glaciers et banquises
+    # inertie_sno : pour la neige (y compris sur les points "ter"
+    #      quand orchidee n'est pas actif)
+    # inertie_lic/sic : pour les "land ice" et "sea ice"
+inertie_sno=350
+    # iflag_sic=1 pour forcer a garder l'inertie des "sea ice"
+    # a inertie_sic meme si la surface est couverte de neige
+iflag_sic=1
+inertie_lic=2000
+inertie_sic=2000
+    # tau_gl = constante de rappel de la temperature a la surface de la glace (en nombre de jours)
+    # Utilise selon le forcage : tau_gl=5 pour Clim, tau_gl=2 pour Amip et tau_gl=1 pour Hadisst2
+tau_gl=2
+
+          # ???
+iflag_albedo=1
+
+          # Frequence appel convection. Nombre appels par jour
+          # (nbapp_cv=0 ==> appel tous les pas de temps)
+nbapp_cv=48
+
+          # Pour pouvoir revenir a un bug purement informatique de duplication
+          # des constantes thermodynamiques entre LMDZ et RRTM
+          # n : FORTEMENT CONSEILLE
+ok_bad_ecmwf_thermo=n
+
+###############################################################
+# Parametres nuages
+###############################################################
+
+          # Activation condensation de grande echelle epuree par E.Vignon : ok_new_lscp=y + ok_icefra_lscp=y
+ok_new_lscp=y
+ok_icefra_lscp=y
+
+
+          #seuils de la precipitation des nuages stratiformes (D: 2.6e-4 2.6e-4)
+cld_lc_lsc=0.00096841923606032
+cld_lc_con=0.00096841923606032
+
+          #constante de temps pour eliminer eau lsc et convective
+          # (D: 3600. 3600.)
+cld_tau_lsc=8151.30533480977
+cld_tau_con=8151.30533480977
+
+          #facteurs correctifs sur la vitesse de chute des cristaux de glace (D: 1 1)
+ffallv_lsc=0.451743342889991
+ffallv_con=0.451743342889991
+
+          # coefficient sur la reevaporation de la pluie (D: 2.e-5 n)
+          # regle a 3.e-5 sur des cas de cumulus en 1D
+coef_eva=0.00490135133370158
+          # coefficient sur la sublimation de la glace (D: 2.e-5 n)
+coef_eva_i=0.00490135133370158
+
+### iflag_autoconversion= 0 avec reste de bogue sur l autoconversion (default)
+### iflag_autoconversion= 1, sans bogue sur l autoconversion 
+### iflag_autoconversion= 2 : 1 + new numerical integration of autoconversion
+###                         : have to add niter_lscp=1
+iflag_autoconversion=1
+
+          # flag pour la reevaporation de la pluie
+          # 0: rien
+          # 1: a l ancienne
+          # 2: prend en compte la surface nuageuse max au dessus pour calculer
+          #    le max de reevaporation, comme ecart a la saturation. CONSEILLE
+          # 4: Ludo original
+          # 5: Ludo + correction Etienne
+iflag_evap_prec=5
+
+          # min local intensity for rainfall before decrease of the associated precipitation fraction (mm/s)
+          # rain_int_min is used if iflag_evap_prec>=4
+rain_int_min=1.44716149682657e-5
+
+          # Modification de la plage de temperature pour la phase mixte
+          # liquide/glace des nuages
+          # Controle par les t_glace_min/max, exposant_glace,
+          # iflag_t_glace (D=0)
+t_glace_min=243.15
+t_glace_max=273.15
+exposant_glace=0.5
+iflag_t_glace=3
+
+          # Activation de la thermodynamique de la glace
+iflag_ice_thermo=1
+
+          # normalement elle ne sert plus a rien dans fisrtilp
+          # sinon elle est un peu toxique:
+          # IF (zt(i) .LT. t_coup.and.reevap_ice) zrfln(i)=0. (Jean-Baptiste)
+reevap_ice=n
+
+          #Bergeron
+iflag_bergeron=2
+
+          # temperature dependency of qsat for iterative procedure
+          # temporary flag to switch between the old (CMIP5, D=0)
+          # and new versions (=1) versions 
+iflag_fisrtilp_qsat=4
+
+
+          #calcul des proprietes des nuages convectifs (D:1 0.375 1.e-4) 
+iflag_cld_th=6
+fact_cldcon=1
+
+          # ?????
+facttemps=0
+
+          #calcul eau condensee et fraction nuageuse a partir des PDF
+          # (D:0, 0:version avec ratqs sinon nouvelles PDFs)  
+iflag_pdf=1
+
+          #calcul epaisseur optique et emmissivite des nuages (D: y 1 0.01 0.3)
+ok_newmicro=y
+
+          # Controle des ratqs grande echelle
+iflag_ratqs=4
+ratqsp0=48254.6771101272
+ratqsdp=16383.2736877474
+ratqsbas=0.002
+ratqshaut=0.342148444633161
+
+          #rayon effectif pour les nuages de glace (D:35)
+rad_froid=35
+
+          # Rayons effectif des gouttes nuageuses
+          # ACTIF SI ok_cdnc=0 dans config.def. (D~9 a 12 microns)
+rad_chau1=12
+rad_chau2=11
+
+	  # Rayon effectif des gouttes si ok_cdnc=0
+          # Parameters in CDNC-maer link (Boucher&Lohmann 1995)
+bl95_b0=1.3
+bl95_b1=0.2
+
+          #  Choix ou non  de  New oliq               
+new_oliq=y
+
+          # Seuils oliqmax, oicemax (en g/kg).
+          # Au dela on precipite directement au sol
+oliqmax=0.0015
+oicemax=0.0015
+
+          #Rayons min et max effectifs des particules des nuages
+rei_min=18.1205
+rei_max=69.4129
+
+          # Distribution vetricale sous maille pour le schema de de nuages.
+          # 0: rien
+          # 1: Premiere version A. Jam,
+          # 2 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5
+          # 3 : J. Jouhaud, parametre cloudth_vert_alpha D=0.5/C= 1.32
+          #                           cloudth_vert_alpha_th C=2:91
+iflag_cloudth_vert=3
+cloudth_vert_alpha=0.5
+cloudth_vert_alpha_th=0.5
+
+
+######################################################################
+# Convection 
+######################################################################
+
+          #Convection scheme switch
+          # (D:2, 1:LMD, 2:Tiedtke, 3:KE New Physics, 30:KE AR4)
+iflag_con=3
+
+          # correct convection bug
+          #(default = y)
+keepbug_ice_frac=n
+
+          #output level of energy conservation diagnostics
+if_ebil=0
+
+          #maximum efficiency of cld water->precipitation conversion (D: 0.993)
+epmax=0.999071
+
+          #dry convective adjustment at beginning of Emanuel convective scheme (D: n 0)
+ok_adj_ema=n
+iflag_clw=0
+
+          # To get back to the bug where tendencies of the dry adjustment
+          # (applied to (w) before calling convection) where added directly
+          # to wake_deltat and wake_deltaq :
+          # ok_bug_ajs_cv=y
+ok_bug_ajs_cv=n
+
+          #Convective closure (D:1, 0:AR4, 2:ALE and ALP)
+iflag_clos=2
+
+          #Convective entrainment mixing law (D:1, 0:AR4=flat PDF, 1=PDF)
+iflag_mix=1
+
+          #weights of the bell shaped and flat PDF (used only if iflag_mix=1) (D: 1 0)
+qqa1=1
+qqa2=0
+
+          #Parameters of the bell shaped PDF; Fmax=center (D= 0.65), Gammas=width (D=0.05)
+cvl_corr=1
+
+          #Multiplicative factor of Emanuel convective scheme precipitation (D: 1.)
+Fmax=0.65
+
+          #altitude, in hPa, above LCL at which buoybase is computed,
+          # def = original: -40
+dpbase=-40
+
+          #reference fractional area of precipitating downdraughts,
+          # def = original: 0.01
+sigdz=0.003
+
+          #unsaturated fraction of precipitating downdraughts, def = original: 0.15
+spfac=0.15
+
+          #relaxation time of cloud base mass flux, def = original: 8000
+tau=8000
+
+          #flag for wb (= vert velocity at LFC);
+          # 0->wb=wbmax, 1->wb=f(plfc) bounded, 2->wb=f(plfc) linear, D=1
+          # Si iflag_wb>=10 : wbeff_min=iflag_wb*0.1
+          # wbmax : assymptotic value 
+flag_wb=50
+wbmax=2.8
+
+          # Divers tests
+ok_convstop=n
+tau_stop=15000
+ok_intermittent=n
+
+          # nouvel epluchage de C.Rio
+          #shedding coefficient (used when iflag_mix_adiab=1)
+iflag_mix_adiab=0
+coef_peel=0.25
+
+          #Conversion cld water->precip: Ep=f(pressure) if 0; Ep=f(temperature)
+          #if 1, def=1
+flag_epKEorig=1
+
+          # Case : flag_epKEorig=1
+          # elcrit : max condensed water concentration in kg/kg, D=0.0003,
+          # tlcrit : temperature in degre Celsius, at which Ep=Epmax, D=-55. (original)
+          # original: 0.0011, used only if flag_epKEorig=1
+elcrit=0.0003
+tlcrit=-55
+
+          # Case : flag_epKEorig=0
+          # pbcrit : pressure (hPa), above LCL below which Ep=0, D=150 (original)
+          # ptcrit : pressure (hPa), above LCL above which Ep=Epmax, D=500 (original)
+pbcrit=150
+ptcrit=500
+
+          # Precipitation ejection from adiabatic ascents
+          # cvflag_prec_eject: when TRUE ejection is activated. Default: False
+          # WARNING ! : cvflag_prec_eject=y is effective only if icvflag_Tpa=1
+          # ejectliq : 0. =>liquid precipitation is carried up in the adiab updraughts
+          #            1. =>liquid precipitation is ejected from adiab updraughts
+          # ejectice   0. =>solid precipitation is carried up in the adiab updraughts
+          #            1. =>solid precipitation is ejected from adiab updraughts
+          # ejectliq may take only 0. and 1. values
+          # ejectice may take any value between 0. and 1.
+          #   Default : cvflag_prec_eject=n, ejectliq=0., ejectice=0.
+cvflag_prec_eject=n
+ejectliq=0.
+ejectice=0.
+
+          # Check whether parcel level temperature and specific humidity
+          # are reasonable and demand that LCL be high enough
+cv_flag_feed=2
+
+          # save 20% CPU 
+ok_optim_yield=y
+
+          # Convective clouds
+iflag_cld_cv=0
+
+          # To get back to bug where convective output variables
+          # were set to 0 at the start of the physics
+ok_bug_cv_trac=n
+
+           # solving of the temperature equation for the adiabatic ascent
+           # 0 -> AJ scheme in two loops ; ice fraction given by envt. temp.
+           # 1 -> JYG scheme in one loop ; ice fraction given by envt. temp.
+           # 2 -> JYG scheme in one loop ; ice fraction given by adiab. asc. temp.
+           # Default: 0
+           # 2 is costly and ill tested; experts only.
+icvflag_Tpa=0
+
+           # taking into account the dependence of qv on qt at saturation
+           # Default: n
+qsat_depends_on_qt=n
+
+##########################################################
+# Fermeture convective et couplages
+##########################################################
+
+          # Decimal code:  + first digit drives the splitting of thermals: 
+          #                                0 = thermals see the average profiles of T and q
+          #                                1 = thermals are only outside wakes
+          #                + second digit drives the splitting of diffuse PBL:
+          #                                0 = The diffuse PBL sees the average profiles of T and q
+          #                                1 = Two BLs, one in (w) and one in (x), with the same
+          #                                    surface temperature and humidity.
+          #                                2 = Two BLs, one in (w) and one in (x), with separate
+          #                                    surface temperature and humidity.
+          #                                3 = same as 1 over ocean; same as 0 everywhere else.
+          # Example: iflag_pbl_split=12 => thermals outside wakes, split BL with Tsurf and qsurf difference
+          # iflag_pbl_split=10 => thermals outside wakes
+iflag_pbl_split=10
+
+          #rendement sur la puissance fournie par les thermiques a la convection
+alp_bl_k=0.5
+
+          #Constant value to be added to Alp (D: 0.)
+alp_offset=0
+
+          #Convection must reach a level at which T<t_top_max to be active
+t_top_max=1000
+
+          ### nrlmd le 10/04/2012: Declenchement ###
+          #Declenchement par couche-limite (0: deterministe, 1: stochastique)
+iflag_trig_bl=1
+
+          #Surface critique au LCL
+s_trig=12000000
+
+          #Temps de decorrelation entre scenes de cumulus
+tau_trig_shallow=1200
+
+          #Temps de decorrelation entre scenes de cumulonimbus
+tau_trig_deep=1200
+
+          # Pour la stabilisation
+random_notrig_max=0.99
+
+          # Temperature a la quelle on fond les descentes glacees
+tmax_fonte_cv=278.15
+
+          #Fermeture par couche-limite (Nicolas Rochetin)
+          #(0: deterministe, 1: statistique, 2: statistique+convergence
+          # grande echelle)
+iflag_clos_bl=1
+
+          #coefficient in Wape->C_star formula: C_star=stark*sqrt(2*Wape),
+          # def: 0.33
+stark=0.33
+
+          #lifting efficiency of wakes; Alp=Alpk*Fip, Fip=power provided
+          # by wakes, def: 0.25
+alpk=0.25
+
+#
+##########################################################
+# Couche limite convective / modele du thermique
+##########################################################
+
+          # Convection seche (D:0, 0:ajustement sec,=>1:modele du thermique) 
+iflag_thermals=18
+
+          # pas de temps du spliting pour les thermiques
+          # TOURNE MAIS POSE PLUS DE PROBLEME QU ELLE N EN RESOUD
+nsplit_thermals=1
+
+          # tau_thermals pour avoir une constante de temps sur les thermiques.
+          # Non valide
+tau_thermals=0
+
+          # Flag controlant le detrainement et l entrainement
+iflag_thermals_ed=8
+
+          # On va chercher l air a z * ( 1+fact_thermals_ed_dz) pour calculer
+          # le detrainement (A. Jam)
+fact_thermals_ed_dz=0.103105610491144
+
+          # Une option de plus sur les thermiques.
+iflag_thermals_optflux=0
+
+          #couplage ou non avec la convection
+iflag_coupl=5
+
+          # Seuil pour desactiver les thermiques en cas de forte inversion
+          # en sommet de couche limite. Actif si iflag_thermals impair.
+seuil_inversion=-0.08
+
+          # Reactivation du calcul d un zmax continu pour les thermiques
+          # reactivation d une fermeture humide des thermiques
+          # iflag_thermals_closure=2 (D=1)
+iflag_thermals_closure=2
+
+          # Flag permettant de garder le bug induisant des wake_deltaT non nuls
+          # a partir des tendances des thermiques, meme en l'absence de poches.
+          # Default : ok_bug_split_th=y ==> le bug est maintenu
+ok_bug_split_th=n
+
+          # contribution of buoyancy to the plume acceleration
+thermals_afact=1.05641487083682
+
+          # drag term in the plume acceleration
+thermals_fact_epsilon=0.00205282352527457
+
+          # scaling factor for entrainment and detrainment
+thermals_betalpha=0.411960019760327
+
+          # width of the environment subgrid scale water distribution
+cloudth_sigma1s_factor=0.673674753185742
+
+##########################################################
+# Flags wakes
+##########################################################
+
+          #Wake scheme switch (D:0, 0:AR4, 1:New Physics)
+iflag_wake=1
+
+          #multiplicative factor of the damping by gravity waves, def: 4.
+coefgw=4
+
+          #wake density = number of wake centers per m2, def: 8.E-12
+          #wdens_ref=8.E-12
+wdens_ref_o=1e-09
+wdens_ref_l=8e-12
+
+          # Ajustement convectif prealable au calcul des poches
+ok_adjwk=y
+
+          # Prevent some crashes
+          # Filter out bad wakes
+flag_wk_check_trgl=n
+iflag_wk_check_trgl=2
+iflag_alp_wk_cond=1
+
+          #Wake population dynamic switch (D:0):
+          # iflag_wk_pop_dyn=0 ==> wake_dens is set to wdens_ref
+          # iflag_wk_pop_dyn=1 ==> wake_dens is a prognostic variable
+iflag_wk_pop_dyn=0
+
+iflag_wk_profile=1
+wk_pupper=0.6
+          #Bug geopotentiel moyen dans pbl_surface (D:y) 
+ok_bug_zg_wk_pbl=n
+####################################################################
+#  Conservation 
+####################################################################
+          # The source terms of the TKE prognostic equation are diagnosed
+          # from tendencies (du, dv, dT) 
+          # associated with subrgrid scaleamotions and treated as an 
+          # additional heat source.
+          # Controled by a key, iflag_ener_conserv (D=-1):
+          # 0 : no conservation, -1 : old adhoc correction for kinetic E only
+          # (used for CMIP5) 
+          # 1 : conservation, 101 : conversion from kinetic to heat only,
+          # 110 conversion from potential to heat only
+iflag_ener_conserv=10
+
+          #Activation of corrections for water conservation (D: n). 
+          ##The default is n for compatibility with earlier versions.
+          ##generally it should be y.
+ok_conserv_q=y
+
+fl_cor_ebil=1
+
+####################################################################
+# Ozone
+####################################################################
+          # new treatment for ozone into the tropopause
+adjust_tropopause=y
+
+####################################################################
+#  Conservation de l'eau
+#####################################################################
+          # conservation de l eau sur les calottes
+ok_lic_cond=y
+####################################################################
+          #min cloud droplet number concentration
+cdnc_min=10.
Index: LMDZ6/branches/cirrus/DefLists/vert_L39.def
===================================================================
--- LMDZ6/branches/cirrus/DefLists/vert_L39.def	(revision 5202)
+++ LMDZ6/branches/cirrus/DefLists/vert_L39.def	(revision 5202)
@@ -0,0 +1,6 @@
+##  Avec ou sans strato
+ok_strato=y 
+#  Couche eponge dans les couches de pression plus faible que 100 fois la pression de la derniere couche
+iflag_top_bound=2
+#  Coefficient pour la couche eponge (valeur derniere couche)
+tau_top_bound=5.e-5
