Fundamental Technologies

Ulysses HISCALE Pages

Appendix 9 Geometric Factor Study for the Deflected and Unscattered Electrons of HISCALE (Buckley MS Thesis continued)

 

*				PROGRAM I.9				      *
*******************************************************************************
*                                DMAINTESTIN1                                 *
* THE MAIN REASON FOR THIS PROGRAM IS ARRAY INITIALIZATION.  THIS PROGRAM     *
* WILL HAVE TO BE MODIFIED TO FIT THE FORMAT OF YOUR SYSTEMS MAGNETIC         *
* FIELD DATA.                                                                 *
*******************************************************************************
        subroutine WWPLOT(DIFF,NORMB)
      PARAMETER (NB=50,NN=1000)
      REAL X(NB,NN),Y(NB,NN),XI,P(NB,NN),Q(NB,NN),XYM(4)
      REAL EXPT(-5:5,-13:4),DIFF(-5:5,-13:4),NORMB(-5:5,-13:4)
      INTEGER NPT(NB),NLN,I,NX,NY
      CHARACTER*72 TITLE,FIG,XLAB,YLAB
        WRITE(6,*) 'ENTER THE NO. OF LINES TO BE PLOTTED FOR Y'
        READ (5,*) NLN
        do I=1,4
         NPT(I) = 16
        end do
        do I=5,7
         NPT(I) = 18
        end do
        do I=8,11
         NPT(I) = 16
        end do
        YI = 4.0         
        DO J = 1,4
          DO I = 1,16
           X(J,I) = YI/10.0            !POINTS TO BE PLOTTED ON THE X-AXIS
d          print *,'i',i,'j',j,'x',x(j,i)
d          print *,' '
           YI = YI - 1.0
          END DO
        YI = 4.0
        END DO
        YI = 4.0
        DO J=5,7        
          DO I = 1,18
           X(J,I) = YI/10.0            !POINTS TO BE PLOTTED ON THE X-AXIS
d          print *,'i',i,'j',j,'x',X(j,i)
d          print *,' '
           YI = YI - 1.0
          END DO
	YI = 4.0
	END DO
        YI = 4.0         
        DO J = 8,11
          DO I = 1,16
           X(J,I) = YI/10.0            !POINTS TO BE PLOTTED ON THE X-AXIS
d          print *,'i',i,'j',j,'x',X(j,i)
d          print *,' '
           YI = YI - 1.0
          END DO
        YI = 4.0
        END DO
        NY = 4
        NX = 5
        do J=1,4
         do I = 1,16  
          Y(J,I) = DIFF(NX,NY)
          Q(J,I) = NORMB(NX,NY)
d         print *,'i',i,'j',j,'y',y(j,i)
d         print *,' '
          NY = NY - 1
         END DO
        NY = 4
        NX = NX - 1
        END DO
        NY = 4
        NX = 1
        do J=5,7
         do I = 1,18
          Y(J,I) = DIFF(NX,NY)
          Q(J,I) = NORMB(NX,NY)
d         print *,'i',i,'j',j,'y',y(j,i)
d         print *,' '
          NY = NY - 1
         END DO
        NY = 4
        NX = NX - 1
        END DO
        NY = 4
        NX = -2
        do J=8,11
         do I = 1,16
          Y(J,I) = DIFF(NX,NY)
          Q(J,I) = NORMB(NX,NY)
d         print *,'i',i,'j',j,'y',y(j,i)
d         print *,' '
          NY = NY - 1
         END DO
        NY = 4
        NX = NX - 1
        END DO
       TITLE = 'RESULTS OF THE OBSERVED AND CALCULATED VALUES, Z = 0.0"'
       FIG = 'FIGURE C.3a'
       XLAB = 'Y VALUES (INCHES) '
       YLAB = 'OBSERVED FIELD VALUES'   
       WRITE(6,*) 'ENTER LNPT,LNK'
       READ(5,*) LNPT,LNK
       NC = 1
       CALL PLOTR4(NC,X,Y,NPT,NLN,TITLE,FIG,XLAB,YLAB,LNPT,LNK)
       XLAB = 'Y VALUES (INCHES)'
       YLAB = 'CALCULATED FIELD VALUES'
       NC = 2
       CALL PLOTR4(NC,X,Q,NPT,NLN,TITLE,FIG,XLAB,YLAB,LNPT,LNK)
C----------------------------------------------------------------------
        WRITE(6,*) 'ENTER THE NO. OF LINES TO BE PLOTTED FOR X'
        READ (5,*) NLN
        do I=1,16
         NPT(I) = 11
        END DO
     
        do I=17,18
         NPT(I) = 3
	END DO
        XI = 5.0         
        DO J = 1,16
          DO I = 1,11
           X(J,I) = XI/10.0            !POINTS TO BE PLOTTED ON THE X-AXIS
           P(J,I) = XI/10.0
d          print *,'i',i,'j',j,'x',x(j,i)
d          print *,' '
           XI = XI - 1.0
          END DO
        XI = 5.0
        END DO
        XI =1.0
        do J = 17,18
          do I = 1,3
           X(J,I) = XI/10.0            !POINTS TO BE PLOTTED ON THE X-AXIS
           P(J,I) = XI/10.0
d          print *,'i',i,'j',j,'x',x(j,i)
d          print *,' '
           XI = XI - 1.0
          end do
        XI = 1.0
        end do
        NY = 4
        NX = 5
        do J=1,16
         do I = 1,11  
          Y(J,I) = DIFF(NX,NY)
          Q(J,I) = NORMB(NX,NY)
d         print *,'i',i,'j',j,'y',y(j,i)
d         print *,' '
          NX = NX - 1
         END DO
        NX = 5
        NY = NY - 1
        END DO
        NY = -12
        NX = 1
        do J=17,18
         do I = 1,3  
          Y(J,I) = DIFF(NX,NY)
          Q(J,I) = NORMB(NX,NY)
d         print *,'i',i,'j',j,'y',y(j,i)
d         print *,' '
          NX = NX - 1
         END DO
        NX = 1
        NY = NY - 1
        END DO
       XLAB = 'X VALUES (INCHES) '
       YLAB = 'OBSERVED FIELD VALUES'   
       WRITE(6,*) 'ENTER LNPT,LNK'
       READ(5,*) LNPT,LNK
       NC = 3
       CALL PLOTR4(NC,X,Y,NPT,NLN,TITLE,FIG,XLAB,YLAB,LNPT,LNK)
       XLAB = 'X VALUES (INCHES)'
       YLAB = 'CALCULATED FIELD VALUES'
       NC = 4
       CALL PLOTR4(NC,X,Q,NPT,NLN,TITLE,FIG,XLAB,YLAB,LNPT,LNK)
       CALL FLUSH
      return
      END