Changeset 3100 for LMDZ6/trunk/libf

Timestamp:

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

Author:

oboucher

Message:

Adding some diags for the StratAer? configuration

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)
• phytrac_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -482 +482 @@
 !
 !---2D budget variables
+      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_ocs
+      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_so2
       REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_h2so4
       REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_dry_part
-!$OMP THREADPRIVATE(budg_dep_dry_h2so4,budg_dep_dry_part)
+!$OMP THREADPRIVATE(budg_dep_dry_ocs,budg_dep_dry_so2,budg_dep_dry_h2so4,budg_dep_dry_part)
+      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_ocs
+      REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_so2
       REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_h2so4
       REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_dep_wet_part
-!$OMP THREADPRIVATE(budg_dep_wet_h2so4,budg_dep_wet_part)
+!$OMP THREADPRIVATE(budg_dep_wet_ocs,budg_dep_wet_so2,budg_dep_wet_h2so4,budg_dep_wet_part)
       REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_ocs
       REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: budg_emi_so2
@@ -771 +775 @@
       ALLOCATE (tau_strat_1020(klon,klev))
       ALLOCATE (tausum_strat(klon,3))
+      ALLOCATE (budg_dep_dry_ocs(klon))
+      ALLOCATE (budg_dep_wet_ocs(klon))
+      ALLOCATE (budg_dep_dry_so2(klon))
+      ALLOCATE (budg_dep_wet_so2(klon))
       ALLOCATE (budg_dep_dry_h2so4(klon))
       ALLOCATE (budg_dep_wet_h2so4(klon))
@@ -1039 +1047 @@
       DEALLOCATE (budg_3D_nucl)
       DEALLOCATE (budg_3D_cond_evap)
+      DEALLOCATE (budg_dep_dry_ocs)
+      DEALLOCATE (budg_dep_wet_ocs)
+      DEALLOCATE (budg_dep_dry_so2)
+      DEALLOCATE (budg_dep_wet_so2)
       DEALLOCATE (budg_dep_dry_h2so4)
       DEALLOCATE (budg_dep_wet_h2so4)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1303 +1303 @@
     'vsed_aer', 'Strat. aerosol sedimentation velocity (mass-weighted)', 'm/s', (/ ('', i=1, 10) /))
 !--budget  2D
+  TYPE(ctrl_out), SAVE :: o_budg_dep_dry_ocs = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
+    'budg_dep_dry_ocs',   'OCS dry deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_budg_dep_wet_ocs = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
+    'budg_dep_wet_ocs',   'OCS wet deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_budg_dep_dry_so2 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
+    'budg_dep_dry_so2',   'SO2 dry deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_budg_dep_wet_so2 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
+    'budg_dep_wet_so2',   'SO2 wet deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_budg_dep_dry_h2so4 = ctrl_out((/ 11, 11, 11, 11, 11, 11, 11, 11, 11, 1/), &
     'budg_dep_dry_h2so4', 'H2SO4 dry deposition flux', 'kg(S)/m2/s', (/ ('', i=1, 10) /))

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -198 +198 @@
          o_budg_sed_part, o_R2SO4, o_OCS_lifetime, o_SO2_lifetime, &
          o_budg_3D_backgr_ocs, o_budg_3D_backgr_so2, & 
+         o_budg_dep_dry_ocs, o_budg_dep_wet_ocs, &
+         o_budg_dep_dry_so2, o_budg_dep_wet_so2, &
          o_budg_dep_dry_h2so4, o_budg_dep_wet_h2so4, &
          o_budg_dep_dry_part, o_budg_dep_wet_part, & 
@@ -303 +305 @@
          budg_sed_part, R2SO4, OCS_lifetime, SO2_lifetime, &
          budg_3D_backgr_ocs, budg_3D_backgr_so2, & 
+         budg_dep_dry_ocs, budg_dep_wet_ocs, &
+         budg_dep_dry_so2, budg_dep_wet_so2, &
          budg_dep_dry_h2so4, budg_dep_wet_h2so4, &
          budg_dep_dry_part, budg_dep_wet_part, &
@@ -1299 +1303 @@
           CALL histwrite_phy(o_budg_3D_nucl,         budg_3D_nucl)
           CALL histwrite_phy(o_budg_3D_cond_evap,    budg_3D_cond_evap)
+          CALL histwrite_phy(o_budg_dep_dry_ocs,     budg_dep_dry_ocs)
+          CALL histwrite_phy(o_budg_dep_wet_ocs,     budg_dep_wet_ocs)
+          CALL histwrite_phy(o_budg_dep_dry_so2,     budg_dep_dry_so2)
+          CALL histwrite_phy(o_budg_dep_wet_so2,     budg_dep_wet_so2)
           CALL histwrite_phy(o_budg_dep_dry_h2so4,   budg_dep_dry_h2so4)
           CALL histwrite_phy(o_budg_dep_wet_h2so4,   budg_dep_wet_h2so4)

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