1. THIS IS THE DATAPLOT PROGRAM FILE     MATTINGL.DP
2. PURPOSE--YOUDEN PLOT FOR INTERLAB STUDY ON FLOW METER CALIBRATION
3. GEORGE MATTINGLY  3/93
4.
5. -----START POINT-----
6.
7DEVICE 2 POSTSCRIPT
8LET STRING MARK = MATTINGL.DP
9.
10. -----------------------------------
11. STEP 1--DEFINE LAB IDENTIFIERS
12. -----------------------------------
13.
14READ STRING S1 TO S18
15B C D F G H I J K L M N O P Q R S *
16.
17. -----------------------------------
18. STEP 2--READ IN THE DATA
19. -----------------------------------
20SKIP 25
21READ MATTINGL.DAT Y1 Y2 LAB
22.
23. -----------------------------------
24. STEP 3--GENERATE A YOUDEN PLOT
25. -----------------------------------
26.
27Y1LABEL S/N 3082
28X1LABEL S/N 37846
29X3LABEL AUTOMATIC
30TITLE AN 8 - 4 DOWNSTREAM METERS - HIGH FLOW
31CHAR B C D F G H I J K L M N O P Q R S *
32LINES BLANK ALL
33YOUDEN PLOT Y1 Y2 LAB
34CALL MARK.DP
35.
36. ------------------------------------------------------
37. STEP 4--GENERATE A BLOCK PLOT OF THE RAW DATA
38. ------------------------------------------------------
39.
40ER
41LET Y = Y1
42EXTEND Y Y2
43LET LAB2 = LAB
44EXTEND LAB2 LAB
45LET METER2 = SEQUENCE 1 18 1 2
46.
47Y1LABEL FLOW
48X1LABEL LABORATORY
49X2LABEL ORDERED ALPHABETICALLY BY LABORATORY ID
50LEGEND 1 1 = METER 3082
51LEGEND 2 2 = METER 37846
52LEGEND 1 COORINDATES 17 87
53LEGEND 2 COORDINATES 17 84
54CHAR BLANK ALL
55LINES SOLID ALL
56CHAR 1 2
57LINES BLANK BLANK
58X1TICS OFF
59X1TIC LABELS OFF
60BLOCK PLOT Y LAB2 METER2
61.
62FEEDBACK OFF
63JUSTIFICATION CENTER
64LOOP FOR K = 1 1 18
65MOVEDATA K 27850
66LET L = LAB2(K)
67TEXT ^S^L
68END LOOP
69FEEDBACK ON
70CALL MARK.DP
71.
72. ------------------------------------------------------
73. STEP 5--GENERATE A BLOCK PLOT OF RESIDUALS FROM LOCAL MEAN
74. ------------------------------------------------------
75.
76ER
77LET M1 = MEAN Y SUBSET METER2 1
78LET M2 = MEAN Y SUBSET METER2 2
79LET RES = Y-M1 SUBSET METER2 1
80LET RES = Y-M2 SUBSET METER2 2
81Y1LABEL RESIDUAL (= FLOW - METER MEAN)
82X2LABEL ORDERED BY LOCATION
83BLOCK PLOT RES LAB2 METER2
84.
85FEEDBACK OFF
86JUSTIFICATION CENTER
87LOOP FOR K = 1 1 18
88MOVEDATA K -110
89LET L = LAB2(K)
90TEXT ^S^L
91END LOOP
92FEEDBACK ON
93CALL MARK.DP
94.
95. ------------------------------------------------------
96. STEP 6--GENERATE AN ORDERED BLOCK PLOT OF RESIDUALS FROM LOCAL MEAN
97.         WITH ORDERING BASED ON LOCATION
98. ------------------------------------------------------
99.
100ER
101LET LAB3 = LAB2
102LET RES3 = SORTC RES LAB3
103LET LAB4 = DISTINCT LAB3
104LET LAB5 = COCODE LAB2 LAB4
105BLOCK PLOT RES LAB5 METER2
106.
107FEEDBACK OFF
108JUSTIFICATION CENTER
109LOOP FOR K = 1 1 18
110MOVEDATA K -110
111LET L = LAB4(K)
112TEXT ^S^L
113END LOOP
114FEEDBACK ON
115CALL MARK.DP
116.
117. ------------------------------------------------------
118. STEP 7--GENERATE AN ORDERED BLOCK PLOT OF RESIDUAL FROM LOCAL MEAN
119.         WITH ORDERING BASED ON SIGNED (= DIRECTIONAL) VARIATION
120. ------------------------------------------------------
121.
122ER
123LET LAB6 = LAB2
124LET DEL = Y1-Y2
125LET DEL2 = DEL
126EXTEND DEL2 DEL
127LET DEL6 = SORTC DEL2 LAB6
128LET LAB7 = DISTINCT LAB6
129LET LAB8 = COCODE LAB2 LAB7
130X2LABEL ORDERED BY SIGNED VARIATION
131BLOCK PLOT RES LAB8 METER2
132.
133FEEDBACK OFF
134JUSTIFICATION CENTER
135LOOP FOR K = 1 1 18
136MOVEDATA K -110
137LET L = LAB7(K)
138TEXT ^S^L
139END LOOP
140FEEDBACK ON
141CALL MARK.DP
142.
143. ------------------------------------------------------
144. STEP 8--GENERATE AN ORDERED BLOCK PLOT OF RESIDUAL FROM LOCAL MEAN
145.         WITH ORDERING BASED ON MAGNITUDE OF VARIATION
146. ------------------------------------------------------
147.
148ER
149LET LAB6 = LAB2
150LET DEL = (Y1-M1)-(Y2-M2)
151LET DEL2 = DEL
152EXTEND DEL2 DEL
153LET ABSDEL2 = ABS(DEL2)
154LET ABSDEL6 = SORTC ABSDEL2 LAB6
155LET LAB7 = DISTINCT LAB6
156LET LAB8 = COCODE LAB2 LAB7
157X2LABEL ORDERED BY MAGNITUDE OF VARIATION
158BLOCK PLOT RES LAB8 METER2
159.
160FEEDBACK OFF
161JUSTIFICATION CENTER
162LOOP FOR K = 1 1 18
163MOVEDATA K -110
164LET L = LAB7(K)
165TEXT ^S^L
166END LOOP
167FEEDBACK ON
168CALL MARK.DP
169.
170