Changeset 5605 for LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90

Timestamp:

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

Author:

lebasn

Message:

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

File:

• LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90 (modified) (3 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