Index: LMDZ6/trunk/libf/phylmd/macv2sp.f90
===================================================================
--- LMDZ6/trunk/libf/phylmd/macv2sp.f90	(revision 6100)
+++ LMDZ6/trunk/libf/phylmd/macv2sp.f90	(revision 6102)
@@ -58,5 +58,6 @@
   !
   REAL :: zlambda, zweight
-  REAL :: year_fr
+!  REAL :: year_fr  ! also a dimension name in SP aerosol file ; renamed more properly "decimal_year" 
+  REAL :: decimal_year 
   !
   INTEGER band, i, k, Nwv
@@ -79,5 +80,5 @@
   !
 
-  ! Call sp_setup to retrieve 'aerosol_SP_forcing_year' value used in year_fr calculation below
+  ! Call sp_setup to retrieve 'aerosol_SP_forcing_year' value used in decimal_year calculation below
   IF (.NOT.sp_initialized) CALL sp_setup
 
@@ -90,10 +91,10 @@
 ! 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)
+     decimal_year= float(year_cur) + float(days_elapsed)/float(year_len)
   ELSE
-     year_fr= float(aerosols_SP_forcing_year) + float(days_elapsed)/float(year_len)
+     decimal_year= 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
+  IF (decimal_year.LT.1850.0.OR.decimal_year.GE.2024.0) THEN
      CALL abort_physic ('macv2sp','year not supported by plume model',1)
   ENDIF
@@ -103,5 +104,5 @@
   CALL sp_aop_profile                                    ( &
        klev     ,klon ,l443 ,oro    ,xlon     ,xlat      , &
-       year_fr  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
+       decimal_year  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
        asy_prof )
   !
@@ -113,5 +114,5 @@
   CALL sp_aop_profile                                    ( &
        klev     ,klon ,l550 ,oro    ,xlon     ,xlat      , &
-       year_fr  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
+       decimal_year  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
        asy_prof )
   !
@@ -142,5 +143,5 @@
   CALL sp_aop_profile                                    ( &
        klev     ,klon ,l865 ,oro    ,xlon     ,xlat      , &
-       year_fr  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
+       decimal_year  ,z    ,dz   ,dNovrN ,aod_prof ,ssa_prof  , &
        asy_prof )
   !
@@ -185,5 +186,5 @@
     CALL sp_aop_profile                                       ( &
          klev     ,klon ,zlambda ,oro    ,xlon     ,xlat      , &
-         year_fr  ,z    ,dz      ,dNovrN ,aod_prof ,ssa_prof  , &
+         decimal_year  ,z    ,dz      ,dNovrN ,aod_prof ,ssa_prof  , &
          asy_prof )
     !
Index: LMDZ6/trunk/libf/phylmd/mo_simple_plumes.f90
===================================================================
--- LMDZ6/trunk/libf/phylmd/mo_simple_plumes.f90	(revision 6100)
+++ LMDZ6/trunk/libf/phylmd/mo_simple_plumes.f90	(revision 6102)
@@ -290,10 +290,10 @@
   ! week in the year and superimposed on the yearly mean value of the weight. 
   !
-  SUBROUTINE set_time_weight(year_fr)
+  SUBROUTINE set_time_weight(decimal_year)
     !
     ! ---------- 
     !
     REAL, INTENT(IN) ::  &
-         year_fr           !< Fractional Year (1850.0 - 2100.99)
+         decimal_year      !< Fractional Year (1850.0 - 2100.99)
 
     INTEGER          ::  &
@@ -304,6 +304,6 @@
     ! ---------- 
     !
-    iyear = FLOOR(year_fr) - 1849
-    iweek = FLOOR((year_fr - FLOOR(year_fr)) * ntimes) + 1
+    iyear = FLOOR(decimal_year) - 1849
+    iweek = FLOOR((decimal_year - FLOOR(decimal_year)) * ntimes) + 1
 
     IF ((iweek > ntimes) .OR. (iweek < 1) .OR. (iyear > nyears) .OR. (iyear < 1)) THEN 
@@ -328,5 +328,5 @@
   SUBROUTINE sp_aop_profile                                                                           ( &
        nlevels        ,ncol           ,lambda         ,oro            ,lon            ,lat            , &
-       year_fr        ,z              ,dz             ,dNovrN         ,aod_prof       ,ssa_prof       , &
+       decimal_year        ,z              ,dz             ,dNovrN         ,aod_prof       ,ssa_prof       , &
        asy_prof       )
     !
@@ -339,5 +339,5 @@
     REAL, INTENT(IN)           :: &
          lambda,                  & !< wavelength
-         year_fr,                 & !< Fractional Year (1903.0 is the 0Z on the first of January 1903, Gregorian)
+         decimal_year,            & !< Fractional Year (1903.0 is the 0Z on the first of January 1903, Gregorian)
          oro(ncol),               & !< orographic height (m)
          lon(ncol),               & !< longitude 
@@ -393,5 +393,5 @@
     ! get time weights
     !
-    CALL set_time_weight(year_fr)
+    CALL set_time_weight(decimal_year)
     !
     ! initialize variables, including output