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