Changeset 6185

Timestamp:

Apr 22, 2026, 8:39:28 PM (12 days ago)

Author:

asima

Message:

Simplification of aod calculations for 443 and 865nm wavelengths
by adding a specific routine for this type of diagnostic

Location:

LMDZ6/trunk/libf/phylmd

Files:

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

LMDZ6/trunk/libf/phylmd/macv2sp.f90

@@ -15 +15 @@
   !ASima : *_allaer arguments and non-SP "aer" diagnostics arrive from rrtm/aeropt_5wv_rrtm.F90 with values for inca1850 aerosols
 
-  USE mo_simple_plumes, ONLY: sp_aod550_profile, sp_aop_profile, sp_setup       ! subroutines
+  USE mo_simple_plumes, ONLY: sp_aod550_profile, sp_aop_profile, sp_aod_diag, sp_setup       ! subroutines
   USE mo_simple_plumes, ONLY: nplumes, sp_initialized, aerosols_SP_forcing_year ! variables
   USE phys_cal_mod, ONLY : year_cur, year_len, days_elapsed
@@ -43 +43 @@
   REAL,DIMENSION(klon,klev,nplumes) :: aod550
   REAL,DIMENSION(klon,klev) :: aod_prof, ssa_prof, asy_prof
+  REAL,DIMENSION(klon)      :: aod_diag
   REAL,DIMENSION(klon,klev) :: z, dz
   REAL,DIMENSION(klon)      :: oro, zrho, zt
@@ -113 +114 @@
   !--AOD550 calculations for diagnostics
   !
-  ! aod550 profile = sum over the 'nplumes' dimension
+  ! AOD profile = sum over the 'nplumes' dimension
   od550SPaer_prof(:,:)=SUM(aod550(:,:,:),dim=3)
   !
@@ -135 +136 @@
   !--call to sp routine -- 443 nm
   !
-  CALL sp_aop_profile                                    ( &
+  CALL sp_aod_diag                    ( &
        klev     ,klon ,  l443         , &
-       aod550   ,aod_prof ,ssa_prof  , &
-       asy_prof )
-
-  od443aer(:)= od443aer(:)+SUM(aod_prof(:,:),dim=2)
+       aod550   ,aod_diag )
+
+  od443aer(:)= od443aer(:)+aod_diag(:)
 
   !--call to sp routine -- 865 nm
   !
-  CALL sp_aop_profile                                    ( &
+  CALL sp_aod_diag                    ( &
        klev     ,klon ,  l865         , &
-       aod550   ,aod_prof ,ssa_prof  , &
-       asy_prof )
-
-  od865aer(:)=od865aer(:)+SUM(aod_prof(:,:),dim=2)
+       aod550   ,aod_diag )
+
+  od865aer(:)= od865aer(:)+aod_diag(:)   
 
   !

LMDZ6/trunk/libf/phylmd/mo_simple_plumes.f90

@@ -80 +80 @@
        time_weight_bg (nfeatures,nplumes)        !< as time_weight but for natural background in Twomey effect 
 
-  PUBLIC sp_aod550_profile, sp_aop_profile, sp_setup
+  PUBLIC sp_aod550_profile, sp_aop_profile, sp_aod_diag, sp_setup
 
 CONTAINS
@@ -622 +622 @@
     RETURN
   END SUBROUTINE sp_aop_profile
- 
-END MODULE mo_simple_plumes
+
+  ! ------------------------------------------------------------------------------------------------------------------------
+  ! sp_aod_diag:  This subroutine for Simple Plume aerosols is a simplified version of sp_aop_prof 
+  !    to be used for diagnostic of aod at a given wavelength "lambda".
+  !    It uses aod550(ncol,nlevels,nplumes) output by the subroutine sp_aod550_profile.
+  ! ------------------------------------------------------------------------------------------------------------------------ 
+
+  SUBROUTINE sp_aod_diag                                    ( &
+       nlevels        ,ncol           ,lambda                , &
+       aod550         ,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
+
+    REAL, INTENT(OUT)          :: &
+         aod_diag(ncol)             !< aerosol optical depth, diagnostic for wavelength "lambda"
+
+    INTEGER                    :: iplume, icol, k
+
+    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
+
+    ! sum contribution from plumes to construct composite profiles of aerosol optical properties
+    !
+    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 = 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
+      ! version 1 compact
+      !aod_diag(:) = aod_diag(:) + SUM(aod550(:,:,iplume),dim=2) * 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)
+      ENDDO ! icol
+    ENDDO ! iplume
+
+    RETURN
+  END SUBROUTINE sp_aod_diag
+                                                                                                                           END MODULE mo_simple_plumes