Subversion Repositories lagranto.ecmwf

date      PROGRAM datelist

c     **********************************************************************

c     * Handling of date lists                                             *

c     * Michael Sprenger                                                   *

c     **********************************************************************

      implicit none

c     ---------------------------------------------------------------------

c     Declaration of variables

c     ---------------------------------------------------------------------

c     Parameters

      character*80              datefile

      character*80              mode 

      character*80              startdate

      character*80              finaldate

      character*80              refdate

      integer                   interval

      real                      randpercent

c     Date list

      integer                   ndates

      integer,allocatable,dimension(:)    :: year,month,day,hour

c     Auxiliary variables

      integer                     i

      integer                     year1,month1,day1,hour1,min1

      integer                     year2,month2,day2,hour2,min2

      character                   direction

      integer                     date1(5),date2(5)

      character*11                datestr,datestr1,datestr2

      character                   ch

      real                        diff

      real                        time

      character*80                datefile1,datefile2

      integer                     timestamp1,timestamp2

      integer                     oldstamp1 ,oldstamp2

      integer                     state

c     ---------------------------------------------------------------------

c     Preparations

c     ---------------------------------------------------------------------

c     Read parameter file

      open(10,file='datelist.param')

       read(10,*) datefile                           ! General parameters

       read(10,*) mode

       if ( mode.eq.'-create' ) then                 ! Create date list

          read(10,*) startdate

          read(10,*) finaldate

          read(10,*) interval

       elseif ( mode.eq.'-totime' ) then             ! Convert to time list

          read(10,*) refdate

       elseif ( mode.eq.'-todate' ) then             ! Convert to date list

          read(10,*) refdate

       elseif ( mode.eq.'-onlyin1' ) then            ! Dates only in file 1

          read(10,*) datefile1

          read(10,*) datefile2

       else                                          ! Invalid mode

          print*,' ERROR: invalid mode for datelist'

          stop

       endif

      close(10)

c     ---------------------------------------------------------------------

c     Create a date list (-create)

c     ---------------------------------------------------------------------

      if ( mode.ne.'-create' ) goto 100

c     Check whether interval is ok

      if ( ( interval.le.0 ).or.(interval.gt.24) ) then

         print*,'Interval must be between 1 h and 24 h... Stop'

         stop

      endif

c     Extract dates and times

      read(startdate( 1: 4),*) year1

      read(startdate( 5: 6),*) month1

      read(startdate( 7: 8),*) day1

      read(startdate(10:11),*) hour1

      read(startdate(12:13),*) min1

      read(finaldate( 1: 4),*) year2

      read(finaldate( 5: 6),*) month2

      read(finaldate( 7: 8),*) day2

      read(finaldate(10:11),*) hour2

      read(finaldate(12:13),*) min2

c     Get direction of the date file

      if (year2.gt.year1) then

         direction = 'f'

         goto 101

      elseif (year2.lt.year1) then

         direction = 'b'

         goto 101

      endif

      if (month2.gt.month1) then

         direction = 'f'

         goto 101

      elseif (month2.lt.month1) then

         direction = 'b'

         goto 101

      endif

      if (day2.gt.day1) then

         direction = 'f'

         goto 101

      elseif (day2.lt.day1) then

         direction = 'b'

         goto 101

      endif

      if (hour2.gt.hour1) then

         direction = 'f'

         goto 101

      elseif (hour2.lt.hour1) then

         direction = 'b'

         goto 101

      endif

      if (min2.gt.min1) then

         direction = 'f'

         goto 101

      elseif (min2.lt.min1) then

         direction = 'b'

         goto 101

      endif

      direction = 'f'

 101  continue

c     Set the interval step depending on the direction

      if ( direction.eq.'b' ) then

         interval = -interval

      endif

c     Save the dates in arrays

      date1(1) = year1

      date1(2) = month1

      date1(3) = day1

      date1(4) = hour1 

      date1(5) = 0

      date2(1) = year2

      date2(2) = month2

      date2(3) = day2

      date2(4) = hour2 

      date2(5) = 0

c     Get starting and ending date for the date list

      if ( direction.eq.'f' ) then

         do while ( mod(date1(4),interval) .ne. 0 )

            date1(4) = date1(4) - 1

         enddo

         if (min2.ne.0) call newdate(date2,1.,date2)

         do while ( mod(date2(4),interval) .ne. 0 )

            date2(4) = date2(4) + 1

         enddo

      else

         if (min1.ne.0) call newdate(date1,1.,date1)

         do while ( mod(date1(4),interval) .ne. 0 )

            date1(4) = date1(4) + 1

         enddo

         do while ( mod(date2(4),interval) .ne. 0 )

            date2(4) = date2(4) - 1

         enddo

      endif

c     Create and write the datefile

      if ( datefile.ne.'/dev/stdout') then

         open(10,file=datefile)

      endif

 102  continue

       call datestring(datestr,date1(1),date1(2),date1(3),date1(4) )

       if ( datefile.ne.'/dev/stdout') then

          write(10,*) datestr

       else

          write(*,*) datestr

       endif

       if ( ( date1(1).ne.date2(1) ).or.

     >      ( date1(2).ne.date2(2) ).or.

     >      ( date1(3).ne.date2(3) ).or.

     >      ( date1(4).ne.date2(4) ) ) 

     > then

          diff = real(interval)

          call newdate(date1,diff,date1)

          goto 102

       endif

      if ( datefile.ne.'/dev/stdout') then

         close(10)

      endif

 100  continue

c     ---------------------------------------------------------------------

c     Convert dates to a list of times

c     ---------------------------------------------------------------------

      if ( mode.ne.'-totime' ) goto 112

c     Extract reference date

      read(refdate( 1: 4),*) year1

      read(refdate( 5: 6),*) month1

      read(refdate( 7: 8),*) day1

      read(refdate(10:11),*) hour1

      read(refdate(12:13),*) min1

c     Loop through the date file

      open(10,file=datefile)

 111  read(10,*,end=110) datestr

c     Extract date

      read(datestr( 1: 4),*) year2

      read(datestr( 5: 6),*) month2

      read(datestr( 7: 8),*) day2

      read(datestr(10:11),*) hour2

      min1 = 0

c     Get the time difference

      date1(1) = year1

      date1(2) = month1

      date1(3) = day1

      date1(4) = hour1 

      date1(5) = 0

      date2(1) = year2

      date2(2) = month2

      date2(3) = day2

      date2(4) = hour2 

      date2(5) = 0

      call timediff(date2,date1,diff)

c     Write it to screen

      write(*,'(i6)') nint(diff)

      goto 111

 110  continue

c     Close datefile

      close(10)

c     Exit point      

  112 continue

c     ---------------------------------------------------------------------

c     Convert times to a list of dates 

c     ---------------------------------------------------------------------

      if ( mode.ne.'-todate' ) goto 122

c     Extract reference date

      read(refdate( 1: 4),*) year1

      read(refdate( 5: 6),*) month1

      read(refdate( 7: 8),*) day1

      read(refdate(10:11),*) hour1

      read(refdate(12:13),*) min1

c     Loop through the date file

      open(10,file=datefile)

 121  read(10,*,end=120) time

c     Calculate the new date

      date1(1) = year1

      date1(2) = month1

      date1(3) = day1

      date1(4) = hour1 

      date1(5) = 0

      call newdate(date1,time,date2)

      call datestring(datestr,date2(1),date2(2),date2(3),date2(4) )

c     Write it to screen

      write(*,'(a11)') trim(datestr)

      goto 121

 120  continue

c     Close datefile

      close(10)

c     Exit point      

 122  continue

c     ---------------------------------------------------------------------

c     Extract all dates which are only in one datefile 

c     ---------------------------------------------------------------------

      if ( mode.ne.'-onlyin1' ) goto 134

c     Set reference date   

      date1(1) = 1979

      date1(2) = 1

      date1(3) = 1

      date1(4) = 0

      date1(5) = 0

c     Open the output file

      if ( datefile.ne.'/dev/stdout') then

         open(30,file=datefile)

      endif

c     Loop through the input date files

      open(10,file=datefile1)

      open(20,file=datefile2)

c     Loop through both date files

      state     = 0

      oldstamp1 = 0

      oldstamp2 = 0

 131  if ( (state.eq.1).or.(state.eq.0) ) then 

         read(10,*,end=133) datestr1

      endif

      if ( (state.eq.2).or.(state.eq.0) ) then 

         read(20,*,end=130) datestr2

      endif

c     Get time stamp for both date strings

      if ( (state.eq.1).or.(state.eq.0) ) then 

         read(datestr1( 1: 4),*) year2

         read(datestr1( 5: 6),*) month2

         read(datestr1( 7: 8),*) day2

         read(datestr1(10:11),*) hour2

         date2(1) = year2

         date2(2) = month2

         date2(3) = day2

         date2(4) = hour2 

         date2(5) = 0

         call timediff(date2,date1,diff)

         timestamp1 =  nint( diff )

         if ( timestamp1.lt.oldstamp1 ) then

            print*,' ERROR: datelist must be ordered ',trim(datefile1)

            stop

         else

            oldstamp1 = timestamp1

         endif

      endif

      if ( (state.eq.1).or.(state.eq.0) ) then 

         read(datestr2( 1: 4),*) year2

         read(datestr2( 5: 6),*) month2

         read(datestr2( 7: 8),*) day2

         read(datestr2(10:11),*) hour2

         date2(1) = year2

         date2(2) = month2

         date2(3) = day2

         date2(4) = hour2 

         date2(5) = 0

         call timediff(date2,date1,diff)

         timestamp2 =  nint( diff )

         if ( timestamp2.lt.oldstamp2 ) then

            print*,' ERROR: datelist must be ordered ',trim(datefile2)

            stop

         else

            oldstamp2 = timestamp2

         endif

      endif

c     Write output and set new state

      if ( timestamp1.gt.timestamp2 ) then

        state = 2

c        print*,trim(datestr1)//'.'//trim(datestr2)//' >'

      else if ( timestamp1.lt.timestamp2 ) then

        state = 1

        if ( datefile.ne.'/dev/stdout') then

           write(30,*) datestr1

        else

           write(*,*) datestr1

        endif

c        print*,trim(datestr1)//'.'//trim(datestr2)//' -> out'

      else if (timestamp1.eq.timestamp2 ) then

        state = 0

c        print*,trim(datestr1)//'.'//trim(datestr2)//' ='

      endif

      goto 131

c     Exit point for parallel reading through files

 130  continue

c     Write remaining part of datefile 1

 132  continue

        read(10,*,end=133) datestr1

        if ( datefile.ne.'/dev/stdout') then

           write(30,*) datestr1

        else

           write(*,*) datestr1

        endif

      goto 132

 133  continue

c     Close datefile

      close(10)

      close(20)

      close(30)

c     Exit point      

 134  continue

c     ---------------------------------------------------------------------

c     End

c     ---------------------------------------------------------------------

      end