Fundamental Technologies Voyager LECP Pages

An Analysis of the Performance of the Magnetic Deflection System in the Voyager Low Energy Charged Particle Experiment

by Sheela Shodhan

E.20 FDLNPLOTBLN5

*********************************************************************
*                PROGRAM FDLNPLOT                                   *
* PURPOSE : THIS ROUTINE PLOTS THE PROJECTION OF THE MAGNETIC FIELD *
*          LINES ON A PLANE DUE TO TWO TILTED MAGNETS.              *
* VARIABLES :                                                       *
* FNAME : NAME OF THE FILE FROM WHICH THE COORDINATES ARE TO BE PLOT*
*        -TED.                                                      *
* Ho,Vo : COORDINATES TO BE PLOTTED ON THE X-AXIS AND THE Y-AXIS RES*
*        -PECTIVELY.                                                *
* NPT : 1-D ARRAY; (I) = NUMBER OF POINTS IN THE LINE (I).          *
* XYZ : DOUBLE PRECISION 1-D ARRAY. EACH LINE OF THE FILE CONTAINS  *
*       X,Y,Z,Bx,By,Bz WHICH ARE READ INTO THIS ARRAY.              *
* XYZ1 : SINGLE PRECISION 1-D ARRAY WITH THE VALUES OF XYZ.         *
* X,Y : 2-D ARRAYS WHOSE FOR EACH ROW I: COL. 1 >>> COL. NPT(I) CONT*
*      -AINS THE POINTS TO BE PLOTTED ON THE X-AXIS AND THE Y-AXIS  *
*       RESPECTIVELY.                                               *
*********************************************************************

      PROGRAM FDLNPLOT

      PARAMETER (NB=50,NN=1000)
      REAL XYZ(6),X(NB,NN),Y(NB,NN)
      REAL*8 XYZ1(6)
      INTEGER NPT(NB)
      CHARACTER*72 FNAME,TITLE,XLAB,YLAB

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

       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
 3     READ(1,19,ERR=2) (XYZ1(M),M=1,6)
       DO M=1,6
        XYZ(M) = SNGL(XYZ1(M))
       END DO
       NPOINT = NPOINT+1
       X(NLN,NPOINT) = XYZ(Ho)
       Y(NLN,NPOINT) = XYZ(Vo)
       NPT(NLN) = NPOINT
       GO TO 3
 70    CLOSE(1)
       PRINT *,'NLN: ',NLN       
       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 FLUSH

 5    FORMAT(X,'ENTER THE NAME OF THE FILE: ',I3)
 10   FORMAT(A72)
 12   FORMAT(X,'X---1',/, ' Y---2',/,' Z---3',/,
     1     'BX---4',/,' BY---5',/,' BZ---6')
 19   FORMAT(X,6(2X,D18.11))

      STOP 
      END
C--------------------------------------------------------------------
C--------------------------------------------------------------------

Return to thesis table of contents.

Return to Voyager LECP Data Analysis Handbook Table of Contents.
Return to Fundamental Technologies Home Page.

Last modified 12/9/02, Tizby Hunt-Ward
tizby@ftecs.com