- Timestamp:
- Jan 20, 2010, 3:27:21 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/readaerosol_interp.F90
r1279 r1299 126 126 IF(mpi_rank == 0 .AND. debug)then 127 127 ! 0.02 is about 0.5/24, namly less than half an hour 128 OLDNEWDAY = (r_day- FLOAT(iday) < 0.02)128 OLDNEWDAY = (r_day-REAL(iday) < 0.02) 129 129 ! Once per day, update aerosol fields 130 130 lmt_pas = NINT(86400./pdtphys) 131 PRINT*,'r_day- FLOAT(iday) =',r_day-FLOAT(iday)131 PRINT*,'r_day-REAL(iday) =',r_day-REAL(iday) 132 132 PRINT*,'itap =',itap 133 133 PRINT*,'pdtphys =',pdtphys … … 233 233 ! 234 234 DO i = 2, 13 235 month_len(i) = float(ioget_mon_len(year_cur, i-1))235 month_len(i) = REAL(ioget_mon_len(year_cur, i-1)) 236 236 CALL ymds2ju(year_cur, i-1, 1, 0.0, month_start(i)) 237 237 ENDDO 238 month_len(1) = float(ioget_mon_len(year_cur-1, 12))238 month_len(1) = REAL(ioget_mon_len(year_cur-1, 12)) 239 239 CALL ymds2ju(year_cur-1, 12, 1, 0.0, month_start(1)) 240 month_len(14) = float(ioget_mon_len(year_cur+1, 1))240 month_len(14) = REAL(ioget_mon_len(year_cur+1, 1)) 241 241 CALL ymds2ju(year_cur+1, 1, 1, 0.0, month_start(14)) 242 242 month_mid(:) = month_start (:) + month_len(:)/2.
Note: See TracChangeset
for help on using the changeset viewer.