Changeset 6092 for LMDZ6

Timestamp:

Mar 5, 2026, 9:26:31 PM (5 weeks ago)

Author:

asima

Message:

For SP aerosols (flag_aerosol=7):

• new diagnostic od550SPaer,
• corrected diagnostics od550lt1aer and dryod550aer

Location:

LMDZ6/trunk/libf/phylmd

Files:

• macv2sp.f90 (modified) (3 diffs)
• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/macv2sp.f90

@@ -17 +17 @@
   USE dimphy
   USE aero_mod
-  USE phys_local_var_mod, ONLY: t_seri, od443aer, od550aer, od865aer, ec550aer, dryod550aer, od550lt1aer, dNovrN
+  USE phys_local_var_mod, ONLY: t_seri, od443aer, od550aer, od865aer, ec550aer, dryod550aer, od550lt1aer, od550SPaer, dNovrN
   !!USE YOMCST, ONLY : RD, RG
   !
@@ -107 +107 @@
   !
   !--AOD calculations for diagnostics
-  od550aer(:)=od550aer(:)+SUM(aod_prof(:,:),dim=2)
+  ! (ASima : new diagnostic od550SPaer; corrected od550lt1aer and dryod550aer)
+  !--a/ AOD of SP aerosols = vertical sum of SP aod profile
+  od550SPaer(:)=SUM(aod_prof(:,:),dim=2)
   !
-  !--dry AOD calculation for diagnostics
-  dryod550aer(:)=dryod550aer(:)+od550aer(:)
+  !--b/ total AOD = Inca1850 + SP
+  od550aer(:)=od550aer(:)+od550SPaer(:)
   !
-  !--fine-mode AOD calculation for diagnostics
-  od550lt1aer(:)=od550lt1aer(:)+od550aer(:)
+  !--c/ fine-mode AOD = Inca1850(fine mode) + SP
+  ! original, bugged : includes od550aer of Inca1850 
+  !od550lt1aer(:)=od550lt1aer(:)+od550aer(:)
+  od550lt1aer(:)=od550lt1aer(:)+od550SPaer(:)
+  !
+  !--d/ dry AOD 
+  ! original, bugged : includes od550aer of Inca1850
+  !dryod550aer(:)=dryod550aer(:)+od550aer(:)
+  dryod550aer(:)=dryod550aer(:)+od550SPaer(:)
+  !
   !
   !--extinction coefficient for diagnostic
@@ -175 +185 @@
   ENDDO
   !
-  !--renpomalizing cg and piz now that MACv2SP increments have been added
+  !--renormalizing cg and piz now that MACv2SP increments have been added
   cg_allaer(:,:,2,:) =cg_allaer(:,:,2,:) /piz_allaer(:,:,2,:)
   piz_allaer(:,:,2,:)=piz_allaer(:,:,2,:)/tau_allaer(:,:,2,:)

LMDZ6/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -261 +261 @@
       !$OMP THREADPRIVATE(ec550aer) 
       REAL, SAVE, ALLOCATABLE :: od550lt1aer(:) 
-      !$OMP THREADPRIVATE(od550lt1aer) 
+      !$OMP THREADPRIVATE(od550lt1aer)
+      REAL, SAVE, ALLOCATABLE :: od550SPaer(:)
+      !$OMP THREADPRIVATE(od550SPaer)
       REAL, SAVE, ALLOCATABLE :: sconcso4(:) 
       !$OMP THREADPRIVATE(sconcso4) 
@@ -1001 +1003 @@
       ALLOCATE(ec550aer(klon,klev))
       ALLOCATE(od550lt1aer(klon))
+      ALLOCATE(od550SPaer(klon))
       ALLOCATE(sconcso4(klon))
       ALLOCATE(sconcno3(klon))
@@ -1471 +1474 @@
       DEALLOCATE(ec550aer)
       DEALLOCATE(od550lt1aer)
+      DEALLOCATE(od550SPaer)
       DEALLOCATE(sconcso4) 
       DEALLOCATE(sconcno3) 

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1756 +1756 @@
   TYPE(ctrl_out), SAVE :: o_od550lt1aer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'od550lt1aer', 'Fine mode optical depth', '-', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_od550SPaer = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'od550SPaer', 'Simple Plume aerosol optical depth', '-', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_sconcso4 = ctrl_out((/ 2, 6, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'sconcso4', 'Surface Concentration of Sulfate ', 'kg/m3', (/ ('', i=1, 10) /))

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -124 +124 @@
          o_solldown, o_dtsvdfo, o_dtsvdft, &
          o_dtsvdfg, o_dtsvdfi, o_z0m, o_z0h,  o_od443aer, o_od550aer, &
-         o_dryod550aer, o_od865aer, o_abs550aer, o_od550lt1aer, &
+         o_dryod550aer, o_od865aer, o_abs550aer, o_od550lt1aer, o_od550SPaer,  &
          o_sconcso4, o_sconcno3, o_sconcoa, o_sconcbc, &
          o_sconcss, o_sconcdust, o_concso4, o_concno3, &
@@ -402 +402 @@
          pmflxr, pmflxs, prfl, psfl,bsfl, re, fl, rh2m, &
          qsat2m, tpote, tpot, d_ts, od443aer, od550aer, dryod550aer, &
-         od865aer, abs550aer, od550lt1aer, sconcso4, sconcno3, &
+         od865aer, abs550aer, od550lt1aer, od550SPaer, sconcso4, sconcno3, &
          sconcoa, sconcbc, sconcss, sconcdust, concso4, concno3, &
          concoa, concbc, concss, concdust, loadso4, &
@@ -2094 +2094 @@
              CALL histwrite_phy(o_loadno3, loadno3)
              CALL histwrite_phy(o_dryod550aer, dryod550aer)
+             IF (flag_aerosol.EQ.7) THEN
+                ! (ASima) od550 for Simple Plume aerosols
+                CALL histwrite_phy(o_od550SPaer, od550SPaer)
+             ENDIF
              DO naero = 1, naero_tot-1
                 CALL histwrite_phy(o_drytausumaero(naero),drytausum_aero(:,naero))