A9.13  Appendix: Computer Programs

*		         PROGRAM I.18				   *

******************************************************************************
*                              FDLNCOORPRO                                   *
*     THIS ROUTINE PLOTS BOTH THE GEOMETRY OF THE SYSTEM AS                  *
*     WELL AS THE TRAJECTORIES OF SELECTED STARTING POINTS                   *
******************************************************************************
 

      PROGRAM FDLNCOORPRO

      PARAMETER (NB=100,NN=1000)
      REAL XYZ(6),X(NB,NN),Y(NB,NN)
      INTEGER NPT(NB),J
      CHARACTER*72 FNAME,TITLE,XLAB,YLAB
      

      WRITE(6,*) 'ENTER THE NUMBER OF TRAJECTORIES TO BE PLOTTED:'
      READ(5,*) J
      DO I=1,(J+1)
       NC = I
       WRITE(6,5) NC
       READ(5,10) FNAME
C       WRITE(6,*) 'ENTER INTEGERS FOR HORIZONTOL AND VERTICAL AXES:' 
C       WRITE(6,12)
C       READ (5,*) Ho,Vo
       Ho=1
       Vo=2
       OPEN (UNIT=1,STATUS='OLD',FILE=FNAME)
       READ (1,*)
 2     NPOINT = 0
       READ(1,*,END=70) NLN
       READ(1,*)
 3     READ(1,19,ERR=2) (XYZ(M),M=1,3)
       NPOINT = NPOINT+1
       X(NLN,NPOINT) = XYZ(Ho)
       Y(NLN,NPOINT) = XYZ(Vo)
D       WRITE(6,*) XYZ(Ho),XYZ(Vo)
       NPT(NLN) = NPOINT
       GO TO 3
 70    CLOSE(1)
      END DO
      

      PRINT *,'NLN: ',NLN

      WRITE(6,*) 'ENTER THE TITLE'
      WRITE(6,'(A,$)') ' --> '
      READ(5,'(A)') TITLE

      WRITE(6,*) 'ENTER THE X-LABEL'
      WRITE(6,'(A,$)') '--> '
      READ (5,'(A)') XLAB
   

      WRITE(6,*) 'ENTER THE Y-LABEL'
      WRITE(6,'(A,$)') '--> '
      READ (5,'(A)') YLAB
    

      WRITE(6,*) 'ENTER LNPT,LNK'
      READ(5,*) LNPT,LNK

      CALL PLOTR5(X,Y,NPT,NLN,TITLE,XLAB,YLAB,LNPT,LNK)

      CALL FRAME

      CALL FLUSH

 5    FORMAT(X,'ENTER THE NAME OF THE FILE: ',I3)
 10   FORMAT(A72)
 12   FORMAT(X,'X---1',/, ' Y---2',/,' Z---3')
 19   FORMAT(3(X,f9.6))
c 19   FORMAT(3(X,F13.6))
      STOP 
      END