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

Timestamp:

Mar 12, 2026, 2:37:31 PM (13 days ago)

Author:

asima

Message:

Three parameters can be now specified via config.def :

1/ aerosols_SP_forcing_year :

-9999 --> "historical" mode : SP forcing year is the simulation year (between 1850 and 2023 included)

YYYY (a given year between 1850 and 2023) : SP forcing is fixed to year YYYY for the entire simulation

2/ aerosols_SP_coef_bN

coef bN in eq (15) of SP paper Stevens et al (2017);

Published values : 1000. (default) and 20.

3/ aerosols_SP_aod_bg_gl

globally constant part of the natural background aod in SP implementation ;

Default : 0.02

File:

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

LMDZ6/trunk/libf/phylmd/macv2sp.f90

@@ -13 +13 @@
   !--dNovrN   = enhancement factor for CDNC
   !
-  USE mo_simple_plumes, ONLY: sp_aop_profile
+  USE mo_simple_plumes, ONLY: sp_aop_profile, sp_setup, sp_initialized, aerosols_SP_forcing_year
   USE phys_cal_mod, ONLY : year_cur, year_len, days_elapsed
   USE dimphy
@@ -78 +78 @@
   ENDDO
   !
+
+  ! Call sp_setup to retrieve 'aerosol_SP_forcing_year' value used in year_fr calculation below
+  IF (.NOT.sp_initialized) CALL sp_setup
+
   !--fractional year
   !
@@ -83 +87 @@
   ! 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)
+!  year_fr= float(year_cur) + float(days_elapsed)/float(year_len)
+! Choice between yearly vs fixed_year forcing, depending on 'aerosols_SP_forcing_year'
+  IF (aerosols_SP_forcing_year.EQ.-9999) THEN
+     year_fr= float(year_cur) + float(days_elapsed)/float(year_len)
+  ELSE
+     year_fr= float(aerosols_SP_forcing_year) + float(days_elapsed)/float(year_len)
+  ENDIF
 
   IF (year_fr.LT.1850.0.OR.year_fr.GE.2024.0) THEN