Blame | Last modification | View Log | Download | RSS feed
date PROGRAM datelistc **********************************************************************c * Handling of date lists *c * Michael Sprenger *c **********************************************************************implicit nonec ---------------------------------------------------------------------c Declaration of variablesc ---------------------------------------------------------------------c Parameterscharacter*80 datefilecharacter*80 modecharacter*80 startdatecharacter*80 finaldatecharacter*80 refdateinteger intervalreal randpercentc Date listinteger ndatesinteger,allocatable,dimension(:) :: year,month,day,hourc Auxiliary variablesinteger iinteger year1,month1,day1,hour1,min1integer year2,month2,day2,hour2,min2character directioninteger date1(5),date2(5)character*11 datestr,datestr1,datestr2character chreal diffreal timecharacter*80 datefile1,datefile2integer timestamp1,timestamp2integer oldstamp1 ,oldstamp2integer statec ---------------------------------------------------------------------c Preparationsc ---------------------------------------------------------------------c Read parameter fileopen(10,file='datelist.param')read(10,*) datefile ! General parametersread(10,*) modeif ( mode.eq.'-create' ) then ! Create date listread(10,*) startdateread(10,*) finaldateread(10,*) intervalelseif ( mode.eq.'-totime' ) then ! Convert to time listread(10,*) refdateelseif ( mode.eq.'-todate' ) then ! Convert to date listread(10,*) refdateelseif ( mode.eq.'-onlyin1' ) then ! Dates only in file 1read(10,*) datefile1read(10,*) datefile2else ! Invalid modeprint*,' ERROR: invalid mode for datelist'stopendifclose(10)c ---------------------------------------------------------------------c Create a date list (-create)c ---------------------------------------------------------------------if ( mode.ne.'-create' ) goto 100c Check whether interval is okif ( ( interval.le.0 ).or.(interval.gt.24) ) thenprint*,'Interval must be between 1 h and 24 h... Stop'stopendifc Extract dates and timesread(startdate( 1: 4),*) year1read(startdate( 5: 6),*) month1read(startdate( 7: 8),*) day1read(startdate(10:11),*) hour1read(startdate(12:13),*) min1read(finaldate( 1: 4),*) year2read(finaldate( 5: 6),*) month2read(finaldate( 7: 8),*) day2read(finaldate(10:11),*) hour2read(finaldate(12:13),*) min2c Get direction of the date fileif (year2.gt.year1) thendirection = 'f'goto 101elseif (year2.lt.year1) thendirection = 'b'goto 101endifif (month2.gt.month1) thendirection = 'f'goto 101elseif (month2.lt.month1) thendirection = 'b'goto 101endifif (day2.gt.day1) thendirection = 'f'goto 101elseif (day2.lt.day1) thendirection = 'b'goto 101endifif (hour2.gt.hour1) thendirection = 'f'goto 101elseif (hour2.lt.hour1) thendirection = 'b'goto 101endifif (min2.gt.min1) thendirection = 'f'goto 101elseif (min2.lt.min1) thendirection = 'b'goto 101endifdirection = 'f'101 continuec Set the interval step depending on the directionif ( direction.eq.'b' ) theninterval = -intervalendifc Save the dates in arraysdate1(1) = year1date1(2) = month1date1(3) = day1date1(4) = hour1date1(5) = 0date2(1) = year2date2(2) = month2date2(3) = day2date2(4) = hour2date2(5) = 0c Get starting and ending date for the date listif ( direction.eq.'f' ) thendo while ( mod(date1(4),interval) .ne. 0 )date1(4) = date1(4) - 1enddoif (min2.ne.0) call newdate(date2,1.,date2)do while ( mod(date2(4),interval) .ne. 0 )date2(4) = date2(4) + 1enddoelseif (min1.ne.0) call newdate(date1,1.,date1)do while ( mod(date1(4),interval) .ne. 0 )date1(4) = date1(4) + 1enddodo while ( mod(date2(4),interval) .ne. 0 )date2(4) = date2(4) - 1enddoendifc Create and write the datefileif ( datefile.ne.'/dev/stdout') thenopen(10,file=datefile)endif102 continuecall datestring(datestr,date1(1),date1(2),date1(3),date1(4) )if ( datefile.ne.'/dev/stdout') thenwrite(10,*) datestrelsewrite(*,*) datestrendifif ( ( date1(1).ne.date2(1) ).or.> ( date1(2).ne.date2(2) ).or.> ( date1(3).ne.date2(3) ).or.> ( date1(4).ne.date2(4) ) )> thendiff = real(interval)call newdate(date1,diff,date1)goto 102endifif ( datefile.ne.'/dev/stdout') thenclose(10)endif100 continuec ---------------------------------------------------------------------c Convert dates to a list of timesc ---------------------------------------------------------------------if ( mode.ne.'-totime' ) goto 110c Extract reference dateread(refdate( 1: 4),*) year1read(refdate( 5: 6),*) month1read(refdate( 7: 8),*) day1read(refdate(10:11),*) hour1read(refdate(12:13),*) min1c Loop through the date fileopen(10,file=datefile)111 read(10,*,end=110) datestrc Extract dateread(datestr( 1: 4),*) year2read(datestr( 5: 6),*) month2read(datestr( 7: 8),*) day2read(datestr(10:11),*) hour2min1 = 0c Get the time differencedate1(1) = year1date1(2) = month1date1(3) = day1date1(4) = hour1date1(5) = 0date2(1) = year2date2(2) = month2date2(3) = day2date2(4) = hour2date2(5) = 0call timediff(date2,date1,diff)c Write it to screenwrite(*,'(i6)') nint(diff)goto 111110 continuec Close datefileclose(10)c ---------------------------------------------------------------------c Convert times to a list of datesc ---------------------------------------------------------------------if ( mode.ne.'-todate' ) goto 120c Extract reference dateread(refdate( 1: 4),*) year1read(refdate( 5: 6),*) month1read(refdate( 7: 8),*) day1read(refdate(10:11),*) hour1read(refdate(12:13),*) min1c Loop through the date fileopen(10,file=datefile)121 read(10,*,end=120) timec Calculate the new datedate1(1) = year1date1(2) = month1date1(3) = day1date1(4) = hour1date1(5) = 0call newdate(date1,time,date2)call datestring(datestr,date2(1),date2(2),date2(3),date2(4) )c Write it to screenwrite(*,'(a11)') trim(datestr)goto 121120 continuec Close datefileclose(10)c ---------------------------------------------------------------------c Extract all dates which are only in one datefilec ---------------------------------------------------------------------if ( mode.ne.'-onlyin1' ) goto 130c Set reference datedate1(1) = 1979date1(2) = 1date1(3) = 1date1(4) = 0date1(5) = 0c Open the output fileif ( datefile.ne.'/dev/stdout') thenopen(30,file=datefile)endifc Loop through the input date filesopen(10,file=datefile1)open(20,file=datefile2)c Loop through both date filesstate = 0oldstamp1 = 0oldstamp2 = 0131 if ( (state.eq.1).or.(state.eq.0) ) thenread(10,*,end=133) datestr1endifif ( (state.eq.2).or.(state.eq.0) ) thenread(20,*,end=130) datestr2endifc Get time stamp for both date stringsif ( (state.eq.1).or.(state.eq.0) ) thenread(datestr1( 1: 4),*) year2read(datestr1( 5: 6),*) month2read(datestr1( 7: 8),*) day2read(datestr1(10:11),*) hour2date2(1) = year2date2(2) = month2date2(3) = day2date2(4) = hour2date2(5) = 0call timediff(date2,date1,diff)timestamp1 = nint( diff )if ( timestamp1.lt.oldstamp1 ) thenprint*,' ERROR: datelist must be ordered ',trim(datefile1)stopelseoldstamp1 = timestamp1endifendifif ( (state.eq.1).or.(state.eq.0) ) thenread(datestr2( 1: 4),*) year2read(datestr2( 5: 6),*) month2read(datestr2( 7: 8),*) day2read(datestr2(10:11),*) hour2date2(1) = year2date2(2) = month2date2(3) = day2date2(4) = hour2date2(5) = 0call timediff(date2,date1,diff)timestamp2 = nint( diff )if ( timestamp2.lt.oldstamp2 ) thenprint*,' ERROR: datelist must be ordered ',trim(datefile2)stopelseoldstamp2 = timestamp2endifendifc Write output and set new stateif ( timestamp1.gt.timestamp2 ) thenstate = 2c print*,trim(datestr1)//'.'//trim(datestr2)//' >'else if ( timestamp1.lt.timestamp2 ) thenstate = 1if ( datefile.ne.'/dev/stdout') thenwrite(30,*) datestr1elsewrite(*,*) datestr1endifc print*,trim(datestr1)//'.'//trim(datestr2)//' -> out'else if (timestamp1.eq.timestamp2 ) thenstate = 0c print*,trim(datestr1)//'.'//trim(datestr2)//' ='endifgoto 131c Exit point for parallel reading through files130 continuec Write remaining part of datefile 1132 continueread(10,*,end=133) datestr1if ( datefile.ne.'/dev/stdout') thenwrite(30,*) datestr1elsewrite(*,*) datestr1endifgoto 132133 continuec Close datefileclose(10)close(20)close(30)c ---------------------------------------------------------------------c Endc ---------------------------------------------------------------------end