Index: LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90	(revision 3087)
+++ LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90	(revision 3094)
@@ -4,5 +4,5 @@
   USE aerophys
   USE infotrac
-  USE phys_local_var_mod, ONLY: mdw, budg_3D_nucl, budg_3D_cond_evap, R2SO4, DENSO4, f_r_wet
+  USE phys_local_var_mod, ONLY: mdw, budg_3D_nucl, budg_3D_cond_evap, budg_h2so4_to_part, R2SO4, DENSO4, f_r_wet
   USE nucleation_tstep_mod
   USE cond_evap_tstep_mod
@@ -64,10 +64,17 @@
 
   DO ilon=1, klon
+!
+!--initialisation of diagnostic
+  budg_h2so4_to_part(ilon)=0.0
+!
   DO ilev=1, klev
+!
+!--initialisation of diagnostic
+  budg_3D_nucl(ilon,ilev)=0.0
+  budg_3D_cond_evap(ilon,ilev)=0.0
+!
   ! only in the stratosphere
   IF (is_strato(ilon,ilev)) THEN
     ! initialize sulfur fluxes
-    budg_3D_nucl(ilon,ilev)=0.0
-    budg_3D_cond_evap(ilon,ilev)=0.0
     H2SO4_init=tr_seri(ilon,ilev,id_H2SO4_strat)
     ! adaptive timestep for nucleation and condensation
@@ -141,4 +148,7 @@
     budg_3D_cond_evap(ilon,ilev)=budg_3D_cond_evap(ilon,ilev)+mSatom/mH2SO4mol &
              & *evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG
+    ! compute vertically integrated flux due to the net effect of nucleation and condensation/evaporation
+    budg_h2so4_to_part(ilon)=budg_h2so4_to_part(ilon)+(H2SO4_init-tr_seri(ilon,ilev,id_H2SO4_strat)) &
+             & *mSatom/mH2SO4mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys
   ENDIF
   ENDDO