Changeset 6128 for LMDZ6/trunk/libf/phylmd/macv2sp.f90

Timestamp:

Mar 26, 2026, 7:09:02 PM (9 days ago)

Author:

asima

Message:

Major changes (don't change results):
Subroutine sp_aop_profile in the mo_simple_plume.f90 module is split in 2 subroutines:

• sp_aod550_profile, for calculations at 550um ; only called in macv2sp once per timestep ;
• sp_aop_profile only keeps the calculations of profiles of optical properties at other wavelengths ;

it is called in macv2sp for aod diagnostics at 443 and 865um,
and for each of the 25 wavelengths filling the 6 RRTM bands.

Minor :

• subroutine name MACv2SP changed in macv2sp for coherence with file name ;
• dNovrN becomes an argument INTENT(OUT) of macv2sp (instead of making USE from phys_local_var_mod) ;
• lfactor variable name was used with 2 different meanings : eq 10 and eq 12 from Stevens et al. (2017); celui from eq(10) was renamed as 'lextinct'
• parameter nmon = 12 : removed because not used ;
• variable aod_lmd : renamed as aod_lmdz.

File:

• LMDZ6/trunk/libf/phylmd/macv2sp.f90 (modified) (9 diffs)

LMDZ6/trunk/libf/phylmd/macv2sp.f90

@@ -1 @@
-SUBROUTINE MACv2SP(pphis,pplay,paprs,xlon,xlat,tau_allaer,piz_allaer,cg_allaer)
+SUBROUTINE macv2sp(pphis,pplay,paprs,xlon,xlat,tau_allaer,piz_allaer,cg_allaer,dNovrN)
   !
   !--routine to read the MACv2SP plume and compute optical properties
@@ -13 +13 @@
   !--dNovrN   = enhancement factor for CDNC
   !
-  USE mo_simple_plumes, ONLY: sp_aop_profile, sp_setup, sp_initialized, aerosols_SP_forcing_year
+  USE mo_simple_plumes, ONLY: sp_aod550_profile, sp_aop_profile, 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
   USE dimphy
   USE aero_mod
-  USE phys_local_var_mod, ONLY: t_seri, od443aer, od550aer, od865aer, ec550aer, dryod550aer, od550lt1aer, od550SPaer, dNovrN
-  !!USE YOMCST, ONLY : RD, RG
+  USE phys_local_var_mod, ONLY: t_seri, od443aer, od550aer, od865aer, ec550aer, dryod550aer, od550lt1aer, od550SPaer
   !
   USE yomcst_mod_h
@@ -35 +35 @@
   REAL, DIMENSION(klon,klev,2,nbands_sw_rrtm), INTENT(INOUT) :: cg_allaer  !  asymmetry parameter aerosol
   !
+  REAL,DIMENSION(klon),INTENT(OUT)       :: dNovrN   !< anthropogenic increase in cloud drop number concentration (factor)
+  !
+  REAL,DIMENSION(klon,klev,nplumes) :: aod550
   REAL,DIMENSION(klon,klev) :: aod_prof, ssa_prof, asy_prof
   REAL,DIMENSION(klon,klev) :: z, dz
   REAL,DIMENSION(klon)      :: oro, zrho, zt
   !
-  INTEGER, PARAMETER :: nmon = 12
-  !
-  REAL, PARAMETER    :: l443 = 443.0, l550 = 550.0, l865 = 865.0 !--wavelengths in nm
+  REAL, PARAMETER    :: l443 = 443.0, l865 = 865.0 !--wavelengths in nm for diagnostics (other than 550 nm calculated by default)
   !
   INTEGER, PARAMETER :: Nwvmax=25
@@ -58 +59 @@
   !
   REAL :: zlambda, zweight
-!  REAL :: year_fr  ! also a dimension name in SP aerosol file ; renamed more properly "decimal_year" 
   REAL :: decimal_year 
   !
@@ -83 +83 @@
   IF (.NOT.sp_initialized) CALL sp_setup
 
-  !--fractional year
-  !
-  ! Original version, bugged :
-  ! year_fr = FLOAT(year_cur) + (FLOAT(day_cur)-0.5) / FLOAT(year_len)
-  ! Correction ASima & FH :
-!  year_fr= float(year_cur) + float(days_elapsed)/float(year_len)
-! Choice between yearly vs fixed_year forcing, depending on 'aerosols_SP_forcing_year'
+  !--fractional year : 
+  !   Original version, bugged :       year_fr = FLOAT(year_cur) + (FLOAT(day_cur)-0.5) / FLOAT(year_len)
+  !   Correction 2026-03, FH & ASima : year_fr= float(year_cur) + float(days_elapsed)/float(year_len)
+  ! Name changed in decimal_year ; choice between yearly vs fixed_year forcing, depending on 'aerosols_SP_forcing_year'
   IF (aerosols_SP_forcing_year.EQ.-9999) THEN
      decimal_year= float(year_cur) + float(days_elapsed)/float(year_len)
@@ -99 +96 @@
      CALL abort_physic ('macv2sp','year not supported by plume model',1)
   ENDIF
-  !
-  !--call to sp routine -- 443 nm
-  !
-  CALL sp_aop_profile                                    ( &
-       klev     ,klon ,l443 ,oro    ,xlon     ,xlat      , &
-       decimal_year  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
-       asy_prof )
-  !
-  !--AOD calculations for diagnostics
-  od443aer(:)= od443aer(:)+SUM(aod_prof(:,:),dim=2)
-  !
-  !--call to sp routine -- 550 nm
-  !
-  CALL sp_aop_profile                                    ( &
-       klev     ,klon ,l550 ,oro    ,xlon     ,xlat      , &
-       decimal_year  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
-       asy_prof )
-  !
-  !--AOD calculations for diagnostics
-  ! (ASima : new diagnostic od550SPaer; corrected od550lt1aer and dryod550aer)
-  !--a/ AOD of SP aerosols = vertical sum of SP aod profile
+
+  !
+  !--call SUBROUTINE sp_aod550_profile once ; all the wavelength-dependent profiles and 2D diagnostics will use aod550 profile
+  !
+  CALL sp_aod550_profile                                                 ( &
+       klev           ,klon       ,oro        ,xlon           ,xlat      , &
+       decimal_year   ,z          ,dz         ,dNovrN         , aod550   )
+
+  !--AOD550 calculations for diagnostics
+  !
+  ! aod550 profile = sum over the 'nplumes' dimension
+  aod_prof(:,:)=SUM(aod550(:,:,:),dim=3)
+  !
+  !--a/ AOD of SP anthropic aerosols = vertical sum of SP aod profile
   od550SPaer(:)=SUM(aod_prof(:,:),dim=2)
   !
@@ -126 +116 @@
   !
   !--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
   ec550aer(:,:)=ec550aer(:,:)+aod_prof(:,:)/dz(:,:)
   !
+
+  ! Diagnostic AOD calculations for other wavelengths : 443 and 865 nm  
+  !--call to sp routine -- 443 nm
+  !
+  CALL sp_aop_profile                                    ( &
+       klev     ,klon ,  l443         , &
+       aod550   ,aod_prof ,ssa_prof  , &
+       asy_prof )
+
+  od443aer(:)= od443aer(:)+SUM(aod_prof(:,:),dim=2)
+
   !--call to sp routine -- 865 nm
   !
   CALL sp_aop_profile                                    ( &
-       klev     ,klon ,l865 ,oro    ,xlon     ,xlat      , &
-       decimal_year  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
+       klev     ,klon ,  l865         , &
+       aod550   ,aod_prof ,ssa_prof  , &
        asy_prof )
-  !
-  !--AOD calculations for diagnostics
+
   od865aer(:)=od865aer(:)+SUM(aod_prof(:,:),dim=2)
+
   !
   !--re-weighting of piz and cg arrays before adding the anthropogenic aerosols
@@ -184 +180 @@
     ENDIF
     !
-    CALL sp_aop_profile                                       ( &
-         klev     ,klon ,zlambda ,oro    ,xlon     ,xlat      , &
-         decimal_year  ,z    ,dz      ,dNovrN ,aod_prof ,ssa_prof  , &
-         asy_prof )
+  CALL sp_aop_profile                                    ( &
+       klev     ,klon ,  zlambda         , &
+       aod550   ,aod_prof ,ssa_prof  , &
+       asy_prof )
+
+
     !
     !--adding up the quantities tau, piz*tau and cg*piz*tau
@@ -200 +198 @@
   piz_allaer(:,:,2,:)=piz_allaer(:,:,2,:)/tau_allaer(:,:,2,:)
   !
-END SUBROUTINE MACv2SP
+END SUBROUTINE macv2sp