Changeset 2803 for LMDZ5/trunk/libf/phylmd/rrtm/read_rsun_rrtm.F90

Timestamp:

Feb 28, 2017, 10:05:52 PM (8 years ago)

Author:

oboucher

Message:

Correcting bug on current day in the year in RSUN forcing

File:

• LMDZ5/trunk/libf/phylmd/rrtm/read_rsun_rrtm.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/rrtm/read_rsun_rrtm.F90

@@ -10 +10 @@
   USE netcdf, ONLY: nf90_get_var, nf90_noerr, nf90_nowrite
 
-  USE phys_cal_mod, ONLY : day_cur, year_len
+  USE phys_cal_mod, ONLY : days_elapsed, year_len
 
   USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root
@@ -23 +23 @@
 
   ! Input arguments
-  LOGICAL, INTENT(IN)          :: debut
+  LOGICAL, INTENT(IN) :: debut
 
 ! Local variables
-  INTEGER               :: ncid, dimid, varid, ncerr, nbday
+  INTEGER :: ncid, dimid, varid, ncerr, nbday
   REAL, POINTER :: wlen(:), time(:)
-  REAL, ALLOCATABLE, SAVE, DIMENSION(:,:)   :: SSI_FRAC
+  REAL, ALLOCATABLE, SAVE, DIMENSION(:,:) :: SSI_FRAC
 !$OMP THREADPRIVATE(SSI_FRAC)
   REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: TSI(:)
@@ -90 +90 @@
     ENDIF
 
-!--only read at beginning of month
-    IF (debut.OR.day_cur.NE.day_pre) THEN
+!--only read at beginning of day
+!--day in year is provided as days_elapsed since the beginning of the year +1
+    IF (debut.OR.days_elapsed+1.NE.day_pre) THEN
 
-!--keep memory of previous month
-      day_pre=day_cur
+!--keep memory of previous day
+      day_pre=days_elapsed+1
 
 !--copy 
-      RSUN(1:NSW)=SSI_FRAC(:,day_cur)
-      solaire=TSI(day_cur)
+      RSUN(1:NSW)=SSI_FRAC(:,days_elapsed+1)
+      solaire=TSI(days_elapsed+1)
 
-      print *,'READ_RSUN_RRTM day=', day_cur,' solaire=', solaire, ' RSUN=', RSUN(1:NSW)
+      print *,'READ_RSUN_RRTM day=', days_elapsed+1,' solaire=', solaire, ' RSUN=', RSUN(1:NSW)
 
     ENDIF !--fin allocation