Changeset 3100 for LMDZ6/trunk/libf/phylmd/phytrac_mod.F90

Timestamp:

Dec 3, 2017, 6:22:35 PM (6 years ago)

Author:

oboucher

Message:

Adding some diags for the StratAer? configuration

File:

• LMDZ6/trunk/libf/phylmd/phytrac_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/phytrac_mod.F90

@@ -104 +104 @@
     USE traccoag_mod
     USE phys_local_var_mod, ONLY: mdw
+    USE phys_local_var_mod, ONLY: budg_dep_dry_ocs,   budg_dep_wet_ocs
+    USE phys_local_var_mod, ONLY: budg_dep_dry_so2,   budg_dep_wet_so2
     USE phys_local_var_mod, ONLY: budg_dep_dry_h2so4, budg_dep_wet_h2so4
     USE phys_local_var_mod, ONLY: budg_dep_dry_part,  budg_dep_wet_part
-    USE infotrac, ONLY: nbtr_sulgas, id_SO2_strat, id_H2SO4_strat
+    USE infotrac, ONLY: nbtr_sulgas, id_OCS_strat, id_SO2_strat, id_H2SO4_strat
     USE aerophys
 #endif
@@ -701 +703 @@
        IF (type_trac=='coag') THEN
          ! initialize wet deposition flux of sulfur
+         budg_dep_wet_ocs(:)=0.0
+         budg_dep_wet_so2(:)=0.0
          budg_dep_wet_h2so4(:)=0.0
          budg_dep_wet_part(:)=0.0
@@ -709 +713 @@
          DO it = 1, nbtr
          !do not include SO2 because most of it comes trom the troposphere
-           IF (it==id_H2SO4_strat) THEN
+           IF (it==id_OCS_strat) THEN
+             budg_dep_wet_ocs(i)=budg_dep_wet_ocs(i)+d_tr_cv(i,k,it)*(mSatom/mOCSmol) &
+                            & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys
+           ELSEIF (it==id_SO2_strat) THEN
+             budg_dep_wet_so2(i)=budg_dep_wet_so2(i)+d_tr_cv(i,k,it)*(mSatom/mSO2mol) &
+                            & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys
+           ELSEIF (it==id_H2SO4_strat) THEN
              budg_dep_wet_h2so4(i)=budg_dep_wet_h2so4(i)+d_tr_cv(i,k,it)*(mSatom/mH2SO4mol) &
                             & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys
@@ -773 +783 @@
 
          ! initialize dry deposition flux of sulfur
+         budg_dep_dry_ocs(:)=0.0
+         budg_dep_dry_so2(:)=0.0
          budg_dep_dry_h2so4(:)=0.0
          budg_dep_dry_part(:)=0.0
@@ -807 +819 @@
              IF (type_trac=='coag') THEN
                ! compute dry deposition flux of sulfur (sum over gases and particles)
-               IF (it==id_H2SO4_strat) THEN
+               IF (it==id_OCS_strat) THEN
+                 budg_dep_dry_ocs(:)=budg_dep_dry_ocs(:)-source(:,it)*(mSatom/mOCSmol)
+               ELSEIF (it==id_SO2_strat) THEN
+                 budg_dep_dry_so2(:)=budg_dep_dry_so2(:)-source(:,it)*(mSatom/mSO2mol)
+               ELSEIF (it==id_H2SO4_strat) THEN
                  budg_dep_dry_h2so4(:)=budg_dep_dry_h2so4(:)-source(:,it)*(mSatom/mH2SO4mol)
                ELSEIF (it.GT.nbtr_sulgas) THEN
@@ -893 +909 @@
            DO k = 1, klev
            DO it = 1, nbtr
-             IF (it==id_H2SO4_strat) THEN
+             IF (it==id_OCS_strat) THEN
+               budg_dep_wet_ocs(i)=budg_dep_wet_ocs(i)+d_tr_ls(i,k,it)*(mSatom/mOCSmol) &
+                              & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys
+             ELSEIF (it==id_SO2_strat) THEN
+               budg_dep_wet_so2(i)=budg_dep_wet_so2(i)+d_tr_ls(i,k,it)*(mSatom/mSO2mol) &
+                              & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys
+             ELSEIF (it==id_H2SO4_strat) THEN
                budg_dep_wet_h2so4(i)=budg_dep_wet_h2so4(i)+d_tr_ls(i,k,it)*(mSatom/mH2SO4mol) &
                               & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys