Changeset 4743 for LMDZ6/branches/Portage_acc/libf/phylmd/StratAer/ocs_to_so2.F90

Timestamp:

Oct 30, 2023, 5:37:00 PM (7 months ago)

Author:

Laurent Fairhead

Message:

Merge of ACC branch with 4740 revision from trunk

Location:

LMDZ6/branches/Portage_acc

Files:

• . (modified) (1 prop)
• libf/phylmd/StratAer/ocs_to_so2.F90 (modified) (3 diffs)

LMDZ6/branches/Portage_acc/libf/phylmd/StratAer/ocs_to_so2.F90

@@ -8 +8 @@
   USE infotrac_phy
   USE YOMCST, ONLY : RG
-  USE phys_local_var_mod, ONLY : OCS_lifetime, budg_3D_ocs_to_so2, budg_ocs_to_so2 
+  USE phys_local_var_mod, ONLY : OCS_lifetime, budg_3D_ocs_to_so2, budg_ocs_to_so2
+  USE strataer_local_var_mod, ONLY : flag_min_rreduce
 
   IMPLICIT NONE
@@ -23 +24 @@
   ! local variables
   INTEGER                                       :: i,j,k,nb,ilon,ilev
-
+  REAL                                          :: rreduce
+  
 !--convert OCS to SO2
   budg_3D_ocs_to_so2(:,:)=0.0
@@ -29 +31 @@
 
   DO ilon=1, klon
-  DO ilev=1, klev
-  !only in the stratosphere
-  IF (is_strato(ilon,ilev)) THEN
-    IF (OCS_lifetime(ilon,ilev).GT.0.0) THEN
-      budg_3D_ocs_to_so2(ilon,ilev)=tr_seri(ilon,ilev,id_OCS_strat)*(1.0-exp(-pdtphys/OCS_lifetime(ilon,ilev)))
-    ENDIF
-    tr_seri(ilon,ilev,id_OCS_strat)=tr_seri(ilon,ilev,id_OCS_strat) - budg_3D_ocs_to_so2(ilon,ilev)
-    tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) + mSO2mol/mOCSmol*budg_3D_ocs_to_so2(ilon,ilev)
-    !convert budget from kg(OCS)/kgA to kg(S)/m2/layer/s for saving as diagnostic
-    budg_3D_ocs_to_so2(ilon,ilev)=budg_3D_ocs_to_so2(ilon,ilev)*mSatom/mOCSmol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys
-    budg_ocs_to_so2(ilon)=budg_ocs_to_so2(ilon)+budg_3D_ocs_to_so2(ilon,ilev) 
-  ENDIF
-  ENDDO
+     DO ilev=1, klev
+        !only in the stratosphere
+        IF (is_strato(ilon,ilev)) THEN
+          IF (OCS_lifetime(ilon,ilev).GT.0.0.AND.OCS_lifetime(ilon,ilev).LT.1.E10) THEN
+            rreduce = OCS_lifetime(ilon,ilev)
+            ! Check lifetime rreduce < timestep*3 (such as H2SO4 loss > 0.28*H2SO4) with exp(-1/3)=0.72
+            IF(flag_min_rreduce) THEN
+               IF (rreduce .LT. (3.*pdtphys)) rreduce = 3.*pdtphys
+            ENDIF
+            budg_3D_ocs_to_so2(ilon,ilev)=tr_seri(ilon,ilev,id_OCS_strat)*(1.0-exp(-pdtphys/rreduce))
+            tr_seri(ilon,ilev,id_OCS_strat)=tr_seri(ilon,ilev,id_OCS_strat) - budg_3D_ocs_to_so2(ilon,ilev)
+            tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) + mSO2mol/mOCSmol*budg_3D_ocs_to_so2(ilon,ilev)
+            !convert budget from kg(OCS)/kgA to kg(S)/m2/layer/s for saving as diagnostic
+            budg_3D_ocs_to_so2(ilon,ilev)=budg_3D_ocs_to_so2(ilon,ilev)*mSatom/mOCSmol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys
+            budg_ocs_to_so2(ilon)=budg_ocs_to_so2(ilon)+budg_3D_ocs_to_so2(ilon,ilev)
+         ENDIF
+         ! END IF(OCS_lifetime(ilon,ilev).GT.0.0.AND.OCS_lifetime(ilon,ilev).LT.1.E10)
+      ENDIF
+      ! END IF(is_strato(ilon,ilev))
+   ENDDO
   ENDDO