program plt
c
c Plot min,max from given file (min,max from each of 241 histories
c during the 10day period 79079-89)
c
      parameter(mxpts=250)
      dimension xx(mxpts),rminmax(mxpts,2),iasf(13),del(mxpts)
      data lu/7/
      character*10 flnm
      character*40 chrec
c     data flnm/'tnzm4.mnmx'/
      data iasf/13*1/
c
      write(6,"('Enter file name: ',$)")
      read(5,"(a)") flnm
      open(lu,file=flnm,status='OLD')
      call opngks
      call gsclip(0)
      call gsasf(iasf)
      call gsfais(1)
c
      read(lu,"(a)") chrec
      write(6,"(a)") chrec
      ir = 0
      avedel = 0.
      do i=1,mxpts
        read(lu,"(8x,2e12.4,12x,2e12.4)",end=900) 
     +    rminmax(i,1),rminmax(i,2),totmin,totmax
        del(i) = rminmax(i,2)-rminmax(i,1)
        avedel = avedel+del(i)
        ir = ir+1
      enddo
      write(6,"('>>> reached ',i4,' records')") mxpts
900   continue
      write(6,"('EOF at i=',i5,' ir=',i5)") i,ir
      avedel = avedel / float(ir)
      write(chrec(15:35),"('AVE DELTA = ',1pe9.3)") avedel
      iday = 79
      uthr = 0.
      do i=1,ir
        uthr = uthr+1. 
        if (uthr.gt.23.) then
          uthr = 0. 
          iday = iday+1
        endif
        xx(i) = float(iday) + uthr/24.
      enddo
c
c Plot min,max:
c
      call agsetr('FRAME.',2.)
      call ezmxy(xx,rminmax,mxpts,2,ir,chrec)
      call getset(vl,vr,vb,vt,wl,wr,wb,wt,lty)
      do i=1,10
        yr = wl+float(i)
        call line(yr,wb,yr,wt)
      enddo
      call frame
c
c Plot deltas:
c
      call ezxy(xx,del,ir,chrec)
      call getset(vl,vr,vb,vt,wl,wr,wb,wt,lty)
      do i=1,10
        yr = wl+float(i)
        call line(yr,wb,yr,wt)
      enddo
      call frame
      call clsgks
      stop
      end