Rev 3 | Rev 9 | Go to most recent revision | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed
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 intervalc 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 datestrcharacter chreal diffreal timec ---------------------------------------------------------------------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,*) refdateelse ! 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(startdate(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 Endc ---------------------------------------------------------------------end