Changeset 5717 for LMDZ6/branches/contrails/libf/phylmd/StratAer

Timestamp:

Jun 18, 2025, 5:12:20 PM (5 weeks ago)

Author:

aborella

Message:

Merge with trunk r5653

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90 (modified) (3 diffs)
• libf/phylmd/StratAer/so2_to_h2so4.f90 (modified) (4 diffs)
• libf/phylmd/StratAer/strataer_local_var_mod.f90 (modified) (1 diff)

LMDZ6/branches/contrails/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/branches/contrails/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

LMDZ6/branches/contrails/libf/phylmd/StratAer/strataer_local_var_mod.f90

@@ -259 +259 @@
     
     !--initialising dry diameters to geometrically spaced mass/volume (see Jacobson 1994)
-    mdw(1)=mdwmin
-    IF (V_rat.LT.1.62) THEN ! compensate for dip in second bin for lower volume ratio
-       mdw(2)=mdw(1)*2.**(1./3.)
-       DO it=3, nbtr_bin
-          mdw(it)=mdw(it-1)*V_rat**(1./3.)
-       ENDDO
+    IF(nbtr_bin < 3) THEN
+       WRITE(lunout,*) 'WARNING: There are less than 3 sulfur aerosol class, it could be a problem for StratAer usage !'
+       WRITE(lunout,*) 'NBTR_BIN=',nbtr_bin
     ELSE
-       DO it=2, nbtr_bin
-          mdw(it)=mdw(it-1)*V_rat**(1./3.)
-       ENDDO
-    ENDIF
-    IF (is_master) WRITE(lunout,*) 'init mdw=', mdw
+       mdw(1)=mdwmin
+       IF (V_rat.LT.1.62) THEN ! compensate for dip in second bin for lower volume ratio
+          mdw(2)=mdw(1)*2.**(1./3.)
+          DO it=3, nbtr_bin
+             mdw(it)=mdw(it-1)*V_rat**(1./3.)
+          ENDDO
+       ELSE
+          DO it=2, nbtr_bin
+             mdw(it)=mdw(it-1)*V_rat**(1./3.)
+          ENDDO
+       ENDIF
+       IF (is_master) WRITE(lunout,*) 'init mdw=', mdw
+    ENDIF
     
     !   compute particle radius RRSI [cm] and volume Vbin [m3] from diameter mdw [m]