A9.13  Appendix: Computer Programs

*			 PROGRAM I.19			  *

******************************************************************************
*                                  FDLNFOUR                                  *
* THIS PROGRAM CAN BE USED TO PLOT BOTH THE GEOMETRY AND TRAJECTORIES        *
* FOUR TO A PAGE.                                                            *
******************************************************************************
 

      PROGRAM FDLNFOUR

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

      WRITE(6,*) 'ENTER NP: THE NO. OF PLOTS ON PAGE'
      READ (5,*) NP

      DO NC=1,NP
      

      WRITE(6,*) 'ENTER THE NUMBER OF TRAJECTORIES TO BE PLOTTED:'
      READ(5,*) J
      DO I=1,(J+1)
       NL = I
       WRITE(6,5) NL
       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 PLOTR4(NC,X,Y,NPT,NLN,TITLE,XLAB,YLAB,LNPT,LNK)
      END DO

      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