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.16 PLOTR5

**************************************************************************
      SUBROUTINE PLOTR5(X,Y,NPT,NLN,TITLE,XLAB,YLAB,LNPT,LNK)
*                                                                        * 
*     THIS SUBROUTINE PLOTS THE DATA IN THE ARRAYS X AND Y               *
*     BY USING SOME NCAR SUBROUTINES.                                    *
*     THIS IS SLIGHTLY MODIFIED FROM THE ORIGINAL TO PLOT THE PROJECTION *
*     OF THE SENSOR SUBSYSTEM ON THE XY PLANE. THIS PLOTS POINTS TO SHOW *
*     THE DIRECTION OF THE MAGNETIC FIELD INSIDE THE SENSOR SUBSYSTEM.   *
*                                                                        *
*     VARIABLES :                                                        *
*     X(I,J),Y(I,J): 2-D ARRAYS CONTAIN THE DATA. THEY SHOULD BE         *
*                    DIMENSIONED BEFORE ENTRY. NB=50,NN=1000             *
*                    THE ROW INDEX,I, REFERS TO LINES,                   *
*                    THE COLUMN INDEX,J, REFERS TO POINTS ON A LINE      *
*     NPT(I)       : 1-D ARRAY CONTAINS THE NUMBER OF POINTS IN EACH LINE*
*     NLN          : AN INTEGER FOR THE NUMBER OF LINES TO BE DRAWN.     *
*     TITLE        : THE TITLE FOR THE PLOT                              *
*     XLAB         : CAPTION IN HORIZONTAL DIRECTION                     *
*     YLAB         : CAPTION IN VERTICAL DIRECTION                       *
*     LNPT         : A FLAG INDICATING LINEAR & LOG MAPPING              *
*                    LNPT   X DIRECTION  Y DIRECTION                     *
*                     1       LINEAR        LINEAR                       *
*                     2       LINEAR        LOG                          *
*                     3       LOG           LINEAR                       *
*                     4       LOG           LOG                          *
*                    11-14 THE SAME SCHEME WILL BE USED, BUT WITH        *
*                          GRID BACKGROUND                               *
*     LNK          : A FLAG INDICATING THE LINE AND POINT SYMBOLS        *
*                    LNK    DESCRIPTION                                  *
*                     -1    SOLID CURVE                                  *
*                      0    DOTS                                         *
*                      1    CIRCLE                                       *
*                      2    SOLID SQUARE                                 *
*                      3    TRIANGLE                                     *
*                      4    SOLID CIRCLE                                 *
*                      5    DIAMOND                                      *
*                      6    SOLID TRIANGLE                               *
*                      7    CROSS                                        *
*                      8    CIRCLE WITH PLUS SIGN                        *
**************************************************************************
      PARAMETER (NB=50,NN=1000)
      CHARACTER*72 TITLE,XLAB,YLAB
      REAL X(NB,NN),Y(NB,NN),X1(NN),Y1(NN),ST(4)
      INTEGER NPT(NB),MTY(NB)
      DATA (ST(I),I=1,4)/0.15,0.90,0.20,0.90/
      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
      CALL CMAXMIN(X,NPT,NLN,XMIN,XMAX)
      CALL CMAXMIN(Y,NPT,NLN,YMIN,YMAX)
      IF (LTYP.EQ.1) THEN
         CALL NEAR(1,XMAX,XMIN,UMAX,UMIN,MJTKX,MNTKX,LPOWX)
         CALL NEAR(1,YMAX,YMIN,VMAX,VMIN,MJTKY,MNTKY,LPOWY)
      END IF
      CALL OPTN('SSIZE',2)
      CALL OPTN('COLOR','WHITE')
      CALL HEADER1(TITLE)
      CALL SET(ST(1),ST(2),ST(3),ST(4),XMIN,XMAX,YMIN,YMAX,LTYP)
      IF (LNPT.LT.10) THEN
         CALL PERIM(MJTKX,MNTKX,MJTKY,MNTKY)
      ELSE
         CALL GRID(MJTKX,MNTKX,MJTKY,MNTKY)
      END IF
      CALL OPTN('COLOR','YELLOW')
      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 (I .EQ. 1) CALL PWRITX(X1(J),Y1(J),'"140"',5,25,0,0)
         END DO
         IF (I.NE.1) CALL CURVE(X1,Y1,NPT(I))
      END DO
      CALL OPTN('COLOR','WHITE')
      CALL LABEL(1,ST,UMIN,UMAX,VMIN,VMAX,MJTKX,MJTKY,16,
     &    LTYP,LPOWX,LPOWY)
      CALL XLABEL1(XLAB,LPOWX)
      CALL YLABEL1(YLAB,LPOWY)
      RETURN
      END
C-------------------------------------------------------------------------
      SUBROUTINE HEADER1(TITLE)
      CHARACTER TITLE*72
C     
C     TO PLOT THE TITLE AND MAXIMUN CHARACTERS ARE 72           
C
      NT=LENCHAR(TITLE)
      CALL SET(0.1,0.9,0.93,1.0,0.0,60.0,0.0,2.0,1)
      CALL PWRITX(30.0,0.5,%REF(TITLE),NT,20,0,0)
      RETURN
      END
C-------------------------------------------------------------------------
      SUBROUTINE XLABEL1(XLAB,LPOWX)
      CHARACTER XLAB*72
C     
C     TO PLOT THE CAPTION X-DIRECTION AND 
C     MAXIMUN CHARACTERS ARE 72           
C
      NXC=LENCHAR(XLAB)
      CALL SET(0.1,0.95,0.07,0.15,0.0,60.0,0.0,2.0,1)
      CALL PWRITX(29.0,0.5,%REF(XLAB),NXC,17,0,0)
      IF (LPOWX.NE.0) THEN
         ENCODE(19,12,POW) LPOWX
         CALL PWRITX(56.0,0.5,%REF(POW),19,17,0,0)
      END IF
  12  FORMAT('  "137"10"S3"',I3,'"N"')
      RETURN
      END
C-------------------------------------------------------------------------
      SUBROUTINE YLABEL1(YLAB,LPOWY)
      CHARACTER YLAB*72
C     
C     TO PLOT THE CAPTION X-DIRECTION AND 
C     MAXIMUN CHARACTERS ARE 72           
C
      NYC=LENCHAR(YLAB)
      CALL SET(0.03,0.08,0.1,0.9,0.0,2.0,0.0,60.0,1)
      CALL PWRITX(0.5,29.0,%REF(YLAB),NYC,17,90,0)
      IF (LPOWY.NE.0) THEN
         ENCODE(19,12,POW) LPOWY
         CALL PWRITX(0.5,55.0,%REF(POW),19,17,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.

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