WebSVN – lagranto.ecmwf – Diff – /trunk/goodies/datelist.f

-      PROGRAM datelist
+date      PROGRAM datelist
 c     **********************************************************************
 c     * Handling of date lists                                             *
 c     * Michael Sprenger                                                   *
 c     **********************************************************************
       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
       integer                     i
       integer                     year1,month1,day1,hour1,min1
       integer                     year2,month2,day2,hour2,min2
       character                   direction
       integer                     date1(5),date2(5)
-      character*11                datestr
+      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     ---------------------------------------------------------------------
        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
       read(finaldate( 1: 4),*) year2
       read(finaldate( 5: 6),*) month2
       read(finaldate( 7: 8),*) day2
       read(finaldate(10:11),*) hour2
-      read(startdate(12:13),*) min2
+      read(finaldate(12:13),*) min2
 c     Get direction of the date file
       if (year2.gt.year1) then
          direction = 'f'
          goto 101
 c     ---------------------------------------------------------------------
 c     Convert dates to a list of times
 c     ---------------------------------------------------------------------
-      if ( mode.ne.'-totime' ) goto 110
+      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
 continue
 c     Close datefile
       close(10)
+c     Exit point
+continue
 c     ---------------------------------------------------------------------
 c     Convert times to a list of dates
 c     ---------------------------------------------------------------------
-      if ( mode.ne.'-todate' ) goto 120
+      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
 continue
 c     Close datefile
       close(10)
+c     Exit point
+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
+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
+continue
+c     Write remaining part of datefile 1
+continue
+        read(10,*,end=133) datestr1
+        if ( datefile.ne.'/dev/stdout') then
+           write(30,*) datestr1
+        else
+           write(*,*) datestr1
+        endif
+      goto 132
+continue
+c     Close datefile
+      close(10)
+      close(20)
+      close(30)
+c     Exit point
+continue
 c     ---------------------------------------------------------------------
 c     End
 c     ---------------------------------------------------------------------
       end

Subversion Repositories lagranto.ecmwf

(root)/trunk/goodies/datelist.f – Rev 5 → 9