PROGRAM PLOTTOD(OUTPUT=/1000)
      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
      CALL SETDEV(7,7)
      OPEN(12,FILE='NPOINT',STATUS='OLD',FORM='FORMATTED')
      REWIND(12)
      READ(12,*)NPOINT
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   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)
      CALL ENDPL(0)
      CALL DONEPL
      STOP
      END