USGS InfoBank program -- window

window

  Program to window a data set based on time and/or positions.
    

  Enter name of input  file:

       (Will accept any file with time and lat/long)
  Enter names of output files:
  INPUT/OUTPUT format must read
     (ida,ihr,imin,isec,iten,rlat,rlon)
  Enter input format:
  (for as many windows as you have...)
  Enter N,S,W,E degree limits for window 
     (<0,0,0,0> = no area window)
  Enter start-end time limits for window 
     (<0> = no time window)
  dddhhmmsst dddhhmmsst
    

  openem      (external shareable)

  percentage  (external shareable)
  showlines   (external shareable)
  whatfmt     (external shareable)
  converttime (external shareable)
  yesno       (external shareable)
  forerror    (external shareable)
    

  Wm. Clinton Steele, Marine Geology, 1/29/81

   
  Carolyn Degnan, 04/02/85  Made structured, added standard goodies,
                              changed to new time/format, and ran MAT.
  Carolyn Degnan, 04/25/86  Changed to unformatted reads for window limits
  Clint Steele,   06/01/87  Added ability to have window crossing 180 
                              longitude.
  Carolyn Degnan, 07/27/88  If no records are written to an output file,
                            that output file is deleted when closed.
  Clint Steele,   10/27/88  Added option to include first point outside
                              of window.
                            Changed i/o to internal reads.
  Clint Steele,    4/24/92  Added call to DEGREES
  Carolyn Degnan, 10/30/98  Added call to YESNO

    

     read (5, 29)
*       iyr , ida , ihr , imin , isec , iten ,
*       jyr , jda , jhr , jmin , jsec , jten 
     read (90, '(q,a)', iostat = ierr, err = 52, end = 99)
*        ALength, ARecord
     read (ARecord, afmt, iostat = ierr, err = 52)
*        kyr, kda ,khr, kmin, ksec, kten, rlat, rlon

    

     write (6, 54) iread
             write (91, '(a)') LastRecord(1:LastLength)
         write (91, '(a)') ARecord(1:ALength)
             write (91, '(a)') ARecord(1:ALength)
         write (92, '(a)') ARecord(1:ALength)
 write (6, 104) iwrite ,
*       North , South , West , East ,
*       iyr , ida , ihr , imin ,
*       isec , iten ,
*       jyr , jda , jhr , jmin ,
*       jsec , jten 
    

 call openem (1, -2)
     call degrees (North , South , West , East )
     call converttime ('COMPRESS', ida , ihr , imin ,
*                 isec , iten , itime1 , ErrorFlag)
     call converttime ('COMPRESS', jda , jhr , jmin ,
*                 jsec , jten , itime2 , ErrorFlag)
         call converttime ('COMPRESS', jda , jhr , 
*                  jmin , jsec , jten , itime2 , 
*                  ErrorFlag)
 call yesno ('Y[N]', AExtend)
     call converttime ('COMPRESS', kda, khr, kmin, ksec, kten, 
*                      ktime, ErrorFlag)
     call forerror (ierr)