program time

  implicit none

  integer :: dsm(12)   ! elapsed Days on Start of Month
  integer :: dpm(12)   ! Days Per Month
  integer :: n
  real    :: t0,dt

  data     dsm  / 0,31,59, 90,120,151, 181,212,243, 273,304,334/
  data     dpm  /31,28,31, 30, 31, 30,  31, 31, 30,  31, 30, 31/

!-------------------------------------------------------------------------------

   t0 = 15.5
   dt = 30.5
   do n=1,12
   !  write(*,*) t0 + (n-1)*dt               ! bill large algorithm
   !  write(*,*) float(dsm(n)) + 15          ! 15th of month
      write(*,*) float(dsm(n)) + dpm(n)/2.0  ! mid month
   end do

stop
end