Changeset 3094

Timestamp:

Nov 30, 2017, 2:57:51 PM (6 years ago)

Author:

oboucher

Message:

Adding calculation of H2SO4 to particles flux diag

File:

• LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90

@@ -4 +4 @@
   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 +64 @@
 
   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 +148 @@
     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