!\begin{verbatim}!===============================================================================! MFM_Date F90 Unit Tests and Examples!===============================================================================      program main       use MFM_TimeMgmt      implicit none      integer, parameter :: START_DATE=20011128, START_SECS=43200      integer, parameter :: STOP_DATE=20041201, STOP_SECS=1200      integer, parameter :: DAY_INC=6, SEC_INC=43200      integer :: retCalDate, retDays, retSecs      logical :: test, isLater      type(MFM_Date) :: startDateG, stopDateG, retDateG      type(MFM_Date) :: startDateN, stopDateN, retDateN      type(MFM_Date) :: copy      type(MFM_Time) :: incTime, retTime      type(MFM_Date) :: calDayTester      real(8) :: floatDay, floatTest      character(60) :: str      print *, "=================================================="      print *, "MFM_Date F90 Unit Tests and Examples"      print *, "=================================================="      startDateG = MFM_DateInit(MFM_GREGORIAN, START_DATE, START_SECS)       stopDateG = MFM_DateInit(MFM_GREGORIAN, STOP_DATE, STOP_SECS)      retDateG = MFM_DateInit()       startDateN = MFM_DateInit(MFM_NO_LEAP, START_DATE, START_SECS)       stopDateN = MFM_DateInit(MFM_NO_LEAP, STOP_DATE, STOP_SECS)      retDateN = MFM_DateInit()       copy = MFM_DateCopyInit(stopDateG)      incTime = MFM_TimeInit(DAY_INC, SEC_INC)       retTime = MFM_TimeInit()      call MFM_DateGet(startDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs       test = ((retCalDate==START_DATE) .AND. (retSecs==START_SECS))      str = "MFM_DateInit, MFM_DateGet:  init date and get attr"      call MFM_ErrorTest(test, str)      call MFM_DateGet(copy, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs       test = ((retCalDate==STOP_DATE) .AND. (retSecs==STOP_SECS))      str = "MFM_CopyInit:  init date and get attr"      call MFM_ErrorTest(test, str)      call MFM_DateCopy(copy, startDateG)      call MFM_DateGet(copy, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs       test = ((retCalDate==START_DATE) .AND. (retSecs==START_SECS))      str = "MFM_DateCopy:  copy date"      call MFM_ErrorTest(test, str)      retDateG = MFM_DateIncrementSec(startDateG, SEC_INC)      call MFM_DateGet(retDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs       test = ((retCalDate==20011129) .AND. (retSecs==0))      str = "MFM_DateIncrementSec:  increment Gregorian date by seconds"        call MFM_ErrorTest(test, str)      retDateG = MFM_DateIncrementDay(startDateG, DAY_INC)      call MFM_DateGet(retDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs       test = ((retCalDate==20011204) .AND. (retSecs==43200))      str = "MFM_DateIncrementDay:  increment Gregorian date by days"        call MFM_ErrorTest(test, str)      retDateG = MFM_DateIncrement(startDateG, incTime)      call MFM_DateGet(retDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20011205) .AND. (retSecs==0))      str = "MFM_DateIncrement:  increment Gregorian date"      call MFM_ErrorTest(test, str)      retDateG = MFM_DateDecrement(retDateG, incTime)      call MFM_DateGet(retDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20011128) .AND. (retSecs==43200))      str = "MFM_DateDecrement:  decrement Gregorian date in place"      call MFM_ErrorTest(test, str)      call MFM_DateDiff(startDateG, stopDateG, retTime, isLater)      call MFM_TimeGet(retTime, retDays, retSecs)      print *, "Ret days ", retDays, " Ret secs ",      &         retSecs, "isLater", isLater       test = ((retDays==1098) .AND. (retSecs==44400)      &       .AND. (isLater))      print *, "test", test      str = "MFM_DateDiff:  diff of two Gregorian dates"      call MFM_ErrorTest(test, str)      call MFM_TimeSet(incTime, 1098, 44400)      retDateG = MFM_DateDecrement(stopDateG, incTime)      call MFM_DateGet(retDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20011128) .AND. (retSecs==43200))      str = "MFM_DateDecrement:  decrement Gregorian date"      call MFM_ErrorTest(test, str)      call MFM_TimeSet(incTime, DAY_INC, SEC_INC)      retDateN = MFM_DateIncrement(startDateN, incTime)      call MFM_DateGet(retDateN, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20011205) .AND. (retSecs==0))      str = "MFM_DateIncrement:  increment no leap date"      call MFM_ErrorTest(test, str)      retDateN = MFM_DateDecrement(retDateN, incTime)      call MFM_DateGet(retDateN, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20011128) .AND. (retSecs==43200))      str = "MFM_DateDecrement:  decrement no leap date in place"      call MFM_ErrorTest(test, str)      call MFM_DateDiff(startDateN, stopDateN, retTime, isLater)      call MFM_TimeGet(retTime, retDays, retSecs)      print *, "Ret days ", retDays, " Ret secs ", retSecs,     &         "isLater", isLater      test = ((retDays==1097) .AND. (retSecs==44400)      &       .AND. (isLater))      str = "MFM_DateDiff:  diff of two no leap dates"      call MFM_ErrorTest(test, str)      call MFM_DateSet(startDateN, MFM_NO_LEAP, 19700101, 0)       call MFM_DateSet(stopDateN, MFM_NO_LEAP, 19710101, 0)      call MFM_DateDiff(startDateN, stopDateN, retTime, isLater)      call MFM_TimeGet(retTime, retDays, retSecs)      print *, "Ret days ", retDays, " Ret secs ", retSecs,     &         "isLater", isLater      test = ((retDays==365) .AND. (retSecs==0)      &       .AND. (isLater))      str = "MFM_DateDiff:  diff of two no leap dates"      call MFM_ErrorTest(test, str)      call MFM_DateSet(startDateN, MFM_NO_LEAP, 20011128, 43200)       call MFM_DateSet(stopDateN, MFM_NO_LEAP, 20041201, 1200)      call MFM_TimeSet(incTime, 1097, 44400)      retDateN = MFM_DateDecrement(stopDateN, incTime)      call MFM_DateGet(retDateN, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20011128) .AND. (retSecs==43200))      str = "MFM_DateDecrement:  decrement no leap date"      call MFM_ErrorTest(test, str)      call MFM_DateSet(retDateN, MFM_NO_LEAP, 20031224, 3500)      call MFM_DateGet(retDateN, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20031224) .AND. (retSecs==3500))      str = "MFM_DateSet:  set a no leap date"      call MFM_ErrorTest(test, str)                  call MFM_DateIsLater(startDateG, stopDateG, isLater)      print *, "isLater", isLater      test = isLater      str = "MFM_DateIsLater:  compare Gregorian dates (true result)"      call MFM_ErrorTest(test, str)      call MFM_DateIsLater(stopDateG, startDateG, isLater)      test = .NOT. isLater      str = "MFM_DateIsLater:  compare Gregorian dates (false result)"      call MFM_ErrorTest(test, str)      call MFM_DateSet(startDateG, MFM_GREGORIAN, 20011031, 0)      retDateG = MFM_DateIncrementMonth(startDateG, 4)      call MFM_DateGet(retDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20020228) .AND. (retSecs==0))      str = "MFM_DateIncrementMonth:  handle invalid day"      call MFM_ErrorTest(test, str)            call MFM_DateSet(startDateG, MFM_GREGORIAN, 20040229, 0)      retDateG = MFM_DateIncrementYear(startDateG, 1)      call MFM_DateGet(retDateG, retCalDate, retSecs)      print *, "Ret date ", retCalDate, " Ret secs ", retSecs      test = ((retCalDate==20050228) .AND. (retSecs==0))      str = "MFM_DateIncrementYear:  handle invalid day"      call MFM_ErrorTest(test, str)            call MFM_DatePrint(startDateG)      calDayTester = MFM_DateInit(MFM_GREGORIAN, START_DATE, START_SECS)      print *, "Calendar Day Test"      call MFM_DatePrint(calDayTester)      floatDay = MFM_DateGetFltDayOfYear(calDayTester)      call MFM_DateGet(calDayTester, retCalDate, retSecs)      retDays = MFM_DateGetDayOfYear(calDayTester)      print *, "float days = ", floatDay      print *, "floor(floatDay) = ", floor(floatDay)      print *, "secs=", (floatDay - floor(floatDay))      floatTest = retDays + (retSecs / 86400.)      print *, "Local calculation got:", floatTest      print *, "Difference: ", (floatTest - floatDay)      test = (abs(floatTest - floatDay) < 0.00001)      str = "MFM_DateFltDayOfYear:  return correct days.seconds"      call MFM_ErrorTest(test, str)            end program main !\end{verbatim}