Changeset 2752 for LMDZ5/trunk/libf/phylmd/StratAer/so2_to_h2so4.F90

Timestamp:

Jan 8, 2017, 4:44:10 PM (8 years ago)

Author:

oboucher

Message:

Changes to StratAer? module, including setting up a new output file
hence dimension of output metadata is increased to 10 items.

File:

• LMDZ5/trunk/libf/phylmd/StratAer/so2_to_h2so4.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/StratAer/so2_to_h2so4.F90

@@ -1 @@
-SUBROUTINE SO2_TO_H2SO4(pdtphys,tr_seri,t_seri,pplay,paprs,sh,is_strato)
+SUBROUTINE SO2_TO_H2SO4(pdtphys,tr_seri,t_seri,pplay,paprs,is_strato)
 
   USE dimphy, ONLY : klon,klev
@@ -5 +5 @@
   USE infotrac
   USE YOMCST, ONLY : RG
-  USE phys_local_var_mod, ONLY : SO2_lifetime, sulf_convert
+  USE phys_local_var_mod, ONLY : SO2_lifetime, budg_3D_so2_to_h2so4, budg_so2_to_h2so4
 
   IMPLICIT NONE
@@ -17 +17 @@
   REAL,DIMENSION(klon,klev),INTENT(IN)          :: pplay   ! pression pour le mileu de chaque couche (en Pa)
   REAL,DIMENSION(klon,klev+1),INTENT(IN)        :: paprs   ! pression pour chaque inter-couche (en Pa)
-  REAL,DIMENSION(klon,klev),INTENT(IN)          :: sh      ! humidite specifique   
   LOGICAL,DIMENSION(klon,klev),INTENT(IN)       :: is_strato ! stratospheric flag
 
@@ -24 +23 @@
 
 !--convert SO2 to H2SO4
-  sulf_convert(:,:)=0.0
+  budg_3D_so2_to_h2so4(:,:)=0.0
+  budg_so2_to_h2so4(:)=0.0
+
   DO ilon=1, klon
   DO ilev=1, klev
@@ -30 +31 @@
   IF (is_strato(ilon,ilev)) THEN
     IF (SO2_lifetime(ilon,ilev).GT.0.0) THEN
-      sulf_convert(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/SO2_lifetime(ilon,ilev)))
+      budg_3D_so2_to_h2so4(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/SO2_lifetime(ilon,ilev)))
     ENDIF
-    tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - sulf_convert(ilon,ilev)
-    tr_seri(ilon,ilev,id_H2SO4_strat)=tr_seri(ilon,ilev,id_H2SO4_strat) + mH2SO4mol/mSO2mol*sulf_convert(ilon,ilev)
-    !convert sulf_convert from kg(SO2)/kgA to kg(S)/m2/layer/s for saving as diagnostic
-    sulf_convert(ilon,ilev)=sulf_convert(ilon,ilev)*mSatom/mSO2mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys
+    tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - budg_3D_so2_to_h2so4(ilon,ilev)
+    tr_seri(ilon,ilev,id_H2SO4_strat)=tr_seri(ilon,ilev,id_H2SO4_strat) + mH2SO4mol/mSO2mol*budg_3D_so2_to_h2so4(ilon,ilev)
+    !convert budget from kg(SO2)/kgA to kg(S)/m2/layer/s for saving as diagnostic
+    budg_3D_so2_to_h2so4(ilon,ilev)=budg_3D_so2_to_h2so4(ilon,ilev)*mSatom/mSO2mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys
+    budg_so2_to_h2so4(ilon)=budg_so2_to_h2so4(ilon)+budg_3D_so2_to_h2so4(ilon,ilev)
   ENDIF
   ENDDO