subroutine geth_newdate (ndate, odate, idts) implicit none !********************************************************************** ! ! purpose - from old date ('YYYY-MM-DD*HH:MM:SS') and time in ! seconds, compute the new date. ! ! on entry - odate - the old hdate. ! idts - the change in time in seconds. ! ! on exit - ndate - the new hdate. ! idts - the change in time in seconds. ! !********************************************************************** integer idts character*(*) ndate, odate integer nlen, olen ! ! Local Variables ! ! yrold - indicates the year associated with "odate" ! moold - indicates the month associated with "odate" ! dyold - indicates the day associated with "odate" ! hrold - indicates the hour associated with "odate" ! miold - indicates the minute associated with "odate" ! scold - indicates the second associated with "odate" ! ! yrnew - indicates the year associated with "ndate" ! monew - indicates the month associated with "ndate" ! dynew - indicates the day associated with "ndate" ! hrnew - indicates the hour associated with "ndate" ! minew - indicates the minute associated with "ndate" ! scnew - indicates the second associated with "ndate" ! ! mday - a list assigning the number of days in each month ! dth - the number of hours represented by "idts" ! i - loop counter ! nday - the integer number of days represented by "idts" ! nhour - the integer number of hours in "idts" after taking out ! all the whole days ! nmin - the integer number of minutes in "idts" after taking out ! all the whole days and whole hours. ! nsec - the integer number of minutes in "idts" after taking out ! all the whole days, whole hours, and whole minutes. ! integer yrnew, monew, dynew, hrnew, minew, scnew integer yrold, moold, dyold, hrold, miold, scold integer mday(12), nday, nhour, nmin, nsec, i real dth logical opass logical noLeapCalendar !************************* Subroutine Begin ************************* ! ! Determine whether the routine should use a no-leap year calendar ! noLeapCalendar = .false. ! ! Assign the number of days in a months ! mday( 1) = 31 mday( 2) = 28 mday( 3) = 31 mday( 4) = 30 mday( 5) = 31 mday( 6) = 30 mday( 7) = 31 mday( 8) = 31 mday( 9) = 30 mday(10) = 31 mday(11) = 30 mday(12) = 31 ! ! Break down old hdate into parts ! hrold = 0 miold = 0 scold = 0 olen = len(odate) read(odate(1:4), '(I4)') yrold read(odate(6:7), '(I2)') moold read(odate(9:10), '(I2)') dyold if (olen.ge.13) then read(odate(12:13),'(I2)') hrold if (olen.ge.16) then read(odate(15:16),'(I2)') miold if (olen.ge.19) then read(odate(18:19),'(I2)') scold endif endif endif ! ! Set the number of days in February for that year. ! mday(2) = 28 if (.not. noLeapCalendar) then if (mod(yrold,4).eq.0) then mday(2) = 29 if (mod(yrold,100).eq.0) then mday(2) = 28 if (mod(yrold,400).eq.0) then mday(2) = 29 endif endif endif endif ! ! Check that ODATE makes sense. ! opass = .TRUE. ! Check that the month of ODATE makes sense. if ((moold.gt.12).or.(moold.lt.1)) then print*, 'GETH_NEWDATE: Month of ODATE = ', moold opass = .FALSE. endif ! Check that the day of ODATE makes sense. if ((dyold.gt.mday(moold)).or.(dyold.lt.1)) then ! ! Send a custom message if any leap-day files if it is a No-Leap Calendar ! if (noLeapCalendar .and. (moold .eq. 2) .and. (dyold .eq. 29)) then print*,'GET_NEWDATE: Using a no-Leap Calendar, but data for 2/29 was found.' else print*, 'GET_NEWDATE: Day of ODATE = ', dyold endif opass = .FALSE. endif ! Check that the hour of ODATE makes sense. if ((hrold.gt.23).or.(hrold.lt.0)) then print*, 'GET_NEWDATE: Hour of ODATE = ', hrold opass = .FALSE. endif ! Check that the minute of ODATE makes sense. if ((miold.gt.59).or.(miold.lt.0)) then print*, 'GET_NEWDATE: Minute of ODATE = ', miold opass = .FALSE. endif ! Check that the second of ODATE makes sense. if ((scold.gt.59).or.(scold.lt.0)) then print*, 'GET_NEWDATE: Second of ODATE = ', scold opass = .FALSE. endif if (.not.opass) then print*, 'Crazy ODATE: ', odate(1:olen), olen STOP 'Error_odate' ! stop endif ! ! Date Checks are completed. Continue. ! ! ! Compute the number of days, hours, minutes, and seconds in idts ! nday = idts/86400 ! Integer number of days in delta-time nhour = mod(idts,86400)/3600 nmin = mod(idts,3600)/60 nsec = mod(idts,60) scnew = scold + nsec if (scnew .ge. 60) then scnew = scnew - 60 nmin = nmin + 1 end if minew = miold + nmin if (minew .ge. 60) then minew = minew - 60 nhour = nhour + 1 end if hrnew = hrold + nhour if (hrnew .ge. 24) then hrnew = hrnew - 24 nday = nday + 1 end if dynew = dyold monew = moold yrnew = yrold do i = 1, nday dynew = dynew + 1 if (dynew.gt.mday(monew)) then dynew = dynew - mday(monew) monew = monew + 1 if (monew .gt. 12) then monew = 1 yrnew = yrnew + 1 mday(2) = 28 if (.not. noLeapCalendar) then if (mod(yrnew,4).eq.0) then mday(2) = 29 if (mod(yrnew,100).eq.0) then mday(2) = 28 if (mod(yrnew,400).eq.0) then mday(2) = 29 endif endif endif endif end if endif enddo ! ! Now construct the new mdate ! nlen = len(ndate) if (nlen.ge.19) then write(ndate,19) yrnew, monew, dynew, hrnew, minew, scnew 19 format(I4,'-',I2.2,'-',I2.2,'_',I2.2,':',I2.2,':',I2.2) else if (nlen.eq.16) then write(ndate,16) yrnew, monew, dynew, hrnew, minew 16 format(I4,'-',I2.2,'-',I2.2,'_',I2.2,':',I2.2) else if (nlen.eq.13) then write(ndate,13) yrnew, monew, dynew, hrnew 13 format(I4,'-',I2.2,'-',I2.2,'_',I2.2) else if (nlen.eq.10) then write(ndate,10) yrnew, monew, dynew 10 format(I4,'-',I2.2,'-',I2.2) endif !************************** Subroutine End ************************** end