Changeset 5605 for LMDZ6/trunk/libf/phylmd/StratAer

Timestamp:

Apr 7, 2025, 4:05:50 PM (4 months ago)

Author:

lebasn

Message:

StratAer?: add new diags (4 new longwaves for ext+AOD and SO2 chem loss) + update field/file_def.

Location:

LMDZ6/trunk/libf/phylmd/StratAer

Files:

• calcaerosolstrato_rrtm.f90 (modified) (3 diffs)
• so2_to_h2so4.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90

@@ -7 +7 @@
   USE iniprint_mod_h
   USE phys_state_var_mod, ONLY: tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, tau_aero_lw_rrtm
-  USE phys_local_var_mod, ONLY: mdw, tausum_aero, tausum_strat, tau_strat_550, tau_strat_1020, stratomask
+  USE phys_local_var_mod, ONLY: mdw, tausum_aero, tausum_strat, tau_strat_443, tau_strat_550, &
+       tau_strat_670, tau_strat_765, tau_strat_1020, tau_strat_10um, stratomask
   USE aero_mod
   USE dimphy
@@ -86 +87 @@
   DO k=1,klev
     IF (stratomask(i,k).GT.0.5) THEN
-      tausum_strat(i,1)=tausum_strat(i,1)+tau_strat_wave(i,k,2)  !--550 nm
-      tausum_strat(i,2)=tausum_strat(i,2)+tau_strat_wave(i,k,5)  !--1020 nm
-      tausum_strat(i,3)=tausum_strat(i,3)+tau_strat_wave(i,k,6)  !--10 um
+      tausum_strat(i,1)=tausum_strat(i,1)+tau_strat_wave(i,k,1)  !--443 nm
+      tausum_strat(i,2)=tausum_strat(i,2)+tau_strat_wave(i,k,2)  !--550 nm
+      tausum_strat(i,3)=tausum_strat(i,3)+tau_strat_wave(i,k,3)  !--670 nm
+      tausum_strat(i,4)=tausum_strat(i,4)+tau_strat_wave(i,k,4)  !--765 nm
+      tausum_strat(i,5)=tausum_strat(i,5)+tau_strat_wave(i,k,5)  !--1020 nm
+      tausum_strat(i,6)=tausum_strat(i,6)+tau_strat_wave(i,k,6)  !--10 um
     ENDIF
   ENDDO
@@ -97 +101 @@
     zrho=pplay(i,k)/t_seri(i,k)/RD            !air density in kg/m3
     zdz=(paprs(i,k)-paprs(i,k+1))/zrho/RG     !thickness of layer in m
-    tau_strat_550(i,k)=tau_strat_wave(i,k,2)/zdz
-    tau_strat_1020(i,k)=tau_strat_wave(i,k,5)/zdz
+        tau_strat_443(i,k)=tau_strat_wave(i,k,1)/zdz
+        tau_strat_550(i,k)=tau_strat_wave(i,k,2)/zdz
+        tau_strat_670(i,k)=tau_strat_wave(i,k,3)/zdz
+        tau_strat_765(i,k)=tau_strat_wave(i,k,4)/zdz
+        tau_strat_1020(i,k)=tau_strat_wave(i,k,5)/zdz
+        tau_strat_10um(i,k)=tau_strat_wave(i,k,6)/zdz
   ENDDO
   ENDDO

LMDZ6/trunk/libf/phylmd/StratAer/so2_to_h2so4.f90

@@ -9 +9 @@
   USE yomcst_mod_h, ONLY : RG, RD
   ! lifetime (sec) et O3_clim (VMR)
-  USE phys_local_var_mod, ONLY : SO2_lifetime, H2SO4_lifetime, O3_clim, budg_3D_so2_to_h2so4, budg_so2_to_h2so4
+  USE phys_local_var_mod, ONLY : SO2_lifetime,H2SO4_lifetime,O3_clim,budg_3D_so2_to_h2so4,budg_so2_to_h2so4,SO2_chlm
   USE strataer_local_var_mod, ONLY : flag_OH_reduced, flag_H2SO4_photolysis, flag_min_rreduce
   
@@ -32 +32 @@
   budg_3D_so2_to_h2so4(:,:)=0.0
   budg_so2_to_h2so4(:)=0.0
-
+  SO2_chlm(:,:)=0.0
+  
   DO ilon=1, klon
      DO ilev=1, klev
@@ -108 +109 @@
            ! IF (SO2_lifetime(ilon,ilev).GT.0.0 .AND. SO2_lifetime(ilon,ilev).LT.1.E10) THEN
            
+           SO2_chlm(ilon,ilev) = tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/rreduce)) * &
+                pplay(ilon,ilev)/(t_seri(ilon,ilev)*1.38e-19)  / pdtphys  !SO2 loss rate [mole.cm-3.s-1]
            
            IF (flag_H2SO4_photolysis) THEN
@@ -145 +148 @@
                 (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)
+        ELSE
+           ! troposphere
+           ! SO2 tropospheric lifetime (in sec) set to 5 days
+           rreduce =  5.0*24.0*60.0*60.0
+           rrate =tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/rreduce))
+           tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - rrate
+           SO2_chlm(ilon,ilev) = rrate * &
+                pplay(ilon,ilev)/(t_seri(ilon,ilev)*1.38e-19)  / pdtphys  !SO2 loss rate [moleccm-3s-1]
         ENDIF
         ! IF (is_strato(ilon,ilev)) THEN