A9.13  Appendix: Computer Programs

*			  PROGRAM I.20			    *

*******************************************************************************
*                                   FDLNCOOR43                                *
* THIS PROGRAM PLOTS THE RESULTS OF *.E*IN AND *.E*OUT FOUR TO A PAGE.        *
*******************************************************************************
 

      PROGRAM FDLNCOOR43

      PARAMETER (NB=50,NN=1500)
      REAL XYZ(6),X(NB,NN),Y(NB,NN)
      INTEGER NC,NPT(NB),NP
      CHARACTER*72 FNAME,TITLE,XLAB,YLAB
     

      WRITE(6,*) 'ENTER NP:,NOS.OR PLOTS ON A PAGE'
      READ (5,*) NP

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

      DO NC=1,NP

       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
	WRITE(6,*) NLN
 3     READ(1,19,ERR=2) (XYZ(M),M=1,2)
D	WRITE (6,*) XYZ(1),XYZ(2)
       NPOINT = NPOINT+1
       X(NLN,NPOINT) = (XYZ(Ho))
       Y(NLN,NPOINT) = (XYZ(Vo))
D	WRITE(6,*) X(NLN,NPOINT),Y(NLN,NPOINT)

       NPT(NLN) = NPOINT
       GO TO 3
 70    CLOSE(1)
       PRINT *,'NLN: ',NLN
D	DO I = 1,NLN
D		WRITE(6,*) NPT(I)
D	END DO       
C       WRITE(6,*) 'ENTER THE X-LABEL'
C       WRITE(6,'(A,$)') '--> '
C       READ (5,'(A)') XLAB
       XLAB = 'THETA'    
C       WRITE(6,*) 'ENTER THE Y-LABEL'
C       WRITE(6,'(A,$)') '--> '
C       READ (5,'(A)') YLAB
       YLAB = 'PHI'   
C       WRITE(6,*) 'ENTER LNPT,LNK'
C       READ(5,*) LNPT,LNK
       LNPT = 1
       LNK = 4 
	WRITE(6,*) NLN
	WRITE(6,*) NLN
	WRITE(6,*) NLN
	WRITE(6,*) NLN
D	DO J = 1,NLN
D		WRITE(6,*) X(NLN,J),Y(NLN,J)
D	END DO
 

       CALL PLOTR4(NC,X,Y,NPT,NLN,TITLE,XLAB,YLAB,LNPT,LNK)
      END DO   
C-----------------------------------------------------------

      CALL FLUSH

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