Voyager LECP Data Analysis Handbook

Instrument Modeling Reports

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.

VOYAGER 1 ELAPSED TIME

VOYAGER 2 ELAPSED TIME

QUICK FACTS

RELATED LINKS