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.19 PLOTR4MYY1

********************************************************************
* THIS ROUTINE IS SUPPLIED WITH THE MINIMUM AND MAXIMUM VALUES     *
* OF X AND Y TO BE PLOTTED ON THE X AND THE Y AXIS.                *
* PRIMARILY, TO PLOT THE EXPT. VALUES AND THE CALCULATED VALUES    *
* OF THE B FIELD.                                                  *
********************************************************************
      SUBROUTINE PLOTR4(NC,X,Y,NPT,NLN,TITLE,XLAB,YLAB,LNPT,LNK)
      PARAMETER (NB=50,NN=1000)
      CHARACTER*72 TITLE,XLAB,YLAB
      INTEGER NPT(NB),MTY(NB)
      REAL X(NB,NN),Y(NB,NN),X1(NN),Y1(NN),SP(4,4),ST(4)
      DATA (SP(1,J),J=1,4)/0.1,0.48,0.55,0.93/
      DATA (SP(2,J),J=1,4)/0.6,0.98,0.55,0.93/
      DATA (SP(3,J),J=1,4)/0.1,0.48,0.08,0.46/
      DATA (SP(4,J),J=1,4)/0.6,0.98,0.08,0.46/
      DATA (MTY(I),I=1,8)/1,2,3,4,5,6,7,8/

      LPOWX=0
      LPOWY=0
      IF (LNPT.GT.10.AND.LNPT.LT.15) THEN
         LTYP=LNPT-10
      ELSE
         LTYP=LNPT
      END IF
C      CALL CMAXMIN(X,NPT,NLN,XMIN,XMAX)
C      CALL CMAXMIN(Y,NPT,NLN,YMIN,YMAX)
       IF ((NC .EQ. 1) .OR. (NC .EQ. 2)) THEN
         XMIN=-1.1
         XMAX=0.5
       ELSE
         XMIN = -0.9
         XMAX = 0.8
       END IF
       YMIN = 20.0
       YMAX = 1765.0
D     WRITE(7,*) ' NC = ',NC
D     WRITE(7,6) XMIN,XMAX,YMIN,YMAX
6     FORMAT(' XMIN=',E11.4,' XMAX=',E11.4,' YMIN=',E11.4,
     X       ' YMAX=',E11.4) 
      IF (LTYP.EQ.1) THEN
         CALL NEAR(2,XMAX,XMIN,UMAX,UMIN,MJTKX,MNTKX,LPOWX)
         CALL NEAR(2,YMAX,YMIN,VMAX,VMIN,MJTKY,MNTKY,LPOWY)
      ELSE
         IF (LTYP.EQ.2) THEN
            CALL NEAR(2,XMAX,XMIN,UMAX,UMIN,MJTKX,MNTKX,LPOWX)
            CALL NEARLOG(YMAX,YMIN,VMAX,VMIN,MJTKY,MNTKY)
         ELSE
            IF (LTYP.EQ.3) THEN
               CALL NEARLOG(XMAX,XMIN,UMAX,UMIN,MJTKX,MNTKX)
               CALL NEAR(2,YMAX,YMIN,VMAX,VMIN,MJTKY,MNTKY,LPOWY)
            ELSE
               CALL NEARLOG(YMAX,YMIN,VMAX,VMIN,MJTKY,MNTKY)
               CALL NEARLOG(XMAX,XMIN,UMAX,UMIN,MJTKX,MNTKX)
            END IF
         END IF
      END IF
D     WRITE(7,6) XMIN,XMAX,YMIN,YMAX
D     WRITE(7,6) UMIN,UMAX,VMIN,VMAX
  

      CALL HEADER4(TITLE)
      CALL SET(SP(NC,1),SP(NC,2),SP(NC,3),SP(NC,4),XMIN,XMAX,
     XYMIN,YMAX,LTYP)
      CALL OPTN('SSIZE',2)
      CALL OPTN('COLOR','WHITE')
      IF (LNPT.LT.10) THEN
         CALL PERIM(MJTKX,MNTKX,MJTKY,MNTKY)
      ELSE
         CALL GRID(MJTKX,MNTKX,MJTKY,MNTKY)
      END IF
      IF (NC.EQ.1) CALL OPTN('COLOR','YELLOW')     
      IF (NC.EQ.2) CALL OPTN('COLOR','RED')
      IF (NC.EQ.3) CALL OPTN('COLOR','GREEN')
      IF (NC.EQ.4) CALL OPTN('COLOR','PINK')
      DO I=1,NLN
         IF (LNK.NE.11) THEN
            NTY=LNK
         ELSE
            NTY=MTY(I)
         END IF 
         DO J=1,NPT(I)
            X1(J)=X(I,J)
            Y1(J)=Y(I,J)
            IF (NTY.EQ.0) CALL PWRITX(X1(J),Y1(J),'"140"',5,18,0,0)
            IF (NTY.EQ.1) CALL PWRITX(X1(J),Y1(J),'"632"',5,15,0,0)
            IF (NTY.EQ.2) CALL PWRITX(X1(J),Y1(J),'"1246"',6,14,0,0)
            IF (NTY.EQ.3) CALL PWRITX(X1(J),Y1(J),'"531"',5,15,0,0)
            IF (NTY.EQ.4) CALL PWRITX(X1(J),Y1(J),'"1245"',6,15,0,0)
            IF (NTY.EQ.5) CALL PWRITX(X1(J),Y1(J),'"532"',5,15,0,0)
            IF (NTY.EQ.6) CALL PWRITX(X1(J),Y1(J),'"1253"',6,15,0,0)
            IF (NTY.EQ.7) CALL PWRITX(X1(J),Y1(J),'"47"',4,15,0,0)
            IF (NTY.EQ.8) CALL PWRITX(X1(J),Y1(J),'"1307"',6,15,0,0)
         END DO
         IF (NTY.EQ.-1.OR.LNK.EQ.11) CALL CURVE(X1,Y1,NPT(I))
      END DO
      CALL OPTN('COLOR','WHITE')
      ST(1)=SP(NC,1)
      ST(2)=SP(NC,2)
      ST(3)=SP(NC,3)
      ST(4)=SP(NC,4)
      CALL LABEL(4,ST,UMIN,UMAX,VMIN,VMAX,MJTKX,MJTKY,13,
     XLTYP,LPOWX,LPOWY)
      CALL XLABEL4(XLAB,NC,LPOWX)
      CALL YLABEL4(YLAB,NC,LPOWY)
      RETURN
      END
C-----------------------------------------------------------
        SUBROUTINE HEADER4(TITLE)
        CHARACTER TITLE*72

        NT=LENCHAR(TITLE)
        CALL SET(0.1,0.9,0.955,0.99,0.0,60.0,0.0,2.0,1)
        CALL PWRITX(30.0,0.5,%REF(TITLE),NT,18,0,0)
        RETURN
        END
C----------------------------------------------------------
        SUBROUTINE XLABEL4(XLAB,NC,LPOWX)
        CHARACTER XLAB*72,POW*25
	REAL XL(4,4)
	DATA (XL(1,J),J=1,4)/0.1,0.5,0.48,0.50/
	DATA (XL(2,J),J=1,4)/0.6,1.0,0.48,0.50/
	DATA (XL(3,J),J=1,4)/0.1,0.5,0.02,0.04/
	DATA (XL(4,J),J=1,4)/0.6,1.0,0.02,0.04/

        NXC=LENCHAR(XLAB)
        CALL SET(XL(NC,1),XL(NC,2),XL(NC,3),XL(NC,4),0.0,60.0,0.0,2.0,1)
        CALL PWRITX(30.0,0.5,%REF(XLAB),NXC,14,0,0)  
        IF (LPOWX.NE.0) THEN
           ENCODE(19,12,POW) LPOWX
           CALL PWRITX(55.0,0.5,%REF(POW),19,14,0,0)
        END IF
  12    FORMAT('  "137"10"S3"',I3,'"N"')
        RETURN
        END
C ----------------------------------------------------------
        SUBROUTINE YLABEL4(YLAB,NC,LPOWY)
        CHARACTER YLAB*72,POW*25

	REAL YL(4,4)
	DATA (YL(1,J),J=1,4)/0.00,0.04,0.58,0.93/
	DATA (YL(2,J),J=1,4)/0.50,0.54,0.58,0.93/
	DATA (YL(3,J),J=1,4)/0.00,0.04,0.10,0.46/
	DATA (YL(4,J),J=1,4)/0.50,0.54,0.10,0.46/

        NYC=LENCHAR(YLAB)
        CALL SET(YL(NC,1),YL(NC,2),YL(NC,3),YL(NC,4),0.0,2.0,0.0,60.0,1)
        CALL PWRITX(0.5,30.0,%REF(YLAB),NYC,14,90,0)
        IF (LPOWY.NE.0) THEN
           ENCODE(19,12,POW) LPOWY
           CALL PWRITX(0.5,50.0,%REF(POW),19,14,90,0)
        END IF
  12    FORMAT('  "137"10"S3"',I3,'"N"')
        RETURN
        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