Changeset 1273 for LMDZ4/branches/LMDZ4-dev/libf/phylmd

Timestamp:

Dec 4, 2009, 2:23:25 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Correction d'un bug sur le calcul de la date courante pour la lecture
des aerosols
Attention, les resultats sont modifies suite a un decalage d'une journee
dans l'ancien calcul de la date
LF

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90 (modified) (4 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol_interp.F90

@@ -105 +105 @@
 
   ! 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 +241 @@
      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 +297 @@
        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 +501 @@
                          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