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

Timestamp:

Mar 5, 2026, 4:37:59 PM (7 days ago)

Author:

asima

Message:

Concerning "year_fr" :

• (related to r6087) : use “float” for “integer” variables “year_cur” and “days_elapsed” ;
• extend the range of available year_weight values to 2023, as in SPv2.1_1850-2023_CMIP7.nc

File:

• LMDZ6/trunk/libf/phylmd/macv2sp.f90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/macv2sp.f90

@@ -80 +80 @@
   !--fractional year
   !
-  ! Adriana Sima et FH : Version bugguée :
+  ! Original version, bugged :
   ! year_fr = FLOAT(year_cur) + (FLOAT(day_cur)-0.5) / FLOAT(year_len)
-  year_fr=year_cur+days_elapsed/float(year_len)
+  ! Correction ASima & FH :
+  year_fr= float(year_cur) + float(days_elapsed)/float(year_len)
 
-  IF (year_fr.LT.1850.0.OR.year_fr.GE.2021.0) THEN
+  IF (year_fr.LT.1850.0.OR.year_fr.GE.2024.0) THEN
      CALL abort_physic ('macv2sp','year not supported by plume model',1)
   ENDIF