|
Fundamental TechnologiesUlysses HISCALE Pages |
Appendix 9 Geometric Factor Study for the Deflected and Unscattered Electrons of HISCALE (Buckley MS Thesis continued)
* 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