*		         PROGRAM I.17				    *

*******************************************************************************
*                                FDLNGEOM                                     *
* THIS IS A PROGRAM THAT IS GOOD FOR PLOTTING THE GEOMETRY OF THE SYSTEM.     *
*******************************************************************************
 

      PROGRAM FDLNGEOM

      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
       WRITE(6,*) 'ENTER INTEGERS FOR HORIZONTOL AND VERTICAL AXES:' 
       WRITE(6,12)
       READ (5,*) Ho,Vo
       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 PLOTR6(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

Return to the Table of Contents for Buckley's Thesis

Return to HISCALE List of Appendices

Return to Ulysses HISCALE Data Analysis Handbook Table of Contents