| 1 $JOB 44-0331-01 LINEAR LEAST-SQUARES CURVE FITTING PROGRAM $EXECUTE IBJOB IBJOB VERSION 5 HAS CONTROL. | $IBJOB GO,MAP,SOURCE,DECK | $IBFTC CF1 $IBFTC CF1 0000 1 12/01/64 PAGE 1 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0| C CF1 0010 | C CF1 0020 C LINEAR LEAST-SQUARES CURVE FITTING PROGRAM CF1 0030 | C CF1 0040 | C CF1 0050 | C PROGRAM WRITTEN AND REVISED AT AMERICAN OIL BY R. A. BREHMER, CF1 0060 | C P. E. PIECHOCKI, W. B. TRAVER, F. M. JACOBSEN, R. J. TOMAN CF1 0070 C AND F. S. WOOD. CF1 0080 | C CF1 0090 | C ROUTINES CF1 THROUGH CF8 ARE WRITTEN IN IBM 7044 FORTRAN IV CF1 0100 | C (32K) AND CF9 IS A MAP ASSEMBLY. CF1 0110 | C CF1 0120 C MAIN ROUTINE... READS CONTROL CARD,(FORMAT AND TRANSFOR- CF1 0130 | C MATION CARDS IF SO INDICATED ON CONTROL CARD,) DATA CARDS AND CF1 0140 | C END CARD CF1 0150 | C CHECKS FOR BZERO, E FORMAT OF PRINTOUT, CONFIDENCE INTERVAL CF1 0160 | C CALCULATION, WEIGHTING FACTORS, INDEX OF CARDS AND SEQUENCE(IF CF1 0170 C DESIRED). CF1 0180 | C LISTS OBSERVED AND TRANSFORMED DATA, CALCULATES DEGREES OF CF1 0190 | C FREEDOM, SUMS OF VARIABLES, RAW AND RESIDUAL SUMS OF SQUARES CF1 0200 | C AND CROSS PRODUCTS, RANGE OF VARIABLES, MEAN AND STANDARD CF1 0210 | C DEVIATION OF VARIABLES, SIMPLE CORRELATION COEFFICIENTS, AND CF1 0220 C PUNCHES CARDS FOR CONFIDENCE INTERVAL AND CP CALCULATIONS. CF1 0230 | C CALCULATES COEFFICIENTS OF EACH INDEPENDENT VARIABLE, CF1 0240 | C CALCULATES AND LISTS STATISTICS OF FIT, CALCULATES Y VALUES AND CF1 0250 | C RESIDUALS, ALSO ORDERS RESIDUALS. CF1 0260 | C CF1 0270 DIMENSION Z(31,31), C(30,30), DATA(51) CF1 0280 | DIMENSION BETA(30), AVG(30), SIGMA(30), STDDEV(30), B(30) CF1 0290 | DIMENSION VAR(30), YRESID(30) CF1 0300 | DIMENSION ITRFM(50), VARTR(50), AVATR(50), ITLOC(50), IOMIT(50) CF1 0310 | DIMENSION ARRAY(21,6),BRRAY(70),CRRAY(56),DRRAY(4),ERRAY(51) CF1 0320 DIMENSION FMT(80),IQ(40),IDZA(1000),IDZB(1000),DELTA(1000) CF1 0330 | DIMENSION YYIN(1000),YYCC(1000),LSORT(1000) CF1 0340 | DIMENSION AMAX(50),AMIN(50),W1(50) CF1 0350 | COMMON YYMIN, YCMIN, YYMAX, YCMAX CF1 0360 | COMMON FMT, IQ, KSL1, KSL2, KSL3, KSL4, KSL5, KSL6 CF1 0370 COMMON IDZA, IDZB, DELTA, YYIN, YYCC, LSORT CF1 0380 | COMMON Z,C, BETA,AVG,SIGMA, B,VAR,YRESID CF1 0390 | COMMON ITRFM,VARTR,AVATR, IOMIT CF1 0400 | COMMON KTIN,KTOU,KTPCH,KTBIN1,KTBIN2,KTBIN3,KTBIN4,KTBIN5 CF1 0410 | COMMON PR1,PR2,PR3 CF1 0420 COMMON BZRO, KNOIND, KNODEP, KNOVAR, KNVP1, JNOIND, JNODEP, JNOVARCF1 0430 | COMMON NOBSP, NOERR, JOBSV, NRDLMT, NOOBSV, MNO, JSEQ, N, M, J, I CF1 0440 | COMMON LNOVAR, NNO, ISEQ, L, KNO, XNOOBS, XNOVAR, XNOIND, DEFR CF1 0450 | COMMON NOVM1, RTDEFR, IPL1, K, RSQD, BZERO, SUMB, ICURY, FACTOR CF1 0460 | COMMON SGBETA, TBETA, PCORK, SIGB, FRATIO, UNBEST, BIEST, SSQ, SSECF1 0470 COMMON IPLOT1, IPLOT2, YCALC, DIFF, IREC, X, IDENT, IOBSV, MP CF1 0480 | COMMON NOVAR, NVP1, IPATH CF1 0490 | EQUIVALENCE (C,DATA), (C(101),STDDEV), (C(201),ITLOC) CF1 0500 | EQUIVALENCE (NOIND,IQ(1)),(NODEP,IQ(2)) CF1 0510 | EQUIVALENCE ( KFM, IQ(22) ) CF1 0520 EQUIVALENCE(ARRAY(71),CRRAY(1)),(ARRAY(1),BRRAY(1)) CF1 0530 | C KSLX = SENSE LITE X. = 1 IF ON, 2 IF OFF CF1 0540 | CALL ATHRUZ(IWDA,6HEND ) CF1 0550 2 | 1 12/01/64 PAGE 2 | CF1 - EFN SOURCE STATEMENT - IFN(S) - 0 CALL ATHRUZ(BLANK,6H ) CF1 0560 4 | DATA BRRAY/4HNONE,4HNATU,4HCOMM,4HANTI,4HANTI,4HEXPO,4HMULT, CF1 0570 | A 4HDIVI,4HADD ,4H ,4H ,4H ,4H ,4HSINE, CF1 0580 | B 4HCOSI,4HC DI,4HLN(L,4HLOG(,4HLOG(,4HANTI,4HANTI, CF1 0590 | C 4H ,4HRAL ,4HON ,4H-NAT,4H-COM,4HNENT,4HIPLY, CF1 0600 D 4HDE B,4HCONS,4HMULT,4HDIVI,4HADD ,4HSUBT,4H, RA, CF1 0610 | E 4HNE, ,4HVIDE,4HN VA,4HLOG ,4HVARI,4H(LN(,4H(LOG, CF1 0620 | F 4H ,4HLOG ,4HLOG ,4HURAL,4HMON ,4HIATI,4H BY , CF1 0630 | G 4HY CO,4HTANT,4HIPLY,4HDE P,4HPOSI,4HRACT,4HDIAN, CF1 0640 | H 4HRADI,4HD BY,4HRIAB,4HVARI,4HABLE,4HLN V,4H(LOG, CF1 0650 I 4H ,4H ,4H ,4H LOG,4HLOG ,4HON B,4HCONS/ CF1 0660 | DATA CRRAY/4HNSTA,4H ,4H POS,4HOSIT,4HTION,4H POS,4HS , CF1 0670 | K 4HANS ,4H VAR,4HLE) ,4HABLE,4H PLU,4HARIA,4H VAR, CF1 0680 | L 4H ,4H ,4H ,4H ,4H ,4HY CO,4HTANT, CF1 0690 | M 4HNT ,4H ,4HITIO,4HIONS,4HS, ,4HITIO,4H , CF1 0700 N 4H ,4HIABL,4H ,4H) ,4HS C),4HBLE),4HIABL, CF1 0710 | O 4H ,4H ,4H ,4H ,4H ,4HNST.,4H , CF1 0720 | P 4H ,4H ,4HNS, ,4H, ,4H ,4HNS, ,4H , CF1 0730 | Q 4H ,4HE ,4H ,4H ,4H ,4H) ,4HE)) / CF1 0740 | C ARRAY PRINTS OUT THE FOLLOWING TRANSFORMATIONS CF1 0750 C 0 NONE CF1 0760 | C 1 NATURAL LOG CF1 0770 | C 2 COMMON LOG CF1 0780 | C 3 ANTI-NATURAL LOG CF1 0790 | C 4 ANTI-COMMON LOG CF1 0800 C 5 EXPONENTIATION BY CONST. CF1 0810 | C 6 MULTIPLY BY CONSTANT CF1 0820 | C 7 DIVIDE BY CONSTANT CF1 0830 | C 8 ADD BY CONSTANT CF1 0840 | C 9 MULTIPLY POSITIONS,--X-- CF1 0850 C 10 DIVIDE POSITIONS, --/-- CF1 0860 | C 11 ADD POSITIONS, --+-- CF1 0870 | C 12 SUBTRACT POSITIONS,----- CF1 0880 | C 13 SINE, RADIANS CF1 0890 | C 14 COSINE, RADIANS CF1 0900 C 15 C DIVIDED BY VARIABLE CF1 0910 | C 16 LN(LN VARIABLE) CF1 0920 | C 17 LOG(LOG VARIABLE) CF1 0930 | C 18 LOG(VARIABLE PLUS C) CF1 0940 | C 19 ANTI(LN(LN VARIABLE)) CF1 0950 C 20 ANTI(LOG(LOG VARIABLE)) CF1 0960 | C CF1 0970 | DATA DRRAY/1HX,1H/,1H+,1H-/ CF1 0980 | C CF1 0990 | C CF1 1000 C ARRAY OF NUMBERS FROM 1 TO 50, 51 IS A BLANK CF1 1010 | C CF1 1020 | DATA ERRAY/2H 1,2H 2,2H 3,2H 4,2H 5,2H 6,2H 7,2H 8,2H 9,2H10, CF1 1030 | 1 2H11,2H12,2H13,2H14,2H15,2H16,2H17,2H18,2H19,2H20, CF1 1040 | 2 2H21,2H22,2H23,2H24,2H25,2H26,2H27,2H28,2H29,2H30, CF1 1050 3 2H31,2H32,2H33,2H34,2H35,2H36,2H37,2H38,2H39,2H40, CF1 1060 | 4 2H41,2H42,2H43,2H44,2H45,2H46,2H47,2H48,2H49,2H50,2H /CF1 1070 | C CF1 1080 | C CF1 1090 | C --- A 2 IN COL 31 OF CONTROL CARD WILL CAUSE PUNCHING OF DATA CF1 1100 C --- REQUIRED FOR CP CALCULATION CF1 1110 | 1 12/01/64 PAGE 3 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0C CF1 1120 | C --- DEFINE BCD INPUT AND OUTPUT FILES AND SCRATCH FILES CF1 1130 C CF1 1140 | KTIN = 5 CF1 1150 | KTOU = 6 CF1 1160 | KTPCH = 7 CF1 1170 | KTBIN2 = 1 CF1 1180 KTBIN3 = 2 CF1 1190 | KTBIN5 = 3 CF1 1200 | C CF1 1210 | C --- READ CONTROL CARD CF1 1220 | C --- AN E IN COL 24 CAUSES LISTING OF RESIDUALS TO BE IN E FORMAT CF1 1230 C CF1 1240 | 10 READ(KTIN,800) PR1,PR2,PR3,IQ(1),IQ(2),XXFRT,(IQ(JZ),JZ=3,40) CF1 1250 11 | 800 FORMAT(3A6,2I2,1X,A1,18I1,I2,28I1) CF1 1260 | C CF1 1270 | C --- CONTROL ON CP PUNCH CF1 1280 C CF1 1290 | KCP=2 CF1 1300 | IF(IQ(9).EQ.2) KCP=1 CF1 1310 | C CF1 1320 | C --- CALCULATE NUMBER OF VARIABLES FROM CONTROL CARD VALUES CF1 1330 C CF1 1340 | NOIND=IQ(1) CF1 1350 | NODEP=IQ(2) CF1 1360 | NOVAR=IQ(1)+IQ(2) CF1 1370 | NVP1=NOVAR+1 CF1 1380 20 KNOIND=NOIND CF1 1390 | KNODEP=NODEP CF1 1400 | KNOVAR=NOVAR CF1 1410 | KNVP1 =NVP1 CF1 1420 | JNOIND=NOIND CF1 1430 JNODEP=NODEP CF1 1440 | JNOVAR=NOVAR CF1 1450 | C CF1 1460 | C --- INITIALIZE INPUT STORAGE AREAS, SWITCHES, AND FILES CF1 1470 | C CF1 1480 DO 15 J=1, 31 CF1 1490 | DO 15 I=1, 31 CF1 1500 | 15 Z(I,J) = 0. CF1 1510 | 18 IF(IQ(20)) 22,21,21 CF1 1520 | 21 REWIND KTBIN5 CF1 1530 44 22 REWIND KTBIN3 CF1 1540 45 | REWIND KTBIN2 CF1 1550 46 | KSL1=2 CF1 1560 | KSL2=2 CF1 1570 | KSL3=2 CF1 1580 KSL4=2 CF1 1590 | KSL5=2 CF1 1600 | KSL6=2 CF1 1610 | C CF1 1620 | C --- WRITE PROBLEM IDENTIFICATION ON OUTPUT FILE CF1 1630 C CF1 1640 | WRITE(KTOU,805) PR1,PR2,PR3 CF1 1650 52 | 805 FORMAT(1H1,32X,42HLINEAR LEAST-SQUARES CURVE FITTING PROGRAM// CF1 1660 | 1 1X, 3A6 ) CF1 1670 | 1 12/01/64 PAGE 4 CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0C CF1 1680 | C --- DETERMINE NUMBER OF EQUATIONS IN PROBLEM AND FORMAT OF INPUT DATA CF1 1690 | C CF1 1700 | IQF=IQ(19)+1 CF1 1710 IQE=IQ(20)+1 CF1 1720 | IF(IQE.GE.1.AND.IQE.LE.3) GO TO 26 CF1 1730 | WRITE(KTOU,24) IQ(20) CF1 1740 57 | 24 FORMAT(1H0,36HBAD VALUE FOR MODEL CONTROL. VALUE =,I2) CF1 1750 | CALL EXIT CF1 1760 58 26 GO TO (28,32,36),IQE CF1 1770 | 28 WRITE(KTOU,30) CF1 1780 60 | 30 FORMAT(1H0,24HPROBLEM HAS ONE EQUATION) CF1 1790 | GO TO 38 CF1 1800 | 32 NOEQ=1 CF1 1810 WRITE(KTOU,34) NOEQ CF1 1820 63 | 34 FORMAT(1H0, 9HEQUATION ,I2,28H OF A MULTI-EQUATION PROBLEM) CF1 1830 | GO TO 38 CF1 1840 | 36 NOEQ=NOEQ+1 CF1 1850 | WRITE(KTOU,34) NOEQ CF1 1860 66 GO TO 80 CF1 1870 | 38 IF(IQF.GE.1.AND.IQF.LE.3) GO TO 42 CF1 1880 | WRITE(KTOU,40) IQ(19) CF1 1890 70 | 40 FORMAT(1H0,37HBAD VALUE FOR FORMAT CONTROL. VALUE =,I2) CF1 1900 | CALL EXIT CF1 1910 71 42 GO TO (44,52,72),IQF CF1 1920 | 44 WRITE(KTOU,50) CF1 1930 73 | 50 FORMAT(1H0,30HDATA READ WITH STANDARD FORMAT) CF1 1940 | GO TO 80 CF1 1950 | 52 WRITE(KTOU,56) CF1 1960 75 56 FORMAT(1H0,29HDATA READ WITH SPECIAL FORMAT) CF1 1970 | IF(KFM.LE.0.OR.KFM.GT.4) KFM=1 CF1 1980 | FMT(80)=BLANK CF1 1990 | DO 62 I=1,79 CF1 2000 | FMT(I)=FMT(80) CF1 2010 62 CONTINUE CF1 2020 | DO 2022 I=1,KFM CF1 2030 | KZ=20*I CF1 2040 | LZ=KZ-19 CF1 2050 | READ(KTIN,1030) (FMT(JZ),JZ=LZ,KZ) CF1 2060 89 1030 FORMAT(20A4) CF1 2070 | WRITE(KTOU,64) I,(FMT(JZ),JZ=LZ,KZ) CF1 2080 94 | 64 FORMAT(1H0,11HFORMAT CARD,I2,1X,20A4) CF1 2090 | 2022 CONTINUE CF1 2100 | GO TO 80 CF1 2110 72 WRITE(KTOU,74) CF1 2120 103 | 74 FORMAT(1H0,31HDATA READ WIH SUBROUTINE REDATA) CF1 2130 | C CF1 2140 | C --- INDICATE ON OUTPUT FILE IF BZERO =0, OR BZERO = CALCULATED VALUE CF1 2150 | C CF1 2160 80 IF(IQ(5)) 17,17,16 CF1 2170 | 16 BZRO=0. CF1 2180 | WRITE(KTOU,980) CF1 2190 107 | 980 FORMAT(10H0BZERO = 0) CF1 2200 | GO TO 86 CF1 2210 17 BZRO=1. CF1 2220 | WRITE(KTOU,985) CF1 2230 110 | 1 12/01/64 PAGE 5 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 985 FORMAT(25H0BZERO = CALCULATED VALUE) CF1 2240 C CF1 2250 | C --- PUNCH CONTROL CARDS FOR CONFIDENCE INTERVAL CALCULATIONS CF1 2260 | C CF1 2270 | 86 IF(IQ(9)-1) 47,46,47 CF1 2280 | 46 WRITE(KTPCH,800) PR1,PR2,PR3, (IQ(JZ), JZ = 1,8) CF1 2290 C CHECK FOR TRANSFORMATION CARD INDICATOR + READ IN CF1 2300 113 | 47 IF(IQ(3) ) 54, 54, 48 CF1 2310 | 48 NOBSP =NOVAR/7 CF1 2320 | IF(MOD(NOVAR,7)) 6001,6001,6003 CF1 2330 | 6001 NOBSP = NOBSP + 1 CF1 2340 6003 J2 = 0 CF1 2350 | 6002 J2 = J2 + 7 CF1 2360 | J1 = J2 - 6 CF1 2370 | IF(J2.GT.NOVAR) J2 = NOVAR CF1 2380 | READ(KTIN,975)(ITRFM(JZ),VARTR(JZ),ITLOC(JZ),IOMIT(JZ),JZ= J1,J2) CF1 2390 131 975 FORMAT(7(I2, F5.0, I2, I1)) CF1 2400 | READ(KTIN,990) (AVATR(JZ),JZ=J1,J2 ) CF1 2410 141 | 990 FORMAT(7(2X,A5,3X)) CF1 2420 | IF(J2.LT.NOVAR) GO TO 6002 CF1 2430 | C PUNCH FOR CONFIDENCE INTERVAL CALCULATIONS CF1 2440 IF(IQ(9) - 1 ) 54, 58, 54 CF1 2450 | 58 WRITE(KTPCH,992) (ITRFM(JZ),AVATR(JZ),ITLOC(JZ),IOMIT(JZ), JZ=1, CF1 2460 | 1 NOVAR ) CF1 2470 | C CARD OUTPUT FOR TRANSFORMATIONS IF CONFIDENCE INTVL PUNCH-OUTCF1 2480 152 | 992 FORMAT(7(I2, A5, I2, I1)) CF1 2490 C WRITE HEADING LINE IF INPUT DATA TO BE LISTED CF1 2500 | 54 IF (IQ(10)) 60, 55, 60 CF1 2510 | 55 WRITE(KTOU,815) IQ(1) , IQ(2) CF1 2520 164 | 815 FORMAT(11H0DATA INPUT , 6X, I5, CF1 2530 | 122H INDEPENDENT VARIABLES,I5,20H DEPENDENT VARIABLES// CF1 2540 2 111H OBSV. SEQ. 1-11-21 2-12-22 3-CF1 2550 | 313-23 4-14-24 5-15-25 6-16-26 7-17-27 8-18-28 9-19-29 CF1 2560 | 4 10-20-30) CF1 2570 | C INITIALIZE VALUES FOR READING OBSERVATIONS CF1 2580 | 60 NOERR= 0 CF1 2590 KSL2 = 1 CF1 2600 | JOBSV= 0 CF1 2610 | C ADJUST NUMBER OF VARIABLES FOR WEIGHTING FACTOR CF1 2620 | IF (IQ(4)) 65, 65, 70 CF1 2630 | 65 NRDLMT=NOVAR CF1 2640 GO TO 75 CF1 2650 | 70 NRDLMT=NVP1 CF1 2660 | C INITIALIZE INSTRUCTIONS + COUNTERS TO READ OBSERVATIONS CF1 2670 | 75 NOOBSV = 0 CF1 2680 | C SET UP INDEX LIMITS TO WRITE FIRST CARD OF DATA GROUP CF1 2690 IF(NOVAR - 10) 76,66,66 CF1 2700 | 76 MNO=NOVAR CF1 2710 | GO TO 77 CF1 2720 | 66 MNO=10 CF1 2730 | 77 JSEQ=0 CF1 2740 N = 1 CF1 2750 | IF(IQ(19)) 67, 68, 67 CF1 2760 | 67 M= IQ(21) CF1 2770 | MP=IQ(21) CF1 2780 | GO TO 69 CF1 2790 1 12/01/64 PAGE 6 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 68 M=10 CF1 2800 | MP=10 CF1 2810 | 69 KSL1 = 1 CF1 2820 C READ OBSERVATIONS- MP VALUES PER DATA CARD CF1 2830 | 71 IF(IQ(20) -2 ) 83, 144, 83 CF1 2840 | 83 IF(IQ(19) - 1 ) 78, 73, 84 CF1 2850 | 84 CALL REDATA(NOVAR, NOOBSV, READIN) CF1 2860 | C VARIABLE READIN IS UNDEFINED IN THIS MAINLINE,MR7207 CF1 2870 194 GO TO 82 CF1 2880 | 73 READ(KTIN,FMT) IDENT,IOBSV,ISEQ, (DATA(JZ), JZ = N,M ) CF1 2890 196 | GO TO 82 CF1 2900 | 78 READ(KTIN,820) IDENT,IOBSV,ISEQ, (DATA(JZ), JZ = N,M ) CF1 2910 | C DATA CARD INPUT CF1 2920 207 820 FORMAT(1A6, I4, I2, 10F6.3) CF1 2930 | 82 IF(IDENT - IWDA) 85,142,85 CF1 2940 | C CHECK OBSERVATION NO. SEQUENCE IF DESIRED CF1 2950 | 85 IF( IQ(6) ) 125, 125, 100 CF1 2960 | C CHECK IF FIRST CARD OF OBSERVATION HAS OBSERVATION NO. CF1 2970 C LARGER THAN PREVIOUS OBSERVATION NO. CF1 2980 | 100 GO TO (105, 115), KSL1 CF1 2990 | 105 KSL1 = 2 CF1 3000 | IF( IOBSV - JOBSV) 110, 110, 125 CF1 3010 | C COUNT ERROR CF1 3020 110 NOERR= NOERR +1 CF1 3030 | C WRITE ERROR MESSAGE CF1 3040 | WRITE(KTOU,840) IOBSV,JOBSV CF1 3050 | C ERROR MSG IF OBSV. NO. NOT IN ASCENDING SEQUENCE CF1 3060 226 | 840 FORMAT(11H0OBSV. NO. I4,41H IS NOT IN ASCENDING SEQUENCE W/OBSV. NCF1 3070 1O I4) CF1 3080 | GO TO 125 CF1 3090 | C CHECK IF CARDS WITHIN OBSERVATION HAVE SAME OBSERVATION CF1 3100 | C NO. CF1 3110 | 115 IF(IOBSV-JOBSV)120,125,120 CF1 3120 120 NOERR=NOERR + 1 CF1 3130 | WRITE(KTOU,835) JOBSV CF1 3140 | C ERROR MSG IF DATA CARDS W/IN OBSV. HAVE INCORRECT OBSV. NO. CF1 3150 231 | 835 FORMAT(23H0DATA CARDS W/IN OBSV. I4, 25H HAVE INCORRECT OBSV. NO.)CF1 3160 | C CHECK IF SEQUENCE NO. WITHIN OBSERVATION IS ONE GREATER THAN CF1 3170 C PREVIOUS CF1 3180 | 125 IF(IQ(7)) 140, 140, 130 CF1 3190 | 130 IF(ISEQ-JSEQ-1)135,140,135 CF1 3200 | 135 NOERR=NOERR +1 CF1 3210 | C WRITE ERROR MESSAGE CF1 3220 WRITE(KTOU,845) IOBSV CF1 3230 | C ERROR MSG IF SEQ. NO. OF OBSV. NOT 1 GREATER THAN PREVIOUS CF1 3240 237 | 845 FORMAT(25H0SEQ. NO. W/IN OBSV. NO. I4,10H INCORRECT) CF1 3250 | C INCREMENT SEQUENCE NO. IF PRESENT ONE IN ERROR CF1 3260 | ISEQ =ISEQ + 1 CF1 3270 C MAKE CURRENT SEQUENCE NO. = PREVIOUS SEQUENCE NO. CF1 3280 | C MAKE CURRENT OBSV. NO = PREVIOUS OBSV. NO CF1 3290 | 140 JSEQ = ISEQ CF1 3300 | JOBSV= IOBSV CF1 3310 | IF(M - NOVAR)141,142,142 CF1 3320 C INCREMENT COUNTERS FOR VARIABLES ON NEXT DATA CARD CF1 3330 | 141 N=N+MP CF1 3340 | M=M+MP CF1 3350 | 1 12/01/64 PAGE 7 | CF1 - EFN SOURCE STATEMENT - IFN(S) - 0 GO TO 83 CF1 3360 | 142 IF(IQ(20)) 145, 145, 143 CF1 3370 | 143 WRITE(KTBIN5) IDENT,IOBSV,ISEQ, (DATA(JZ), JZ = 1, NOVAR ) CF1 3380 247 | GO TO 145 CF1 3390 | 144 READ (KTBIN5) IDENT,IOBSV,ISEQ, (DATA(JZ), JZ = 1, NOVAR ) CF1 3400 C CHECK FOR CARD WITH WORD END IN FIRST 6 COLUMNS FOR LAST OBSVCF1 3410 255 | 145 IF(IDENT - IWDA)146,250,146 CF1 3420 | C LIST DATA IF DESIRED CF1 3430 | 146 IF(IQ(10)) 149, 147, 149 CF1 3440 | C INCLUDE OBSERVATION NO. IF FIRST CARD OF OBSERVATION CF1 3450 147 ISEQ=1 CF1 3460 | WRITE(KTOU,825) IOBSV,ISEQ, (DATA(JZ), JZ = 1,MNO ) CF1 3470 | C DATA CARD LISTING-CARD 1 OF THE OBSERVATION CF1 3480 270 | 825 FORMAT(1H I4, I6, 10F10.3) CF1 3490 | KKK=10 CF1 3500 170 IF(NOVAR-KKK)149,149,148 CF1 3510 | 148 KONE=KKK+1 CF1 3520 | IF(NOVAR-(KKK+10))166,166,165 CF1 3530 | 165 KTWO=KKK+10 CF1 3540 | GO TO 167 CF1 3550 166 KTWO=NOVAR CF1 3560 | 167 ISEQ=ISEQ+1 CF1 3570 | WRITE(KTOU,830) ISEQ, (DATA(JZ), JZ = KONE, KTWO ) CF1 3580 | C DATA CARD LISTING-CARDS OTHER THAN 1ST OF THE OBSERVATION CF1 3590 286 | 830 FORMAT(5X, I6, 10F10.3) CF1 3600 KKK=KKK+10 CF1 3610 | GO TO 170 CF1 3620 | C END OF LOOP WHICH READS ONE OBSERVATION CF1 3630 | C INCREMENT NO. OF OBSERVATIONS COUNTER CF1 3640 | 149 NOOBSV = NOOBSV + 1 CF1 3650 C AFTER READING + CHECKING ALL VARIABLES, CHECK IF ANY SEQUEN- CF1 3660 | C CING ERRORS OCCURRED. IF NONE OCCUR, CONTINUE PROCESSING THE CF1 3670 | C OBSERVATION. IF THE NO. OF ERRORS IS LESS THAN FIVE, CONTINUECF1 3680 | C READING DATA + CHECKING THE SEQUENCE W/O PROCESSING. WHEN THECF1 3690 | C NO. OF ERRORS REACHES FIVE, CONTINUE READING DATA BUT DISCON-CF1 3700 C TINUE CHECKING SEQUENCING + PROCESSING. PRINTER MESSAGES CF1 3710 | C NOTIFY USER OF THESE ALTERNATIVES. CF1 3720 | IF(NOERR)150, 185, 150 CF1 3730 | 150 IF(NOERR-5)247,155, 155 CF1 3740 | 155 GO TO (160,77), KSL2 CF1 3750 160 KSL2 = 2 CF1 3760 | INSQOB = 0 CF1 3770 | INSQCD = 0 CF1 3780 | INTRFM = 0 CF1 3790 | WRITE (KTOU,855) PR1,PR2,PR3 CF1 3800 C ERROR MSG IF MORE THAN 5 ERRORS FOUND IN SEQ CHECKING DATA CF1 3810 304 | 855 FORMAT(9H0PROBLEM 3A6,91H HAS 5 ERRORS IN SEQUENCING OF DATA. PROBCF1 3820 | 1LEM SKIPPED AT THIS POINT. NO FURTHER DATA CHECKS.) CF1 3830 | GO TO 77 CF1 3840 | C DO ALGEBRAIC TRANSFORMATIONS ON SCALED DATA CF1 3850 185 IF(IQ(3) ) 230,230,190 CF1 3860 | 190 CALL AE04 CF1 3870 | C CHECK IF ANY SPILLS IN DOING TRANSFORMATIONS CF1 3880 | C COMPUTE RAW SUMS + CROSS PRODUCTS, WEIGHTED OR UNWEIGHTED CF1 3890 309 | 230 IF(IQ(4) ) 235,235,240 CF1 3900 235 DATA(NVP1) = 1.0 CF1 3910 | 1 12/01/64 PAGE 8 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 240 DO 245 J=1, KNOVAR CF1 3920 | Z(J, KNVP1) = Z(J, KNVP1) + DATA (J) * DATA(NVP1) CF1 3930 DO 245 M=J, KNOVAR CF1 3940 | 245 Z(J,M) = Z(J, M) + DATA(J) * DATA (M) * DATA (NVP1) CF1 3950 | C CHECK FOR SPILLS CF1 3960 | KNVP1 = KNVP1 CF1 3970 | Z(KNVP1, KNVP1) = Z(KNVP1, KNVP1) + DATA(NVP1) CF1 3980 C WRITE OBSERVATION ON BINARY TAPE -FOR LATER USE- CF1 3990 | WRITE(KTBIN3) IDENT,IOBSV, (DATA(JZ), JZ = 1, KNOVAR ) CF1 4000 337 | 247 JOBSV=IOBSV CF1 4010 | C END OF LOOP FOR READING ALL OBSERVATIONS CF1 4020 | GO TO 77 CF1 4030 250 IF(NOERR) 255, 265, 255 CF1 4040 | 255 WRITE(KTOU,860) PR1,PR2,PR3, NOERR CF1 4050 | C ERROR MSG IF LESS THAN 5 ERRORS FOUND AFTER ALL SEQ CHECKS CF1 4060 | C COMPLETED. CF1 4070 348 | 860 FORMAT(9H0PROBLEM 3A6,10H SKIPPED. I2,93H ERRORS. IF THE NO. OF ERCF1 4080 1RORS IS LESS THAN 5, ALL THE OBSERVATIONS HAVE BEEN SEQUENCE CHECKCF1 4090 | 2ED) CF1 4100 | GO TO 10 CF1 4110 | C FINALIZE THE NO. VARIABLE QUANTITIES CF1 4120 | 265 NOVAR = KNOVAR CF1 4130 NOIND = KNOIND CF1 4140 | NODEP = KNODEP CF1 4150 | NVP1 = KNVP1 CF1 4160 | C CALCULATE DEGREES OF FREEDOM + CHECK CF1 4170 | XNOOBS = Z(NVP1, NVP1) CF1 4180 XNOVAR = NOVAR CF1 4190 | XNOIND = NOIND CF1 4200 | DEFR = XNOOBS - XNOIND - BZRO CF1 4210 | IF(DEFR) 266, 266, 267 CF1 4220 | 266 WRITE(KTOU,935) DEFR CF1 4230 C PRINTER + TAPE MSG IF DEGREES OF FREEDOM ZERO OR NEGATIVE CF1 4240 357 | 935 FORMAT(40H0THE CALCULATED DEGREES OF FREEDOM ARE F5.0,39H NO FURTCF1 4250 | 1HER CALCULATIONS. NEXT PROBLEM.) CF1 4260 | GO TO 10 CF1 4270 | C PUNCH XNOBS, DEGREES OF FREEDOM, KNOIND + KNODEP FOR CONFDNCECF1 4280 C INTERVAL CALCULATIONS. CF1 4290 | 267 IF(IQ(9) - 1 ) 269, 268, 269 CF1 4300 | 268 WRITE (KTPCH,1020) XNOOBS, DEFR, KNOIND, KNODEP CF1 4310 | C OCTAL PUNCH OUT FOR CONFIDENCE INTERVAL DATA CF1 4320 361 | 1020 FORMAT(6O12) CF1 4330 C LIST TRANSFORMED DATA IF DESIRED CF1 4340 | 269 REWIND KTBIN3 CF1 4350 362 | IF(IQ(3)) 3032,3032,264 CF1 4360 | 264 IF(IQ(11)) 3032,270,3032 CF1 4370 | C CF1 4380 C TRANSFORMED DATA LISTING CF1 4390 | 270 WRITE(KTOU,863) CF1 4400 366 | 863 FORMAT(21H0DATA TRANSFORMATIONS) CF1 4410 | C LIST TRANSFORMATION CARDS CF1 4420 | WRITE(KTOU,995) CF1 4430 C TITLE LINE FOR TRANSFORMATION CARD LISTING CF1 4440 367 | 995 FORMAT(1H ,7X,8HPOSITION,3X,4HCODE,3X,9HOPERATION,17X, CF1 4450 | 1 8HCONSTANT,2X,8HLOCATION,2X,4HOMIT,2X,8HVARIABLE) CF1 4460 | MCOUNT = 0 CF1 4470 | 1 12/01/64 PAGE 9 CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 DO 1010 I=1, JNOVAR CF1 4480 | C CF1 4490 | C --- SIMPLIFY TRANSFORMATION SYMBOLS FOR PRINTOUT CF1 4500 | C CF1 4510 272 IOPCDE=ITRFM(I)+1 CF1 4520 | C CF1 4530 | C --- PUT CODE SYMBOL IN F ARRAY VIA E ARRAY CF1 4540 | C CF1 4550 | NAAA = ITRFM(I) CF1 4560 IF(ITRFM(I).EQ.0) GO TO 1001 CF1 4570 | FRRAY= ERRAY(NAAA) CF1 4580 | GO TO 1002 CF1 4590 | 1001 FRRAY= ERRAY(51) CF1 4600 | 1002 NAAB = ITLOC(I) CF1 4610 C CF1 4620 | C --- PUT LOCATION SYMBOL IN G ARRAY VIA E ARRAY CF1 4630 | C CF1 4640 | IF(ITLOC(I).EQ.0) GO TO 1003 CF1 4650 | GRRAY=ERRAY(NAAB) CF1 4660 GO TO 1004 CF1 4670 | 1003 GRRAY= ERRAY(I) CF1 4680 | C CF1 4690 | C --- PUT VARIABLE IN H ARRAY VIA E ARRAY CF1 4700 | C CF1 4710 1004 IF(IOMIT(I).EQ.1) GO TO 1005 CF1 4720 | MCOUNT = MCOUNT + 1 CF1 4730 | HRRAY = ERRAY(MCOUNT) CF1 4740 | GO TO 1006 CF1 4750 | 1005 HRRAY = ERRAY(51) CF1 4760 1006 IF(ITRFM(I).GE.9.AND.ITRFM(I).LE.12) GO TO 1008 CF1 4770 | C CF1 4780 | C --- PRINTOUT OPERATIONS THAT USE CONSTANTS CF1 4790 | C CF1 4800 | WRITE(KTOU,1007)I,FRRAY,(ARRAY(IOPCDE,JZ),JZ=1,6),AVATR(I), CF1 4810 1 GRRAY, IOMIT(I), HRRAY CF1 4820 | C OUTPUT FORMAT OF TRANSFORMATION CARD CF1 4830 407 | 1007 FORMAT(1H ,10X,I2,7X,A2,4X,6A4,3X,A5,7X,A2,6X,I1,7X,A2) CF1 4840 | GO TO 1010 CF1 4850 | C CF1 4860 C --- PRINTOUT OPERATIONS THAT UTILIZE POSITIONS CF1 4870 | C CF1 4880 | C CF1 4890 | C --- PUT POSITION CHANGES IN D ARRAY CF1 4900 | C CF1 4910 1008 JAA = ITRFM(I) - 8 CF1 4920 | MA = INT(0.01*VARTR(I)) CF1 4930 | NA = INT(AMOD(VARTR(I),100.) ) CF1 4940 | WRITE(KTOU,1009) I,FRRAY,(ARRAY(IOPCDE,JZ),JZ=1,6), MA,DRRAY(JAA),CF1 4950 | 1 NA, GRRAY, IOMIT(I), HRRAY CF1 4960 C OUTPUT FORMAT OF TRANSFORMATION CARD CF1 4970 421 | 1009 FORMAT(1H ,10X,I2,7X,A2, 6A4, I2,A1,I2,14X,A2,6X,I1,7X,A2) CF1 4980 | 1010 CONTINUE CF1 4990 | C CF1 5000 | IF(KNOIND+KNODEP-30)299,299,296 CF1 5010 296 WRITE(KTOU,297) CF1 5020 432 | 297 FORMAT(1H0///96H0AFTER PERFORMING TRANSFORMATIONS, THE NUMBER OF CF1 5030 | 1 12/01/64 PAGE 10 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 1INDEPENDENT PLUS DEPENDENT VARIABLES EXCEEDS 30/17H0PROBLEM SKIPPECF1 5040 2D.) CF1 5050 | GO TO 10 CF1 5060 | 299 IF(KNOIND)282,282,281 CF1 5070 | 281 IF(KNODEP)283,283,284 CF1 5080 | 282 WRITE(KTOU,286) CF1 5090 438 286 FORMAT(1H0///74H0AFTER PERFORMING THE TRANSFORMATIONS, THERE ARE NCF1 5100 | 1O INDEPENDENT VARIABLES./17H0PROBLEM SKIPPED.) CF1 5110 | GO TO 10 CF1 5120 | 283 WRITE(KTOU,287) CF1 5130 440 | 287 FORMAT(1H0///72H0AFTER PERFORMING THE TRANSFORMATIONS, THERE ARE NCF1 5140 1O DEPENDENT VARIABLES./17H0PROBLEM SKIPPED.) CF1 5150 | GO TO 10 CF1 5160 | C CF1 5170 | 284 WRITE(KTOU,864) KNOIND, KNODEP CF1 5180 442 | 864 FORMAT(27H0DATA AFTER TRANSFORMATIONS, 6X, CF1 5190 1 24HTHE FITTED EQUATION HAS ,I2, CF1 5200 | 224H INDEPENDENT VARIABLES, I2,20H DEPENDENT VARIABLES / ) CF1 5210 | WRITE (KTOU,862) CF1 5220 443 | 862 FORMAT(120H OBSV. 1-10-19 2-11-20 3-12-21 4-13-22 CF1 5230 | 1 5-14-23 6-15-24 7-16-25 8-17-26 9-18-27 )CF1 5240 C AMAX(50) AND AMIN(50) ARE LOCATIONS FOR MAX + MIN OF VARIABLES CF1 5250 | 3032 DO 3000 I=1,NOVAR CF1 5260 | AMIN(I)=1.0E32 CF1 5270 | AMAX(I)=-AMIN(I) CF1 5280 | 3000 CONTINUE CF1 5290 DO 275 I=1, NOOBSV CF1 5300 | READ (KTBIN3) IDENT,IOBSV, (DATA(JZ), JZ = 1, NOVAR ) CF1 5310 456 | DO 3008 JZ=1,NOVAR CF1 5320 | IF(DATA(JZ)-AMIN(JZ)) 3002,3004,3004 CF1 5330 | 3002 AMIN(JZ)=DATA(JZ) CF1 5340 3004 IF(DATA(JZ)-AMAX(JZ)) 3008,3008,3006 CF1 5350 | 3006 AMAX(JZ)=DATA(JZ) CF1 5360 | 3008 CONTINUE CF1 5370 | IF(IQ(3)) 275,275,3034 CF1 5380 | 3034 IF(IQ(11)) 275,3036,275 CF1 5390 3036 WRITE(KTOU,865)IOBSV,(DATA(JZ),JZ=1,NOVAR) CF1 5400 483 | 865 FORMAT(1X,I4,1X, 1P9E12.5/(6X,1P9E12.5)) CF1 5410 | 275 CONTINUE CF1 5420 | C W1(50) IS LOCATION FOR RANGE OF VARIABLES CF1 5430 | DO 3010 JZ=1,NOVAR CF1 5440 W1(JZ)=AMAX(JZ)-AMIN(JZ) CF1 5450 | 3010 CONTINUE CF1 5460 | C LIST SUMS OF VARIABLES CF1 5470 | 280 IF(IQ(12) ) 290,285,290 CF1 5480 | 285 WRITE(KTOU,870) CF1 5490 C SUM OF VARIABLES LISTING CF1 5500 502 | 870 FORMAT(18H0SUMS OF VARIABLES) CF1 5510 | WRITE(KTOU,875) (Z(I,NVP1), I =1,NOVAR) CF1 5520 503 | 875 FORMAT(6X, 1P9E12.5) CF1 5530 | C LIST RAW SUMS + CROSS PRODUCTS IF DESIRED CF1 5540 290 IF(IQ(13)) 305,295,305 CF1 5550 | 295 WRITE(KTOU,885) CF1 5560 | C HEADING LINE FOR RAW SUMS OF SQUARES + CROSS PRODUCTS CF1 5570 512 | 885 FORMAT(37H0RAW SUMS OF SQUARES + CROSS PRODUCTS/120H 1-1CF1 5580 | 10-19 2-11-20 3-12-21 4-13-22 5-14-23 6-15-24 CF1 5590 1 12/01/64 PAGE 11 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 2 7-16-25 8-17-26 9-18-27 ) CF1 5600 | DO 300J=1, NOVAR CF1 5610 | WRITE(KTOU,890) (Z(JZ,J), JZ = 1,J ) CF1 5620 515 890 FORMAT(6X, 1P9E12.5) CF1 5630 | 300 CONTINUE CF1 5640 | C PUNCH RAW SUMS + CROSS PRODUCTS, TRIANGULAR, IF CONFIDENCE CF1 5650 | C INTERVAL DATA NEEDED + BZERO = 0. CF1 5660 | C SAVE Z(1,1) IF JOINT CONFIDENCE INTERVAL PUNCH OUT RQSTD CF1 5670 305 B(30) = Z(1,1) CF1 5680 | IF(IQ(5)) 307,307,306 CF1 5690 | 306 IF(IQ(9) - 1 ) 325,308,325 CF1 5700 | 308 WRITE(KTPCH,1020)((Z(JZA,JZB),JZB=1,KNOIND),JZA=1,KNOIND) CF1 5710 527 | GO TO 325 CF1 5720 C COMPUTE RESIDUAL SUMS CF1 5730 | 307 DO 310 I=1, NOVAR CF1 5740 | DO 310 J=I, NOVAR CF1 5750 | 310 Z(I,J) = Z(I,J)-(Z(I,NVP1)*Z(J,NVP1)/Z(NVP1,NVP1)) CF1 5760 | C CHECK FOR SPILLS CF1 5770 C PUNCH RESIDUALS IF CONFIDENCE INTERVAL DATA NEEDED + BZERO CF1 5780 | C NOT EQUAL TO 0. CF1 5790 | IF(IQ(9) - 1 ) 314, 312, 314 CF1 5800 | 312 WRITE(KTPCH,1020)((Z(JZA,JZB),JZB=1,KNOIND),JZA=1,KNOIND) CF1 5810 | C LIST RESIDUALS IF DESIRED CF1 5820 555 314 IF(IQ(14) ) 325,315,325 CF1 5830 | 315 WRITE(KTOU,895) CF1 5840 | C RESIDUAL SUMS CF1 5850 567 | 895 FORMAT(42H0RESIDUAL SUMS OF SQUARES + CROSS PRODUCTS) CF1 5860 | DO 320 J=1, NOVAR CF1 5870 320 WRITE(KTOU,890) (Z(JZ,J), JZ = 1,J ) CF1 5880 | C CHECK DIAGONALS OF RESIDUALS CF1 5890 571 | 325 NOVM1 =NOVAR-1 CF1 5900 | DO 350 I=1, NOVAR CF1 5910 | IF(Z(I,I))330,335,340 CF1 5920 330 WRITE(KTOU,900) I CF1 5930 585 | 900 FORMAT(34H0RESIDUAL DIFFERENCE FOR VARIABLE I2,30H IS NEGATIVE. PRCF1 5940 | 1OBLEM SKIPPED.) CF1 5950 | GO TO 10 CF1 5960 | 335 WRITE(KTOU,905) I CF1 5970 587 905 FORMAT(34H0RESIDUAL DIFFERENCE FOR VARIABLE I2,31H IS ZERO. VARIABCF1 5980 | 1LE IS CONSTANT.) CF1 5990 | C PROGRAM CHANGE MADE HERE ON 2-13-62. PROGRAM NOW GIVES UP IF CF1 6000 | C Z(I,I) = 0 CF1 6010 | WRITE (KTOU,906) CF1 6020 588 906 FORMAT(16H PROBLEM SKIPPED) CF1 6030 | GO TO 10 CF1 6040 | 340 SIGMA(I)= SQRT(Z(I,I)) CF1 6050 592 | 350 Z(I,I) = 1. CF1 6060 | C SAVE DEPENDENT VARIABLE PORTION ABOVE DIAGONAL OF THE RESIDU-CF1 6070 C AL MATRIX FOR STATISTICAL COMPUTATIONS TO COME. CF1 6080 | L=NOIND+1 CF1 6090 | DO 353J=L, NOVAR CF1 6100 | WRITE(KTBIN2) (Z(JZ,J), JZ = 1,NOIND ) CF1 6110 602 | 353 CONTINUE CF1 6120 C CALCULATE MEAN + STD. DEV. OF VARIABLES CF1 6130 | 355 RTDEFR = SQRT(Z(NVP1,NVP1)-1.0) CF1 6140 612 | DO 360 I=1, NOVAR CF1 6150 | 1 12/01/64 PAGE 12 | CF1 - EFN SOURCE STATEMENT - IFN(S) - 0 AVG(I)= Z(I,NVP1)/Z(NVP1,NVP1) CF1 6160 | 360 STDDEV(I)=SIGMA(I)/RTDEFR CF1 6170 | C LIST CALCULATED MEANS + STD. DEVS. IF DESIRED CF1 6180 | IF(IQ(15) ) 365,362,365 CF1 6190 | 362 WRITE(KTOU,910) CF1 6200 625 910 FORMAT(19H0MEANS OF VARIABLES) CF1 6210 | WRITE(KTOU,890) (AVG(JZ), JZ = 1, NOVAR ) CF1 6220 626 | WRITE(KTOU,915) CF1 6230 633 | 915 FORMAT(33H0STANDARD DEVIATIONS OF VARIABLES) CF1 6240 | WRITE(KTOU,890) (STDDEV(JZ), JZ = 1,NOVAR ) CF1 6250 C CALCULATE SIMPLE CORRELATION COEFFICIENTS, R(I,J) CF1 6260 634 | 365 GO TO(3305,369), KCP CF1 6270 | C CP PUNCH CF1 6280 | 3305 WRITE(KTPCH,3307) PR1,PR2,PR3,NOVAR,NOOBSV CF1 6290 642 | 3307 FORMAT( 6HCPCPCP, 5O12 ) CF1 6300 WRITE(KTPCH,1020) (STDDEV(JZ), JZ = 1,NOVAR) CF1 6310 643 | 369 DO 370 I = 1, NOVM1 CF1 6320 | IPL1 =I+1 CF1 6330 | DO 370 J=IPL1, NOVAR CF1 6340 | Z(I,J)=Z(I,J)/(SIGMA(I)*SIGMA(J)) CF1 6350 C SQUARE THE TRIANGULAR R(I,J) MATRIX CF1 6360 | 370 Z(J,I)= Z(I,J) CF1 6370 | C CHECK FOR SPILLS CF1 6380 | C PUNCH R(I,J) VALUES, TRIANGULAR, IF CONFIDENCE INTERVAL DATA CF1 6390 | C NEEDED. CF1 6400 IF(IQ(9) -1) 373, 372,373 CF1 6410 | 372 WRITE(KTPCH,1020)((Z(JZA,JZB),JZB=1,KNOIND),JZA=1,KNOIND) CF1 6420 | C LIST THE SIMPLE CORRELATION COEFFICIENTS IF DESIRED CF1 6430 668 | 373 IF(IQ(16) ) 385,375,385 CF1 6440 | 375 WRITE(KTOU,920) CF1 6450 680 920 FORMAT(1H0,45HSIMPLE CORRELATION COEFFICIENTS, R(I,I PRIME)) CF1 6460 | DO 380 J=1, NOVAR CF1 6470 | 380 WRITE(KTOU,3038) (Z(JZ,J),JZ=1,J) CF1 6480 684 | 3038 FORMAT(1H ,5X,9F12.3) CF1 6490 | C INVERT INDEPENDENT VARIABLE PORTION OF THE SIMPLE CORRELATIONCF1 6500 C COEFFICIENT MATRIX, R(I,J), I=J=NO. IND. VARIABLES CF1 6510 | 385 GO TO(3315,386), KCP CF1 6520 | C CP PUNCH CF1 6530 | 3315 DO 3317 J = 1, NOVAR CF1 6540 | 3317 WRITE(KTPCH,1020) (Z(J,JZ), JZ = J,NOVAR) CF1 6550 696 386 CALL AE02(NOIND) CF1 6560 | C CHECK FOR SPILLS CF1 6570 | C AFTER INVERTING, THE INVERSE IS C(I,J) CF1 6580 | C LIST INVERSE IF DESIRED CF1 6590 704 | IF(IQ(17) ) 400,390,400 CF1 6600 390 WRITE(KTOU,925) CF1 6610 | C INVERSE LISTING CF1 6620 707 | 925 FORMAT(1H0,21HINVERSE, C(I,I PRIME)) CF1 6630 | DO 395 I=1, NOIND CF1 6640 | 395 WRITE(KTOU,3040) (C(I,JZ),JZ=1,NOIND) CF1 6650 711 3040 FORMAT(1H ,5X,9F12.3) CF1 6660 | C BEGIN CALCULATIONS FOR EACH DEPENDENT VARIABLE CF1 6670 | 400 REWIND KTBIN2 CF1 6680 | C PUNCH REMAINDER OF NON-REPITITIVE DATA FOR CONFIDENCE INTER- CF1 6690 | C VAL CALCULATIONS. CF1 6700 718 IF(IQ(9) -1 ) 401,510,401 CF1 6710 | 1 12/01/64 PAGE 13 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0C PUNCH MEANS OF THE INDEPENDENT VARIABLES CF1 6720 | 510 WRITE(KTPCH,1020) (AVG(JZ), JZ = 1,KNOIND ) CF1 6730 C PUNCH INVERSE OF R(I,J) MATRIX CF1 6740 720 | WRITE(KTPCH,1020)((C(JZA,JZB),JZB=1,KNOIND),JZA=1,KNOIND) CF1 6750 | C REARRANGE TRFMTN INFO IF BACK TRANSFORMING DEPENDENT VARIABLECF1 6760 727 | 401 IF(IQ(3) ) 410,410,402 CF1 6770 | 402 IF(IQ(8) ) 410,410,403 CF1 6780 403 K=1 CF1 6790 | DO 405I=1, JNOVAR CF1 6800 | IF(IOMIT(I)) 404, 404, 405 CF1 6810 | 404 ITRFM(K) = ITRFM(I) CF1 6820 | VARTR(K) = VARTR(I) CF1 6830 K =K+1 CF1 6840 | 405 CONTINUE CF1 6850 | 410 K = L CF1 6870 | C SET LOOP TO CALCULATE STATISTICS FOR EACH DEPENDENT VARIABLE CF1 6860 | 411 RSQD = 0.0 CF1 6880 BZERO =0.0 CF1 6890 | SUMB =0.0 CF1 6900 | REWIND KTBIN3 CF1 6910 758 | YYMIN = 1.0E32 CF1 6920 | YCMIN =1.0E32 CF1 6930 YYMAX = -YYMIN CF1 6940 | YCMAX = -YCMIN CF1 6950 | C INITIALIZE BACK-TRANSFORMATION ROUTINE CF1 6960 | IF(IQ(8) ) 412,412,413 CF1 6970 | 412 M = 2 CF1 6980 GO TO 414 CF1 6990 | 413 M = 1 CF1 7000 | C WRITE HEADING LINE FOR STATISTICS CF1 7010 | 414 ICURY = K - NOIND CF1 7020 | WRITE(KTOU,5102)PR1,PR2,PR3, ICURY CF1 7030 C HEADING LINE FOR STATISTICAL CALCULATIONS CF1 7040 764 | 5102 FORMAT(1H1,3A6,10H DEP VAR ,I2) CF1 7050 | WRITE(KTOU,945) CF1 7060 765 | 945 FORMAT(1H0,10HIND.VAR(I),2X,9HCOEF.B(I),6X,10HS.E. COEF.,3X,7HT-VACF1 7070 | 1LUE,4X,8HINV. DG.,6X,8HMIN X(I),7X,8HMAX X(I),5X,10HRANGE X(I),2X,CF1 7080 212HREL.INF.X(I)) CF1 7090 | C READ RESIDUAL VALUES FOR DEPENDENT VARIABLES CF1 7100 | READ(KTBIN2) (YRESID(JZ), JZ =1,NOIND ) CF1 7110 | C CALCULATE BETA VALUES CF1 7120 766 | DO 420 I=1, NOIND CF1 7130 BETA(I) = 0.0 CF1 7140 | DO 415 J=1, NOIND CF1 7150 | 415 BETA(I)=BETA(I)+C(I,J)*Z(J,K) CF1 7160 | C CALCULATE B(I), MULT. REGRN. COEFFICIENT CF1 7170 | B(I) =BETA(I)*SIGMA(K)/SIGMA(I) CF1 7180 BZERO = BZERO +(B(I)*AVG(I)) CF1 7190 | C CALCULATE R ZERO SQ, OVER-ALL MULT. CORREL. COEFF. CF1 7200 | 420 RSQD=RSQD + BETA(I)*Z(I,K) CF1 7210 | C CALC BZERO CF1 7220 | BZERO =AVG(K)-BZERO CF1 7230 C CALCULATE CONSTANT TERM FOR SIGMA BETA CALC CF1 7240 | FACTOR = SQRT((1.0 - RSQD)/DEFR) CF1 7250 | C CALCULATE STATISTICS + LIST FOR EACH IND. VARIABLE CF1 7260 797 | WRITE(KTOU,3012) BZERO CF1 7270 798 | 1 12/01/64 PAGE 14 CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 3012 FORMAT(1H ,5X,1H0,4X,1PE12.5) CF1 7280 | DO 425 I=1, NOIND CF1 7290 | C CALCULATE STANDARD DEV. OF BETA CF1 7300 | SGBETA = SQRT(C(I,I)) * FACTOR CF1 7310 C PERFORM T TEST ON BETA CF1 7320 803 | TBETA = ABS(BETA(I)/SGBETA) CF1 7330 | C CALCULATE PARTIAL CORRELATION COEFFICIENT CF1 7340 | PCORK =TBETA**2/(TBETA**2+DEFR) CF1 7350 | C CALCULATE SIGMA B CF1 7360 SIGB =SGBETA*SIGMA(K)/SIGMA(I) CF1 7370 | C CALCULATE B(I)*X(I) FOR BZERO CALCULATIONS CF1 7380 | BWY=B(I)*W1(I)/W1(K) CF1 7390 | C WRITE THESE VALUES ON TP 6 CF1 7400 | WRITE(KTOU,950) I,B(I),SIGB,TBETA,C(I,I),AMIN(I),AMAX(I),W1(I), CF1 7410 1BWY CF1 7420 810 | 950 FORMAT(1H ,4X,I2,4X,1PE12.5,5X,1PE9.2,4X,0PF5.1,5X,1PE9.2,1P2E15.3CF1 7430 | 1,4X,1PE10.3,4X,0PF6.2) CF1 7440 | 425 SUMB=SUMB+(B(I)*YRESID(I)) CF1 7450 | C CHECK FOR SPILLS CF1 7460 C PERFORM F RATIO TEST ON OVERALL REGRESSION CF1 7470 | FRATIO =RSQD*DEFR/(XNOIND*(1.0-RSQD)) CF1 7480 | C CALCULATE THE SIGMA OF THE REGRN. EQN - UNBIASED ESTIMATE CF1 7490 | UNBEST = SQRT((SIGMA(K)**2 - SUMB)/DEFR) CF1 7500 | C CALCULATE THE BIASED ESIMATE CF1 7510 823 BIEST = UNBEST* SQRT(DEFR/XNOOBS) CF1 7520 | C CALCULATE S SQD CF1 7530 824 | SSQ =UNBEST**2 CF1 7540 | C CALCULATE SSE CF1 7550 | SSE =SIGMA(K)**2-SUMB CF1 7560 SSYT=SIGMA(K)**2 CF1 7570 | C LIST THE STATISTICS COMPUTED CF1 7580 | WRITE(KTOU,3014) NOOBSV CF1 7590 829 | 3014 FORMAT(1H0,19HNO. OF OBSERVATIONS,11X,I5) CF1 7600 | WRITE(KTOU,3016) NOIND CF1 7610 830 3016 FORMAT(1H ,21HNO. OF IND. VARIABLES,10X,I4) CF1 7620 | NDEFR=DEFR+0.01 CF1 7630 | WRITE(KTOU,3018) NDEFR CF1 7640 831 | 3018 FORMAT(1H , 27HRESIDUAL DEGREES OF FREEDOM,4X,I4) CF1 7650 | WRITE(KTOU,3020) FRATIO CF1 7660 832 3020 FORMAT(1H ,7HF-VALUE,24X,F6.1) CF1 7670 | WRITE(KTOU,3022) UNBEST CF1 7680 833 | 3022 FORMAT(1H ,25HRESIDUAL ROOT MEAN SQUARE,3X,F16.8) CF1 7690 | WRITE(KTOU,3024) SSQ CF1 7700 834 | 3024 FORMAT(1H ,20HRESIDUAL MEAN SQUARE,8X,F16.8) CF1 7710 WRITE(KTOU,3026) SSE CF1 7720 835 | 3026 FORMAT(1H ,23HRESIDUAL SUM OF SQUARES,5X,F16.8) CF1 7730 | WRITE(KTOU,3028) SSYT CF1 7740 836 | 3028 FORMAT(1H ,20HTOTAL SUM OF SQUARES,8X,F16.8) CF1 7750 | WRITE(KTOU,3030) RSQD CF1 7760 837 3030 FORMAT(1H ,27HMULT. CORREL. COEF. SQUARED,8X,F5.4) CF1 7770 | C CHECK FOR SPILLS CF1 7780 | C CALCULATE FITTED Y AND LIST WITH RESIDUALS CF1 7790 | C INITIALIZE B ZERO WITH ZERO IF YCALC HAS BZERO = 0 CF1 7800 | IF(IQ(5)) 438,438,427 CF1 7810 427 BZERO = 0.0 CF1 7820 | 438 DO 475 IB = 1, NOOBSV CF1 7830 | 1 12/01/64 PAGE 15 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 READ(KTBIN3) IDZA(IB),IDZB(IB), (VAR(JZ), JZ = 1,NOVAR ) CF1 7840 843 YCALC = BZERO CF1 7850 | DO 440 J=1, NOIND CF1 7860 | 440 YCALC = YCALC + B(J) * VAR(J) CF1 7870 | C DO BACK TRANSFORMATION--DECISION MADE BY XAE12F CF1 7880 | GO TO(445, 470), M CF1 7890 445 CALL YBACK(K,YCALC) CF1 7900 858 | CALL YBACK(K,VAR(K) ) CF1 7910 860 | 470 DIFF = VAR(K)- YCALC CF1 7920 | DELTA(IB) = DIFF CF1 7930 | YYIN(IB) = VAR(K) CF1 7940 YYCC(IB) = YCALC CF1 7950 | CALL MINMAX(YYMIN, YYMAX, VAR(K) ) CF1 7960 868 | CALL MINMAX(YCMIN, YCMAX, YCALC ) CF1 7970 871 | 475 CONTINUE CF1 7980 | C CHECK FOR SPILLS IN COMPUTING Y CALC + DIFF CF1 7990 LSORT(1) = 1 CF1 8000 | DO 1306 ISA = 2, NOOBSV CF1 8010 | TEMP= DELTA(ISA) CF1 8020 | ISAX= ISA - 1 CF1 8030 | DO 1301 ISB = 1, ISAX CF1 8040 IBB = LSORT(ISB) CF1 8050 | IF(TEMP - DELTA(IBB) ) 1302,1302,1301 CF1 8060 | 1301 CONTINUE CF1 8070 | LSORT(ISA) = ISA CF1 8080 | GO TO 1306 CF1 8090 1302 KSA = ISA CF1 8100 | 1303 KSA = KSA - 1 CF1 8110 | LSORT(KSA + 1) = LSORT(KSA) CF1 8120 | IF(KSA - ISB ) 1305, 1305, 1303 CF1 8130 | 1305 LSORT(ISB) = ISA CF1 8140 1306 CONTINUE CF1 8150 | WRITE(KTOU,1015) CF1 8160 | C HEADING IF LISTING ORDERED RESIDUALS WITH PLOTS CF1 8170 900 | 1015 FORMAT(1H0,6HIDENT.,3X,8HOBS. NO.,4X,6HOBS. Y,5X,8HFITTED Y,5X,12HCF1 8180 | 1RESIDUAL ***,2X,8HOBS. NO.,4X,6HOBS. Y,5X,8HFITTED Y,2X,14HORDEREDCF1 8190 2 RESID.,5H SEQ) CF1 8200 | DO 1220 IC = 1, NOOBSV CF1 8210 | JC = LSORT(IC) CF1 8220 | WRITE(KTOU,1219) IDZA(IC),IDZB(IC),YYIN(IC),YYCC(IC),DELTA(IC), CF1 8230 | 1IDZB(JC),YYIN(JC),YYCC(JC),DELTA(JC),IC CF1 8240 C CALCULATED Y + DIFFERENCE CF1 8250 905 | 1219 FORMAT(1H ,A6,4X,I4,3F13.3,7X,I4,3F13.3,3X,I4) CF1 8260 | 1220 CONTINUE CF1 8270 | IF(IQ(18) ) 486, 480, 486 CF1 8280 | C DO PLOTTING CF1 8290 480 CALL PITCHA CF1 8300 | C PUNCH REPETITIVE DATA FOR CONFIDENCE INTERVAL CALCULATION CF1 8310 920 | 486 IF(IQ(9)- 1) 485, 482, 483 CF1 8320 | 482 WRITE(KTPCH,1020) SSQ,BZERO, (B(JZ), JZ=1,KNOIND ) CF1 8330 923 | GO TO 485 CF1 8340 C PUNCH CARD FOR ELLIPTIC CONF. INTVL. CF1 8350 | 483 ME = KNVP1 CF1 8360 | WRITE(KTPCH,1025) PR1,NOOBSV,B(30),Z(1,ME),BZERO,B(1),SSQ CF1 8370 932 | 1025 FORMAT(A6, 3X, I3, 5O12) CF1 8380 | 485 K = K + 1 CF1 8390 1 12/01/64 PAGE 16 | CF1 - EFN SOURCE STATEMENT - IFN(S) - | 0 IF( K - NOVAR) 411, 411, 10 CF1 8400 | END CF1 8410 | 1 12/01/64 PAGE 17 | 0| -| $IBFTC CF2 $IBFTC CF2 0000 | 1 12/01/64 PAGE 18 CF2 - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE REDATA(KAAAAA,JAAAAA,XXXXX) CF2 0010 | C CF2 0020 | C SUBROUTINE TO READ DATA. CAN BE REVISED TO SUIT SPECIAL NEEDSCF2 0030 C CF2 0040 | DIMENSION Z(31,31), C(30,30), DATA(51) CF2 0050 | DIMENSION BETA(30), AVG(30), SIGMA(30), STDDEV(30), B(30) CF2 0060 | DIMENSION VAR(30), YRESID(30) CF2 0070 | DIMENSION ITRFM(50), VARTR(50), AVATR(50), ITLOC(50), IOMIT(50) CF2 0080 DIMENSION FMT(80), IQ(40) CF2 0090 | DIMENSION IDZA(1000),IDZB(1000),DELTA(1000) CF2 0100 | DIMENSION YYIN(1000),YYCC(1000),LSORT(1000) CF2 0110 | COMMON YYMIN, YCMIN, YYMAX, YCMAX CF2 0120 | COMMON FMT, IQ, KSL1, KSL2, KSL3, KSL4, KSL5, KSL6 CF2 0130 COMMON IDZA, IDZB, DELTA, YYIN, YYCC, LSORT CF2 0140 | COMMON Z,C, BETA,AVG,SIGMA, B,VAR,YRESID CF2 0150 | COMMON ITRFM,VARTR,AVATR, IOMIT CF2 0160 | COMMON KTIN,KTOU,KTPCH,KTBIN1,KTBIN2,KTBIN3,KTBIN4,KTBIN5 CF2 0170 | COMMON PR1,PR2,PR3 CF2 0180 COMMON BZRO, KNOIND, KNODEP, KNOVAR, KNVP1, JNOIND, JNODEP, JNOVARCF2 0190 | COMMON NOBSP, NOERR, JOBSV, NRDLMT, NOOBSV, MNO, JSEQ, N, M, J, I CF2 0200 | COMMON LNOVAR, NNO, ISEQ, L, KNO, XNOOBS, XNOVAR, XNOIND, DEFR CF2 0210 | COMMON NOVM1, RTDEFR, IPL1, K, RSQD, BZERO, SUMB, ICURY, FACTOR CF2 0220 | COMMON SGBETA, TBETA, PCORK, SIGB, FRATIO, UNBEST, BIEST, SSQ, SSECF2 0230 COMMON IPLOT1, IPLOT2, YCALC, DIFF, IREC, X, IDENT, IOBSV, MP CF2 0240 | COMMON NOVAR, NVP1, IPATH CF2 0250 | EQUIVALENCE (C,DATA), (C(101),STDDEV), (C(201),ITLOC) CF2 0260 | EQUIVALENCE (NOIND,IQ(1)),(NODEP,IQ(2)) CF2 0270 | EQUIVALENCE (KSW45,IQ(22)), (KSW46,IQ(23)) CF2 0280 C 45 = PARAM , 46 = NEW DATA PUNCH CF2 0290 | DIMENSION P(6), XA(10), XC(6),XB(10) CF2 0300 | GO TO (20,2), KSL6 CF2 0310 | 2 KSL6 = 1 CF2 0320 | CALL ATHRUZ(IEND,6HEND ) CF2 0330 5 CALL ATHRUZ(JTR, 6HTRIAL ) CF2 0340 7 | 60 IF(KSW45.NE.1) KSW45 = 2 CF2 0350 | IF(KSW46.NE.1) KSW46 = 2 CF2 0360 | GO TO (62,20), KSW45 CF2 0370 | 62 READ(KTIN,4) JJTR,JIG,KSWA,P CF2 0380 14 4 FORMAT(A6,I3,I1,1X,1PE12.5,2X,0P5F10.5) CF2 0390 | IF(JJTR. EQ. JTR) GO TO 10 CF2 0400 | 6 WRITE(KTOU,8) CF2 0410 22 | 8 FORMAT(24H0TRIAL CARD ERROR. BYE. ) CF2 0420 | CALL EXIT CF2 0430 23 10 WRITE(KTOU,12) JJTR,JIG, P CF2 0440 24 | 12 FORMAT(1H0 A6,I4,1PE20.5,0P5F20.6) CF2 0450 | 14 WRITE(KTOU,16) CF2 0460 26 | 16 FORMAT(1H0 17X,1HT,14X,1HZ,14X,3HTZZ,12X,1HU,14X,1HV,14X,4HTZZZ,/)CF2 0470 | 20 GO TO (22,30), KSW46 CF2 0480 22 READ(KTIN,24) IDENT,IOBSV, XA CF2 0490 28 | 24 FORMAT(A6,I4, 10F5.0 ) CF2 0500 | IF(IDENT .NE. IEND ) GO TO 26 CF2 0510 | 90 WRITE(KTPCH,28) IEND CF2 0520 35 | 99 RETURN CF2 0530 26 XB(1)=XA(8)/(XA(3)*(1.0-XA(4)+(XA(4)+XA(5))*XA(6)/XA(7))) CF2 0540 | XB(2)= ALOG(XA(10)) CF2 0550 38 | 1 12/01/64 PAGE 19 | CF2 - EFN SOURCE STATEMENT - IFN(S) - | 0 WRITE(KTPCH,28) IDENT, IOBSV, XA(9),XB(1),XA(2),XA(1),XB(2) CF2 0560 39 28 FORMAT(A6,I4,1P5E14.7 ) CF2 0570 | GO TO 40 CF2 0580 | 30 READ(KTIN,28) IDENT, IOBSV, XA(9),XB(1),XA(2),XA(1),XB(2) CF2 0590 41 | IF(IDENT.NE.IEND) GO TO 40 CF2 0600 | GO TO 99 CF2 0610 40 XC(4) = XA(2) * EXP(-P(4)/XA(9) ) CF2 0620 49 | XC(5) = XA(1) * EXP(-P(6)/XA(9) ) CF2 0630 50 | XC(1) = P(1) * XB(1) * EXP(-P(2)/XA(9) ) CF2 0640 51 | XC(2) =1.0/(1.0 + P(3)* XC(4) + P(5)* XC(5) ) CF2 0650 | XC(3) = XC(1) * XC(2) * XC(2) CF2 0660 XC(6) = XC(3) * XC(2) CF2 0670 | GO TO(42,50), KSW45 CF2 0680 | 42 WRITE(KTOU, 44) IDENT, IOBSV, XC CF2 0690 53 | 44 FORMAT(1H A6,I6,1P6E15.4 ) CF2 0700 | 50 DATA( 7) = 1.0 / XA(9) CF2 0710 DATA( 8) = 2.0 * XC(4) * XC(2) CF2 0720 | DATA( 9) = 2.0 * XC(5) * XC(2) CF2 0730 | DATA(10) = - DATA(8) * P(3) / XA(9) CF2 0740 | DATA(11) = - DATA(9) * P(5) / XA(9) CF2 0750 | DO 52 IA = 1,5 CF2 0760 52 DATA(IA) = DATA(IA + 6) * XC(3) CF2 0770 | DATA(6) = XC(3) CF2 0780 | DATA(12) = EXP(XB(2)) CF2 0790 64 | DATA(13) = XB(2) CF2 0800 | DATA(14) = XB(2)/XC(3) - 1.0 CF2 0810 GO TO 99 CF2 0820 | END CF2 0830 | 1 12/01/64 PAGE 20 | | 0 -| $IBFTC CF3 $IBFTC CF3 0000 | 1 12/01/64 PAGE 21 | CF3 - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE AE04 CF3 0010 | C CF3 0020 | C SUBROUTINE TO PERFORM TRANSFORMATION OF VARIABLES CF3 0030 | C CF3 0040 | DIMENSION Z(31,31), C(30,30), DATA(51) CF3 0050 DIMENSION BETA(30), AVG(30), SIGMA(30), STDDEV(30), B(30) CF3 0060 | DIMENSION VAR(30), YRESID(30) CF3 0070 | DIMENSION ITRFM(50), VARTR(50), AVATR(50), ITLOC(50), IOMIT(50) CF3 0080 | DIMENSION FMT(80), IQ(40) CF3 0090 | DIMENSION IDZA(1000),IDZB(1000),DELTA(1000) CF3 0100 DIMENSION YYIN(1000),YYCC(1000),LSORT(1000) CF3 0110 | COMMON YYMIN, YCMIN, YYMAX, YCMAX CF3 0120 | COMMON FMT, IQ, KSL1, KSL2, KSL3, KSL4, KSL5, KSL6 CF3 0130 | COMMON IDZA, IDZB, DELTA, YYIN, YYCC, LSORT CF3 0140 | COMMON Z,C, BETA,AVG,SIGMA, B,VAR,YRESID CF3 0150 COMMON ITRFM,VARTR,AVATR, IOMIT CF3 0160 | COMMON KTIN,KTOU,KTPCH,KTBIN1,KTBIN2,KTBIN3,KTBIN4,KTBIN5 CF3 0170 | COMMON PR1,PR2,PR3 CF3 0180 | COMMON BZRO, KNOIND, KNODEP, KNOVAR, KNVP1, JNOIND, JNODEP, JNOVARCF3 0190 | COMMON NOBSP, NOERR, JOBSV, NRDLMT, NOOBSV, MNO, JSEQ, N, M, J, I CF3 0200 COMMON LNOVAR, NNO, ISEQ, L, KNO, XNOOBS, XNOVAR, XNOIND, DEFR CF3 0210 | COMMON NOVM1, RTDEFR, IPL1, K, RSQD, BZERO, SUMB, ICURY, FACTOR CF3 0220 | COMMON SGBETA, TBETA, PCORK, SIGB, FRATIO, UNBEST, BIEST, SSQ, SSECF3 0230 | COMMON IPLOT1, IPLOT2, YCALC, DIFF, IREC, X, IDENT, IOBSV, MP CF3 0240 | COMMON NOVAR, NVP1, IPATH CF3 0250 EQUIVALENCE (C,DATA), (C(101),STDDEV), (C(201),ITLOC) CF3 0260 | EQUIVALENCE (NOIND,IQ(1)),(NODEP,IQ(2)) CF3 0270 | C THE FOLLOWING TRANSFORMATION CODES ARE REPRESENTED CF3 0280 | C BLANK = NO TRANSFORMATION CF3 0290 | C 0 = NO TRANSFORMATION CF3 0300 C 1 = NATURAL LOG (BASE E) CF3 0310 | C 2 = COMMON LOG (BASE 10) CF3 0320 | C 3 = ANTI-LOG OF NATURAL LOG CF3 0330 | C 4 = ANTI-LOG OF COMMON LOG CF3 0340 | C 5 = POWER RAISING CF3 0350 C 6 = PRODUCT OF VARIABLE AND ARGUMENT (CON) CF3 0360 | C 7 = QUOTIENT OF VARIABLE AND ARGUMENT (CON) CF3 0370 | C 8 = ALGEBRAIC SUM OF VARIABLE AND ARGUMENT (CON) CF3 0380 | C 9 = PRODUCT OF TWO VARIABLES CF3 0390 | C 10 = QUOTIENT OF TWO VARIABLES CF3 0400 C 11 = SUM OF TWO VARIABLES CF3 0410 | C 12 = DIFFERENCE OF TWO VARIABLES CF3 0420 | C 13 = SIN OF VAR IN RADIANS CF3 0430 | C 14 = COS OF VAR IN RADIANS CF3 0440 | C 15 = QUOTIENT OF ARGUMENT DIVIDED BY VARIABLE CF3 0450 C 16 = LN(LN(X)) WHERE LN = NATURAL LOG CF3 0460 | C 17 = LOG(LOG(X)) CF3 0470 | C 18 = LOG(X+K) = BASE 10 LOG OF (VARIABLE + .OR. - ARGUMENT) CF3 0480 | C 19 = E**(E**X) = INVERSE OF NUMBER 16 CF3 0490 | C 20 = 10**(10**X)=INVERSE OF NUMBER 17 CF3 0500 DO 65 IA= 1, NOVAR CF3 0510 | VX = DATA(IA) CF3 0520 | VY = VARTR(IA) CF3 0530 | K = ITRFM(IA) CF3 0540 | IF( (K.GT.0) .AND. (K.LT.21) ) GO TO 50 CF3 0550 1 12/01/64 PAGE 22 | CF3 - EFN SOURCE STATEMENT - IFN(S) - | 0 VALUE = VX CF3 0560 | GO TO 60 CF3 0570 | C DETERMINE TRANSFORMATION CF3 0580 50 GO TO(1,2,3,4,5,6,7,8,9, 9, 9, 9,13,14,15,16,17,18,19,20), K CF3 0590 | C NATURAL LOG CF3 0600 | 1 VALUE = ALOG(VX) CF3 0610 15 | GO TO 60 CF3 0620 | C COMMON LOG CF3 0630 2 VALUE = ALOG10(VX) CF3 0640 18 | GO TO 60 CF3 0650 | C ANTI-LOG OF NATURAL LOG CF3 0660 | 3 VALUE = EXP(VX) CF3 0670 21 | GO TO 60 CF3 0680 C ANTI-LOG OF COMMON LOG CF3 0690 | 4 VALUE = EXP(2.302585093 * VX ) CF3 0700 24 | GO TO 60 CF3 0710 | C POWER RAISING CF3 0720 | 5 VALUE = VX ** VY CF3 0730 27 GO TO 60 CF3 0740 | C PRODUCT OF VAR + CON CF3 0750 | 6 VALUE = VX * VY CF3 0760 | GO TO 60 CF3 0770 | C QUOTIENT OF VAR + CON CF3 0780 7 VALUE = VX / VY CF3 0790 | GO TO 60 CF3 0800 | C SUM OF VAR + CON CF3 0810 | 8 VALUE = VX + VY CF3 0820 | GO TO 60 CF3 0830 9 MA = INT( 0.01 * VY) CF3 0840 | NA = INT(AMOD(VY,100.)) CF3 0850 | K = K - 8 CF3 0860 | GO TO (91,10,11,12), K CF3 0870 | C PRODUCT OF TWO VARIABLES CF3 0880 91 VALUE = DATA(MA) * DATA(NA) CF3 0890 | GO TO 60 CF3 0900 | C QUOTIENT OF TWO VARIABLES CF3 0910 | 10 VALUE = DATA(MA) / DATA(NA) CF3 0920 | GO TO 60 CF3 0930 C SUM OF TWO VARIABLES CF3 0940 | 11 VALUE = DATA(MA) + DATA(NA) CF3 0950 | GO TO 60 CF3 0960 | C DIFFERENCE OF TWO VARIABLES CF3 0970 | 12 VALUE = DATA(MA) - DATA(NA) CF3 0980 GO TO 60 CF3 0990 | C SIN OF VAR IN RADIANS CF3 1000 | 13 VALUE = SIN(VX) CF3 1010 56 | GO TO 60 CF3 1020 | C COS OF VAR IN RADIANS CF3 1030 14 VALUE = COS(VX) CF3 1040 59 | GO TO 60 CF3 1050 | 15 VALUE = VY / VX CF3 1060 | GO TO 60 CF3 1070 | 16 VALUE = ALOG(ALOG(VX)) CF3 1080 64 65 GO TO 60 CF3 1090 | 17 VALUE = ALOG10(ALOG10(VX)) CF3 1100 68 69 | GO TO 60 CF3 1110 | 1 12/01/64 PAGE 23 | CF3 - EFN SOURCE STATEMENT - IFN(S) - 0 18 VALUE = ALOG10(VX + VY) CF3 1120 72 | GO TO 60 CF3 1130 | 19 VALUE = EXP(EXP(VX)) CF3 1140 75 76 | GO TO 60 CF3 1150 | 20 VALUE = EXP(2.302585093 * EXP(2.302585093 * VX)) CF3 1160 79 80 GO TO 60 CF3 1170 | C DETERMINE LOCATION TO PLACE VARIABLE CF3 1180 | 60 IF( ITLOC(IA) .GT. 0) GO TO 62 CF3 1190 | DATA(IA) = VALUE CF3 1200 | GO TO 65 CF3 1210 62 NA = ITLOC(IA) CF3 1220 | DATA(NA) = VALUE CF3 1230 | 65 CONTINUE CF3 1240 | C OMIT VARIABLES SPECIFIED, RECALCULATING NET NO. LEFT CF3 1250 | L=NOIND +1 CF3 1260 KNO = 1 CF3 1270 | DO 115 IA = 1, NOIND CF3 1280 | IF( IOMIT(IA) .NE. 0) GO TO 115 CF3 1290 | DATA(KNO) = DATA(IA) CF3 1300 | KNO = KNO + 1 CF3 1310 115 CONTINUE CF3 1320 | KNOIND = KNO -1 CF3 1330 | DO 125 IA = L, NOVAR CF3 1340 | IF( IOMIT(IA) .NE. 0) GO TO 125 CF3 1350 | DATA(KNO) = DATA(IA) CF3 1360 KNO = KNO +1 CF3 1370 | 125 CONTINUE CF3 1380 | KNVP1 = KNO CF3 1390 | KNOVAR= KNO -1 CF3 1400 | KNODEP= KNOVAR -KNOIND CF3 1410 RETURN CF3 1420 | END CF3 1430 | 1 12/01/64 PAGE 24 | | 0 -| $IBFTC CF4 $IBFTC CF4 0000 | 1 12/01/64 PAGE 25 | CF4 - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE AE02(NUT) CF4 0010 | C CF4 0020 | C MATRIX INVERSION-UNIT MATRIX METHOD SUBROUTINE AE020000CF4 0030 | C CF4 0040 | DIMENSION Z(31,31), C(30,30), DATA(51) CF4 0050 DIMENSION BETA(30), AVG(30), SIGMA(30), STDDEV(30), B(30) CF4 0060 | DIMENSION VAR(30), YRESID(30) CF4 0070 | DIMENSION ITRFM(50), VARTR(50), AVATR(50), ITLOC(50), IOMIT(50) CF4 0080 | DIMENSION FMT(80), IQ(40) CF4 0090 | DIMENSION IDZA(1000),IDZB(1000),DELTA(1000) CF4 0100 DIMENSION YYIN(1000),YYCC(1000),LSORT(1000) CF4 0110 | COMMON YYMIN, YCMIN, YYMAX, YCMAX CF4 0120 | COMMON FMT, IQ, KSL1, KSL2, KSL3, KSL4, KSL5, KSL6 CF4 0130 | COMMON IDZA, IDZB, DELTA, YYIN, YYCC, LSORT CF4 0140 | COMMON Z,C, BETA,AVG,SIGMA, B,VAR,YRESID CF4 0150 COMMON ITRFM,VARTR,AVATR, IOMIT CF4 0160 | COMMON KTIN,KTOU,KTPCH,KTBIN1,KTBIN2,KTBIN3,KTBIN4,KTBIN5 CF4 0170 | COMMON PR1,PR2,PR3 CF4 0180 | COMMON BZRO, KNOIND, KNODEP, KNOVAR, KNVP1, JNOIND, JNODEP, JNOVARCF4 0190 | COMMON NOBSP, NOERR, JOBSV, NRDLMT, NOOBSV, MNO, JSEQ, N, M, J, I CF4 0200 COMMON LNOVAR, NNO, ISEQ, L, KNO, XNOOBS, XNOVAR, XNOIND, DEFR CF4 0210 | COMMON NOVM1, RTDEFR, IPL1, K, RSQD, BZERO, SUMB, ICURY, FACTOR CF4 0220 | COMMON SGBETA, TBETA, PCORK, SIGB, FRATIO, UNBEST, BIEST, SSQ, SSECF4 0230 | COMMON IPLOT1, IPLOT2, YCALC, DIFF, IREC, X, IDENT, IOBSV, MP CF4 0240 | COMMON NOVAR, NVP1, IPATH CF4 0250 EQUIVALENCE (C,DATA), (C(101),STDDEV), (C(201),ITLOC) CF4 0260 | EQUIVALENCE (NOIND,IQ(1)),(NODEP,IQ(2)) CF4 0270 | C CONSTRUCT THE UNIT MATRIX CF4 0280 | DO 15 I=1,NUT CF4 0290 | DO 15 J=1,NUT CF4 0300 IF (I-J) 10, 5, 10 CF4 0310 | 5 C(I,J)= 1.0 CF4 0320 | GO TO 15 CF4 0330 | 10 C(I,J)= 0.0 CF4 0340 | 15 CONTINUE CF4 0350 DO 45 I=1,NUT CF4 0360 | C DIVIDE ROW(I) BY PIVOT ELEMENT CF4 0370 | PIVOT = 1.0/Z(I,I) CF4 0380 | DO 20 J=I,NUT CF4 0390 | 20 Z(I,J)=Z(I,J)*PIVOT CF4 0400 DO 25 J=1,NUT CF4 0410 | 25 C(I,J)=C(I,J)*PIVOT CF4 0420 | C ELEMINATE KTH VARIABLE FROM ALL ROWS BUT THE I(TH) CF4 0430 | DO 45 K=1,NUT CF4 0440 | IF(K-I) 30, 45, 30 CF4 0450 30 EMPYER = -1.0*Z(K,I) CF4 0460 | DO 35 J=I,NUT CF4 0470 | 35 Z(K,J) =Z(K,J)+Z(I,J)*EMPYER CF4 0480 | DO 40 J=1,NUT CF4 0490 | 40 C(K,J) =C(K,J)+C(I,J)*EMPYER CF4 0500 45 CONTINUE CF4 0510 | RETURN CF4 0520 | END CF4 0530 | 1 12/01/64 PAGE 26 | 0| -| $IBFTC CF5 $IBFTC CF5 0000 | 1 12/01/64 PAGE 27 | CF5 - EFN SOURCE STATEMENT - IFN(S) - 0| SUBROUTINE MINMAX(TMIN,TMAX,T ) CF5 0010 | C CF5 0020 | C SUBROUTINE TO FIND MINIMUM AND MAXIMUM OF ANY GIVEN ARRAY CF5 0030 | C CF5 0040 DIMENSION Z(31,31), C(30,30), DATA(51) CF5 0050 | DIMENSION BETA(30), AVG(30), SIGMA(30), STDDEV(30), B(30) CF5 0060 | DIMENSION VAR(30), YRESID(30) CF5 0070 | DIMENSION ITRFM(50), VARTR(50), AVATR(50), ITLOC(50), IOMIT(50) CF5 0080 | DIMENSION FMT(80), IQ(40) CF5 0090 DIMENSION IDZA(1000),IDZB(1000),DELTA(1000) CF5 0100 | DIMENSION YYIN(1000),YYCC(1000),LSORT(1000) CF5 0110 | COMMON YYMIN, YCMIN, YYMAX, YCMAX CF5 0120 | COMMON FMT, IQ, KSL1, KSL2, KSL3, KSL4, KSL5, KSL6 CF5 0130 | COMMON IDZA, IDZB, DELTA, YYIN, YYCC, LSORT CF5 0140 COMMON Z,C, BETA,AVG,SIGMA, B,VAR,YRESID CF5 0150 | COMMON ITRFM,VARTR,AVATR, IOMIT CF5 0160 | COMMON KTIN,KTOU,KTPCH,KTBIN1,KTBIN2,KTBIN3,KTBIN4,KTBIN5 CF5 0170 | COMMON PR1,PR2,PR3 CF5 0180 | COMMON BZRO, KNOIND, KNODEP, KNOVAR, KNVP1, JNOIND, JNODEP, JNOVARCF5 0190 COMMON NOBSP, NOERR, JOBSV, NRDLMT, NOOBSV, MNO, JSEQ, N, M, J, I CF5 0200 | COMMON LNOVAR, NNO, ISEQ, L, KNO, XNOOBS, XNOVAR, XNOIND, DEFR CF5 0210 | COMMON NOVM1, RTDEFR, IPL1, K, RSQD, BZERO, SUMB, ICURY, FACTOR CF5 0220 | COMMON SGBETA, TBETA, PCORK, SIGB, FRATIO, UNBEST, BIEST, SSQ, SSECF5 0230 | COMMON IPLOT1, IPLOT2, YCALC, DIFF, IREC, X, IDENT, IOBSV, MP CF5 0240 COMMON NOVAR, NVP1, IPATH CF5 0250 | EQUIVALENCE (C,DATA), (C(101),STDDEV), (C(201),ITLOC) CF5 0260 | EQUIVALENCE (NOIND,IQ(1)),(NODEP,IQ(2)) CF5 0270 | IF (T-TMIN) 2,3,3 CF5 0280 | 2 TMIN = T CF5 0290 3 IF (T-TMAX) 5,5,4 CF5 0300 | 4 TMAX = T CF5 0310 | 5 RETURN CF5 0320 | END CF5 0330 | 1 12/01/64 PAGE 28 | 0| -| $IBFTC CF6 $IBFTC CF6 0000 | 1 12/01/64 PAGE 29 CF6 - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE YBACK (KQ, ZBAK) CF6 0010 | C CF6 0020 | C SUBROUTINE TO PERFORM BACK TRANSFORMATION OF VARIABLES CF6 0030 C CF6 0040 | DIMENSION Z(31,31), C(30,30), DATA(51) CF6 0050 | DIMENSION BETA(30), AVG(30), SIGMA(30), STDDEV(30), B(30) CF6 0060 | DIMENSION VAR(30), YRESID(30) CF6 0070 | DIMENSION ITRFM(50), VARTR(50), AVATR(50), ITLOC(50), IOMIT(50) CF6 0080 DIMENSION FMT(80), IQ(40) CF6 0090 | DIMENSION IDZA(1000),IDZB(1000),DELTA(1000) CF6 0100 | DIMENSION YYIN(1000),YYCC(1000),LSORT(1000) CF6 0110 | COMMON YYMIN, YCMIN, YYMAX, YCMAX CF6 0120 | COMMON FMT, IQ, KSL1, KSL2, KSL3, KSL4, KSL5, KSL6 CF6 0130 COMMON IDZA, IDZB, DELTA, YYIN, YYCC, LSORT CF6 0140 | COMMON Z,C, BETA,AVG,SIGMA, B,VAR,YRESID CF6 0150 | COMMON ITRFM,VARTR,AVATR, IOMIT CF6 0160 | COMMON KTIN,KTOU,KTPCH,KTBIN1,KTBIN2,KTBIN3,KTBIN4,KTBIN5 CF6 0170 | COMMON PR1,PR2,PR3 CF6 0180 COMMON BZRO, KNOIND, KNODEP, KNOVAR, KNVP1, JNOIND, JNODEP, JNOVARCF6 0190 | COMMON NOBSP, NOERR, JOBSV, NRDLMT, NOOBSV, MNO, JSEQ, N, M, J, I CF6 0200 | COMMON LNOVAR, NNO, ISEQ, L, KNO, XNOOBS, XNOVAR, XNOIND, DEFR CF6 0210 | COMMON NOVM1, RTDEFR, IPL1, K, RSQD, BZERO, SUMB, ICURY, FACTOR CF6 0220 | COMMON SGBETA, TBETA, PCORK, SIGB, FRATIO, UNBEST, BIEST, SSQ, SSECF6 0230 COMMON IPLOT1, IPLOT2, YCALC, DIFF, IREC, X, IDENT, IOBSV, MP CF6 0240 | COMMON NOVAR, NVP1, IPATH CF6 0250 | EQUIVALENCE (C,DATA), (C(101),STDDEV), (C(201),ITLOC) CF6 0260 | EQUIVALENCE (NOIND,IQ(1)),(NODEP,IQ(2)) CF6 0270 | LBJ = ITRFM(KQ) CF6 0280 KV = 1 CF6 0290 | IF( KQ .EQ. 0 ) KV = 2 CF6 0300 | VK = VARTR(KQ) CF6 0310 | IF( LBJ .LT.1) GO TO 9 CF6 0320 | IF( LBJ .GT.8) GO TO 10 CF6 0330 GO TO(1,2,3,4,5,6,7,8 ), LBJ CF6 0340 | 1 ZBAK = EXP(ZBAK) CF6 0350 16 | GO TO 9 CF6 0360 | 2 ZBAK = EXP(2.302585093 * ZBAK) CF6 0370 19 | GO TO 9 CF6 0380 3 ZBAK = ALOG(ZBAK) CF6 0390 22 | GO TO 9 CF6 0400 | 4 ZBAK = ALOG10(ZBAK) CF6 0410 25 | GO TO 9 CF6 0420 | 5 GO TO(51, 9), KV CF6 0430 51 ZBAK = EXP( ALOG(ZBAK) / VK) CF6 0440 29 30 | GO TO 9 CF6 0450 | 6 GO TO(61, 9), KV CF6 0460 | 61 ZBAK = ZBAK / VK CF6 0470 | GO TO 9 CF6 0480 7 GO TO(71, 9), KV CF6 0490 | 71 ZBAK = ZBAK * VK CF6 0500 | GO TO 9 CF6 0510 | 8 GO TO(81, 9), KV CF6 0520 | 81 ZBAK = ZBAK - VK CF6 0530 9 RETURN CF6 0540 | 10 IF( LBJ .LT. 13) GO TO 9 CF6 0550 | 1 12/01/64 PAGE 30 | CF6 - EFN SOURCE STATEMENT - IFN(S) - | 0 KBJ = LBJ - 12 CF6 0560 GO TO (13,14,15,16,17,18,19,20), KBJ CF6 0570 | 13 ZBAK = ARSIN(ZBAK) CF6 0580 46 | GO TO 9 CF6 0590 | 14 ZBAK = ARCOS(ZBAK) CF6 0600 49 | GO TO 9 CF6 0610 15 GO TO(151, 9), KV CF6 0620 | 151 ZBAK = VK / ZBAK CF6 0630 | GO TO 9 CF6 0640 | 16 ZBAK = EXP(EXP(ZBAK)) CF6 0650 55 56 | GO TO 9 CF6 0660 17 ZBAK = EXP(2.302585093 * ZBAK) CF6 0670 59 | GO TO 2 CF6 0680 | 18 GO TO(181, 9), KV CF6 0690 | 181 ZBAK = EXP(2.302585093 * ZBAK) - VK CF6 0700 63 | GO TO 9 CF6 0710 19 ZBAK = ALOG(ZBAK) CF6 0720 66 | GO TO 3 CF6 0730 | 20 ZBAK = ALOG10(ZBAK) CF6 0740 69 | GO TO 4 CF6 0750 | END CF6 0760 1 12/01/64 PAGE 31 | | 0| -| $IBFTC CF7 $IBFTC CF7 0000 1 12/01/64 PAGE 32 | CF7 - EFN SOURCE STATEMENT - IFN(S) - | 0| SUBROUTINE PITCHA CF7 0010 | C CF7 0020 C SUBROUTINE PLOTS CUMULATIVE FREQUENCY OF RESIDUALS AND CF7 0030 | C RESIDUALS VS. FITTED Y CF7 0040 | C CF7 0050 | DIMENSION GRIDA(53,17),GRIDB(53,17),NSW(6),PRAX(51) CF7 0060 | DIMENSION XAX(11),RESID(9),POS(8),NEG(8) CF7 0070 DIMENSION Z(31,31), C(30,30), DATA(51) CF7 0080 | DIMENSION BETA(30), AVG(30), SIGMA(30), STDDEV(30), B(30) CF7 0090 | DIMENSION VAR(30), YRESID(30) CF7 0100 | DIMENSION ITRFM(50), VARTR(50), AVATR(50), ITLOC(50), IOMIT(50) CF7 0110 | DIMENSION FMT(80), IQ(40) CF7 0120 DIMENSION IDZA(1000),IDZB(1000),DELTA(1000) CF7 0130 | DIMENSION YYIN(1000),YYCC(1000),LSORT(1000) CF7 0140 | COMMON YYMIN, YCMIN, YYMAX, YCMAX CF7 0150 | COMMON FMT, IQ, KSL1, KSL2, KSL3, KSL4, KSL5, KSL6 CF7 0160 | COMMON IDZA, IDZB, DELTA, YYIN, YYCC, LSORT CF7 0170 COMMON Z,C, BETA,AVG,SIGMA, B,VAR,YRESID CF7 0180 | COMMON ITRFM,VARTR,AVATR, IOMIT CF7 0190 | COMMON KTIN,KTOU,KTPCH,KTBIN1,KTBIN2,KTBIN3,KTBIN4,KTBIN5 CF7 0200 | COMMON PR1,PR2,PR3 CF7 0210 | COMMON BZRO, KNOIND, KNODEP, KNOVAR, KNVP1, JNOIND, JNODEP, JNOVARCF7 0220 COMMON NOBSP, NOERR, JOBSV, NRDLMT, NOOBSV, MNO, JSEQ, N, M, J, I CF7 0230 | COMMON LNOVAR, NNO, ISEQ, L, KNO, XNOOBS, XNOVAR, XNOIND, DEFR CF7 0240 | COMMON NOVM1, RTDEFR, IPL1, K, RSQD, BZERO, SUMB, ICURY, FACTOR CF7 0250 | COMMON SGBETA, TBETA, PCORK, SIGB, FRATIO, UNBEST, BIEST, SSQ, SSECF7 0260 | COMMON IPLOT1, IPLOT2, YCALC, DIFF, IREC, X, IDENT, IOBSV, MP CF7 0270 COMMON NOVAR, NVP1, IPATH CF7 0280 | EQUIVALENCE (C,DATA), (C(101),STDDEV), (C(201),ITLOC) CF7 0290 | EQUIVALENCE (NOIND,IQ(1)),(NODEP,IQ(2)) CF7 0300 | DATA POS/1HP,1HO,1HS,1HI,1HT,1HI, 1HV,1HE/,NEG/1HN,1HE,1HG,1HA,1HTCF7 0310 | 1,1HI,1HV,1HE/ CF7 0320 PRAX(01) = 0.000233 CF7 0330 | PRAX(02) = 0.000302 CF7 0340 | PRAX(03) = 0.000390 CF7 0350 | PRAX(04) = 0.000501 CF7 0360 | PRAX(05) = 0.000641 CF7 0370 PRAX(06) = 0.000816 CF7 0380 | PRAX(07) = 0.001035 CF7 0390 | PRAX(08) = 0.001306 CF7 0400 | PRAX(09) = 0.001641 CF7 0410 | PRAX(10) = 0.002052 CF7 0420 PRAX(11) = 0.002555 CF7 0430 | PRAX(12) = 0.003167 CF7 0440 | PRAX(13) = 0.003907 CF7 0450 | PRAX(14) = 0.004797 CF7 0460 | PRAX(15) = 0.005868 CF7 0470 PRAX(16) = 0.007143 CF7 0480 | PRAX(17) = 0.008656 CF7 0490 | PRAX(18) = 0.010444 CF7 0500 | PRAX(19) = 0.012545 CF7 0510 | PRAX(20) = 0.015003 CF7 0520 PRAX(21) = 0.017864 CF7 0530 | PRAX(22) = 0.021178 CF7 0540 | PRAX(23) = 0.025588 CF7 0550 | 1 12/01/64 PAGE 33 | CF7 - EFN SOURCE STATEMENT - IFN(S) - 0 PRAX(24) = 0.029379 CF7 0560 | PRAX(25) = 0.034380 CF7 0570 | PRAX(26) = 0.040059 CF7 0580 | PRAX(27) = 0.046479 CF7 0590 | PRAX(28) = 0.053699 CF7 0600 PRAX(29) = 0.061780 CF7 0610 | PRAX(30) = 0.070781 CF7 0620 | PRAX(31) = 0.080757 CF7 0630 | PRAX(32) = 0.091759 CF7 0640 | PRAX(33) = 0.103835 CF7 0650 PRAX(34) = 0.117023 CF7 0660 | PRAX(35) = 0.131357 CF7 0670 | PRAX(36) = 0.146859 CF7 0680 | PRAX(37) = 0.163543 CF7 0690 | PRAX(38) = 0.181411 CF7 0700 PRAX(39) = 0.200454 CF7 0710 | PRAX(40) = 0.220650 CF7 0720 | PRAX(41) = 0.241964 CF7 0730 | PRAX(42) = 0.264347 CF7 0740 | PRAX(43) = 0.287740 CF7 0750 PRAX(44) = 0.312067 CF7 0760 | PRAX(45) = 0.337243 CF7 0770 | PRAX(46) = 0.363169 CF7 0780 | PRAX(47) = 0.389739 CF7 0790 | PRAX(48) = 0.416834 CF7 0800 PRAX(49) = 0.444330 CF7 0810 | PRAX(50) = 0.472097 CF7 0820 | PRAX(51) = 0.500000 CF7 0830 | XNP = NOOBSV CF7 0840 | CALL ATHRUZ(NSW(1), 6H ) CF7 0850 2 CALL ATHRUZ(NSW(2), 6H / ) CF7 0860 4 | CALL ATHRUZ(NSW(3), 6H S ) CF7 0870 6 | CALL ATHRUZ(NSW(4), 6H T ) CF7 0880 8 | CALL ATHRUZ(NSW(5), 6H U ) CF7 0890 10 | CALL ATHRUZ(NSW(6), 6H V ) CF7 0900 12 JJC = LSORT(NOOBSV) CF7 0910 | DELP = DELTA(JJC) CF7 0920 | IIC = LSORT(1) CF7 0930 | YIW = (1.005 * DELP - DELTA(IIC) ) / 51.0 CF7 0940 | XIW = (1.005 * YCMAX - YCMIN )/101.0 CF7 0950 DO 200 IA=1,NOOBSV CF7 0960 | JAC=LSORT(IA) CF7 0970 | IF(DELTA(JAC).LT.0.) GO TO 200 CF7 0980 | L1=INT((DELP-DELTA(JAC))/YIW)+2 CF7 0990 | GO TO 202 CF7 1000 200 CONTINUE CF7 1010 | 202 CALL GRID(GRIDA,GRIDB,L1) CF7 1020 34 | DO 30 IA = 1, NOOBSV CF7 1030 | JAC = LSORT(IA) CF7 1040 | LINE=INT((DELP-DELTA(JAC))/YIW)+2 CF7 1050 IF(LINE.LT.2.OR.LINE.GT.52) GO TO 30 CF7 1060 | LOCX1 = INT((YYCC(JAC) - YCMIN)/XIW) + 1 CF7 1070 | IF((LOCX1.LT.1).OR.(LOCX1.GT.101)) GO TO 10 CF7 1080 | IF(LOCX1.EQ.1) LOCX1=2 CF7 1090 | LOCX2= LOCX1 + 6 CF7 1100 LCHAR = MOD(LOCX2,6) CF7 1110 | 1 12/01/64 PAGE 34 | CF7 - EFN SOURCE STATEMENT - IFN(S) - | 0 LWORD = LOCX2/6 CF7 1120 | IF(LCHAR) 4,2,4 CF7 1130 2 LCHAR = 6 CF7 1140 | LWORD = LWORD - 1 CF7 1150 | 4 NCHAR = NSW(LCHAR) CF7 1160 | CALL PACK(GRIDA(LINE,LWORD), NCHAR, 6H . ) CF7 1170 62 | 10 PRLOC = (FLOAT(IA) - 0.5)/XNP CF7 1180 IF(PRLOC - 0.5) 12,12,18 CF7 1190 | 12 DO 14 IB = 1,51 CF7 1200 | IF(PRLOC - PRAX(IB)) 16,16,14 CF7 1210 | 14 CONTINUE CF7 1220 | IB = 51 CF7 1230 16 LOCP1 = IB CF7 1240 | GO TO 24 CF7 1250 | 18 DO 20 IB = 1, 51 CF7 1260 | IC = 52 - IB CF7 1270 | IF(PRLOC + PRAX(IC) - 1.0) 22,22,20 CF7 1280 20 CONTINUE CF7 1290 | 22 LOCP1 = IB + 50 CF7 1300 | 24 LOCP2 = LOCP1 + 6 CF7 1310 | LCHAR = MOD(LOCP2,6) CF7 1320 | LWORD = LOCP2/6 CF7 1330 IF(LCHAR) 28,26,28 CF7 1340 | 26 LCHAR = 6 CF7 1350 | LWORD = LWORD - 1 CF7 1360 | 28 NCHAR = NSW(LCHAR) CF7 1370 | CALL PACK(GRIDB(LINE,LWORD), NCHAR, 6H + ) CF7 1380 94 30 CONTINUE CF7 1390 | CALL ATHRUZ(RESID(1), 6H R) CF7 1400 99 | CALL ATHRUZ(RESID(2), 6H E) CF7 1410 101 | CALL ATHRUZ(RESID(3), 6H S) CF7 1420 103 | CALL ATHRUZ(RESID(4), 6H I) CF7 1430 105 CALL ATHRUZ(RESID(5), 6H D) CF7 1440 107 | CALL ATHRUZ(RESID(6), 6H U) CF7 1450 109 | CALL ATHRUZ(RESID(7), 6H A) CF7 1460 111 | CALL ATHRUZ(RESID(8), 6H L) CF7 1470 113 | CALL ATHRUZ(RESID(9),6H S) CF7 1480 115 XAX(1) = YCMIN + 0.5 * XIW CF7 1490 | XIW = XIW * 10.0 CF7 1500 | DO 34 IA = 2,11 CF7 1510 | 34 XAX(IA) = XAX(IA - 1) + XIW CF7 1520 | WRITE(KTOU,102) PR1,PR2,PR3, ICURY CF7 1530 124 102 FORMAT(1H1,3A6,10H DEP VAR ,I2,23X, CF7 1540 | 136HCUMULATIVE DISTRIBUTION OF RESIDUALS) CF7 1550 | WRITE(KTOU,104) CF7 1560 125 | 104 FORMAT(1H ) CF7 1570 | WRITE(KTOU,106) CF7 1580 126 106 FORMAT(1H ,17X,103H.0002 .001 .005 .01 .02 .05 .1 .2 .3CF7 1590 | 1 .4 .5 .6 .7 .8 .9 .95 .98 .99 .995 .999 ) CF7 1600 | DO 40 IA=1,3 CF7 1610 | 40 WRITE(KTOU,108) (GRIDB(IA,JZ),JZ=1,17) CF7 1620 130 | 108 FORMAT(1H ,19X,17A6) CF7 1630 DO 42 IA=4,11 CF7 1640 | IB=IA-3 CF7 1650 | 42 WRITE(KTOU,44) POS(IB),(GRIDB(IA,JZ),JZ=1,17) CF7 1660 140 | 44 FORMAT(1H ,17X,A1,1X,17A6) CF7 1670 | 1 12/01/64 PAGE 35 CF7 - EFN SOURCE STATEMENT - IFN(S) - | 0 DO 46 IA=12,22 CF7 1680 | IF(IA.EQ.L1) GO TO 300 CF7 1690 | WRITE(KTOU,108) (GRIDB(IA,JZ),JZ=1,17) CF7 1700 153 | GO TO 46 CF7 1710 300 WRITE(KTOU,302) (GRIDB(IA,JZ),JZ=1,17) CF7 1720 159 | 302 FORMAT(1H ,17X,1H0,1X,17A6) CF7 1730 | 46 CONTINUE CF7 1740 | DO 48 IA=23,31 CF7 1750 | IB=IA-22 CF7 1760 IF(IA.EQ.L1) GO TO 304 CF7 1770 | WRITE(KTOU,110) RESID(IB),(GRIDB(IA,JZ),JZ=1,17) CF7 1780 173 | 110 FORMAT(1H ,A6,13X,17A6) CF7 1790 | GO TO 48 CF7 1800 | 304 WRITE(KTOU,306) RESID(IB),(GRIDB(IA,JZ),JZ=1,17) CF7 1810 180 306 FORMAT(1H ,A6,11X,1H0,1X,17A6) CF7 1820 | 48 CONTINUE CF7 1830 | DO 308 IA=32,42 CF7 1840 | IF(IA.EQ.L1) GO TO 310 CF7 1850 | WRITE(KTOU,108) (GRIDB(IA,JZ),JZ=1,17) CF7 1860 194 GO TO 308 CF7 1870 | 310 WRITE(KTOU,302) (GRIDB(IA,JZ),JZ=1,17) CF7 1880 200 | 308 CONTINUE CF7 1890 | DO 312 IA=43,50 CF7 1900 | IB=IA-42 CF7 1910 312 WRITE(KTOU,44) NEG(IB),(GRIDB(IA,JZ),JZ=1,17) CF7 1920 211 | DO 314 IA=51,53 CF7 1930 | 314 WRITE(KTOU,108) (GRIDB(IA,JZ),JZ=1,17) CF7 1940 222 | WRITE(KTOU,106) CF7 1950 228 | WRITE(KTOU,112) CF7 1960 229 112 FORMAT(1H0,54X,33HCUMULATIVE FREQUENCY, NORMAL GRID) CF7 1970 | WRITE(KTOU,122) PR1,PR2,PR3,ICURY CF7 1980 230 | 122 FORMAT(1H1,3A6,10H DEP VAR ,I2,29X,21HRESIDUAL VS. FITTED Y ) CF7 1990 | WRITE(KTOU,104) CF7 2000 231 | WRITE(KTOU,126) (XAX(JZ), JZ = 1,11 ) CF7 2010 232 126 FORMAT(1H ,13X,11F10.3 ) CF7 2020 | DO 50 IA=1,3 CF7 2030 | 50 WRITE(KTOU,108) (GRIDA(IA,JZ),JZ=1,17) CF7 2040 243 | DO 52 IA=4,11 CF7 2050 | IB=IA-3 CF7 2060 52 WRITE(KTOU,44) POS(IB),(GRIDA(IA,JZ),JZ=1,17) CF7 2070 253 | DO 54 IA=12,22 CF7 2080 | IF(IA.EQ.L1) GO TO 56 CF7 2090 | WRITE(KTOU,108) (GRIDA(IA,JZ),JZ=1,17) CF7 2100 266 | GO TO 54 CF7 2110 56 WRITE(KTOU,302) (GRIDA(IA,JZ),JZ=1,17) CF7 2120 272 | 54 CONTINUE CF7 2130 | DO 58 IA=23,31 CF7 2140 | IB=IA-22 CF7 2150 | IF(IA.EQ.L1) GO TO 60 CF7 2160 WRITE(KTOU,110) RESID(IB),(GRIDA(IA,JZ),JZ=1,17) CF7 2170 286 | GO TO 58 CF7 2180 | 60 WRITE(KTOU,306) RESID(IB),(GRIDA(IA,JZ),JZ=1,17) CF7 2190 293 | 58 CONTINUE CF7 2200 | DO 62 IA=32,42 CF7 2210 IF(IA.EQ.L1) GO TO 64 CF7 2220 | WRITE(KTOU,108) (GRIDA(IA,JZ),JZ=1,17) CF7 2230 307 | 1 12/01/64 PAGE 36 | CF7 - EFN SOURCE STATEMENT - IFN(S) - | 0 GO TO 62 CF7 2240 64 WRITE(KTOU,302) (GRIDA(IA,JZ),JZ=1,17) CF7 2250 313 | 62 CONTINUE CF7 2260 | DO 66 IA=43,50 CF7 2270 | IB=IA-42 CF7 2280 | 66 WRITE(KTOU,44) NEG(IB),(GRIDA(IA,JZ),JZ=1,17) CF7 2290 324 DO 68 IA=50,53 CF7 2300 | 68 WRITE(KTOU,108) (GRIDA(IA,JZ),JZ=1,17) CF7 2310 335 | WRITE(KTOU,126) (XAX(JZ), JZ = 1,11 ) CF7 2320 341 | WRITE(KTOU,132) CF7 2330 348 | 132 FORMAT(1H0,66X, 8HFITTED Y ) CF7 2340 RETURN CF7 2350 | END CF7 2360 | 1 12/01/64 PAGE 37 | | 0 -| $IBFTC CF8 $IBFTC CF8 0000 | 1 12/01/64 PAGE 38 | CF8 - EFN SOURCE STATEMENT - IFN(S) - | 0 SUBROUTINE GRID(GRIDA,GRIDB,L1) CF8 0010 | C CF8 0020 | C SUBROUTINE MAKES THE GRID FOR THE PLOTS CF8 0030 | C CF8 0040 | DIMENSION GRIDA(53,17),GRIDB(53,17) CF8 0050 CALL ATHRUZ(BLANK,6H ) CF8 0060 2 | CALL ATHRUZ(BLINK,6H ' ) CF8 0070 4 | DO 42 IZ=2,52 CF8 0080 | DO 40 JZ = 2,16 CF8 0090 | GRIDA(IZ,JZ) = BLANK CF8 0100 40 GRIDB(IZ,JZ) = BLANK CF8 0110 | 42 GRIDB(IZ,9)=BLINK CF8 0120 | CALL ATHRUZ( GRIDA( 1, 1) ,6H*----*) CF8 0130 19 | CALL ATHRUZ( GRIDA( 1, 2) ,6H----*-) CF8 0140 21 | CALL ATHRUZ( GRIDA( 1, 3) ,6H---*--) CF8 0150 23 CALL ATHRUZ( GRIDA( 1, 4) ,6H--*---) CF8 0160 25 | CALL ATHRUZ( GRIDA( 1, 5) ,6H-*----) CF8 0170 27 | GRIDA( 1,16) = GRIDA( 1, 1) CF8 0180 | GRIDA( 1,17) = GRIDA( 1, 2) CF8 0190 | DO 10 JA = 1,5 CF8 0200 GRIDA(1, JA + 5) = GRIDA(1,JA) CF8 0210 | 10 GRIDA(1, JA +10) = GRIDA(1,JA) CF8 0220 | CALL ATHRUZ( GRIDB( 1, 3) ,6H--*--*) CF8 0230 39 | CALL ATHRUZ( GRIDB( 1, 8) ,6H-*---*) CF8 0240 41 | CALL ATHRUZ( GRIDB( 1,15) ,6H-*--*-) CF8 0250 43 CALL ATHRUZ( GRIDB( 1,16) ,6H------) CF8 0260 45 | CALL ATHRUZ( GRIDB( 1, 1), 6H*-----) CF8 0270 47 | GRIDB( 1, 2) = GRIDB( 1, 1) CF8 0280 | GRIDB( 1, 4) = GRIDA( 1, 3) CF8 0290 | GRIDB( 1, 5) = GRIDA( 1, 3) CF8 0300 GRIDB( 1, 6) = GRIDA( 1, 4) CF8 0310 | GRIDB( 1, 7) = GRIDA( 1, 4) CF8 0320 | GRIDB( 1, 9) = GRIDA( 1, 3) CF8 0330 | GRIDB( 1,10) = GRIDB( 1, 8) CF8 0340 | GRIDB( 1,11) = GRIDA( 1, 2) CF8 0350 GRIDB( 1,12) = GRIDA( 1, 2) CF8 0360 | GRIDB( 1,13) = GRIDA( 1, 3) CF8 0370 | GRIDB( 1,14) = GRIDA( 1, 3) CF8 0380 | GRIDB( 1,17) = GRIDB( 1, 1) CF8 0390 | DO 12 IA = 1,17 CF8 0400 GRIDA(53,IA)=GRIDA(1,IA) CF8 0410 | 12 GRIDB(53,IA)=GRIDB(1,IA) CF8 0420 | CALL ATHRUZ(GRIDA( 2, 1), 6HI ) CF8 0430 59 | CALL ATHRUZ(GRIDA( 2,17), 6H I ) CF8 0440 61 | DO 14 IA=3,52 CF8 0450 GRIDA(IA, 1) = GRIDA(IA-1, 1) CF8 0460 | 14 GRIDA(IA,17) = GRIDA(IA-1,17) CF8 0470 | DO 30 IA=2,16 CF8 0480 | GRIDA(L1,IA)=GRIDB(1,16) CF8 0490 | GRIDB(L1,IA)=GRIDB(1,16) CF8 0500 30 CONTINUE CF8 0510 | CALL ATHRUZ(HOR1,6HI-----) CF8 0520 81 | CALL ATHRUZ(HOR2,6H-----I) CF8 0530 83 | GRIDA(L1,1)=HOR1 CF8 0540 | GRIDA(L1,17)=HOR2 CF8 0550 1 12/01/64 PAGE 39 | CF8 - EFN SOURCE STATEMENT - IFN(S) - | 0 19 DO 20 IA=2,52 CF8 0560 | GRIDB(IA, 1) = GRIDA(IA, 1) CF8 0570 | 20 GRIDB(IA,17) = GRIDA(IA,17) CF8 0580 RETURN CF8 0590 | END CF8 0600 | 1 12/01/64 PAGE 40 | CF1 STORAGE MAP | 0 MAIN PROGRAM - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 20370 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | YYMIN 00000 R YCMIN 00001 R YYMAX 00002 R | YCMAX 00003 R FMT 00004 R IQ 00124 I KSL1 00174 I KSL2 00175 I KSL3 00176 I | KSL4 00177 I KSL5 00200 I KSL6 00201 I | IDZA 00202 I IDZB 02152 I DELTA 04122 R | YYIN 06072 R YYCC 10042 R LSORT 12012 I | Z 13762 R C 15663 R BETA 17467 R AVG 17525 R SIGMA 17563 R B 17621 R | VAR 17657 R YRESID 17715 R ITRFM 17753 I | VARTR 20035 R AVATR 20117 R IOMIT 20201 I | KTIN 20263 I KTOU 20264 I KTPCH 20265 I | KTBIN1 20266 I KTBIN2 20267 I KTBIN3 20270 I KTBIN4 20271 I KTBIN5 20272 I PR1 20273 R | PR2 20274 R PR3 20275 R BZRO 20276 R | KNOIND 20277 I KNODEP 20300 I KNOVAR 20301 I | KNVP1 20302 I JNOIND 20303 I JNODEP 20304 I | JNOVAR 20305 I NOBSP 20306 I NOERR 20307 I JOBSV 20310 I NRDLMT 20311 I NOOBSV 20312 I | MNO 20313 I JSEQ 20314 I N 20315 I | M 20316 I J 20317 I I 20320 I | LNOVAR 20321 I NNO 20322 I ISEQ 20323 I | L 20324 I KNO 20325 I XNOOBS 20326 R XNOVAR 20327 R XNOIND 20330 R DEFR 20331 R | NOVM1 20332 I RTDEFR 20333 R IPL1 20334 I | K 20335 I RSQD 20336 R BZERO 20337 R | SUMB 20340 R ICURY 20341 I FACTOR 20342 R | SGBETA 20343 R TBETA 20344 R PCORK 20345 R SIGB 20346 R FRATIO 20347 R UNBEST 20350 R | BIEST 20351 R SSQ 20352 R SSE 20353 R | IPLOT1 20354 I IPLOT2 20355 I YCALC 20356 R | DIFF 20357 R IREC 20360 I X 20361 R | IDENT 20362 I IOBSV 20363 I MP 20364 I NOVAR 20365 I NVP1 20366 I IPATH 20367 I | DATA 15663 R STDDEV 16027 R ITLOC 16173 I | NOIND 00124 I NODEP 00125 I KFM 00151 I | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE ARRAY 20371 R BRRAY 20371 R CRRAY 20477 R | DRRAY 20567 R ERRAY 20573 R AMAX 20656 R | AMIN 20740 R W1 21022 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE 1 12/01/64 PAGE 41 | CF1 STORAGE MAP | 0IWDA 21104 I BLANK 21105 R XXFRT 21106 R | KCP 21107 I IQF 21110 I IQE 21111 I | NOEQ 21112 I KZ 21113 I LZ 21114 I J2 21115 I J1 21116 I READIN 21117 R | KKK 21120 I KONE 21121 I KTWO 21122 I | INSQOB 21123 I INSQCD 21124 I INTRFM 21125 I | MCOUNT 21126 I IOPCDE 21127 I NAAA 21130 I | FRRAY 21131 R NAAB 21132 I GRRAY 21133 R HRRAY 21134 R JAA 21135 I MA 21136 I | NA 21137 I JZA 21140 I BWY 21141 R | SSYT 21142 R NDEFR 21143 I IB 21144 I | ISA 21145 I TEMP 21146 R ISAX 21147 I | ISB 21150 I IBB 21151 I KSA 21152 I IC 21153 I JC 21154 I ME 21155 I | - ENTRY POINTS | 0 ...... SECTION 3 | - SUBROUTINES CALLED | 0 ATHRUZ SECTION 4 .FVIO. SECTION 5 .FRDD. SECTION 6 .FRWT. SECTION 7 .FWRD. SECTION 8 EXIT SECTION 9 | REDATA SECTION 10 .FWRB. SECTION 11 .FRDB. SECTION 12 | AE04 SECTION 13 SQRT SECTION 14 AE02 SECTION 15 | YBACK SECTION 16 MINMAX SECTION 17 PITCHA SECTION 18 | .FXEM. SECTION 19 .FRTN. SECTION 20 .FCNV. SECTION 21 .FFIL. SECTION 22 .FWLR. SECTION 23 .FBLT. SECTION 24 | .FBDT. SECTION 25 .FRLR. SECTION 26 E.1 SECTION 27 | E.2 SECTION 28 E.3 SECTION 29 E.4 SECTION 30 | CC.1 SECTION 31 CC.2 SECTION 32 CC.3 SECTION 33 | CC.4 SECTION 34 SYSLOC SECTION 35 - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 10 11A 22463 800 FORMAT 21220 20 24A 22540 | 15 37A 22571 18 42A 22576 22 45A 22612 | 21 44A 22601 805 FORMAT 21225 26 59A 22730 24 FORMAT 21242 28 60A 22740 32 62A 22754 | 36 65A 22774 30 FORMAT 21252 38 68A 23015 | 34 FORMAT 21260 80 104A 23232 42 72A 23046 | 40 FORMAT 21271 44 73A 23056 52 75A 23072 | 72 103A 23217 50 FORMAT 21302 56 FORMAT 21311 62 82A 23124 2022 99A 23214 1030 FORMAT 21320 | 64 FORMAT 21321 74 FORMAT 21327 17 109A 23252 | 16 106A 23235 980 FORMAT 21336 86 111A 23267 | 985 FORMAT 21341 47 120A 23320 46 113A 23272 | 54 162A 23472 48 122A 23323 6001 125A 23335 6003 126A 23340 6002 127A 23341 975 FORMAT 21346 | 990 FORMAT 21352 58 152A 23441 992 FORMAT 21355 | 60 165A 23513 55 164A 23474 815 FORMAT 21361 | 65 170A 23522 70 172A 23525 75 173A 23527 | 1 12/01/64 PAGE 42 CF1 STORAGE MAP | 076 176A 23534 66 178A 23537 77 179A 23541 | 67 183A 23546 68 186A 23553 69 188A 23557 | 71 189A 23561 83 191A 23564 144 255A 24060 | 78 207A 23631 73 196A 23577 84 193A 23570 82 217A 23662 820 FORMAT 21436 85 219A 23665 | 142 245A 24024 125 232A 23753 100 221A 23670 | 105 222A 23677 115 228A 23730 110 225A 23705 | 840 FORMAT 21442 120 230A 23733 835 FORMAT 21455 | 140 238A 24005 130 234A 23756 135 236A 23762 845 FORMAT 21467 141 242A 24015 145 265A 24110 | 143 247A 24027 146 267A 24113 250 346A 24434 | 149 294A 24227 147 269A 24115 825 FORMAT 21477 | 170 277A 24150 148 279A 24154 166 284A 24171 | 165 282A 24165 167 285A 24173 830 FORMAT 21503 150 297A 24234 185 306A 24276 247 344A 24431 | 155 299A 24240 160 300A 24247 855 FORMAT 21506 | 230 310A 24304 190 308A 24301 235 312A 24307 | 240 314A 24312 245 324A 24356 255 348A 24436 | 265 350A 24462 860 FORMAT 21531 266 357A 24526 267 359A 24544 935 FORMAT 21557 269 362A 24572 | 268 361A 24547 1020 FORMAT 21577 3032 444A 25226 | 264 364A 24606 270 366A 24610 863 FORMAT 21600 | 995 FORMAT 21605 1010 428A 25112 272 372A 24647 | 1001 383A 24666 1002 384A 24670 1003 393A 24704 1004 395A 24706 1005 402A 24723 1006 403A 24725 | 1008 415A 25001 1007 FORMAT 21631 1009 FORMAT 21642 | 299 434A 25135 296 432A 25121 297 FORMAT 21654 | 282 438A 25144 281 436A 25140 283 440A 25160 | 284 442A 25174 286 FORMAT 21702 287 FORMAT 21724 864 FORMAT 21746 862 FORMAT 22005 3000 451A 25236 | 275 488A 25352 3008 478A 25321 3002 468A 25311 | 3004 471A 25313 3006 475A 25317 3034 481A 25326 | 3036 483A 25330 865 FORMAT 22032 3010 497A 25363 | 280 500A 25365 290 510A 25432 285 502A 25367 870 FORMAT 22040 875 FORMAT 22044 305 523A 25511 | 295 512A 25434 885 FORMAT 22047 300 520A 25507 | 890 FORMAT 22103 307 538A 25565 306 525A 25516 | 325 578A 25772 308 527A 25521 310 545A 25630 | 314 565A 25713 312 555A 25650 315 567A 25715 895 FORMAT 22106 320 571A 25740 350 594A 26070 | 330 585A 26011 335 587A 26027 340 590A 26060 | 900 FORMAT 22116 905 FORMAT 22133 906 FORMAT 22150 | 353 607A 26137 355 610A 26144 360 619A 26175 | 365 641A 26300 362 625A 26204 910 FORMAT 22154 915 FORMAT 22160 3305 642A 26307 369 650A 26357 | 3307 FORMAT 22167 370 661A 26424 373 678A 26502 | 372 668A 26437 385 691A 26561 375 680A 26504 | 920 FORMAT 22172 380 684A 26527 3038 FORMAT 22203 | 3315 692A 26570 386 703A 26635 3317 696A 26603 400 718A 26721 390 707A 26643 925 FORMAT 22206 | 395 711A 26671 3040 FORMAT 22213 401 737A 27023 | 510 720A 26735 410 756A 27061 402 739A 27026 | 403 741A 27031 405 753A 27057 404 748A 27042 | 411 757A 27063 412 760A 27112 413 762A 27115 414 763A 27117 5102 FORMAT 22216 945 FORMAT 22223 | 420 791A 27260 415 779A 27240 3012 FORMAT 22253 | 1 12/01/64 PAGE 43 | CF1 STORAGE MAP | 0425 816A 27434 950 FORMAT 22257 3014 FORMAT 22275 3016 FORMAT 22303 3018 FORMAT 22312 3020 FORMAT 22322 | 3022 FORMAT 22326 3024 FORMAT 22336 3026 FORMAT 22345 | 3028 FORMAT 22354 3030 FORMAT 22363 438 840A 27712 | 427 839A 27711 475 872A 30042 440 852A 27756 | 445 857A 27773 470 862A 30010 1306 898A 30136 1301 886A 30104 1302 890A 30111 1303 891A 30113 | 1305 896A 30131 1015 FORMAT 22373 1220 915A 30240 | 1219 FORMAT 22423 486 921A 30247 480 919A 30244 | 485 934A 30342 482 923A 30254 483 931A 30304 | 1025 FORMAT 22432 THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 30377.| 1 12/01/64 PAGE 44 | CF2 STORAGE MAP | 0 SUBROUTINE REDATA | - COMMON VARIABLES 0 COMMON BLOCK // ORIGIN 00001 LENGTH 20370 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | YYMIN 00000 R YCMIN 00001 R YYMAX 00002 R | YCMAX 00003 R FMT 00004 R IQ 00124 I | KSL1 00174 I KSL2 00175 I KSL3 00176 I KSL4 00177 I KSL5 00200 I KSL6 00201 I | IDZA 00202 I IDZB 02152 I DELTA 04122 R | YYIN 06072 R YYCC 10042 R LSORT 12012 I | Z 13762 R C 15663 R BETA 17467 R | AVG 17525 R SIGMA 17563 R B 17621 R VAR 17657 R YRESID 17715 R ITRFM 17753 I | VARTR 20035 R AVATR 20117 R IOMIT 20201 I | KTIN 20263 I KTOU 20264 I KTPCH 20265 I | KTBIN1 20266 I KTBIN2 20267 I KTBIN3 20270 I | KTBIN4 20271 I KTBIN5 20272 I PR1 20273 R PR2 20274 R PR3 20275 R BZRO 20276 R | KNOIND 20277 I KNODEP 20300 I KNOVAR 20301 I | KNVP1 20302 I JNOIND 20303 I JNODEP 20304 I | JNOVAR 20305 I NOBSP 20306 I NOERR 20307 I | JOBSV 20310 I NRDLMT 20311 I NOOBSV 20312 I MNO 20313 I JSEQ 20314 I N 20315 I | M 20316 I J 20317 I I 20320 I | LNOVAR 20321 I NNO 20322 I ISEQ 20323 I | L 20324 I KNO 20325 I XNOOBS 20326 R | XNOVAR 20327 R XNOIND 20330 R DEFR 20331 R NOVM1 20332 I RTDEFR 20333 R IPL1 20334 I | K 20335 I RSQD 20336 R BZERO 20337 R | SUMB 20340 R ICURY 20341 I FACTOR 20342 R | SGBETA 20343 R TBETA 20344 R PCORK 20345 R | SIGB 20346 R FRATIO 20347 R UNBEST 20350 R BIEST 20351 R SSQ 20352 R SSE 20353 R | IPLOT1 20354 I IPLOT2 20355 I YCALC 20356 R | DIFF 20357 R IREC 20360 I X 20361 R | IDENT 20362 I IOBSV 20363 I MP 20364 I | NOVAR 20365 I NVP1 20366 I IPATH 20367 I DATA 15663 R STDDEV 16027 R ITLOC 16173 I | NOIND 00124 I NODEP 00125 I KSW45 00151 I | KSW46 00152 I | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE P 20371 R XA 20377 R XC 20411 R | XB 20417 R | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | 1 12/01/64 PAGE 45 CF2 STORAGE MAP | 0IEND 20431 I JTR 20432 I JJTR 20433 I | JIG 20434 I KSWA 20435 I | - ENTRY POINTS | 0 REDATA SECTION 3 - SUBROUTINES CALLED | 0 ATHRUZ SECTION 4 .FVIO. SECTION 5 .FRDD. SECTION 6 | .FSLI. SECTION 7 .FWRD. SECTION 8 EXIT SECTION 9 | .FSLO. SECTION 10 ALOG SECTION 11 EXP SECTION 12 | .FXEM. SECTION 13 .FRTN. SECTION 14 .FCNV. SECTION 15 .FFIL. SECTION 16 SYSLOC SECTION 17 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 20 27A 20677 2 3A 20533 60 8A 20547 | 62 14A 20570 4 FORMAT 20460 10 24A 20640 6 22A 20622 8 FORMAT 20466 12 FORMAT 20473 | 14 26A 20664 16 FORMAT 20500 22 28A 20706 | 30 41A 21031 24 FORMAT 20512 26 37A 20754 | 90 35A 20736 28 FORMAT 20515 99 36A 20753 | 40 48A 21067 42 53A 21170 50 55A 21214 44 FORMAT 20520 52 59A 21246 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 21315.| 1 12/01/64 PAGE 46 | CF3 STORAGE MAP | 0 SUBROUTINE AE04 - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 20370 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | YYMIN 00000 R YCMIN 00001 R YYMAX 00002 R | YCMAX 00003 R FMT 00004 R IQ 00124 I KSL1 00174 I KSL2 00175 I KSL3 00176 I | KSL4 00177 I KSL5 00200 I KSL6 00201 I | IDZA 00202 I IDZB 02152 I DELTA 04122 R | YYIN 06072 R YYCC 10042 R LSORT 12012 I | Z 13762 R C 15663 R BETA 17467 R AVG 17525 R SIGMA 17563 R B 17621 R | VAR 17657 R YRESID 17715 R ITRFM 17753 I | VARTR 20035 R AVATR 20117 R IOMIT 20201 I | KTIN 20263 I KTOU 20264 I KTPCH 20265 I | KTBIN1 20266 I KTBIN2 20267 I KTBIN3 20270 I KTBIN4 20271 I KTBIN5 20272 I PR1 20273 R | PR2 20274 R PR3 20275 R BZRO 20276 R | KNOIND 20277 I KNODEP 20300 I KNOVAR 20301 I | KNVP1 20302 I JNOIND 20303 I JNODEP 20304 I | JNOVAR 20305 I NOBSP 20306 I NOERR 20307 I JOBSV 20310 I NRDLMT 20311 I NOOBSV 20312 I | MNO 20313 I JSEQ 20314 I N 20315 I | M 20316 I J 20317 I I 20320 I | LNOVAR 20321 I NNO 20322 I ISEQ 20323 I | L 20324 I KNO 20325 I XNOOBS 20326 R XNOVAR 20327 R XNOIND 20330 R DEFR 20331 R | NOVM1 20332 I RTDEFR 20333 R IPL1 20334 I | K 20335 I RSQD 20336 R BZERO 20337 R | SUMB 20340 R ICURY 20341 I FACTOR 20342 R | SGBETA 20343 R TBETA 20344 R PCORK 20345 R SIGB 20346 R FRATIO 20347 R UNBEST 20350 R | BIEST 20351 R SSQ 20352 R SSE 20353 R | IPLOT1 20354 I IPLOT2 20355 I YCALC 20356 R | DIFF 20357 R IREC 20360 I X 20361 R | IDENT 20362 I IOBSV 20363 I MP 20364 I NOVAR 20365 I NVP1 20366 I IPATH 20367 I | DATA 15663 R STDDEV 16027 R ITLOC 16173 I | NOIND 00124 I NODEP 00125 I | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE VX 20371 R VY 20372 R VALUE 20373 R | MA 20374 I NA 20375 I | - ENTRY POINTS | 0 AE04 SECTION 3 | 1 12/01/64 PAGE 47 CF3 STORAGE MAP | 0 SUBROUTINES CALLED | 0 ALOG SECTION 4 ALOG10 SECTION 5 EXP SECTION 6 | .XP3. SECTION 7 SIN SECTION 8 COS SECTION 9 | .FXEM. SECTION 10 E.1 SECTION 11 E.2 SECTION 12 E.3 SECTION 13 E.4 SECTION 14 CC.1 SECTION 15 | CC.2 SECTION 16 CC.3 SECTION 17 CC.4 SECTION 18 | SYSLOC SECTION 19 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 65 92A 21070 50 13A 20443 60 82A 21051 | 1 14A 20474 2 17A 20502 3 20A 20510 | 4 23A 20516 5 26A 20544 6 29A 20553 | 7 31A 20557 8 33A 20563 9 35A 20567 | 13 55A 20706 14 58A 20714 15 61A 20722 16 63A 20726 17 67A 20741 18 71A 20754 | 19 74A 20765 20 78A 21000 91 39A 20636 | 10 43A 20650 11 47A 20662 12 51A 20674 | 62 88A 21061 115 106A 21117 125 119A 21144 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 21202.| 1 12/01/64 PAGE 48 | CF4 STORAGE MAP | 0 SUBROUTINE AE02 | - COMMON VARIABLES 0 COMMON BLOCK // ORIGIN 00001 LENGTH 20370 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | YYMIN 00000 R YCMIN 00001 R YYMAX 00002 R | YCMAX 00003 R FMT 00004 R IQ 00124 I | KSL1 00174 I KSL2 00175 I KSL3 00176 I KSL4 00177 I KSL5 00200 I KSL6 00201 I | IDZA 00202 I IDZB 02152 I DELTA 04122 R | YYIN 06072 R YYCC 10042 R LSORT 12012 I | Z 13762 R C 15663 R BETA 17467 R | AVG 17525 R SIGMA 17563 R B 17621 R VAR 17657 R YRESID 17715 R ITRFM 17753 I | VARTR 20035 R AVATR 20117 R IOMIT 20201 I | KTIN 20263 I KTOU 20264 I KTPCH 20265 I | KTBIN1 20266 I KTBIN2 20267 I KTBIN3 20270 I | KTBIN4 20271 I KTBIN5 20272 I PR1 20273 R PR2 20274 R PR3 20275 R BZRO 20276 R | KNOIND 20277 I KNODEP 20300 I KNOVAR 20301 I | KNVP1 20302 I JNOIND 20303 I JNODEP 20304 I | JNOVAR 20305 I NOBSP 20306 I NOERR 20307 I | JOBSV 20310 I NRDLMT 20311 I NOOBSV 20312 I MNO 20313 I JSEQ 20314 I N 20315 I | M 20316 I J 20317 I I 20320 I | LNOVAR 20321 I NNO 20322 I ISEQ 20323 I | L 20324 I KNO 20325 I XNOOBS 20326 R | XNOVAR 20327 R XNOIND 20330 R DEFR 20331 R NOVM1 20332 I RTDEFR 20333 R IPL1 20334 I | K 20335 I RSQD 20336 R BZERO 20337 R | SUMB 20340 R ICURY 20341 I FACTOR 20342 R | SGBETA 20343 R TBETA 20344 R PCORK 20345 R | SIGB 20346 R FRATIO 20347 R UNBEST 20350 R BIEST 20351 R SSQ 20352 R SSE 20353 R | IPLOT1 20354 I IPLOT2 20355 I YCALC 20356 R | DIFF 20357 R IREC 20360 I X 20361 R | IDENT 20362 I IOBSV 20363 I MP 20364 I | NOVAR 20365 I NVP1 20366 I IPATH 20367 I DATA 15663 R STDDEV 16027 R ITLOC 16173 I | NOIND 00124 I NODEP 00125 I | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | PIVOT 20371 R EMPYER 20372 R - ENTRY POINTS | 0 AE02 SECTION 3 | 1 12/01/64 PAGE 49 | CF4 STORAGE MAP | 0 SUBROUTINES CALLED 0 SYSLOC SECTION 4 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 15 13A 20435 10 11A 20434 5 8A 20431 | 45 52A 20625 20 23A 20507 25 29A 20526 30 36A 20565 35 40A 20577 40 47A 20615 | | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 20666.| 1 12/01/64 PAGE 50 | CF5 STORAGE MAP 0 SUBROUTINE MINMAX | - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 20370 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | YYMIN 00000 R YCMIN 00001 R YYMAX 00002 R YCMAX 00003 R FMT 00004 R IQ 00124 I | KSL1 00174 I KSL2 00175 I KSL3 00176 I | KSL4 00177 I KSL5 00200 I KSL6 00201 I | IDZA 00202 I IDZB 02152 I DELTA 04122 R | YYIN 06072 R YYCC 10042 R LSORT 12012 I Z 13762 R C 15663 R BETA 17467 R | AVG 17525 R SIGMA 17563 R B 17621 R | VAR 17657 R YRESID 17715 R ITRFM 17753 I | VARTR 20035 R AVATR 20117 R IOMIT 20201 I | KTIN 20263 I KTOU 20264 I KTPCH 20265 I KTBIN1 20266 I KTBIN2 20267 I KTBIN3 20270 I | KTBIN4 20271 I KTBIN5 20272 I PR1 20273 R | PR2 20274 R PR3 20275 R BZRO 20276 R | KNOIND 20277 I KNODEP 20300 I KNOVAR 20301 I | KNVP1 20302 I JNOIND 20303 I JNODEP 20304 I JNOVAR 20305 I NOBSP 20306 I NOERR 20307 I | JOBSV 20310 I NRDLMT 20311 I NOOBSV 20312 I | MNO 20313 I JSEQ 20314 I N 20315 I | M 20316 I J 20317 I I 20320 I | LNOVAR 20321 I NNO 20322 I ISEQ 20323 I L 20324 I KNO 20325 I XNOOBS 20326 R | XNOVAR 20327 R XNOIND 20330 R DEFR 20331 R | NOVM1 20332 I RTDEFR 20333 R IPL1 20334 I | K 20335 I RSQD 20336 R BZERO 20337 R | SUMB 20340 R ICURY 20341 I FACTOR 20342 R SGBETA 20343 R TBETA 20344 R PCORK 20345 R | SIGB 20346 R FRATIO 20347 R UNBEST 20350 R | BIEST 20351 R SSQ 20352 R SSE 20353 R | IPLOT1 20354 I IPLOT2 20355 I YCALC 20356 R | DIFF 20357 R IREC 20360 I X 20361 R IDENT 20362 I IOBSV 20363 I MP 20364 I | NOVAR 20365 I NVP1 20366 I IPATH 20367 I | DATA 15663 R STDDEV 16027 R ITLOC 16173 I | NOIND 00124 I NODEP 00125 I | - ENTRY POINTS 0 MINMAX SECTION 3 | - SUBROUTINES CALLED | 0 SYSLOC SECTION 4 | - EFN IFN CORRESPONDENCE | 1 12/01/64 PAGE 51 CF5 STORAGE MAP | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 2 3A 20401 3 4A 20403 5 7A 20411 | 4 6A 20407 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 20441. 1 12/01/64 PAGE 52 | CF6 STORAGE MAP | 0 SUBROUTINE YBACK | - COMMON VARIABLES | 0 COMMON BLOCK // ORIGIN 00001 LENGTH 20370 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | YYMIN 00000 R YCMIN 00001 R YYMAX 00002 R | YCMAX 00003 R FMT 00004 R IQ 00124 I | KSL1 00174 I KSL2 00175 I KSL3 00176 I | KSL4 00177 I KSL5 00200 I KSL6 00201 I IDZA 00202 I IDZB 02152 I DELTA 04122 R | YYIN 06072 R YYCC 10042 R LSORT 12012 I | Z 13762 R C 15663 R BETA 17467 R | AVG 17525 R SIGMA 17563 R B 17621 R | VAR 17657 R YRESID 17715 R ITRFM 17753 I VARTR 20035 R AVATR 20117 R IOMIT 20201 I | KTIN 20263 I KTOU 20264 I KTPCH 20265 I | KTBIN1 20266 I KTBIN2 20267 I KTBIN3 20270 I | KTBIN4 20271 I KTBIN5 20272 I PR1 20273 R | PR2 20274 R PR3 20275 R BZRO 20276 R KNOIND 20277 I KNODEP 20300 I KNOVAR 20301 I | KNVP1 20302 I JNOIND 20303 I JNODEP 20304 I | JNOVAR 20305 I NOBSP 20306 I NOERR 20307 I | JOBSV 20310 I NRDLMT 20311 I NOOBSV 20312 I | MNO 20313 I JSEQ 20314 I N 20315 I M 20316 I J 20317 I I 20320 I | LNOVAR 20321 I NNO 20322 I ISEQ 20323 I | L 20324 I KNO 20325 I XNOOBS 20326 R | XNOVAR 20327 R XNOIND 20330 R DEFR 20331 R | NOVM1 20332 I RTDEFR 20333 R IPL1 20334 I K 20335 I RSQD 20336 R BZERO 20337 R | SUMB 20340 R ICURY 20341 I FACTOR 20342 R | SGBETA 20343 R TBETA 20344 R PCORK 20345 R | SIGB 20346 R FRATIO 20347 R UNBEST 20350 R | BIEST 20351 R SSQ 20352 R SSE 20353 R IPLOT1 20354 I IPLOT2 20355 I YCALC 20356 R | DIFF 20357 R IREC 20360 I X 20361 R | IDENT 20362 I IOBSV 20363 I MP 20364 I | NOVAR 20365 I NVP1 20366 I IPATH 20367 I | DATA 15663 R STDDEV 16027 R ITLOC 16173 I NOIND 00124 I NODEP 00125 I | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | LBJ 20371 I KV 20372 I VK 20373 R | KBJ 20374 I - ENTRY POINTS | 0 YBACK SECTION 3 | 1 12/01/64 PAGE 53 | CF6 STORAGE MAP | 0 SUBROUTINES CALLED 0 EXP SECTION 4 ALOG SECTION 5 ALOG10 SECTION 6 | ARSIN SECTION 7 ARCOS SECTION 8 .FXEM. SECTION 9 | E.1 SECTION 10 E.2 SECTION 11 E.3 SECTION 12 | E.4 SECTION 13 SYSLOC SECTION 14 | - EFN IFN CORRESPONDENCE 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 9 40A 20613 10 41A 20614 1 15A 20460 | 2 18A 20466 3 21A 20514 4 24A 20522 | 5 27A 20530 6 32A 20553 7 35A 20566 | 8 38A 20601 51 28A 20537 61 33A 20562 71 36A 20575 81 39A 20610 13 45A 20641 | 14 48A 20647 15 51A 20655 16 54A 20670 | 17 58A 20703 18 61A 20731 19 65A 20767 | 20 68A 20775 151 52A 20664 181 62A 20740 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 21073.| 1 12/01/64 PAGE 54 | CF7 STORAGE MAP | 0 SUBROUTINE PITCHA | - COMMON VARIABLES 0 COMMON BLOCK // ORIGIN 00001 LENGTH 20370 | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | YYMIN 00000 R YCMIN 00001 R YYMAX 00002 R | YCMAX 00003 R FMT 00004 R IQ 00124 I | KSL1 00174 I KSL2 00175 I KSL3 00176 I KSL4 00177 I KSL5 00200 I KSL6 00201 I | IDZA 00202 I IDZB 02152 I DELTA 04122 R | YYIN 06072 R YYCC 10042 R LSORT 12012 I | Z 13762 R C 15663 R BETA 17467 R | AVG 17525 R SIGMA 17563 R B 17621 R VAR 17657 R YRESID 17715 R ITRFM 17753 I | VARTR 20035 R AVATR 20117 R IOMIT 20201 I | KTIN 20263 I KTOU 20264 I KTPCH 20265 I | KTBIN1 20266 I KTBIN2 20267 I KTBIN3 20270 I | KTBIN4 20271 I KTBIN5 20272 I PR1 20273 R PR2 20274 R PR3 20275 R BZRO 20276 R | KNOIND 20277 I KNODEP 20300 I KNOVAR 20301 I | KNVP1 20302 I JNOIND 20303 I JNODEP 20304 I | JNOVAR 20305 I NOBSP 20306 I NOERR 20307 I | JOBSV 20310 I NRDLMT 20311 I NOOBSV 20312 I MNO 20313 I JSEQ 20314 I N 20315 I | M 20316 I J 20317 I I 20320 I | LNOVAR 20321 I NNO 20322 I ISEQ 20323 I | L 20324 I KNO 20325 I XNOOBS 20326 R | XNOVAR 20327 R XNOIND 20330 R DEFR 20331 R NOVM1 20332 I RTDEFR 20333 R IPL1 20334 I | K 20335 I RSQD 20336 R BZERO 20337 R | SUMB 20340 R ICURY 20341 I FACTOR 20342 R | SGBETA 20343 R TBETA 20344 R PCORK 20345 R | SIGB 20346 R FRATIO 20347 R UNBEST 20350 R BIEST 20351 R SSQ 20352 R SSE 20353 R | IPLOT1 20354 I IPLOT2 20355 I YCALC 20356 R | DIFF 20357 R IREC 20360 I X 20361 R | IDENT 20362 I IOBSV 20363 I MP 20364 I | NOVAR 20365 I NVP1 20366 I IPATH 20367 I DATA 15663 R STDDEV 16027 R ITLOC 16173 I | NOIND 00124 I NODEP 00125 I | - DIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | GRIDA 20371 R GRIDB 22176 R NSW 24003 I PRAX 24011 R XAX 24074 R RESID 24107 R | POS 24120 R NEG 24130 I | - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | 1 12/01/64 PAGE 55 CF7 STORAGE MAP | 0XNP 24140 R JJC 24141 I DELP 24142 R | IIC 24143 I YIW 24144 R XIW 24145 R | IA 24146 I JAC 24147 I L1 24150 I | LINE 24151 I LOCX1 24152 I LOCX2 24153 I LCHAR 24154 I LWORD 24155 I NCHAR 24156 I | PRLOC 24157 R IB 24160 I LOCP1 24161 I | IC 24162 I LOCP2 24163 I | - ENTRY POINTS | 0 PITCHA SECTION 3 - SUBROUTINES CALLED | 0 ATHRUZ SECTION 4 GRID SECTION 5 PACK SECTION 6 | .FVIO. SECTION 7 .FWRD. SECTION 8 .FFIL. SECTION 9 | .FCNV. SECTION 10 E.1 SECTION 11 E.2 SECTION 12 | E.3 SECTION 13 E.4 SECTION 14 CC.1 SECTION 15 CC.2 SECTION 16 CC.3 SECTION 17 CC.4 SECTION 18 | SYSLOC SECTION 19 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION | 200 30A 24760 202 33A 24762 30 96A 25260 10 64A 25131 4 59A 25110 2 57A 25103 | 12 66A 25147 18 75A 25167 14 71A 25160 | 16 73A 25164 24 84A 25214 20 81A 25207 | 22 83A 25211 28 91A 25237 26 89A 25232 | 34 120A 25347 102 FORMAT 24337 104 FORMAT 24360 106 FORMAT 24361 40 130A 25432 108 FORMAT 24405 | 42 140A 25472 44 FORMAT 24410 46 164A 25615 | 300 159A 25567 302 FORMAT 24414 48 186A 25722 | 304 180A 25667 110 FORMAT 24420 306 FORMAT 24423 | 308 205A 26012 310 200A 25764 312 211A 26024 314 222A 26066 112 FORMAT 24427 122 FORMAT 24437 | 126 FORMAT 24451 50 243A 26227 52 253A 26267 | 54 277A 26412 56 272A 26364 58 299A 26517 | 60 293A 26464 62 318A 26607 64 313A 26561 | 66 324A 26621 68 335A 26663 132 FORMAT 24454 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 26767.| 1 12/01/64 PAGE 56 | CF8 STORAGE MAP | 0 SUBROUTINE GRID - UNDIMENSIONED PROGRAM VARIABLES | 0SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE SYMBOL LOCATION TYPE | BLANK 00001 R BLINK 00002 R IZ 00003 I | HOR1 00004 R HOR2 00005 R | - ENTRY POINTS 0 GRID SECTION 2 | - SUBROUTINES CALLED | 0 ATHRUZ SECTION 3 SYSLOC SECTION 4 | - EFN IFN CORRESPONDENCE | 0EFN IFN LOCATION EFN IFN LOCATION EFN IFN LOCATION 42 15A 00111 40 12A 00103 10 34A 00156 | 12 54A 00250 14 68A 00272 30 78A 00312 | 19 87A 00334 20 93A 00337 | THE FIRST LOCATION NOT USED BY THIS PROGRAM IS 00541. 1 12/01/64 PAGE 57 | | 0$IBMAP CF9 LIST,DECK,RELMOD,REF $IBMAP CF9 0000 1 CF9 12/01/64 PAGE 58 | 7090 RELMOD ASSEMBLY. | 0| $IBLDR CF9 12/01/64 CF9 0000| 1 CF9 12/01/64 PAGE 59 | ASSEMBLED TEXT. | 0 $TEXT CF9 CF9 0001| | * CF9 0010 * SUBROUTINE PLACES THE APPROPRIATE SYMBOLS ON THE PLOTS CF9 0020| * CF9 0030| 00000 ENTRY PACK FORTRAN CALL PACK(X, N, 1HC) PACKA 01 CF9 0040| 00000 PACK SAVE 1 | 0BINARY CARD ID. CF9 0002 00000 1 00000 0 00004 10001 | 00001 0774 00 1 00000 10000 | 00002 0774 00 4 00000 10000 | 00003 0020 00 4 00001 10000 | 00004 0634 00 4 04000 10011 00005 0634 00 4 00041 10001 | 00006 0634 00 4 00002 10001 | 00007 0634 00 1 00001 10001 | 00010 4500 60 4 00005 10000 CAL* 5,4 CHAR | 00011 4320 00 0 00043 10001 ANA =77 CLEAN UP 00012 4130 00 0 00000 10000 XCL TO MQ | 00013 0500 60 4 00004 10000 CLA* 4,4 LOAD COUNT | 00014 0734 00 1 00000 10000 PAX 0,1 TO XRA | 00015 0500 00 1 00033 10001 CLA SCNT+1,1 LOAD SHIFT COUNT | 00016 0621 00 0 01001 10011 STA *+1 MODIFY SHIFT COUNT 00017 4773 00 0 00000 10000 RQL ** SHIFT CHAR INTO PLACE | 00020 4500 00 1 00041 10001 CAL MSK+1,1 LOAD MASK | 00021 0320 60 4 00003 10000 ANS* 3,4 MASK OUT CHAR. | 00022 4130 00 0 00000 10000 XCL CHAR TO AC | 0BINARY CARD ID. CF9 0003 00023 4602 60 4 00003 10000 ORS* 3,4 PUT CHAR IN PLACE | 00024 RETURN PACK | 00025 0 00000 0 00000 10000 PZE 0*6 6 | 00026 0 00000 0 00006 10000 PZE 1*6 5 | 00027 0 00000 0 00014 10000 PZE 2*6 4 00030 0 00000 0 00022 10000 PZE 3*6 3 | 00031 0 00000 0 00030 10000 PZE 4*6 2 | 00032 0 00000 0 00036 10000 SCNT PZE 5*6 1 | 00033 777777777700 10000 OCT 777777777700 6 | 00034 777777770077 10000 OCT 777777770077 5 00035 777777007777 10000 OCT 777777007777 4 | 00036 777700777777 10000 OCT 777700777777 3 | 00037 770077777777 10000 OCT 770077777777 2 | 00040 007777777777 10000 MSK OCT 007777777777 1 | 00041 000000000000 10000 *LDIR 00042 232611606060 10000 | 00043 000000000115 10000 *LORG| 00000 01111 END | 1 CF9 12/01/64 PAGE 60 | CONTROL DICTIONARY 0 $CDICT CF9 CF9 0004| | 0BINARY CARD ID. CF9 0005| 000044000000 PREFACE START=0,LENGTH=36,TYPE=7090,CMPLX=5| 000000000005 232611606060 CF9 DECK | 000044000000 | 472123426060 PACK REAL SECT. 2,LOC=0,LENGTH=0| 000000000000 | 472123426060 PACK REAL SECT. 3,LOC=0,LENGTH=0 000000000000 | 627062434623 SYSLOC VIRTUAL SECT. 4| 200000000000 | - $DKEND CF9 CF9 0006| -NO MESSAGES FOR THIS ASSEMBLY | 1 CF9 12/01/64 PAGE 61 | SYMBOL REFERENCE DATA | 0| REFERENCES TO DEFINED SYMBOLS. 0 CLASS SYMBOL VALUE REFERENCES| | ..0001 00002 6,7| ..0002 00003 | ..0003 00004 0 MSK 00040 20| PACK 00000 0,0,24| LCTR BLCTR | QUAL UNQS | LCTR // SCNT 00032 15| -REFERENCES TO VIRTUAL SYMBOLS.| | SYSLOC 4 4 1 12/01/64 PAGE 62 | | 0$IBMAP CF10 LIST,DECK,RELMOD,REF $IBMAP CF9 0000 1 CF10 12/01/64 PAGE 63 | 7090 RELMOD ASSEMBLY. | 0| $IBLDR CF10 12/01/64 CF100000| 1 CF10 12/01/64 PAGE 64 | ASSEMBLED TEXT. | 0 $TEXT CF10 CF100001| | * IBMAP ROUTINE FOR PROGRAMS SIFTED THRU SIFT 40 BOOL0100 00000 ENTRY ATHRUZ BOOL0110| BOOL0120| ALPHAMERIC CONVERSION BOOL0130| ***CALL ATHRUZ(Q000HL,6H EOF )*** BOOL0140| BOOL0150 00000 ATHRUZ SAVE 4 BOOL0160| 0BINARY CARD ID. CF100002| 00000 1 00000 0 00003 10001 | 00001 0774 00 4 00000 10000 | 00002 0020 00 4 00001 10000 00003 0634 00 4 04000 10011 | 00004 0634 00 4 00011 10001 | 00005 0634 00 4 00001 10001 | 00006 4500 60 4 00004 10000 CAL* 4,4 GET ALPHAMERIC INFO BOOL0170| 00007 0602 60 4 00003 10000 SLW* 3,4 AND STORE IT BOOL0180 00010 RETURN ATHRUZ BOOL0190| 00011 000000000000 10000 *LDIR| 00012 232601006060 10000 | 00000 01111 END BOOL0200| 1 CF10 12/01/64 PAGE 65 CONTROL DICTIONARY | 0 $CDICT CF10 CF100003| | 0BINARY CARD ID. CF100004| 000013000000 PREFACE START=0,LENGTH=11,TYPE=7090,CMPLX=5 000000000005 | 232601006060 CF10 DECK | 000013000000 | 216330516471 ATHRUZ REAL SECT. 2,LOC=0,LENGTH=0| 000000000000 216330516471 ATHRUZ REAL SECT. 3,LOC=0,LENGTH=0| 000000000000 | 627062434623 SYSLOC VIRTUAL SECT. 4| 200000000000 | - $DKEND CF10 CF100005 | -NO MESSAGES FOR THIS ASSEMBLY | 1 CF10 12/01/64 PAGE 66 | SYMBOL REFERENCE DATA | 0 REFERENCES TO DEFINED SYMBOLS.| 0 CLASS SYMBOL VALUE REFERENCES| | ATHRUZ 00000 0,0,10| ..0001 00001 5 ..0002 00002 | ..0003 00003 0| LCTR BLCTR | QUAL UNQS | LCTR // -REFERENCES TO VIRTUAL SYMBOLS.| | SYSLOC 4 3 1 12/01/64 PAGE 67 | | 0$IBMAP CF11 BOOL0240 1 CF11 12/01/64 PAGE 68 | 7090 RELMOD ASSEMBLY. | 0| $IBLDR CF11 12/01/64 CF110000| 1 CF11 12/01/64 PAGE 69 | ASSEMBLED TEXT. | 0 $TEXT CF11 CF110001| | * IBMAP ROUTINE FOR PROGRAMS SIFTED THRU SIFT 40 BOOL0250 00000 ENTRY BOOL ZERO VS NON'ZERO TEST BOOL0260| BOOL0270| ZERO VS NON'ZERO TEST BOOL0280| ***BOOL(ARG.)*** BOOL0290| BOOL0300 00000 BOOL SAVE 4 BOOL0310| 0BINARY CARD ID. CF110002| 00000 1 00000 0 00003 10001 | 00001 0774 00 4 00000 10000 | 00002 0020 00 4 00001 10000 00003 0634 00 4 04000 10011 | 00004 0634 00 4 00010 10001 | 00005 0634 00 4 00001 10001 | 00006 4500 60 4 00003 10000 CAL* 3,4 BOOL0320| 00007 RETURN BOOL BOOL0330 00010 000000000000 10000 *LDIR| 00011 232601016060 10000 | 00000 01111 END BOOL0340| 1 CF11 12/01/64 PAGE 70 | CONTROL DICTIONARY 0 $CDICT CF11 CF110003| | 0BINARY CARD ID. CF110004| 000012000000 PREFACE START=0,LENGTH=10,TYPE=7090,CMPLX=5| 000000000005 232601016060 CF11 DECK | 000012000000 | 224646436060 BOOL REAL SECT. 2,LOC=0,LENGTH=0| 000000000000 | 224646436060 BOOL REAL SECT. 3,LOC=0,LENGTH=0 000000000000 | 627062434623 SYSLOC VIRTUAL SECT. 4| 200000000000 | - $DKEND CF11 CF110005| -NO MESSAGES FOR THIS ASSEMBLY | 1 CF11 12/01/64 PAGE 71 | SYMBOL REFERENCE DATA | 0| REFERENCES TO DEFINED SYMBOLS. 0 CLASS SYMBOL VALUE REFERENCES| | BOOL 00000 0,0,7| ..0001 00001 5| ..0002 00002 ..0003 00003 0| LCTR BLCTR | QUAL UNQS | LCTR // | -REFERENCES TO VIRTUAL SYMBOLS. | SYSLOC 4 3 1 12/01/64 PAGE 72 | | 0$IBMAP HLCT BOOL0380 1 HLCT 12/01/64 PAGE 73 | 7090 RELMOD ASSEMBLY. | 0| $IBLDR HLCT 12/01/64 HLCT0000| 1 HLCT 12/01/64 PAGE 74 | ASSEMBLED TEXT. | 0 $TEXT HLCT HLCT0001| | * IBMAP ROUTINE FOR PROGRAMS SIFTED THRU SIFT 40 BOOL0390 00000 ENTRY HOLOCT OCTAL CONVERSION BOOL0400| BOOL0410| OCTAL CONVERSION BOOL0420| ***CALL HOLOCT(Q000CT,12H77777 )*** BOOL0430| BOOL0440 00000 HOLOCT SAVE 1,2 BOOL0450| 0BINARY CARD ID. HLCT0002| 00000 1 00000 0 00005 10001 | 00001 0774 00 2 00000 10000 | 00002 0774 00 1 00000 10000 00003 0774 00 4 00000 10000 | 00004 0020 00 4 00001 10000 | 00005 0634 00 4 04000 10011 | 00006 0634 00 4 00040 10001 | 00007 0634 00 4 00003 10001 00010 0634 00 1 00002 10001 | 00011 0634 00 2 00001 10001 | 00012 0500 00 4 00004 10000 CLA 4,4 ADDRESS OF OCTAL NUMBER BOOL0460| 00013 0621 00 0 00017 10001 STA LOAD BOOL0470| 00014 0600 60 4 00003 10000 STZ* 3,4 BOOL0480 00015 0774 00 2 00000 10000 AXT 0,2 BOOL0490| 00016 0774 00 1 00006 10000 AXT 6,1 SIX CHARACTERS BOOL0500| 00017 0560 00 2 00000 10000 LOAD LDQ **,2 BOOL0510| 00020 4754 00 0 00000 10000 ZAC BOOL0520| 00021 4763 00 0 00006 10000 LGL 6 BOOL0530 00022 4340 00 0 00042 10001 LAS =O60 IS IT A BLANK BOOL0540| 00023 RET RETURN HOLOCT TERM. CONV., COULD CHECK FOR ERROR BOOL0550| 0BINARY CARD ID. HLCT0003| 00024 0020 00 0 00035 10001 TRA SKIP A BLANK. SKIP TO NEXT NUMBER BOOL0560| 00025 4340 00 0 00043 10001 LAS =8 IS IT AN OCTAL NUMBER BOOL0570 00026 0020 00 0 00023 10001 TRA ERROR NO, ILLEGAL CHARACTER BOOL0580| 00027 0020 00 0 00023 10001 TRA ERROR SAME BOOL0590| 00030 0602 00 0 00037 10001 SLW TEMP YES. PROCESS IT BOOL0600| 00031 4500 60 4 00003 10000 CAL* 3,4 GET THE OCTAL BEING FORMED BOOL0610| 00032 0767 00 0 00003 10000 ALS 3 MAKE ROOM FOR NEXT NUMBER BOOL0620 00033 4501 00 0 00037 10001 ORA TEMP ADD IT BOOL0630| 00034 0602 60 4 00003 10000 SLW* 3,4 AND STORE BACK BOOL0640| 00035 2 00001 1 00020 10001 SKIP TIX LOAD+1,1,1 SKIP TO THE NEXT NUMBER BOOL0650| 00036 1 77777 2 00016 10001 TXI LOAD-1,2,-1 GET NEXT WORD BOOL0660| 00037 0 00000 0 00000 10000 TEMP PZE TEMPORARY STORAGE BOOL0670 00023 ERROR SYN RET BOOL0680| 00040 000000000000 10000 *LDIR| 00041 304323636060 10000 | 00042 000000000060 10000 *LORG| 00043 000000000010 10000 00000 01111 END BOOL0690| 1 HLCT 12/01/64 PAGE 75 | CONTROL DICTIONARY | 0 $CDICT HLCT HLCT0004| 0BINARY CARD ID. HLCT0005| 000044000000 PREFACE START=0,LENGTH=36,TYPE=7090,CMPLX=5| 000000000005 | 304323636060 HLCT DECK | 000044000000 304643462363 HOLOCT REAL SECT. 2,LOC=0,LENGTH=0| 000000000000 | 304643462363 HOLOCT REAL SECT. 3,LOC=0,LENGTH=0| 000000000000 | 627062434623 SYSLOC VIRTUAL SECT. 4 200000000000 | - $DKEND HLCT HLCT0006| | -NO MESSAGES FOR THIS ASSEMBLY | 1 HLCT 12/01/64 PAGE 76 SYMBOL REFERENCE DATA | 0| REFERENCES TO DEFINED SYMBOLS.| 0 CLASS SYMBOL VALUE REFERENCES| ERROR 00023 26,27| HOLOCT 00000 0,0,23| ..0001 00003 7,10,11| ..0002 00004 | ..0003 00005 0 LOAD 00017 13,35,36| RET 00023 40| LCTR BLCTR | QUAL UNQS | LCTR // SKIP 00035 24| TEMP 00037 30,33| -REFERENCES TO VIRTUAL SYMBOLS.| | SYSLOC 4 5 1 12/01/64 PAGE 77 | | 0$ENTRY | $DATA 0LEVEL=2 UNRECOGNIZABLE PARAMETER ON $IBJOB CD... DECK IS IGNORED.| 1 IBLDR 12/01/64 PAGE 78 | * MEMORY MAP * | 0 | 0SYSTEM 00000 THRU 02717 FILE BLOCK ORIGIN 02720| FILES 1. UNIT01 | 2. UNIT02 | 3. UNIT03 | 4. UNIT04 5. UNIT05 | 6. UNIT06 | 7. UNIT07 | 8. UNIT08 | FILE LIST ORIGIN 03060 PRE-EXECUTION INITIALIZATION 03100| CALL ON OBJECT PROGRAM 03137| OBJECT PROGRAM 03144 THRU 42110| 0 DECK ORIGIN CONTROL SECTIONS (/NAME/=NON 0 LENGTH, (LOC)=DELETED, 'LOC'=MOVED, *=NOT REFERENCED)| 1. CF1 03144 /// /(57410) ...... 13135 * | 2. CF2 13153 /// /(57410) REDATA 14064 | 3. CF3 14100 /// /(57410) AE04 14676 | 4. CF4 14712 /// /(57410) AE02 15161 | 5. CF5 15210 /// /(57410) MINMAX 15234 6. CF6 15261 /// /(57410) YBACK 15703 | 7. CF7 15764 /// /(57410) PITCHA 24345 | 8. CF8 24363 GRID 24731 | 9. CF9 25124 PACK 25124 PACK (25124) | 10. CF10 25170 ATHRUZ 25170 ATHRUZ (25170) 11. CF11 25203 BOOL 25203 BOOL (25203) | 12. HLCT 25215 HOLOCT 25215 HOLOCT (25215) | 13. .LXCON 25261 .LXSTR 25261 * .LXSTP 25264 .LXOUT 25332 * .LXERR 25341 .LXCAL 25344 * | .LXRTN 25344 IBEXIT 25344 * .DBCLS 25526 * .LXARG 25675 * .L0 25720 * | .CLSE 25726 .LFBL 25727 * .LUNB 25730 .DFOUT 25731 14. .IODEF 25735 .DEFIN 25735 .ATTAC 25741 * .CLOSE 25743 .OPEN 25745 .READ 25747 | .WRITE 25751 .BSR 25761 * .READR 25771 .RELES 25773 * .LAREA 26004 | .LFBLK 26022 .LTSX 26025 * .RLHLD 26035 * .AREA1 26040 .LUNBL 26046 | .ENTRY 26052 .GOA 26107 .GO 26113 .DERR 26127 .NOPXI 26130 | .COMXI 26132 .EX34 26154 15. .LXSL 26161 .LXSLB 26161 * .LXSEL 26162 .LXSL1 26163 .LXTST 26166 * .LXOVL 26231 * | .LXMOD 26273 * .LXIND 26317 * .LXDIS 26322 * .LXFLG 26323 * .LTCH 26324 | 16. .FPTRP 26333 .FFPT. 26333 * .FPOUT 26462 * .FPARG 26470 * /.COUNT/ 26472 * OVFLOW 26536 * | 17. .ERAS. 26543 E.1 26543 E.2 26544 E.3 26545 E.4 26546 | 18. .XCC. 26547 CC.1 26547 CC.2 26550 CC.3 26551 CC.4 26552 19. XIT 26553 EXIT 26553 .EXIT. 26553 * | 20. FXEM 26554 .FXEM. 26554 TRACE 26562 * /.OPTW./ 27117 * | 21. FCOM 27133 .FCOM. 27133 /FORCOM/ 27236 * | 22. FRDB 27250 .FRDB. 27250 .FWRB. 27252 | 23. FRDD 27302 .FRDD. 27302 .FWRD. 27305 24. FRWT 27417 .FRWT. 27417 | 25. FIOS 27505 .FIOS. 27505 .FSEL. 27674 .FRTB. 27711 .FRTD. 27726 /.FBF1./ 30160 * | 26. FIOB 30164 .FIOB. 30164 .FBLT. 30226 .FBDT. 30254 .FRLR. 30325 .FWLR. 30325 | 27. FIOH 30621 .FIOH. 30621 .FCNV. 31400 .FFIL. 31475 .FRTN. 31475 | 28. FCNV 31614 .FCON. 31614 29. FVIO 34702 .FVIO. 34702 ..FVIO 34702 * | 30. FIOT 34705 .FIOT. 34705 .FIOE. 34776 | 1 IBLDR 12/01/64 PAGE 79 | * MEMORY MAP * | 0 31. UN01 35024 .UN01. 35024 32. UN02 35025 .UN02. 35025 | 33. UN03 35026 .UN03. 35026 | 34. UN04 35027 .UN04. 35027 | 35. UN05 35030 .UN05. 35030 | 36. UN06 35031 .UN06. 35031 ..UN06 35031 * 37. UN07 35032 .UN07. 35032 | 38. UN08 35033 .UN08. 35033 | 39. FLOG 35034 ALOG10 35034 ALOG 35035 | 40. FXPF 35240 EXP 35240 | 41. FSCN 35361 COS 35361 SIN 35362 42. FSQR 35555 SQRT 35555 | 43. FXP3 35630 .XP3. 35630 | 44. FSLDI 35755 .FSLI. 35773 .FSDI. 36001 * | 45. FSLI 36012 .SLI. 36012 .SLI1. 36017 .SDI. 36025 .SDI1. 36033 | 46. FSLDO 36051 .FSLO. 36067 .FSDO. 36075 * 47. FSLO 36106 .SLO. 36106 .SLO2. 36114 .SDO. 36121 .SDO2. 36133 | 48. FASC 36145 ARCOS 36145 ARSIN 36146 | 49. .IOCS 36276 .L(0) 36276 .MONSW 36316 .TEOR 36365 .DEFI. 36445 .JOINX 36511 * | .CLOS. 36530 .ATTC. 36543 .SH1 36755 * .SH9 37017 * .OPEN. 37040 | .OP4 37066 * .OP7 37117 * .OP9.2 37133 * .RLSE. 37205 .RER2. 37205 .READ. 37206 .RER1. 37231 .WRIT. 37233 .MNT1A 37423 * .EOFEX 37504 * | .FEEIT 37554 .GTIOX 37575 .RW7 37713 * .RE7 40332 * .ENDTR 40773 | .SEL59 40775 * .BSR. 41414 .EOTOF 41541 .ETOF3 41547 * .SWITC 41576 | .TCHEX 42103 .BASIO 42106 * | 50. .IOCSM 42111 51. // 57410 | 0I/O BUFFERS 42111 THRU 57372| 0UNUSED CORE 57373 THRU 57407| 0 1 LINEAR LEAST-SQUARES CURVE FITTING PROGRAM EXAMPLE BROWNLEE 0EQUATION 1 OF A MULTI-EQUATION PROBLEM 0DATA READ WITH STANDARD FORMAT 0BZERO = CALCULATED VALUE 0DATA INPUT 3 INDEPENDENT VARIABLES 3 DEPENDENT VARIABLES OBSV. SEQ. 1-11-21 2-12-22 3-13-23 4-14-24 5-15-25 6-16-26 7-17-27 8-18-28 9-19-29 10-20-30 2 1 49.000 31.000 6.000 355.000 -0. -0. 3 1 181.000 38.000 8.000 631.000 -0. -0. 4 1 4.000 17.000 2.000 174.000 -0. -0. 5 1 22.000 20.000 4.000 603.000 -0. -0. 6 1 152.000 39.000 6.000 331.000 -0. -0. 7 1 75.000 30.000 7.000 759.000 -0. -0. 8 1 54.000 29.000 7.000 282.000 -0. -0. 9 1 43.000 35.000 6.000 316.000 -0. -0. 10 1 41.000 31.000 5.000 490.000 -0. -0. 11 1 17.000 23.000 4.000 457.000 -0. -0. 12 1 22.000 21.000 3.000 282.000 -0. -0. 13 1 16.000 8.000 3.000 174.000 -0. -0. 14 1 10.000 23.000 3.000 269.000 -0. -0. 15 1 63.000 37.000 6.000 240.000 -0. -0. 16 1 170.000 40.000 8.000 525.000 -0. -0. 17 1 125.000 38.000 6.000 257.000 -0. -0. 18 1 15.000 25.000 4.000 240.000 -0. -0. 19 1 221.000 39.000 7.000 331.000 -0. -0. 20 1 171.000 33.000 7.000 331.000 -0. -0. 21 1 97.000 38.000 6.000 457.000 -0. -0. 22 1 254.000 39.000 8.000 776.000 -0. -0. 0DATA TRANSFORMATIONS POSITION CODE OPERATION CONSTANT LOCATION OMIT VARIABLE 1 NONE STRL 1 0 1 2 NONE 124 2 0 2 3 NONE 8 3 0 3 4 2 COMMON LOG 6 5 1 5 8 ADD CONSTANT 6 1 6 6 MULTIPLY BY CONSTANT 6 0 4 0DATA AFTER TRANSFORMATIONS THE FITTED EQUATION HAS 3 INDEPENDENT VARIABLES, 1 DEPENDENT VARIABLES OBSV. 1-10-19 2-11-20 3-12-21 4-13-22 5-14-23 6-15-24 7-16-25 8-17-26 9-18-27 2 4.90000E 01 3.10000E 01 6.00000E 00 5.50228E 01 3 1.81000E 02 3.80000E 01 8.00000E 00 8.00029E 01 4 4.00000E 00 1.70000E 01 2.00000E 00 2.40549E 01 5 2.20000E 01 2.00000E 01 4.00000E 00 7.80317E 01 6 1.52000E 02 3.90000E 01 6.00000E 00 5.19828E 01 7 7.50000E 01 3.00000E 01 7.00000E 00 8.80242E 01 8 5.40000E 01 2.90000E 01 7.00000E 00 4.50249E 01 9 4.30000E 01 3.50000E 01 6.00000E 00 4.99687E 01 10 4.10000E 01 3.10000E 01 5.00000E 00 6.90196E 01 11 1.70000E 01 2.30000E 01 4.00000E 00 6.59916E 01 12 2.20000E 01 2.10000E 01 3.00000E 00 4.50249E 01 13 1.60000E 01 8.00000E 00 3.00000E 00 2.40549E 01 14 1.00000E 01 2.30000E 01 3.00000E 00 4.29752E 01 15 6.30000E 01 3.70000E 01 6.00000E 00 3.80211E 01 16 1.70000E 02 4.00000E 01 8.00000E 00 7.20159E 01 17 1.25000E 02 3.80000E 01 6.00000E 00 4.09933E 01 18 1.50000E 01 2.50000E 01 4.00000E 00 3.80211E 01 19 2.21000E 02 3.90000E 01 7.00000E 00 5.19828E 01 20 1.71000E 02 3.30000E 01 7.00000E 00 5.19828E 01 21 9.70000E 01 3.80000E 01 6.00000E 00 6.59916E 01 22 2.54000E 02 3.90000E 01 8.00000E 00 8.89862E 01 0SUMS OF VARIABLES 1.80200E 03 6.34000E 02 1.16000E 02 1.16717E 03 0RAW SUMS OF SQUARES + CROSS PRODUCTS 1-10-19 2-11-20 3-12-21 4-13-22 5-14-23 6-15-24 7-16-25 8-17-26 9-18-27 2.72692E 05 6.47460E 04 2.07180E 04 1.22340E 04 3.77800E 03 7.08000E 02 1.13543E 05 3.67385E 04 6.85658E 03 7.20773E 04 0RESIDUAL SUMS OF SQUARES + CROSS PRODUCTS 1.18063E 05 1.03428E 04 1.57724E 03 2.28010E 03 2.75905E 02 6.72381E 01 1.33884E 04 1.50093E 03 4.09329E 02 7.20610E 03 0MEANS OF VARIABLES 8.58095E 01 3.01905E 01 5.52381E 00 5.55797E 01 0STANDARD DEVIATIONS OF VARIABLES 7.68320E 01 8.88042E 00 1.83355E 00 1.89817E 01 0SIMPLE CORRELATION COEFFICIENTS, R(I,I PRIME) 1.000 0.758 1.000 0.809 0.847 1.000 0.459 0.445 0.588 1.000 0INVERSE, C(I,I PRIME) 3.062 -0.785 -1.813 -0.785 3.745 -2.538 -1.813 -2.538 4.617 1EXAMPLE BROWNLEE DEP VAR 1 0IND.VAR(I) COEF.B(I) S.E. COEF. T-VALUE INV. DG. MIN X(I) MAX X(I) RANGE X(I) REL.INF.X(I) 0 2.46786E 01 1 -2.49264E-03 8.42E-02 0.0 3.06E 00 4.000E 00 2.540E 02 2.500E 02 -0.01 2 -3.95989E-01 8.05E-01 0.5 3.74E 00 8.000E 00 4.000E 01 3.200E 01 -0.20 3 7.79718E 00 4.33E 00 1.8 4.62E 00 2.000E 00 8.000E 00 6.000E 00 0.72 0NO. OF OBSERVATIONS 21 NO. OF IND. VARIABLES 3 RESIDUAL DEGREES OF FREEDOM 17 F-VALUE 3.1 RESIDUAL ROOT MEAN SQUARE 16.52488446 RESIDUAL MEAN SQUARE 273.07180405 RESIDUAL SUM OF SQUARES 4642.22070313 TOTAL SUM OF SQUARES 7206.10437012 MULT. CORREL. COEF. SQUARED .3558 0IDENT. OBS. NO. OBS. Y FITTED Y RESIDUAL *** OBS. NO. OBS. Y FITTED Y ORDERED RESID. SEQ ASTRIA 2 55.023 59.064 -4.041 8 45.025 67.641 -22.616 1 CANADA 3 80.003 71.557 8.446 13 24.055 44.862 -20.807 2 CEYLON 4 24.055 33.531 -9.476 15 38.021 56.653 -18.632 3 CHILE 5 78.032 47.893 30.139 17 40.993 56.102 -15.109 4 DENMK 6 51.983 55.639 -3.656 20 51.983 65.765 -13.782 5 FINLND 7 88.024 67.192 20.832 19 51.983 63.264 -11.282 6 FRANCE 8 45.025 67.641 -22.616 4 24.055 33.531 -9.476 7 GMANY 9 49.969 57.495 -7.526 18 38.021 45.930 -7.909 8 IRELND 10 69.020 51.287 17.733 9 49.969 57.495 -7.526 9 ISRAEL 11 65.992 46.717 19.274 2 55.023 59.064 -4.041 10 ITALY 12 45.025 39.700 5.325 6 51.983 55.639 -3.656 11 JAPAN 13 24.055 44.862 -20.807 16 72.016 70.793 1.223 12 MEXICO 14 42.975 38.937 4.038 14 42.975 38.937 4.038 13 NTHLDS 15 38.021 56.653 -18.632 12 45.025 39.700 5.325 14 NWZEAL 16 72.016 70.793 1.223 3 80.003 71.557 8.446 15 NORWAY 17 40.993 56.102 -15.109 21 65.992 56.172 9.819 16 PORTGL 18 38.021 45.930 -7.909 10 69.020 51.287 17.733 17 SWEDEN 19 51.983 63.264 -11.282 22 88.986 70.979 18.007 18 SWITZD 20 51.983 65.765 -13.782 11 65.992 46.717 19.274 19 UNKING 21 65.992 56.172 9.819 7 88.024 67.192 20.832 20 U.S.A. 22 88.986 70.979 18.007 5 78.032 47.893 30.139 21 1EXAMPLE BROWNLEE DEP VAR 1 CUMULATIVE DISTRIBUTION OF RESIDUALS .0002 .001 .005 .01 .02 .05 .1 .2 .3 .4 .5 .6 .7 .8 .9 .95 .98 .99 .995 .999 *-----*-------*--*---*-----*----*-----*----*---*---*---*---*----*-----*----*-----*---*--*-------*----- I ' I I ' I P I ' I O I ' I S I ' I I I ' I T I ' I I I ' I V I ' I E I ' I I ' I I ' I I ' I I ' I I ' I I ' I I ' I I ' I I ' I I ' I I ' I R I ' I E I ' I S I ' I I I ' I D I ' I U I ' I A 0 I----------------------------------------------------------------------------------------------------I L I ' I S I ' I I ' I I ' I I 0 I I ' I I ' I I ' I I ' I I ' I I ' I I ' I I ' I N I ' I E I ' I G I ' I A I ' I T I ' I I I ' I V I ' I E I ' I I ' I I ' I *-----*-------*--*---*-----*----*-----*----*---*---*---*---*----*-----*----*-----*---*--*-------*----- .0002 .001 .005 .01 .02 .05 .1 .2 .3 .4 .5 .6 .7 .8 .9 .95 .98 .99 .995 .999 0 CUMULATIVE FREQUENCY, NORMAL GRID 1EXAMPLE BROWNLEE DEP VAR 1 RESIDUAL VS. FITTED Y 33.721 37.522 41.322 45.122 48.923 52.723 56.523 60.324 64.124 67.925 71.725 *----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*- I Z I I I P I I O I I S I I I I I T I I I I I V I U I E I I I U I I U Z + I I I I I I I I I I I I I I I Z I I U+ R I I E I I S I U I I I I D I Z I U I I A 0 I--------------------------------------------------------------------------------------------------R-+ L I I S I I I I I I I U Z I I I I I I I I Z U I I I + Z I I Z I I I N I I E I Z I G I Z I A I I T I I I I I V I Z I E I I I I I Z I I Z I *----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*----*- 33.721 37.522 41.322 45.122 48.923 52.723 56.523 60.324 64.124 67.925 71.725 0 FITTED Y 1 2874 LINES OUTPUT. $IBSYS $STOP | 1 PERIPHERAL UNIT POSITIONS AT END OF JOBS | 1 SYSPP1 IS B2 REC. 00514, FILE 00000 SYSOU1 IS B1 REC. 00779, FILE 00000 SYSIN1 IS A2 REC. 00002, FILE 00001| 1 END OF JOBS | 1