PROGRAM PLOTTOD(OUTPUT=/1000)
C     PROGRAM CON(OUTPUT=/1000,TEKPL=/1000)
C     DIMENSION W(40,40)
      DIMENSION XT(2000),YT(2000)
      COMMON WORK(5000)
      DATA XMIN,XMAX,YMIN,YMAX/-1.0,2.,-1.0,2./
      DATA M,N/40,20/
      DATA ZSCALE/.50/
      CALL HP2623
C     CALL ZETA(3600,11,8)
C     CALL ZETA
C     CALL TK4014(960,1)
C     CALL HWSPEC('TEKPL','FILE')
      CALL SETDEV(7,7)
      OPEN(12,FILE='NPOINT',STATUS='OLD',FORM='FORMATTED')
      REWIND(12)
      READ(12,*)NPOINT
C
C   CREATE POTENTIAL SURFACE ARRAY
C     CLOSE(UNIT=7,STATUS='KEEP')
C     OPEN(7,FILE='DATA',STATUS='OLD',FORM='FORMATTED')
C     REWIND(7)
C     READ(7,*)DX,DY,AX,AY,XD,YD,WX,WY
C     DO 1 J=1,N
C     Y=(YMAX-YMIN)*(FLOAT(J-1)/FLOAT(N-1))+YMIN
C     DO 1 I=1,M
C     X=(XMAX-XMIN)*(FLOAT(I-1)/FLOAT(M-1))+YMIN
C     W(I,J)=VX+VY
C1    CONTINUE
C
C   CREATE THE CONTOURS
      CALL PAGE(11,8.5)
      CALL AREA2D(6.00,6.00)
      CALL HEIGHT(.25)
      CALL YNAME('Y',1)
      CALL XNAME('X',1)
      CALL GRAF(XMIN,1.0,XMAX,YMIN,1.0,YMAX)
      CALL RLMESS('PYIN=-0.30',10,0.00,1.50)
      CALL FRAME
      CALL BCOMON(10000)
C     CALL CONMAK(W,M,N,ZSCALE)
C     C
C   PLOT CONTOURS
C     CALL CONLIN(0,'SOLID','NOLABELS',2,5)
C     CALL CONANG(90.)
C     CALL RASPLN(0.25)
C
C   PLOT TRAJECTORIES
      OPEN(11,FILE='TRAJ',STATUS='OLD',FORM='FORMATTED')
      REWIND(11)
      DO 20 IP=1,NPOINT
      READ(11,*)XT(IP),YT(IP)
 20   CONTINUE
      CALL CURVE(XT,YT,NPOINT,0)
C     CALL CONTUR(1,'NOLABELS','DRAW')
      CALL ENDPL(0)
      CALL DONEPL
      STOP
      END