Changeset 6205

Timestamp:

Apr 24, 2026, 12:15:21 PM (10 days ago)

Author:

asima

Message:

Follow-up of r6184 : further simplification of diagnostic aod calculations at 443 and 865nm

Location:

LMDZ6/trunk/libf/phylmd

Files:

• macv2sp.f90 (modified) (2 diffs)
• mo_simple_plumes.f90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/macv2sp.f90

@@ -42 +42 @@
   !
   REAL,DIMENSION(klon,klev,nplumes) :: aod550
+  REAL,DIMENSION(klon,nplumes) :: aod550_pl
   REAL,DIMENSION(klon,klev) :: aod_prof, ssa_prof, asy_prof
   REAL,DIMENSION(klon)      :: aod_diag
@@ -132 +133 @@
   ec550aer(:,:)=ec550aer(:,:)+od550SPaer_prof(:,:)/dz(:,:)
   !
-
-  ! Diagnostic AOD calculations for other wavelengths : 443 and 865 nm  
+  !
+  ! AOD calculations for diagnostics at other wavelengths : 443 and 865 nm  
+  !   using AOD550(icol,iplume) : 
+  aod550_pl(:,:)=SUM(aod550(:,:,:),dim=2)
+  !
   !--call to sp routine -- 443 nm
   !
   CALL sp_aod_diag                    ( &
-       klev     ,klon ,  l443         , &
-       aod550   ,aod_diag )
+       klon ,  l443         , &
+       aod550_pl   ,aod_diag )
 
   od443aer(:)= od443aer(:)+aod_diag(:)
-
+  !
   !--call to sp routine -- 865 nm
   !
   CALL sp_aod_diag                    ( &
-       klev     ,klon ,  l865         , &
-       aod550   ,aod_diag )
+       klon ,  l865         , &
+       aod550_pl   ,aod_diag )
 
   od865aer(:)= od865aer(:)+aod_diag(:)   

LMDZ6/trunk/libf/phylmd/mo_simple_plumes.f90

@@ -630 +630 @@
 
   SUBROUTINE sp_aod_diag                                    ( &
-       nlevels        ,ncol           ,lambda                , &
-       aod550         ,aod_diag       )
+       ncol           ,lambda                , &
+       aod550_pl      ,aod_diag       )
     !
     ! ---------- 
     !
     INTEGER, INTENT(IN)        :: &
-         nlevels,                 & !< number of levels
          ncol                       !< number of columns
 
     REAL, INTENT(IN)           :: &
          lambda,                  & !< wavelength
-         aod550(ncol,nlevels,nplumes) !< anthropogenic aod550 profiles by individual plumes
+         aod550_pl(ncol,nplumes) !< anthropogenic aod550 map by individual plumes
 
     REAL, INTENT(OUT)          :: &
          aod_diag(ncol)             !< aerosol optical depth, diagnostic for wavelength "lambda"
 
-    INTEGER                    :: iplume, icol, k
+    INTEGER                    :: iplume, icol
 
     REAL                       ::  &
          lextinct                    !< anthropogenic aerosol extinction (function of wavelenth and aerosol type/plume)
 
-    ! initialize variables, including output
-    !
-      DO icol=1,ncol
-        aod_diag(icol) = 0.0
-      ENDDO
+    ! initialize output
+    DO icol=1,ncol
+      aod_diag(icol) = 0.0
+    ENDDO
 
     ! sum contribution from plumes to construct composite profiles of aerosol optical properties
@@ -661 +659 @@
     DO iplume=1,nplumes
       !
-      ! distribute plume optical properties across its vertical profile weighting by optical depth and scaling for
-      ! wavelength using the angstrom parameter. 
-      !
-      !   lextinct = factor for aerosol extiction, eq(10) in Stevens et al 2017
-      !    Depends on 'iplume' via 'angstrom' (Note : angstrom=2. is prescribed for all plumes)
+      ! lextinct = factor for aerosol extiction, eq(10) in Stevens et al 2017
+      !   Depends on 'iplume' via 'angstrom' (Note : angstrom=2. is prescribed for all plumes)
       lextinct = EXP(-angstrom(iplume) * LOG(lambda/550.0))
 
-      ! Vertical sum of aod550(icol,k,,iplume) profile to 2D (icol,plume)
-      !  multipled by lextinct(iplume) to get aod for "lambda" wavelength from aod550
+      ! aod550_pl(icol,iplume) multiplied by lextinct(iplume) to get aod for "lambda" wavelength
       ! version 1 compact
-      !aod_diag(:) = aod_diag(:) + SUM(aod550(:,:,iplume),dim=2) * lextinct
+      !aod_diag(:) = aod_diag(:) + aod550_pl(:,iplume) * lextinct
       ! version 2, loops
-      DO k=1,nlevels
-        DO icol = 1,ncol
-           aod_diag(icol)         = aod_diag(icol)+ aod550(icol,k,iplume) * lextinct
-        ENDDO ! k (levels)
+      DO icol = 1,ncol
+         aod_diag(icol) = aod_diag(icol)+ aod550_pl(icol,iplume) * lextinct
       ENDDO ! icol
     ENDDO ! iplume