REM read the Orbit.dat file and create a plttek.data file
REM to be processed by Tek2plot file for display by the
REM program WPLOT
REM GE Copeland
REM Dept of Physics
REM Old Dominion University
REM Jan 2001

OPEN "Orbit.dat" FOR INPUT AS #2
OPEN "plttek.dat" FOR OUTPUT AS #3

REM read the header of Orbit.dat
  FOR i = 1 TO 9
LINE INPUT #2, dum$(i)
PRINT dum$(i)
         NEXT i
         
         b$ = " "
Save1$ = RTRIM$(LTRIM$(dum$(3))) + b$ + RTRIM$(LTRIM$(dum$(4))) + b$ + RTRIM$(LTRIM$(dum$(5)) + b$ + LTRIM$(dum$(6)))
    PRINT Save1$

    DIM x(400), Y(400)
c = 0
DO UNTIL EOF(2)
      INPUT #2, dum, x(c), Y(c)
      c = c + 1
   REM   PRINT c
LOOP
REM got data
PRINT "Count= "; c
REM now write out to plttek.dat
CLOSE #2
WRITE #3, 1
WRITE #3, 1
WRITE #3, c
WRITE #3, "Y  in AU"
WRITE #3, "X in AU"
WRITE #3, "Plot of data from Orbit.exe - Phys 313"
WRITE #3, Save1$
FOR i = 1 TO c
WRITE #3, x(i), Y(i)
NEXT i
CLOSE #3
CHAIN "Tek2wplt.bas"


         END