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

Timestamp:

Dec 7, 2016, 7:27:29 PM (8 years ago)

Author:

oboucher

Message:

Correction to be more flexible on the rsun solar forcing in case of noleap calendar

File:

• LMDZ5/trunk/libf/phylmd/rrtm/read_rsun_rrtm.F90 (modified) (1 diff)

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

@@ -56 +56 @@
        ENDIF
 
-       IF (size(time).NE.year_len) THEN 
+!--test if time is different from year_len but allow a mismatch of 1 day
+       IF (size(time).NE.year_len.AND.size(time).NE.year_len+1) THEN
          PRINT *,'read_rsun_rrtm time <> year_len = ', size(time), year_len
          CALL abort_physic('read_rsun_rrtm','time dim should be the number of days in year',1)
+       ENDIF
+!--warning only if forcing file has 366 days but year_len has only 365
+       IF (size(time).EQ.year_len+1) THEN 
+         PRINT *,'Warning read_rsun_rrtm uses a leap year rsun for a noleap year'
        ENDIF