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