Trailing-Edge
-
PDP-10 Archives
-
decuslib10-02
-
43,50247/genplt.f4
There are no other files named genplt.f4 in the archive.
SUBROUTINE GENPLT (MODE1 ,MODE2 ,MODE3 ,XPOINT,YPOINT,KRDBGN,
1KRDEND,XFIRST,YFIRST,XFINAL,YFINAL,FIRSTX,FIRSTY,FINALX,FINALY,
2XWIDTH,YWIDTH,XCORNR,YCORNR,RELHIT,LETTER,LASTX ,LASTY )
C 09/18/67
C BASIC DATA ARRAY PLOTTING SUBROUTINE FOR THE GENPLT-II PACKAGE
C
C DONALD BARTH, C/O K.B. WIBERG, DEPT. OF CHEMISTRY, YALE UNIVERSITY
C
C SUBROUTINE GENPLT ARGUMENT LIST DEFINITIONS.
C
C MODE1 = 1 DOES NOT CHECK IF POINTS ARE IN REQUESTED PLOT AREA.
C MODE1 = 2 SETS SCALES TO FIT COORDINATES OF PLOTTED POINTS.
C MODE1 = 3 DOES NOT PLOT POINT COORDINATES EXCEEDING LIMITS.
C
C MODE2 = 1 REQUIRES BOTH INPUT X AND Y COORDINATE ARRAYS.
C MODE2 = 2 SIMULATES X ARRAY LEFT TO RIGHT AT EQUAL INTERVALS.
C
C MODE3 = 1 GIVES PLOT WITH LINES CONNECTING AREAS AROUND POINTS.
C MODE3 = 2 GIVES PLOT WITH LINE CONNECTING POINTS.
C MODE3 = 3 GIVES PLOT WITH DOT AT EACH POINT.
C MODE3 = 4 GIVES PLOT WITH ALPHANUMERIC CHARACTER AT EACH POINT.
C MODE3 = 5 GIVES PLOT WITH SPECIAL SYMBOL AT EACH POINT.
C
C XPOINT = ARRAY OF X COORDINATES OF POINTS TO BE PLOTTED.
C YPOINT = ARRAY OF Y COORDINATES OF POINTS TO BE PLOTTED.
C KRDBGN = SUBSCRIPT OF FIRST POINT TO BE PLOTTED.
C KRDEND = SUBSCRIPT OF FINAL POINT TO BE PLOTTED.
C
C XFIRST = X COORDINATE TO BE PLACED AT LEFT BORDER OF PLOT AREA.
C YFIRST = Y COORDINATE TO BE PLACED AT LOWER BORDER OF PLOT AREA.
C XFINAL = X COORDINATE TO BE PLACED AT RIGHT BORDER OF PLOT AREA.
C YFINAL = Y COORDINATE TO BE PLACED AT UPPER BORDER OF PLOT AREA.
C
C SCALES MAY EITHER INCREASE OR DECREASE IF MODE1 = 1 OR 3.
C SCALES ARE SELECTED TO INCREASE TOWARD RIGHT AND TOP IF MODE1 = 2.
C INPUT LIMITS XFIRST,YFIRST,XFINAL,YFINAL ARE IGNORED IF MODE1 = 2.
C
C FIRSTX = X COORDINATE PLACED AT LEFT BORDER OF PLOT AREA.
C FIRSTY = Y COORDINATE PLACED AT LOWER BORDER OF PLOT AREA.
C FINALX = X COORDINATE PLACED AT RIGHT BORDER OF PLOT AREA.
C FINALY = Y COORDINATE PLACED AT UPPER BORDER OF PLOT AREA.
C FIRSTX, FIRSTY, FINALX AND FINALY ARE OUTPUT ARGUMENTS.
C THESE VARIABLES MAY BE USED AS LIMITING COORDINATES FOR SCLPLT.
C THESE WILL VARY FROM THE INPUT VARIABLES ONLY IF MODE1 = 2.
C
C XWIDTH = HORIZONTAL WIDTH OF PLOT AREA.
C YWIDTH = VERTICAL WIDTH OF PLOT AREA.
C XCORNR = X DISTANCE FROM PLOTTING TABLE ORIGIN TO LEFT PLOT EDGE.
C YCORNR = Y DISTANCE FROM PLOTTING TABLE ORIGIN TO LOWER PLOT EDGE.
C
C RELHIT = HEIGHT IN UNITS OF 0.02*YWIDTH OF THE LETTERS, SYMBOLS,
C OR BLANK AREAS CENTERED AT EACH POINT FOR MODE3=1, 4 OR 5
C RELHIT = 1.0 GIVES LETTERS YWIDTH/50 HIGH.
C
C LETTER = ALPHANUMERIC CHARACTER CENTERED AT POINTS IF MODE3 = 4.
C THIS IS FIRST LETTER READ BY AN A FORMAT OR IN HOLLERITH FIELD.
C LETTER = INTEGER SELECTING SPECIAL SYMBOL IF MODE3 = 5.
C SYMBOLS ARE SPECIFIED BY LETTER VALUES 1,2=TRIANGLE 3,4=SQUARE
C 5,6=PENTAGON 7,8=HEXAGON 9,10=OCTAGON AND 11,12=CROSS.
C
C LASTX = LAST X GRID COORDINATE PLOTTED.
C LASTY = LAST Y GRID COORDINATE PLOTTED.
C
DIMENSION XPOINT(1000),YPOINT(1000),LCT(12),KRD(76),KRDX(12),
1KRDY(12)
COMMON/PPARM/FACTOR,OFSETX,OFSETY,IERR,IPEN,NTAPE,MODE,IPOINT,
1IFREER,ILINE
DATA (LCT(I),I=1,12)/0103,0403,0704,1104,1505,2005,2506,3106,3708,
1 4508,5312,6512/
DATA (KRD(I),I=1,76)/00,36,60,06,66,30,00,06,66,60,03,36,63,30,04,
136,64,50,10,16,56,62,30,02,03,26,46,63,40,20,04,36,64,62,30,02,04,
226,46,64,62,40,20,02,03,15,36,55,63,51,30,11,04,24,26,46,44,64,62,
342,40,20,22,02,05,16,34,56,65,43,61,50,32,10,01,23/
C
C ************CHECK FOR LIMITS OF INPUT COORDINATE DATA*************
IF(MODE1-2)11,1,11
1 YLOWER=YPOINT(KRDBGN)
YUPPER=YPOINT(KRDBGN)
DO 5 I=KRDBGN,KRDEND
IF(YLOWER-YPOINT(I))3,3,2
2 YLOWER=YPOINT(I)
3 IF(YUPPER-YPOINT(I))4,5,5
4 YUPPER=YPOINT(I)
5 CONTINUE
IF(MODE2-2)6,12,91
6 XLEFT=XPOINT(KRDBGN)
XRIGHT=XPOINT(KRDBGN)
DO 10 I=KRDBGN,KRDEND
IF(XLEFT-XPOINT(I))8,8,7
7 XLEFT=XPOINT(I)
8 IF(XRIGHT-XPOINT(I))9,10,10
9 XRIGHT=XPOINT(I)
10 CONTINUE
GO TO 13
11 YLOWER=YFIRST
YUPPER=YFINAL
12 XLEFT=XFIRST
XRIGHT=XFINAL
13 FIRSTX=XLEFT
FINALX=XRIGHT
FIRSTY=YLOWER
FINALY=YUPPER
C
C **************CALCULATE SCALING FACTORS AND OFFSETS***************
XSCALE=FACTOR*XWIDTH/(XRIGHT-XLEFT)
YSCALE=FACTOR*YWIDTH/(YUPPER-YLOWER)
XMIN=XCORNR*FACTOR+OFSETX
YMIN=YCORNR*FACTOR+OFSETY
XMAX=XMIN+(FACTOR*XWIDTH)
YMAX=YMIN+(FACTOR*YWIDTH)
DIF=FACTOR*XWIDTH/(FLOAT(KRDEND-KRDBGN))
IF(MODE3-2)14,15,70
C
C ****************************LINE PLOT*****************************
14 HALF=RELHIT*FACTOR*YWIDTH/100.0
WIDTH=1.0+(2.0*HALF)
15 LINE=-1
MODE=IFREER
IF(MODE2-2)17,16,91
16 THISX=XMIN
GO TO 18
17 THISX=XMIN+(XSCALE*(XPOINT(KRDBGN)-XLEFT))
18 THISY=YMIN+(YSCALE*(YPOINT(KRDBGN)-YLOWER))
INDEX=KRDBGN
19 INDEX=INDEX+1
IF(INDEX-KRDEND)20,20,91
20 THATX=THISX
THATY=THISY
IF(MODE2-2)22,21,91
21 THISX=XMIN+(DIF*(FLOAT(INDEX-KRDBGN)))
GO TO 23
22 THISX=XMIN+(XSCALE*(XPOINT(INDEX)-XLEFT))
23 THISY=YMIN+(YSCALE*(YPOINT(INDEX)-YLOWER))
DISTX=THISX-THATX
DISTY=THISY-THATY
C
C ********DETERMINE COORDINATES IF LINES CONNECT BLANK AREAS********
IF(MODE3-1)91,25,24
24 STARTX=THATX
STARTY=THATY
HALTX=THISX
HALTY=THISY
IF(MODE1-3)43,46,91
25 IF(DISTX)27,26,26
26 IF(DISTY)31,28,28
27 IF(DISTY)29,30,30
28 IF(DISTX-DISTY)33,32,32
29 IF(DISTX-DISTY)34,35,35
30 IF(DISTX+DISTY)34,33,33
31 IF(DISTX+DISTY)35,32,32
32 IF(DISTX-WIDTH)19,36,36
33 IF(DISTY-WIDTH)19,39,39
34 IF(DISTX+WIDTH)37,37,19
35 IF(DISTY+WIDTH)40,40,19
36 STARTX=THATX+HALF
HALTX=THISX-HALF
GO TO 38
37 STARTX=THATX-HALF
HALTX=THISX+HALF
38 STARTY=THATY+((DISTY*(STARTX-THATX))/DISTX)
HALTY=THISY+((DISTY*(HALTX-THISX))/DISTX)
GO TO 42
39 STARTY=THATY+HALF
HALTY=THISY-HALF
GO TO 41
40 STARTY=THATY-HALF
HALTY=THISY+HALF
41 STARTX=THATX+((DISTX*(STARTY-THATY))/DISTY)
HALTX=THISX+((DISTX*(HALTY-THISY))/DISTY)
42 LINE=-1
IF(MODE1-3)43,46,91
C
C *******************PLOT THE LINE WITHOUT TESTS********************
43 IF(LINE)44,44,45
44 LINE=1
NEWX=STARTX
NEWY=STARTY
CALL PENUP (LASTX,LASTY,NEWX,NEWY)
45 NEWX=HALTX
NEWY=HALTY
CALL PENDWN(LASTX,LASTY,NEWX,NEWY)
GO TO 19
C
C ***************TEST HORIZONTAL COORDINATES OF LINE****************
46 IF(MODE2-2)47,57,91
47 IF(STARTX-XMIN)48,49,49
48 STARTX=XMIN
LINE=-1
IF(HALTX-XMIN)19,53,53
49 IF(STARTX-XMAX)51,51,50
50 STARTX=XMAX
LINE=-1
IF(HALTX-XMAX)51,51,19
51 IF(HALTX-XMIN)52,53,53
52 HALTX=XMIN
GO TO 55
53 IF(HALTX-XMAX)55,55,54
54 HALTX=XMAX
55 IF(ABS(DISTX)-0.1)57,56,56
56 STARTY=THATY+((DISTY*(STARTX-THATX))/DISTX)
HALTY=THISY-((DISTY*(THISX-HALTX))/DISTX)
C
C ****************TEST VERTICAL COORDINATES OF LINE*****************
57 IF(STARTY-YMIN)58,59,59
58 STARTY=YMIN
LINE=-1
IF(HALTY-YMIN)19,63,63
59 IF(STARTY-YMAX)61,61,60
60 STARTY=YMAX
LINE=-1
IF(HALTY-YMAX)61,61,19
61 IF(HALTY-YMIN)62,63,63
62 HALTY=YMIN
GO TO 65
63 IF(HALTY-YMAX)65,65,64
64 HALTY=YMAX
65 IF(ABS(DISTY)-0.1)67,66,66
66 STARTX=THATX+((DISTX*(STARTY-THATY))/DISTY)
HALTX=THISX-((DISTX*(THISY-HALTY))/DISTY)
C
C *******************PLOT THE TESTED LINE SEGMENT*******************
67 IF(LINE)68,68,69
68 LINE=1
NEWX=STARTX
NEWY=STARTY
CALL PENUP (LASTX,LASTY,NEWX,NEWY)
69 NEWX=HALTX
NEWY=HALTY
CALL PENDWN(LASTX,LASTY,NEWX,NEWY)
GO TO 19
C
C ************PREPARE FOR POINT, LETTER OR SYMBOL PLOT**************
70 KSCALE=RELHIT*FACTOR*YWIDTH/300.0
MINX=XMIN
MAXX=XMAX
MINY=YMIN
MAXY=YMAX
IF(MODE3-4)71,72,73
71 MODE=IPOINT
GO TO 77
72 NEWX=0
CALL LTRPLT (LETTER,1,1,0,KSCALE,LASTX,LASTY,NEWX,NEWY)
NEWX=(NEWX-2*KSCALE)/2
XMIN=XMIN-FLOAT(NEWX)
YMIN=YMIN-FLOAT(3*KSCALE)
MINX=XMIN
MINY=YMIN
MAXX=MAXX-NEWX
MAXY=MAXY-3*KSCALE
MODE=ILINE
GO TO 77
73 IF(LETTER)91,91,74
74 IF(LETTER-12)75,75,91
75 INDEX=LCT(LETTER)/100
NUM=LCT(LETTER)-100*INDEX
DO 76 I=1,NUM
KRDX(I)=KRD(INDEX)/10
KRDY(I)=KRD(INDEX)-10*KRDX(I)
KRDX(I)=KSCALE*KRDX(I)
KRDY(I)=KSCALE*KRDY(I)
76 INDEX=INDEX+1
XMIN=XMIN-FLOAT(3*KSCALE)
YMIN=YMIN-FLOAT(3*KSCALE)
MINX=XMIN
MINY=YMIN
MAXX=MAXX-3*KSCALE
MAXY=MAXY-3*KSCALE
MODE=ILINE
C
C **************FORM THE POINT, LETTER OR SYMBOL PLOT***************
77 INDEX=KRDBGN-1
78 INDEX=INDEX+1
IF(INDEX-KRDEND)79,79,91
79 IF(MODE2-2)80,81,91
80 NEWX=XMIN+(XSCALE*(XPOINT(INDEX)-XLEFT))
NEWY=YMIN+(YSCALE*(YPOINT(INDEX)-YLOWER))
IF(MODE1-3)86,82,91
81 NEWX=XMIN+(DIF*(FLOAT(INDEX-KRDBGN)))
NEWY=YMIN+(YSCALE*(YPOINT(INDEX)-YLOWER))
IF(MODE1-3)86,84,91
82 IF(NEWX-MINX)78,84,83
83 IF(NEWX-MAXX)84,84,78
84 IF(NEWY-MINY)78,86,85
85 IF(NEWY-MAXY)86,86,78
86 IF(MODE3-4)87,88,89
87 CALL PENUP (LASTX,LASTY,NEWX,NEWY)
GO TO 78
88 CALL LTRPLT (LETTER,1,1,1,KSCALE,LASTX,LASTY,NEWX,NEWY)
GO TO 78
89 NOWX=NEWX+KRDX(NUM)
NOWY=NEWY+KRDY(NUM)
CALL PENUP (LASTX,LASTY,NOWX,NOWY)
DO 90 I=1,NUM
NOWX=NEWX+KRDX(I)
NOWY=NEWY+KRDY(I)
90 CALL PENDWN(LASTX,LASTY,NOWX,NOWY)
GO TO 78
91 RETURN
END