*  	 	 	 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

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