Index: LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90
===================================================================
--- LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90	(revision 1272)
+++ LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90	(revision 1273)
@@ -105,14 +105,16 @@
 
   ! Use phys_cal_mod
-  !iday= day_cur
-  !iyr = year_cur
-  !im  = mth_cur
-
-  iday = INT(r_day)
-  iyr  = iday/360
-  iday = iday-iyr*360         ! day of the actual year
-  iyr  = iyr + annee_ref      ! year of the run   
-  im   = iday/30 +1           ! the actual month
+  iday= day_cur
+  iyr = year_cur
+  im  = mth_cur
+
+!  iday = INT(r_day)
+!  iyr  = iday/360
+!  iday = iday-iyr*360         ! day of the actual year
+!  iyr  = iyr + annee_ref      ! year of the run   
+!  im   = iday/30 +1           ! the actual month
   CALL ymds2ju(iyr, im, iday, 0., jDay)
+!   CALL ymds2ju(iyr, im, iday-(im-1)*30, 0., jDay)
+
 
   IF(MOD(itap-1,NINT(86400./pdtphys)) == 0)THEN
@@ -239,4 +241,9 @@
      CALL ymds2ju(year_cur+1, 1, 1, 0.0, month_start(14))
      month_mid(:) = month_start (:) + month_len(:)/2.
+     
+     if (debug) then
+       write(lunout,*)' month_len = ',month_len
+       write(lunout,*)' month_mid = ',month_mid
+     endif
 
   END IF  ! IF ( (first .OR. iday==0) .AND. lnewday ) THEN 
@@ -290,4 +297,7 @@
        CALL abort_gcm('readaerosol_interp', 'number of months undefined',1)
      ENDIF
+     if (debug) then
+       write(lunout,*)' jDay, day1, day2, im, im2 = ', jDay, day1, day2, im, im2
+     endif
 
  
@@ -491,6 +501,6 @@
                          var_year(i,k,im2,id_aero) - var_year(i,k,im,id_aero)
                  END IF
-                 
                  WRITE(lunout,*) 'stop for aerosol : ',name_aero(id_aero)
+                 WRITE(lunout,*) 'day1, day2, jDay = ', day1, day2, jDay 
                  CALL abort_gcm('readaerosol_interp','Error in interpolation 1',1)
               END IF