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

Timestamp:

Apr 1, 2026, 6:36:00 PM (7 days ago)

Author:

asima

Message:

Bug fix for Simple Plume aerosols:
the orography was incorrectly accounted for in the aod vertical profile (since the original CMIP6 implementation)
ASima, FHourdin

File:

• LMDZ6/trunk/libf/phylmd/macv2sp.f90 (modified) (3 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_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
@@ -21 +21 @@
   USE aero_mod
   USE phys_local_var_mod, ONLY: t_seri, od443aer, od550aer, od865aer, ec550aer, dryod550aer, od550lt1aer
-  ! SP-specific diagnostics ; an alternative to od550SPaer_prof would be ec550SPaer = od550SPaer_prof/dz
+  ! SP-specific diagnostics (an SP extinction coefficient could be also defined: ec550SPaer = od550SPaer_prof/dz)
   USE phys_local_var_mod, ONLY: od550SPaer, od550SPaer_prof 
   !
@@ -107 +107 @@
   !--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   )
+  CALL sp_aod550_profile                                            ( &
+       klev           ,klon       ,xlon       ,xlat      , &
+       decimal_year   ,z          ,dz         ,dNovrN    , aod550   )
 
   !--AOD550 calculations for diagnostics