Urban and Regional Report No. 31-24 _ DOCUMENTATION FOR THE LDC URBAN HOUSING MODEL (in two volumes) VOLUME II Raymond J. Struyk The Urban Institute and The World Bank August 1976 These materials are for internal use only and are circulated to stimulate discussicn and critical comment. Views are those of the author and should not be interpreted as reflecting the views of the World Bank. References in publications to Reports should be cleared with the author to protect the tentative character of these papers. These papers were prepared under Research Project 671-37. ACKNOWLEDGEMENT The research reported in this paper was performed under contract with the International Bank for Reconstruction and* Development. Douglas Keare and Bertrand Renaud of the Bank's Urban and Regional Division, Development Economics Department have provided extremely valuable guidance throughout the work. Very useful comments were. obtained form an Urban Projects Department seminar at the Bank. An additional acknowledgement is given to.Larry Ozanne for general discussions about this project. The views expressed in this paper are the authors' and not neces- sarily those of the World Bank or the Urban Institute. TABLE OF CONTENTS Page Introduction ..................... 1 SummaryofSubroutines ..................... .......... 2 Model Inputs ..................... . .............. 5 Program Listing . . . . . . . . . . . . . . . . . . . . . . . . . 19 SampleProgram ............................ .......... ... 52 ii 1. INTRODUCTION This manual provides various descriptive material necessary for the LDC model computer program to be run. The presentation is composed of four parts: (a) brief description of the purpose of each of the subroutines;- (b) definition of the inputs to the model; (c) listing of the full program; and (d) the listing of the inputs and outputs for a sample problem. The prior papers contain material essential for understanding the model which is not included here. The "basic reference" for the general theoretical structure of the model is, "A Simulation Model of Urban Housing Markets in Developing Countries." 1 The working of the solution algorithm is described in some detail in The Web of Urban Housing.2 The LDC program is written in FORTRAN IV, in a version compatible with a CDC 6600 machine. Its core requirement is about 60,000 CDC words (10 bytes). 'R. Struyk, Urban Institute Working Paper 5062-1, 1976. 2Chapter 3 and Appendix A of F. de Leeuw and R. Struyk, The Urban Institute, 1976. 2. SUMARKY DESCRIPTION OF SUBROUTINES Table 1 LIST OF SUBROUTINES SUBROUTINE PURPOSE SUBROUTINE CALL ARGUMENTS, IF ANY ADJUST It sets the occupancy status of each The Jth dwelling dwelling based on number of households choosing each dwelling. PRINT Prints out after N iterations the status Controlled by of the household and dwelling arrays MPRINT FUNCTU Governs the calculation of the utility The Ith household function for new and existing dwellings The Jth dwelling and selects the utility maximizing new dwelling for each household. SUPPLY Modifies the price and quantity of ser- The Jth dwelling. vices of each dwelling depending on occu- pancy status and using the supply function. ZONE It calculates for each zone a profile of The Jth dwelling households and dwellings. MAXIM Matches households to dwellings based The Mth dwelling on the utility each.household derives from each dwelling. SNPSHOT Prints out the contents of the utility The Jth dwelling matrix and dwelling iriformation. It is controlled by the parameter X2. ERRTB Calculates the error of fit of the current solution compared to the actual end-of- simulation period conditions in the market. Controlled by TFLG parameter. UCAL Subroutines is called by FUNCTU and calcu- The Ith household lates the utility fraction for a given The Jth dwelling dwelling and household The Kth zone The Price PSTAR The Quantity of Service QSTAR 4 Table 1 (continued) SUBROUTINE PURPOSE SUBROUTINE CALL ARGU1MENS, IF ANY SETPRC Modifies the price per zone PRC(XZ) based on 1) The input arrays AP, BP; 2) the maximum number of dwellings allowed for each new zone, MDWL. Controlled by STPFLG. CROWD Aggregates subdwellings to the master dwelling level. The identity of the master dwelling is st.ored in DWL (ll,J). Controlled by NC,ROWD. POLICY Calculates the change between the current solution and a prior solution in housing quantities, prices, expenditures, and household utility for specified groups of households. Controlled by MPLCY. 5 3. MODEL INPUTS There are basically two kinds of card inputs. (The program does not use tapes as inputs or outputs.) The first cards follow the CDC procedure of a NAMELIST FORTRAN input. These cards contain parameters and control flags needed to activate several routines called only under special conditions; these inputs are listed in Table 2. The rest of the cards, which are data cards, are documented in Table 3. These data are the input arrays for model households and dwellings, travel times by income class, and the results of a "base case" solution; the last are used in performing selected, initial analysis of the results of a policy simulation. 6 Table 2 INPUT STATEMENTS CARD #1 $INPUT CARD #2 VARIABLE TYPE MEANING ACCEPTABLE EXPLANATION XAM VALUES NHH Integer Number of households 1-55 Constant MDL Integer Number of dwellings 1-90 Constant N=HYPS Integer # Types of household Constant NZONES Integer Total # of zones 1-10 NZONES NZE + Constant NZN NZE Integer Totl # of existing zones 1-5 Constant NZN Integer Total # of new zones 1-5 Constant NINC Integer Total # of income classes 1-5 Constant for transport costs NITER Integer # of minor iterations per No Standard range Constant major iterationu limits is 90 - 200. MIT Integer Maximum # of major no Standard value Constant iterations allowed before limits is 80. execution stops 7 Table 2 (continued) CARD #3 VARIABLE TYPE MEANING ACCEPTABLE EXPLANATION NAME VALUES MPRINT Integar # of the iterations be- No Standard range Constant tween printing preliminary limits 5-20 output MODP Integer Option for grouping exist- 1, 2-3 Set to one Constant ing dwellings for no grouping MPLCY Integer Flag to call policy 0, 1. Set to one Constant subroutine to call policy KGPRHH Integer Option flag to group 0, 1 Set to one-to Constant households for policy group house- analysis hold If one,input cards type F must follow the cards type E NGRPH Integer # of groups desired when 0, 1-4 Constant grouping households CARD #4 VARIABLE TYPE MEANING ACCErABLE EXPLANATION NAMLE VALUES TITLE. ALPHA Title of run ALPHA (80 Char- acter) 8 Table 2 (continued) CARD $INPUT 1 CARD #6 VARIABLE TYPE NEANING ACCEPTABLE EXPLANATION NA~ VALUES GAMMA I Real Utility function parameter 0-1.0 See basic Constant reference. GAMMA 2 Real Utility function parameter Positi7e See basic Constant reference. GAMMA 3 Real Utility function parameter Positive See basic Constant reference. GAMMA 4 Real Utility function parameter Positive See basic Constant reference. Xl Real Size of price adjustment .0001 - Constant step for existing dwell- .03 ings X2 Real Flag to call SNAPSHOT 0, 1 Set to 1 to Constant subroutine. get printout of U matrix. X3 Real Price increment for Used only in Constant Zone 1. special situa- tions where opportunity cost of land might not be reflected in market prices, e.g., urban renewal. WT Real Fraction of hourly wage 0-2.0 Usually set Constant rate at which commutation at 0.5. time is value 9 Table 2 (cbntinued) CARD #7 $INPUT 2 CARD #8 VARIABLE TYPE MEANING ACCEPIABLE EXPLANATION NA1IE VALUES ALPHA Real Utility function parameter 0-1.0 One value for array each household (10 type i, See values) basic reference. 3'THAl Real Supply function parameter 0-1.0 One value for array each zone of (10 existing values) housing (NZE); See basic reference. BETHA2 Real Supply function parameter Positive One value for array each zone of (10 existing values) housing (NZE); See basic reference. CARD #9 VARIABLE TYPE MEANING ACCEPTABLE EXPLA.ATION NAbE VALUES BETHA3 Real Supply function parameter Positive One value for array each zorne of (10 existing values) housing (NZE); See basic reference. BEX Real Fraction of residential O-1.0 One value for array sites in zone which are each zone (10 serviced (NZOTES) ; termed values) Bj in basic reference; NZE followed by NZN. 10 Tab-le 2 (continued) CARD #9 (cont.) VARIABLE TYPE MEANING ACCE CABLE EXFLANATION NAME VALUES PZ Real Price per unit of operat- Positive One value for array ing units (existing-new) each zone (10 (NZONES); the values) new market sec- tor must be the first of the new construction zones. NZE fol- lowed by NZN. CARD #10 VARIABLE TYPE MEANING ACCEPTABLE EXPLANATION NA1tE VALUES PC Real Price per unit of Positive One value for array capital inputs each zone (10 (NZONES) , the values) new market sec- tor must be the first of the new construction zones. NZE fol- fowed by NZN. Al Real Parameter of cash 0 or One value for array transfer equation Positive each household (10 type. values) A2 Real Parameter of cash 0 br One value for array transfer equation Positive each household (10 type. values) Table 2 (continued) CARD #11 VARIABLE -MEANING ACCEPTABLE EXPLANATION NAME VALUES B1 Real Parameter for 0 or One value for array tax equation Positive each household (10 type values) B2 Real Parameter for 0 or One value for array tax equation Positive each household (10 type values) QMT Real Minimum quantity of 0 or One value for array services which mist be* Positive each household (10 consumed to quality for type. values) some government housing programs. RDT Real Race or social group 0 or 1 One value for array indicator each household (10 type values) DWLM Real Multiple dwelling Positive Set to zero if array indicator each represents (10 only itself; values) one value for each zone CARD #12 $INPUT 3 CARD #13 VARIABLE TYPE MEANING ACCEPTABLE EXPLANATION NAME VALUES QOMIN Real Minimum quantity of hous- 0, Positive Onea value for arLay ing services taquired each new con- (5. struction zone vslues) (NZN) QMAX Real Maximum quantity of Positive One value Lor array housing services each new con- (5 allawed. struction zone values) (NZN) 12 Table 2 (continued) CARD #12 (cont.) VARIABLE TYPE MEANING ACCEPABLE EAATION NAME VALUES CIDA Real Hourly travel cost by 0, Positive One value for array income class each income (5 class (NINC) values) CARD #14 VARABLE TYPE MEANING ACCEPTABLE EXAMINATION NAI,E VALUES CLWYINC Real Lower bound for income 0, Positive One value for array classes for travel costs each income (5 class (NINC) values) CHYINC Real Upper bound for ineome Positive One value for array classes for travel each income (5 costs class (NINC) values) CARD #15 $INPUT 4 CARD #16 VARIABLE YPE MEANING ACCEPTABLE EXAMINATION NAIAE VALUES Z Real Proportion of minority 0-1.0 One value.for array social or racial group each zone; (11 in zone in base period See deLeeuw- values) Struyk, Chp. 3. NZE followed by NZN. VAC Real Vacancy indicator; used Positive Variable value array only for special testing integer is dwelling ID (5 purposes values) 13 Table 2 (continued) CARD #17 $INPUT 5 CARD #18 VARIABLE TYPE MEANING ACCEPTABLE EXPLANATION NAME VALUES IPUNCH -Integer Flag to call PUNCH 0, 1 Set to one to constant output option call PUNCH YCL Real Lower income bound to Positive array group households for (5 analysis of solution values) output YCU Real Upper income bound to Positive array group households for (5 analysis of solution values) output CARD #19 $INPUT 6 CARD #20 VARIABLE TYE MEANING ACCEPTABLE EXPLANATION NAME VALUES TFLG Real Flag to call error sub- 0, 1, 2 1 - 1960 constant routine 2 = 1970 To call ERRTS subroutine, used only when fitting model to actual city EX?LM Real Low expenditure limit Positive Maximum month- constant ly rent for low-expendi-. ture unit 14 Table 2:(continued) CARD #21 VARIABLE TYPE M4EANING ACCEPTABLE EXPLANATION NAME VALUES ERRMAT Real Actual values of error Positive 1960 matrix criteria by zone (except 1. Income (10 x 5) low expend-70) 2. Percent Zone Variable white 3. Rel. price 4. Prop. low- expend. 1970 1. Income 2. Percent white 3. Expenditure 4. Utilization 5. Prop. low expenditure units CARD #22 VARIABLE TYPE MEANING ACCEPTABLE EXPLANATION NAME VALUES ERNEAN Real Mean values of error Positive array criteria; i.e., the actual (7 end-of-period values for values) each criteria 15 Table 2 (continued) CARD #23 $INPUT CARD #24 VARIABLE TYPE MEANING ACCEPTABLE EXPLANATION NAAME VALUES NCROWD Integer # of crowded dwellings 0, Positive If positive constant subroutine CROWD is called; See basic refer- ence QMST Real Quantity of services pro- Positive See basic array vided by "full" dw-elling reference (10 (i.e. sum of subdwellings) values) in base period STPFLG Real Flag to call SETPRC 0, 1 Set to 1 to constant subroutine call SETPRC DELTPR Real Price increment for new 0-1.0 Used by SETPRC constant units under most conditions to modify price of new dwell-' ings M) WL Real Maximum # of dwellings for 0, Positive Set to 99 if array each new zone no restriction (5 to construc- values) tion present; one value for each new zone (NZN). Fill in array with 99's, if less than 5 new construction zones. CARD #25 - VARIABLE TYPE MEANING ACCEPTANCE EXPLANATI(G1 NAME VALUE AP Real Intercept in price deter- Used by SETPRC; array mination equation See basic (5 reference values) BP Real Slope parameter in price Used by SETPRC; array determination equation See basic (5 reference values) 16 Table 3 DATA INPUTS CARD TYPE A TRAVEL TIME MATRIX Field No. Field Name Cols. Values Explanation 1 Travel Time for Zone 1 1-5 Positive Z Travel Time for' Zone 2 6-10 Positive 3 Travel Time for Zone 3 11-15 Positive 4. Travel Time for Zone 4 16-20 Positive 5 Travel Tine for Zone 5 21-25 Positive 6 Travel Time for Zone 6 26-30 Positive 7 Travel Time for Zone 7 31-35 Positive 8 Travel Time for Zone 8 36-40 Positive 9 Travel Time for Zone 9 41-45 Positive 10 Travel Time for Zone 10 46-50 Positive CARD TYPE B HOUSEHOLD ARRAY Field No. Field Name Cols. Values Explanation 1 Household ID 1-5 2 Monthy Income 6-10 3 Household TYPE 15 1-9 4 Assistance Eligibility 20 0, 1 1 = Eligible for particular housing program 5 Model Income 46-50 17 Table 3 (continued) CARD TYPE C: DWELLING ARRAY Field No. Field Name Cols. Values Explanation 1 Dwelling ID 1-5 Positive 2 Initial quantity of 6-10 Positive service 3 Zone 14-15 1-10 4 Crowding ID of master 54-55 0, 1-11 Only if crowding is dwelling used CARD TYPE D: TITLES FOR POLICY RUN (Only needed if MPLCY = 1) Field No. Field Name Cols. Values Explanation 1 Title of Policy 1-40 ALPHA Simulation 2 Title of Base 41-80 ALPHA Simulation 18 Table 3 (contiaued) CARD TYPE E: BASE HOUSEHOLD MATRIX (Only needed if MPLCY 1) Field No. Field-Name Cols./Decimal Explanatian place 1 Household ID 2-5/0 2 Monthly Income 6-15/3 3 Household Type 16-20/0 4 Assistance Eligibility 21-25/0 5 Income Shift 26-33/0 Transfer or tax 6 ID of Dwelling Selected 34-38/0 7 Price of Dwelling 39-48/5 Selected 8 Quantity of services of 49-58/3 selected dwelling 9 Zone of dwelling 59-63/0 10 Model income of 64-73/3 household 11 Utility of household from dwelling selected CARD TYPE F: HOUSEHOLD GROUPS FOR POLICY ANALYSIS EXpLANtATION (Only needed if MPLDY = 1) Field No. Field Name Cols. Explanation 1 Household ID 1-2 2 Household ID 3-4 3 Household ID 5-6 39 Household ID 77-78 40 Group Number 79-80 6T RJOEE2U FTN V3.0P380 OPT=l 07/28/76 08.56,00, PAGE 1 C dASIC LECK - - - - PITTX1 'MAIN 2 PROGRAII HOUSt2( INPUT,OUTPUT,PUNUH,TAPE5ziNPUT,TAPE6-OUTPUT) MAIN 3 C COMMON /1TT/ITER MAIN 4 COMMON /HOSHL0/ HH(10,55) MAIN 5 5 CUrI11014/DWELNG/ DWL(ilogG) MAIN 6 C CO11il0N /UF/ U(55p41).USTAR(55)pMUL1 MAIN 7 GCQM10ON /ALPH AS,'ALPHA-1 o) MAIN ColIIOrN /PRM2/NHHDMUL MAIN 9 COlCII JX/KKNT/ KOUNT(5),KITER MAIN 10 1 0 COII01ON /liULT/ DWLM(10).TDWL(9o) MAIN 11 COilMON /SOLV/LSOLV#KSOLV MAIN 12 C COilMON /DEIUG/ V(21h)MODP MAIN 13 CO`MMfiOfl /WAGES/lT MAIN 14 CO11MON /ASSIST/A1(10) ,A2(10),B1(10),B2(10),OMT(1O,RODT(10) MAIN 1c 15 COfI140N /FIX/Z(11.'.VACC5) MAIN 16 CO8MON /PRTABX/TABX(6,1j) ,SMXX MAIN 17 COIIMON/PUNCC/IPUNCH. yCL(10). yC (10) TITLE(8),HHY(55) MAIN 18 C COMMON /JPARMi/ TFLG EXPLM,ERMATV(6,7).RMEAN(7) MAIN 19 CO:llI1OtJ/il'Z/i4HTYPS ,NONeS,NZE.NZN MAIN 20 20 COH.IIU1NPARAM/ KFLGoXl,X2,X3eX4 MAIN 21 C COili1lIJ/HVCTORIHINC(55) MAIN 22 COMMON/HETHAS/ dETHAI(10),BETHA2(10)gBETHA3C10) MAIN 23 COIMMON/GAMMAS/GAMMA1,GAMMA2#GAMMA3#GAMMA4 HAIN 24 C COM1lON/PRICE/PZ(10),PC(1O),PRC1lO) MAIN 25 25 CO,.,MON/DIST/CLMSOA(5),TT(5.10) MAIN 26 COfiH4ON/QINMXI OMIN(5)oQMAXC5) MAIN 27 C COtiMuN / EXES/ ByXCO) MAIN 28 COMMON/SETP/ AP5hBPA5),MXDWL(5) ,DLTPR MAIN 29 COMMoN'/CRWD/ NCROWD,oI4ST(10) MAIN 30 C 30 COiM11ON /GRPH/ NGRPH,NHG(4),IXHH(4O,4) MAIN 31 CO*1IMUN /BASEHH/ HHHS(10,55),U8S(55) MAIN 32 COlIMOiN /TTLE/ TITLEP(4),TITLEB(4) MAIN 33 Q DIMEtNSION CLWYNC(5)#CHIYNC(5)tKKXT(10) MAIN 34 NAMELIST /INPUT/NIHHDMDLINHTYPSDNZONES,NZE,NZI.NINC.NITERfMAXIT. - MAIN 35 3) 1 14PRINT.MODP,MPLCY,KGRPHH,NGRPH hAIN 36 c NiAMELIS T/INPUTI/GAMMA1aGAMMA2,GAMMA3,GAMMA4,XlX2#X3,WT MAIN 37 NAMELIST/INPUT2/ALPHAeBETHA1,BETHA2DBETHA3,BEXePZePC,AIeA2q MAIN 38 1 Bl*B2vQ11T#RUT#Dl4LM MAIN 3 NAMELST/AICNPUT3/QLN,QMAX,CL LNC.WYNCCCHIYNC MAIN 40 40 NAMELIST/INPUT4/Z#VAC MAIN 41 NAMELIST/INPUT5/IPUNCH,YCL,YCU MAIN 42 C NAMELIST/INPUT6/TFLGDEXPLMDERMATIRMEAN MAIN 43 NAMELIST/INPUT7/NCROWD,QMST,STPFLG,DELTPR,MXDWLIAP,SP MAIN 44 C MAIN 45 0 4! C MAIN 46 C MAIN 47 C HOUSEHOLD DATA INDICES ARE UEFINEU AS FOLLOWS, MAIN 48 C 1. IDENTIFICATION NUMBER MAIN 49 C 2. MONTHLY INCOME MAIN 50 50 C 3. HH TYPE (1;WHITE FAMILY, 2^WHITE SINGLE. 3aNONWHITE FAMILYf MAIN 51 C 4. ASSIS.TANCE ELIGIBILITY dI.ELIGIBLEt 0;-NOT) MAIN 52 C 5, INCOM1E SHIFT (NiEGATIVE FOR TAX, POSITIVE FOR TRANSFER) MAIN 53 C *6, IDENTIFICATION NUMBER OF DWELLING SELECTED MAIN 54 C 7, PRICE OF DWELLING SELECTED MAIN 55 55 C B. UUANTITY OF SERVICES OF DWELLING SELECTED MAIN 56 PROGRAl M uUsE2 RJO1E2O FTN V3,O-p38 Opt-1 01/28/76 08,56,00, PAGE 2 C 9, ZO E OF OWEI. I ,yG SELECTED MAIN 57 c 10. IODEL' IIOMES MAIN 58 C MAIN 59 C [J4LLLING DATA 1NDICES ARE DEFINED AS FOLLOWS, MAIN 60 6U C 1, IDENTIFICATION NUMBER MAIN 61 C 2. lINITIAL QUANTITY OF SERVICES MAIN 62 C 3. ZO tE MAIN 63 C i, CURRENT QUANTITY OF SERVICES MAIN 64 C 5, AVERAGE NET RENT (P*Q-PZ*O) IN ZONE MAIN 65 65 C 6. PRI C.E MAIN 66 C 7, WHITE PROPORTION OF TOTAL HOUSEHOLDS IN ZONE MAIN 67 C 8. OcCUPANCY INDICATOR (1;OVER, 2-JUSF, 3mVACAN?j? GT PZ)D MAIN 68 C 4zVACANT(P=PZ)o 50VACANT(P LT PZ1) MAIN 69 C 9. NUMBER OF HOUSEHOLDS CHOOSING DWELLING MAIN 70 7U C 10. SUPpLy MOVEMENT INDICATOR CO-NO MOVEMENT THIS ITERATON# MAIN 71 C -1iDOWNWARD MOVEM.ENT, ,42'UPWARD MOVEMENT) MADtI 72 C MAIN 73 C MAIN 74 C MAIN 75 7> KASE=O MAIN 76 9999 CONTINUE MAIN 77 LSOLV=o MAIN 78 READ(5Pl NPuT) MAIN 79 WRITE(6,.NPUT) MAIIN 80 80 I IF(NHH ,EQ. 0) STOP MAIN 81 READt5 907) TITLE MAIN 82 C vi1EN Xi -O. ZoNEl DOES NOT RECEIVE SPECIAL TREATMENT, MAiN 83 - 1Y=5 MAIN 84 IF(KASE ,GT, 0) GO TO 9998 MAIN 85 READ(5,INPUT1) HAIN 86 WRITEC6, iNPUT1) MAIN 87 READ(5, INPUT2) MAIN 88 WRITE(6. INJPUT2) MAIN 89 READ(5, INPUT3) MAIN 90 9U READ(5s INPUT4) MAIN 91 907 FORMAT(8A10) MAIN 92 READ(5, INPUT5) MAIN 93 READ(5, INPUT6) MAIN 94 READ(5, INPUT7) MAIN 95 9 WRITE(6.INPUT3) MAIN 96 C) WRITE(6,INPUT4) MAIN 97 WRITE(6,INPUT5) MAIN 98 WRITEC6, ltJPUT6) MAIN 99 CV;n' I TE; IJPIU 7; MAIN c0e 100 WRITE(6,197) MAIN 101 97 FORMATl1H1,2X,*TT VALUES*//1X,*INCLASS*/) MAIN 102 ;) D0 102 IM=1,NINC MAIN 103 READ(5,103)(TTCIM,MZ),MZz1,10) MAIN 104 WRITEC6.195)lt,(TT(IM,MZ),MZ;1q0) MAIN 105 C) 10' 195 FORMAT(2X,13,4X,10F5.O) MAIN 106 102 COINTINUE MAIN 107 103 FORMAT(10F5.0) MAIN 108 0 DO 202 I=1,NHH MAIN 109 READ(IN,199)(HH(J J,JII1,0) MAIN 110 HH(2.I1)-,HH (2,I)/12. MAIN 111 C ') - . . .-....... PROGRAM' HOUSE2 RJOEE2O FTN V3.0-P380 OPT-1 07/28/76 08,56,00, PAGE HHY() -H (2# 1) MAIN 112 HH(1O,I)=Hnt(10,1 )/12. MAIN 113 C BtOVE IS CONVElHSION TO MONTHLY INCOMES. MAIN 114 J=Hi(3. I) MAIN 115 1 1! I F(HH( ) 4 EO, U. ) GO TO 401 MAIN 116E HHi(. I)=AlCJ)-A2(J)#HH 2 t I MAIN 117 GO TO 202 MAIN 118 401 IF(HH(2,I) ,LT. 82(J)) GO TO 202 MAIN 119 HH(5,I)=8l(J)*(82(J)-HH(2sI)) - MAIN 120 20U C INCOME SHIFT FORMULAS, MAIN 121 202 CONTINUIE MAIN 122 WRITEC6,4) MAIN 123 4 FORMAT(lHO#*INCOME DATA*/lX.2HlP,4Xq*INCOME*,2Xe,TRANSFER*g?X, MAIN 124 X *SUM*,4X,*MOUDL INcOME*,2X.*TRAVEL COSTt,2X,eINC CLASS*/) MAIN 125 12! DO 207 I1,NHH MAIN 126 Y=HH(2I)+HH(5, 5) * MAIN 127 JgHH(1 I) MAIN 128 C INCOME CLAfSIFICATION FOR EACH HOUSEHOLD MAIN 129 FINC=HH(2D I) MAIN 130 130 DO 205 IW91,5 MAIN 131 IF(FINC,GE,CLWyNC(IW),AND,FINC,LE,CHIyNC(lW)I HWIN (I)lIW MAIN 132 ?05 CONTINUE MAIN 133 ITX=HINC(I) MAIN 134 WRITE(6,7)J,H (2: )IHHC5s )VYHHHlIIlI)#CIMBDA(ITX).iTX MAIN 135 1~3 HH(2.1)-,5*(H(2, )+H (HIO#I)' MAIN 136 207 CONTINU4 MAIN 137 7 FORMAT1X, 12,!F1O,2, 4X, I2) MAIN 138 C MAIN 139 C MODEL INCOMES ARE DEFINED AS YM ; (Yet,6) * (YAVG**.4) WHERE YAVO MAIN 140 140 C IS THE AVERAGE INCOME FOR THAT FAMILY TYP9. IT IS USED TO SIMPLlj MAIN 141 C FY THE UTILITY FUNCTION BECAUSE OF THE FACT THAT ON A NATIONWIDE MAIN 142 C BASIS, THERE IS APPRQXIMATELY UNIT ELASTICITY BETWEEN RENT AND MAIN 143 C MODEL INCOME. MAIN 144 C MAIN 145 * WRITE(6,5) MAIN 146 5 FORMAT(1HO#*DWELLING DATA*/1Xe2HIDf?Xf.SERVICESe,2X,*ADJUSTED*# MAIN 147 . 6X*Z04NE.,3XfoCR0WDED*/) MAIN 148 SUMNo MAIN 149 00 201 I11,MDL MAIN 150 154 READ IN.199)(DWL(J!I)J;1,14) MAIN 151 i99 FORMAT(11F5.O) MAIN 152 C MAIN 153 C MAIN 154 C NOTE -- FOR THE ALGORITHM TO WORK PROPERLY. IT IS NEC-cSSARY FOR MAIN 155 c THE INPUT CARDS REPRESENTING EXISTING GTOCK TO ENTERE4 IN DEw MAIN 156 C SCENDING ORDER OF INITIAL SERVICES, THAT ISt THE FIRST SUCH PARD MAIN 157 o READ MUST HAVE THE HIGHEST 00. IT IS ALSO z TRONGLY ADVISABLE TO MAIN 158 C ENTER ALL HOUSEHOLDS IN DESCENDING ORDER OF MODEL INCOME A5 MAIN 159 C WELL (FOR EVEN MORE PRECISION, THIS SHOULD BE ALPHA*INCOME), MAIN 160 18P C MAIN 161 C MAIN 162 K=DHL(ioI) MAIN 163 DW-DWL(2,I) MAIN 164 WRITE(6,3) K,#DL(2,I),DW#DWL(3#1),DWL(I.1eI)? MAIN 165 16F 3 FORMAT(1X, I2s5F1O,2) MAIN 166 * PR)GHAPI HOLJSE2 RJOEE20 FTN V3.0-P380 OPT= 07/28/76 08,56,00, PAGE 4 6 PWL(2, I )D W MAIN 167 o LzPWL(3,I) MAIN 168 |UlII Dl&M(L)+SUM MAIN 169 TDWL( I )=SUM MAIN 170 170 201 CONTINJUE MAIN 171 C WRITE(6*909) TITLE HAIN 172 909 FORMAT(//jOX.*M0DEL VERSION*#5X,8A1O/) MAIN 173 C) CC IF MPLCY ; 1 HEAD HH BASE ARRAY FOR POLICY MODULE MAIN 174 C If(MPLCY.NE.1) GO TO 200 MAIN 175 I READ(INoY2 o) TITLEP,TITLEB MAIN 176 920 FoRliAT(4AI0I4AIO) MAIN 177 uO 50 IW,1 ,tIHH MAIN 178 READ(IN,921)(IIHBS(lV,I}w)IV-t1.1tl,u8S(IW) MAIN 179 50 CONTINUE MAIN 180 w 180 921 FURMAT(1X,F4.QoF1Uo3,2F5,QuF1O,3#F5,O#f1F,5,F10,4.F5,o#F6O,§3F$,?) MAIN 181 IF(KG PHH.NE.1) GO TO 200 MAIN 182 WRITE(6,924) MAIN 183 - 1 I924 FORMAT(/ 2X ;iHOUSEHnLD GROUPS FOR POL#cY MUAINOU //;; q HAl N 184 I 115X. NHG/) MAIN 185 IV 9?3 FORMAT(2X,I0#,9(1XI2)#?X#Ij? MAIN 186 DO 55 IW=1jNGRPH MAIN 187 READ(IN,922)(IXHH(IXIIW),.IX;139 ),NIHG( W) MAIN 18O WRITE(6,924)IWf(IXHH(-IX,IW)qIXe2539)NHG IGW) MAIN 169 55 CONTINuE MAIN 190 Igo 922 FOR1IAT(4012) MAIN 191 200 CONTINUE MAIN 192 CA DO 250 IZzl#NZONES MAIN 193 N PRC( IZ)rPZ( IZ)tFPC(Z) MAIN 194 Q_ 250 CONTINUE MAIN 195 19! 9998 CONTINUE MAIN 196 KFLG-o MAIN 197 K K ITERzO MAIN 196 300 KSOLV-O MAIN 199 LSOLV=O MAIN 200 O ?00 299 CONTINUE MAIN 201 C KFLG=KFLGt1 MAIN 202 IF(KFLG *EQ. 1) 1020 MAIN 2032 O C INITIALIZE UTILITY MATRIX FOR EXISTING DWELLING$ MAIN 204 1U 00 12 l14,NHH MAIN 205 20$ DO 12 J=1IMIL MAIN 206 O 12 U(I,J);-1.Et6 MAIN 207 C REMEI1HER# UTILITY HERE IS TH9 LOGARITHM OF TRUE YTIUTY, MAIN 208 MDL1-MDL+I . MAIN 209 O C INITIALIZE UTILITY MATRIX FOR NOW HOUSINQ MAIN 210 210 DO 15 I;1,NHH MAIN 211 15 CALL FNCTU(IJ,MDL1) MAIN Z12 O C INITIALIZE OCCUPANCY -" EVERYONE IN NEW HMUSINQ MAIN 213 LZN=NZE+1 MAIN 214 DO 16 J.1,MDL , MAIN 215 215 DWL(6#J)zPZ(LZN)+PC(LZN)t,25 MAIN 210 DWL( 8 J)=3, MAIN - 217 16 DWL(9#J)zO. MAIN 218 *ZJ I I;MUL+NHH MAIN 219 jDo 17 J;MD L1,II MAIN ?20 2?U 17 OWL(9,J);1,0 IlAIN I N,. ,, . , a PROOGRiAtM HOUSE2 RJOIE2O FTN y3,0-P380 OPT-i 47/26/76 0O.56,00, PAtE c INITIALIZE USIAR VECTOR AND HOUSEHOLD CHOICES MAIN 222 o c0 16 I,rwNHH MAI N 223 USTARt I zI=U( I ,MDL1j) MAIN R24 18 HH( 6, I);IiDL+l MAIN 225- G GO TO 22 1AIN 226 20 CALL ZONE MAIN 227 IF(KFLG.LT.4' GO TO 22 MAI1 228 lF(STPFLG,EG.1.) CALL SETPRC. tiAIN 229 22 DO 23 Kz1.5 MAIN 230 230 23 KOuriT ;K IMAIN 231 START of MAIN ALGORI TiM MAIN 232 8 SWEEP SERIALLY OVER ALl DwELLJNQS MAIN 233 25 DO 1DO J-1,14DLi MAI N 234 IF(J ,EQ. MDL1) Go TO 31 "AIN 235 2$> CALL A.JUST(J) MAIN 236 ITER=0 MAIN zDjWL(3#J) MAIN 238 D-DWLMlK) MAIN 239 30 CALL SUPPLYlJ) MAIN 240 I 240 ITER4ITER+1 MAIN 241 IF(KFLG,LT, 3 )ITER6NIWR MAIN 242 31 00 32 l 1 ,tNHH MAIN 243 32 CALL FNCTU(I J} MAIN 244 M33 CALL HAXIM(J) MAIN 240 0 NOTE THAT WE CALL MAXIti EVEN IF JTH DWELLI1G JU$T OCCUPIgp qEcAUSa MAIN 246 C T MAY NOT PE THE SAME PQL6LNq AO OCCURED 14 TH9 PREY 9US PTRf MAIN 24? IF(J.EO,MDL1,AND,IX2EQ41,0) cAI4 sNPS1TIJw) MfAN 248 Ir(4 .EQ. MDLI1 SQ TO 100 MAIN 249 CALL ADJUST1J0 MAIN 250 250 IF(X2 ,EQ, i.0) CAk4 SNP4HT14) MAIN 251 a IuJ Qj 1) OR, DWIL.Qf4)-,EQ* It)) ,ANP, $jTgR ILTg MAIN 252 t00 CONT14 UE MAIN 259 DO 110 Jz1#MOL MAIN 255 CALL. ADJUSTWJ) MAIN 256 NN=DilL(8# J) MAIN 257 IF t LSOLV,E, 0) 0DWLi 8J)i NN MAIN 258 C ABOVE TRUNCATION CHANGES PWL8 OF ?V5 TO 2.0 AT THE END OF EACH MAIN 259 C MAJOR ITERATION SO THAT PRICES WILL BE INCREASED SLIGHTLY UPON RE' MAIN 260 ?60 C EtTERING DWELLING ON FIRST MINOR ITERATION OF THE NEXT TRIA A§ MAIN 261 O C PREScRIBED IN SUBROUTINE SUPPLYt MAIN 262. KOUNT(NN)=KOUNT(NN)t1 MAIN 263 .10 GCO1TINUE MAIN 264 o PO 115 KZ=I,NZONES MAIN 265 1.0 1 KKXT(KZ)70 MAIN 266 00 120 KZ-,NHH MAIN 267 O IF(HH(9#KZ)hLE,lZE) GQO To 1a MAIN 268 KTB%HH(9DKZ) MAIN 269 $KXTCKTH);KKXTCKTB)+l MAIN 270 2 ?70 120 CONTINUE M, A IN 271 KI TER-KI TERti MAIN 27? IF'(MOD(KlTER,HPRINT).EQ,O) CALL PRiNT "AIN 273 IF(KITER,GE1tIAXIT) GO TO 902 MAIN 274 2FUKOUNT(1).GE.l.OR.KOUNT(3).GE.1) MP TO 30 AIN 275 275 C FALL THROUGH HERE IS CANJIDATE SQLUTION, MAIN 27q *0 | PROUHAII HOUSE2 RJOEE2O FTN V3,o-P380 OPT.1 07/28/76 08,56,00, pAGE / / KSOLV=I(SO)LV+1 MAIN 277 IF(KSOLV.Lr(,3) GO TO 299 MAIN 278 Do 130 KZ=I,NZN MAIN 279 KTBGNZE+KZ MAIN 280 * Z28 IF(HXDWL(KZ).GQ.99) GO TO 130 MAIN 281 0 Ml.;MMXDWL (KZ)tl MAIN 282 tMJ=MXDWL( KZ)-1 MAIN 283 4 IF(KKXT(KTR).GT.MI) GO TO 299 MAIN 284 IF(KKXT(KTB),LT.MJ) 0O TO 299 MAIN 285 265 130 CONTItNUE MAIN 286 C KSOLV l 0,1,2 IN PHASEI (LSOLVe0), MAIN 287 IF(LSOLV.EQ,0) GO TO 898 MAIN 288 IF(MOD(KITER,MPRINT).NE,0) CALL PRINT MAIN 289 * IFCKSOLV.LT.4) GO TO 299 MAIN 290 @ g90 C KSOLV a 0,1,2.3 IN PHASE2 (LSOLVY)# .WITH PRINTOUT ON BITH FINAL MAIN 291 C AND PENULTIMATE CYCLES, MAIN 292 GO TO 902 MAIN 293 898 LSOLV-1 MAIN 294 KSOLV;O MAIN 295 * GO TO 299 MAIN 296 0 902 CONTINUE MAIN 297. CALL CROWD MAIN ?9g * CALL ZONE MAIN 299 0 908 FORMAT(1Hl,f CO1TENTS OF HH ARRAY AFTER FINAL, ITRATIONt OUT WIT]i MAIN 300 - 30QJ iiHH(2I) AS ACTUAL INCOME AND LOG UTILITYf t/)M MAIN 301 C TEST IF HOUSEHOLP ARRAY JS To BE PUNCHEp MAIN 302 C PRINTS ARRAY EACH RUN MAIN 303 WRITEl6t908) MAIN 301 DO 1907 IJ. NHH MAIN 305 0 40 KFMTP HH'(,) MAIN 306 HH¢2,1) v HHY(I) MAIN 3Q7 -0 HH(4#I)ZHH(4*l)+1 MAIN 308 C SET ASSISTANCE ELIGI8ILITY TO 4 FOfl INNOMES OUST ABOYV THE MAIN 309 C CUTOFF (I,F, THE RANGE YCL TO YCU 1, AND TO 4 FOR MAIN .31Q 310 C HOUSEHOLDS ELFGIBLF BUT NOT ACCEPTING AS$ISTANCE MAIN 3 1 IF((HCH2eI),GE.YCL(FMTP) IAND,II5H(2, l) ,LE,YCVCKFMTP)I) HH(4,#)a;-is MAIN 412 IF((HHit2I),GE,YCLIKFMTP)) ,AN (HH(8IJ).LE,QMT(KFMTP)fl HHj4fl);49 MAIN 313 WRITE (6,906) (HH(J.I),Jj40 lrUsTA'(I)t MAIN 311 . -107 CONTINUE MAIN 315 CALL ERRTB MAIN 316 IF(IPUNCH,EQ.0) GO TO 903 MAIN 317 0 PUNCH 907*TlTLE MAIN 31 DO 905 I;i,NHH MAIN 319 WRITE(6,906) (HH(J, I )#J;1l.Q0)USTAR( j) MAIN 320 32g PUNCH 9o6,UIH(J,I)#J;l#l0)#USTAR(I? MAIN 321 906 FORMAT(lX#F4.O0F,10O,32F5OtF8,3iFfQ@Fl FlOM3fF5 MAIN 322 905 CONTINUE *MAIN 323 903 KASE KASE t1 . MAIN 324 IF(MPLCY.EQO, ,) CALL PQLICY(KQRPHH) MAIN 325 GO TO 9999 MAIN 3.26 iNU MAIN 3?7 9, 5 * -l -.;> r ~~. 0 J>;i sU8ROuTlIE Ai]JuST RJOEE20 FTN V3,"-P363 OPT.i 07/28/76 08,56,00t PAGE SUbROUTINE AUJUST(J) ADJUST 2 COMMOr'N/UWELNu/ DWL(1,#90) ADJUST 3 COMIMON/PARAM/ KFLG#X1,X2#X3,X4 ADJUST 4 COItIMON/PRIcE/PZ(10),PC(10) ,PRC(10) ADJUST 5 COMtMON /MULT DIWLM(4) ADJUST 6 C DWL( 5, WHEN UNDEROCCUPIE0* P ,LT. PZ (ERROR CONDITION) ADJUST 7 C DWL - 4. WHEr.j UNDEROCCUPIED, P ,EQ, pZ (ABANDONED--NO SEARCH ADJUST 8 C DWL( 3, WHEN UINDEROCCUPIED, P ,GT PZ (EMPTY--SEARCHING) ADJUST 9 C DWL( 2.t WHE1N JUST OCCUPIED AND PREVIOUSLY UNDEROCCUPIED ADJUST 10 lo C UDW( t 2. WHEN JUST 0CCUpIED AND PREVIoUSLY oVEROCCUPIEP ADJUST 11 C OWL( 1, WHEN O#EROCCUPIED, ADJUST 12 C ADJUST 13 C NOTE THAT STATUS OF 2,5 MEAN5 NO PRICE INCREASE IN SUPPLY, ADJUST 14 C NOTE THAT STATUS OF 2,0 MEANS DEFINITE PRICE INCREASE IN SUPPLYj ADJUST 15 15 K;DWL(43J) ADJUST 16 -DaDWLM(K) ADJUST 17 PXgPZ(K)+XS ADJUST 18 -FCDwL(9,J)-D ) 35u39,40 ADJUST 19 35 lFC (JL.(6,1J)-PZ(K))1 437?t3 ADJUST go 20. 36 DWL(8UJ)-5, ADJUST 21 00 TO 45 ADJUST 2Z 37 0WL(8,J);A, ADJUST 23 G0 TO 45 ADJUST 24 38 IF((nWL(3#J) ,EQ, 1,) ,AND, lDOW iq ,9,4 PX)1 QQ T9 07 ADJUST 25 DWL(8,#J)3. ADJUST Z6 00 45 ADJUST 27 39 IF (DWL(8,tJ),GE.3) DWLf8#4)i2,5 ADJUST 28 . N IF (DWL(4i,J!,LE,) 0 (10J), ADJUST 29 GO TO 45 ADJUST 30 30 40 RWL(8oJ);1, ADJUST 31 45 ETURN ADJUST 32 END ApJUST 43 * t SUBROuTINE PRINT RJOEE20 FTN V3.0oP380 OPT=1 07/28/76 08.56,00, PA9E i., SUGkOUTINE PRINT PRINT 2 o0MMON /HIOSHLD/..HH(10.55) PRINT 3 COMMON /UF/ U(55.4I)#USTAR(55)#MLll PRINT 4 C0oMM0tJ/UWtLNG/ DWL(I1,90) PRINT 5 COMMON /ALPHAS/AL PHA(1Q) PRINT 6 COMf1oNi/PRICE/PZL1O)DPC(lo),PRC(1O) PRINT 7 C01MON/EBETHAS/ 8THA1(10),BETHA? 10Q} ETHA3(10) PRINT 4 CO0MON/PARAM/ KFLG,9j#X2,X3,X4 PRINT 9 COMMON /PRM2/ NHH,MDL PRINT $0 COMMON /KKNT/ KOUNT(5),KITER PRINT 11 COMMON /DEBUG/ V(21)#MODP PRINT $2 CO.MMoN /MULT/ DWLM(10OTDWL(90) PRINT 13 COMMON/NHZ/NHTYPS NZONES.NZE,NZN PRINT 14 COMMON /ASSIST/Ai10lfA2(*O)oP,%$0)!g?(*Q?fQMT(Q?,RPT($9 PRINT to 1 COMMON /PRhA8X/TABX 6,14,SMXX FA 1NT 16 COMMON/PUNCC/ IPUNCH#YCL(101 YCU(49) T Tlg(§ HY>5? PRINT $7 DIMENSION TA,(6161} TA8711 PRINT DIMENSION HEAD(?,5) PRINT 19 DIMENSIQN QQ(4,$$(4o4)vSU4 ^6),SAYAj0) PRINT - 0 ; 20 DATA HEADO5HUVER ,5HOOCI)P#5JUST o5HOCCUPt PRINT 2* . 9 5 5HVACNTr5H PGPO05HVAcNT,SH PEPO,5HVAvNTI5H PWPQ PRINT 22 N 1T6 PfIlNT .. * MT~0 PRINT 2 Q120, PRINIT 25 25Q;0, PRINT 26 : --P WRITEINfl00Q PRINT 27. FQQ fORMAT (jtHi#fgS.x** 1NQ HU 0IQDI.1/0 - PRINT ?a NRITE(N, u 7i KITER PRINT 2 .. 107 VORMATC 4 * NHa oF JT ATIQN~ ,1l# PRINT3. 3U WRITE (N' 1d. PRINT 106 FCRMAT(4 § Y QCQUPAN9Y. T#TVft0l PRNT DO 20 K 1L5 PRINT 33 .0 WRITE(N,109)814EAP(I J¢K) Ji;ffK0YNT(KI PRINT ;4 20 CONTINUE PRINT 35 -l 0 i 9 FQR11AT(4X,RA5,5XiI*00 PRINT * MTZMT+1 PRINT 37 777 FORMAT(1HO) PRINT § WRITE(N.777. PRINT 3. WRITE(No, 1.1) PRINT 40 II 111 FORMAT(27X,t BY HQUS§HOIp5f/ PRINt 41 WRITE(N,112) PRINT 42 $12 FORMAT(5X1tHTID#4X#*0N§L# 90Pt 9f* e@q PRINT i;. .jQtXPRINT 44 p tCURiNC *5X f PiiM;NC#15X,#TRANSFfRJ q3X,UT1 1TffJ PRINT 4t 00O 25 1$eSPRNT 4 45 25 SUM(I)z0 PRINT 46 Do 30 o1;1,NHH1 PRIlNT 47 . .... UUoUSTAR( I PRINT 40 PF(ViU ,LE. 1.,Et6) UusOt PRINT 4. KKHH(6# I PRINT 50 50 HHH-DWLC1.qKK ) PRINT 51 . 0 WRITE(N,15j HHPI,ft1HHHNw7,IhHH #I Y(I),WH9,IhHI1I 1h PRINT 52 9 Uu PRINT 53 SUM(1)=SUM(1)t+HH(7I1) PIRINT 54 0 SUH( 2) ;SUM (2 )tHH (8, I PRINT ' 52 SUM(3J)SUM(3)+ HHY(I) PRINT P6 *1 'J .E RSUBQUT1NE PRINT RJOEE20 FTN V3,U-P380 OPT-l 07/?8/76 08,5&§00t pGE 5LJiI4) -SUM (4.)+ H,(lo, ) PRINT 57 SUM(5)-SUM1(5)+ HH(5,.) PRINT 58 SUM(6)zSUM(6)* UU PRINT 59 30 COQJTIWNUE PRINT 60 * DO 31 Kzl6 PRINT 61 SUM(K)*SUM(K)/NNHH PRINT 62 31 CONTINUE PRINT 63 SHXXZSUMK(I) PRINT .64 WRITE(hlj116) SUII PRINT 65 115 FOR14AT(1X ,9(F9j2,2X) PRINT 6 - i16 FORMAT(/3X, MEANtS*,15X,6(Fq,?f?X;l PRINT 67 Q TABLE 2 PRINT i WRITEJN,777) PRINt 6970 WRITEMC I0 PNtNT 70i, j 7I~ ;18F0 MATc~5XDBY D#Et.L INQS#//PX.IJWL.lDP,*XeOZONetu9XIftP#iRN 79110 FQRMA7( 5X,#|YPEkNSi/5@U;l#xzN§9*t PRINT 7$ . ' X 9K,lQo,7Xq*OCCUPA C XjqPRQFIT,3XqqCRQWDING*/) PRINT 72 *TOT;F0 PRINT 73 - VA C;; O. PRINT 74- CNEW PRINT 75 * 7 LI. L .iDL 4NHH PRINT 7* 4 TP R 0F 0~ PRINT 77 pO 40 I;i'muL PRINT 78 . KZN-DWL(3#I) PRINT 79 DELPiDWL 6,I)4ZKZN? PRINT 80 PR;PRtI NT 8 EL PF1L; PRINT 82 - PRF2 DELP DL(2 l U*THAWA4%KZN) PRINT 8 PRoFIT;PRF1*PRF' PflNT 3' TPROF;TPROFatPR0FIT PRINT 84 PRINT 05 : 0'DWLM(tZN) PRINT 6 CTOTzCTOTD PR INT 87 ;IF(WL(9,l).LT,D )CVAC*CVACtD 'W.1f1) PRINT -8 WRITE%N,115)DW(LUl, IDil,#pW9 PRINT 09 X .PROFlTsDWL(1j#1) PRINT 90 49 9 4 a CONTINUE PRINT 91 L10Di.+1 PRINT 92 PO 41 IZLI#LL PRINT 93 - IFDWL91),NE41t) GQ TO 41 PRINT 94 CMEW;CNEW+1 PRINT 95 * : CTOTzCTOT+1' PRINT 96 - KzDWL(3#1 PRINT 97 D;DWLM4K) PRINT 98 IFcUWL)L,I ,GT, 3,} CVAC;CvActU90WIL(91I PRINT .99 WRITE(N 115)PWI(4n1fDWL.I1,PW;"6.I),DWLt4,1}DWIi8,I1 PRINt 100 41 CONTINUd PRINT 102 RAT1;CVAC/UTOT PRINT 102 RAT? tCNEW/ jCT OT PRINT 103 WRITE(N,120)RATl,RAT? PRINT 040 * ~ ~~~ X 0FOH"AT ./OXs. VACAN T/ TO AL t ,F'S a q* NF 1 TQ TAL * tFP, ?I* °TQA Pflp PR|N T tvb . R; < X,TPROF~ PRINT jOB MT ,M T+ 1 PRINT :0 DO 53 J1 NZg NES PRINT 110t TA87(J)o. POINT 11- 0. 0 5 NUBROUTIIE PRINT RJOEE2U FTN V3.0-P380 OPT-l 07/28/76 08,56,00, PAGE 0 DO 53 1;,6 PRINT 1*2 0 TAB(I,J)=O PfRlNT 113 53 CONT1t4UE PRINT 11-4 WRITEMN.777) PRINT 115 0 WRITEM(N.1'1) PRINT 116 1?1 FORMAT (25X,l*Y ZONE*//5X,*ZONE. 6X,*CURINC*, 6X, PRMINCP.6x, PRINT *17 .1*IEAN Q*#7X,*MEAIj P*.5X,iktiEAN .EXF##.5Xt*PCT WHITE#/) PRINT 118 0 DO 50 ;1,N1NHH PRINT 119 K.DST .; HH(9,I) PRINT 120 i , 0 KDw;-HH(6 I PRINT *21 * T TB(1.KDST)-T BCl KD$T)+HHY(I) PRINT 122 T B2(2etKDST)=TtB(2g#KDST)tHHj(6a PRNT $2 TkBA3,KDST)-T 8(3 KDST?tHH(7#I) PRINT 124 'T B(4,KD8T)"T6{(4uKDST)+DWL(7,KUW) PRINT 10 TA6B(5KDST);;T A(5,KDST?+HH(QqI)#ItH7f1) PRINT 126 T '7(KDST)zTA 7I(DST)+HH(10fI? PRINT 12? TAB(6 KDSTjTAB(#4KDPST)+1 PRINT 126 04 C&NTINUE PRI NT 129 Do 5? Kz1,NZONES PRINT *30 * IF(TfA(6lK)!LEP,?)QO To 54 PRINT 1 *31 . G 5 JJ;al5 PRINT 132 TA6(jJ,K);TA8CJJ#K(/TAB(PK) PRINT 43; 5; CONTINUE PRINT 134 -Tsa7(K)zTAB7(K)/TAp4IK. PRINT 1 54 CONTINUE PRINT 136 w hRPTER#12NK!TAQ(I#K)tTA87(KT^!8f")ITAP(31K{TA@zfit 4 PRINT 137 Do 5.6 IX;j,5 PRINT 135 TABX( X#Ks)TAB( Xi,K, PRINT *39 o 56 CONTINIUE PRINT 140 194 52 cONTINUE PRINT 141 * ^2 FORHAT.l6X?I2,5Xl6f$0,.4f?X;) PRINT 142 * 00 64 Kt(i,2 PRINT 143 NHT1- 4HTYPSt. PRINT 144 DO 60 J;1 ,NHT1 PRINT s45 * Do 6a .1;1e6 PRINT 146 60 TAB(liJ);-0. PRINT 147 WpITE(N,777) pRINT *48 *I ( KK,EO,?)GO TO 65 PRINT 149 WRITE(N,123) PRINT 150 150 *23 FORMATs20X,*BY HOUSEtOLI TYPl9//5X,#TYPE.Y#5 pCURINP#,l(IPRKt4NJ9, PRINT 151 X 6 6X6M AN Q#,#7X.MEAN P*#7X#SUPST */ PRINT .52 GO To 66 PRINT 153 65 WRITE,(Ni?23) PRINT 154 ?3 FORMAT(15X,*E.EGIB.E HOUSEHOLD QY TY E*//5g.TYP $io§X! UE;uC.I!q6Xt PRINT 155 1 OPRMINCi,6Xt*MEAN Q*t7Xt*MEAN P ,7X,Tt*' PRINT 156 66 DO 61 I;1eNHH PRINT 157 KTY-HH(3,I1).(1.tl-IK)*(1,-HH(4,1)}) PRINT 158 K(Ty IS THE HOUSEHOLV.TyPe WHEN KK;l OR OHEN KK;? ANP EPIQJBL9 PRINT 159 lF(KTYjEQ 0)GO TO 61 PRINT 160 TA6(1,ITY$rTAB(11KTY)+ HHYJI) PRINT 16* TA8(6,#TY)zTAe(6#KTY) t HH(10,. PRINT 162 * TA6(2#KTY);TAB(2oKTY)tHH(6ti? PRINT 163 TAB(3#KTY)4TAB 3#KTY)+HH(17fI) PRINT 16i ' TA8(4*KTY);TA8l(4KTY)+t, PRINT 165 .65 QM=QMT(KTy) PRINT 166 _ 0' . bUBROUTINE PHINT RJOEE2U FTN V3.0-P380 OPT-1 07/28/76 08,56,00, PAGg 9 IF(HH(C8PI)UE.OM)GO TO 61 PRINT 167 6 TA8(5,KTY)mTA35,IKTY)j1, PRINT 165 *61 CONT INuE PRI NT 469 *DO 7 1:1,5 PRINT 170 * j70 DO 67 Jr-1i,HTYPS PRINT 171 67 TA8(I,NHTj)zTA6I ,NHTI)+TAB(I,J) PRINT 172 DO 62 KQ1,NHTYPS PRINT 173 IF(TASC4#K),LE.O,) GO TO 68 PRINT 174 IF(TABC4,K).LE,0.)GO TO 63 PRINT 175 17 TAB(JJe K0TAB (JJ,K)/TAB i0 K) PRINT 176 63 C'NTIN E PRINT TAB(6,k)QTA8(6#K)/TA8fi,K) PRINT 178 68 CONT NUE PRI NT t79 WRlITE(6 124 RK ,TA { ,K TAB ¢K)TAO(? N t TAP( I664 fT g I8 pRINT 18 0 62 COTINUT r PRINT $81 * 64 CONTINUE PRINT 182 .;24 FORMAT(0XqI 2Xt51F1Q,t42X)? PRINT 18l ENV PRINT 18 *. . f I A * , ** * * 0 ,* .' ~*... .9 ,; *0' ! - * 3l *. 0 , b .. t i . - SUBROUTINE FNCTU RJOLE2O FTN V3,O-P38O OPTFI 07/?8/76 0ti5 1.56,O pAqg FNCTU SUDRoUTINE FNCTUiItJ# tCT CNT oTHIS IS A NEW VERSION OF FNCTU OF TAKE INTO ACCOUNT NT4 UP TO 4 NJEW ZONES OF C1NS TRUT ION AND SOME NEW PARAMETER FNCTU C SPECS RAY STRUYK(FCT c PROGRAIiER J, GUERONFNU7 C THE URBAN INSTITUTE APRIL 76 FNCTU8 CC THE ACTUAL CALCULATION OF U IS DONE BY ACALFNT9 CC UCAL(I#JPK#P*Q#U) 'TU0- Q0 WHERE I;;HOESHULD J;DWELINO K;ZQNE P;FFRICE Q;QUAIATY UoUT,IOTY FNCTU 1 UNCHANGED COMMONS fNCTU 4 COMMON/HOSHLID/ HH(10!55) FNCTU 13 COIMMON/DWELNG/ DWL(11.,00 FNCTU l1 COMMIONIALPHAS/ A4PHAC10) FNCTU . COtlMON/PRM2/ NKHHtiDL FN4CTU 16 CO!tlMON/UF/ U(5!ij41),VSTARC50j,MV4 . FNCTU 27 COMMON/WAG ES/WT FNCTU 18 COMMN /ASISTA1CO)eA(1O)B~9I0)IMT%~)v~P~9?FNCTU 1 CCHANGED COMMONS OR NEW cOMMONS FNCTU 2 COIO/H/HYSNOE,?#Z FNCTU 21. OMQtFNCTU 22 CO.MnMON/BEXES/P7E(I0.C(O R01FNcTU 24 *COMMON/QMINMX/ 9M1N(M),MAXQ) FNCTOJ 94 .., OMMONIGAMMAS/GAMIMA1,GGAMMAAMAe A FN4cTU 26 Omo/ETA/ETA(F4CTU 27 ~.OMMON/HVCTORMAN C.~)FCT. ? TG , nmMON/DIST/ CL 80A(5),TTC5,D FNCTU 20 EMMON/UCNE D/yNCIXALPHA,XIALPHtPFRC,914 DIIXFNT 29 ImENSI 0N 0g)Ip) UFN cTU 31 NZNWNNZN FNCTU 32 *RATIOmHH(5,I)/hI(2#f? FNCTU 33 KFMTP;HH(3#11 FNCTU 3 * M*QMT(KFMTP) FNCTV 36 KO IFMTP) FNCTU 5 X(AL HA;ALPH AWKMTP~ FNCTU 38 * I(ALPHI;.'XALPHA FNVTU 39 INX=HINC() F NCTU 40 NZNWIaNZE+l FNCTU 41 *PRClzPRC(NZNWIJ) FNCTL4 4? IF(KFLGeEQ~. PRCClipRCq, FNCTU 43 lIaMDL+l FNCTU 44 * TE$T FOR NEW DWE~-Nq$ FNCTU 46 JIIFNCTU 4 IF(J;EQ.t1DLI) 'iP TO ?09 FNCTU 47 1.00 PSTAR7;DWL(6,J) VNOTU 4 * . KZONE;DWL(3,J) FNCTU 50 *CALL UCALCiJKOE1 A,OTRiX NCTU 51 UtIuj);;UX FNCTU 52 k b GO TO 999 FNICTU 54 '3*200 CONTiNU VNCTV C FOR NEW IJWELINGS IN NEW 40N~$ FNCTU 55 DO 800 L'Z;XuINzpjw -.f N9TU SUBHOUTINE FNCtU RJOEE2Q FTN V3,0-PS80 OPT-I 07/28/76 08,5o,QQ0 PAE -KZNNZEFLZ FNCTU 57 5v. STAR- PRC(KZN 1 FNCTU 50 . IF(LZ.GT.I) GO TO 300 FNCTU 59 U STAR-XALPHAO(YNC-CLMBDA(INX)#TT(INX,KZNI)/PRCI FNCTU 60 GO TO 350 61 c cAICULATION OF 0 NEW ZONE I FNCTU 62 3 00 OsTl-xALPHA*(YNC-CCEM8DA(INX)OTT(INX#KZN)) FNCTU .63 *ST2IGAMMAI#XALPHA*XiALPHHH(loi I FNCTU 65. QST'3-1. -PRC (KZN )/PRCl FNCTU 65 - :USTARZ-(UsTi-OST2*QST3)/PRCtKZN) fNCTU 66 CTHIS SECTloN CUMPARES OSTAR WITH INPUT BOUNDARIES FNCTU 67 00 IF(OSTAR,LT QMINtLz))GO TO 500 fNCTV 6z t- . .1QSTAR#T.Q GO TO 450 FNCTU 69 I:tUSTARoLEIGMAXoLz)U GO To 0oo F.NCTU 70 * QSTARaQmAxt(LZ) FNCTU 75 * G T0 6U00 FNC tU 7Z :.450 9sTARzQM FNCTU 73 s Q TO 600 pFNCTV 7'4 - 50 QSTAR-=QMIN(LZ) .FNcTU 75 *- 78 - IF(t0STAR,LT,QMt QSTARPQM FNCTU 7 -6n C014TINUE FNCTU 77 fAL4 UcAL;l1tJ4,KZN#P§TAR#Q$TARly4I - FNCTU 7B D( Q7 QLZ ;S^ FNCTU 79 ULo (LZ);TaUZPNT 84 boo CONT INUE -FNOTU 8 CC CHOSE MAXIMUM U FNCTO . 3 *F. FNZNW.EQ,3] GQ TQ 900 FNCTU 84 .- - Lml1BNZNW- . FNCT4 6 DO 850 Z1;14ml, FNCTU .86 NL;iZl+1 FNCTU 87 J, , DO 855 IZ?RNIfNZNW . FNCTU 89 IFU(DIZI .QF,UUC1Z2?l) q9 TQ 85W FNCTU 09 J VI(Z]FNVTU 90 . ;, .V;UU(IZ1) . fNCtU 91 FNCTU 92 ID(IZI)rID(i?2. FNCTO 93 ' ID(IZ2)zlSV' FNCTU 9j 855 UU( iZ2)gSv FNCTU 95 855 CONTI4UE FNCTU 96 050 QONTINUE FNCTU 97 C AFTER THIS SORT TH5 R9ATER U 1§ IN THE FIRT 9LEI1NTvQF UU ANP FNCTU 98 9 ITS ZONE IN 1S(t1 FNCTU 99 900 CONTINUE FNCTU $00 9 ' U( IbJ U V1 I FNCTU 301 DWL(i,oI)mlr F'NCTU to? * NZONiID(1) FNCTU $03 DwL(2oII);z, FNCTU 104 DWL(4. I I )}oZJLNZON) FNCTU $05 Q - NZZN;N ZE+NZON F NCTU $06 DWL(C6I0I)PRC(NZZN) FNCTU 107 DWL(3.I1);FLOAT(NZZN) FNCTU / $00 < 999 CONTINUE FNCTU 109 RETURN FNCTU 1Li E . ENU . .FCTO 11. * 0i:, SURfHUTINJ1 -SUPpLy RJOEE20 FTN V3,0-P380 OpT-l 07/28/76 08.56,00, PAGE I SUbiROUTIPJE :iJPPLY(J) SUPPLY 2 COMt1i0 /(ITT/ITER SUPPLY 3 4 COMMON/DWELNG/ DWL(11,90) * SUPPLY 4 COtlMOUN/PARA;l/ KFLG.X1,X2,X3oX4 SUPPLY 5 C0tlMOrj/PrMl2NHlHiIDL SUPPLY 6 COMOJON/PRICE/PZ(10O),PC(1O,PRC(IO) SUPPLY 7 COI.IMOlN/GA91AS/GAMMA1.GAMMA2,GAMMA3, GAMMA4 SUPPLY 8 CDONIrlON/ TiiAS/ ULETHAt1(O)#BETHA2(50),BETUIA310) SUPPLY 9 COUN40N/rvKNT/K0UNT(5).KlTER SUPPLY lO lU COtIMoN /FIX/Z(11).VAC(5) SUPPLY 11 COMIMON /SOLV/LSOLV#KSOLV SUPPLY 12 o)WL6;)WL(6,J) SUPPLY 13 IrtITER,GT.O) GO TO 8 SUPPLY 14 X4 -Xl SUPPLY 15 0 15 IF(KITER,GT,MUL)X4;XI./(1+8*KSOLVt25SLSOLV) SUPPLY 16 X6=1 , SUPPLY 17 8 XSO,. SUPPLY 18 IF((X6,EQ.1).ANI,(DWL(8,J),LE.2))X5;1. SUPPLY 19 IF((X6. EQ.O,ANDg(DWL(8,J),GE,3)!X5;., SUPPLY 20 zU IF(X5,EQ,l)X42,5qX4 SUPPLY 21 C IF OVERoCCUPI ED OR JUST OCCUPIED WITH PREVIOUS OVEROCCIPATJQNe SUPPLY .22 C X4 IS COhNTINUA LLY HALVED -- THAT IS# IN ITS USE AS THE AMOUNT SUPPLY 23 C BY WHICH PRICES ARE INCREASED, THIS PROCEDURE HAS SEEN ADOPTED SO SUPPLY 24 C THAT IF UNDEROCCUPIED, PRICES ARE DROPPED AT THE RATE XX (RgLATIVE SUPPLY 25 2) c LY LARGE), SO AS TO REACH A SOLUTION PLATFORM QUICKLY, UPONtATTR. SUPPLY 26 C ACTING ENOUGH HOUSEHOLDS, HOWEVER! THE X4 PARAMETER WMLL STOP SUPPLY 27 C SHRINKING. SUPPLY 28 C ANOTHER REASON FOR THIS PROCEDURE IS TO EtJABLE THE ALGORITHM TO SUPPLY 29 C CONVERGE IN THE NEIGHIBORHOOD OF AN EXTREMELY NARROW PAND OF POTEN. SUPPLY 30 30 C TIALLY STABLE SOLUTIONS, SUPPLY 31 SUPPLY 32 KZNv-j0WL(3sJ) SUPPLY 33 C0NlBETHA2(KZN)*BETHA3WKZN) SUPPLY 34 PZZ;PZ(KZN) SUPPLY 35 35 PCC=PC(KZN) SUPPLY 36 0 PX-PZZ+Xs SUPPLY 37 IF%LSOLV ,GT- 0) GO TO 2 SUPPLY 38 la;0 SUPPLY 39 I MI+l: SUPPLY 40 40 IF(I,GT.5) Go TO 2 SUPPLY 41 V-VACtl) SUPPLY 42 0 IF(V,EQ.O.) GO TO 2 SUPPLY 43 C V(l) 0 O IS THE SIGNAL THAT NO DWELLINGS ARE TO FROZEN AT MINIMUM, SUPPLY 44 CC OTHERWISE UP TO FIVE DWELLINGS CAN BE FROZEN AT PRICE P? DURING SUPPLY 45 45 C PHASEl OF THE ALGORITH FOR POSSIBLY FASTgR CONVERGENCE, SUPPLY 46 jF(DWL(l,J),NE,V) GO TO 1 SUPPLY 47 DWL6;pZZ SUPPLY '48 GQ1 TO 6 SUPPLY 49 2 CONTINUE SUPPLY 50 5u IF(DWL(8,J)hNE,3.)GO TO 4 SUPPLY 51 DWL6zDWL6-X4 SUPPLY 52 Xo.1. SUPPLY 53 GO TO 5 SUPPLY 54 4 IF(DWL(8,JJ,GT,2) GO TO 5 SUPPLY 55 5) O L6;DWL6+X4 SUPPLY 56 )Q SUHROUTINE SUPPLY RJOEE20 FTN V3,0P360 OPT-1 07/28/76 08,56,00, PAGE 0 X 0, TOSUPPLY 57 5 IF(DWL(3#J),,IE,I)GO TO 7 SUPPLY 58 IF(OWL6,GT,PX)GO TO 7 SIJPPLY 59 OWL6zpX SUPPLY 60 6U GO TO 6 SUPPLY 61 7 IF(DWL6,GT,PZZ) GO TO 6 SUPPLY 62 OWL6=PZZ SUPPLY 63 IF(DWL(8#J),EQ.2.5) DWL0aPZZ+X4 zUPPLY 64 6 CONTINUE SUPPLY 65 6S PWLi6, IQDWL6 SUPPLY 66 Al- iDWL 6#J)"PZZ)/PCC SUPPLY 67 Al;AtAXI( A1, . ,i SUPPLY 68 Al-AlCON14 SUPPLY 9 OWL (4 J CBETHAI(KZN) +A1) *DW. (2 .# ) SUPPLY 70 70 RExURN SUPPLY 71 SUPPLY 79 $s0 *e - - SUBROUTINE ZONE RJOEE2O FTN V3,0-P36O OPT=1 07/28/76 08,56,00, PAGg SUiRROUTINE ZONE ZONE 2 CUMMoN/HOSHLi/HHC1o. 55) ZONE 3 CUiMUON!DWELNG/ OWC11,90) ZONE 4 COMMON/KKNT/KOUNTS5)#,KITER ZONE 5 ; COMfION /ALPHAS/ALPHA(1°) ZONE 6 COMMoN/PARAM/ KFLG,Xl*X2oX3fX4 ZONE 7 COMMION/pRM2/NHH1 .ML ZONE a COv1MON/PRl6CE/PzClO)pCclo),PRC1O0) ZONE 9 0 COMIION/fETHAS/ 1ETHA3(C10),8ETiHAX10)uBETHA4.10) ZONE 10 10 COMMoN/aAMmAS/qAM"MA3.GAMMA2.GAMMA35GAMMA4 ZONE 3. CUlHM0N/fJ14Z/NHITYPS NOL)NES,NZvEpNZN ZONE 12 D COpIMON /ASSIsT/AllO),A210Q)B{1$O),iT3.4O)utii(to)#RPT(3.Q ZONE 63 COMMON/DE8UG/VC 2 ),MODP ZONE 14 COMMON /FI X/Z(C31)# VAC5 ZONE15 *COMMON /SOLv/LSOLV#SSLV ZONE 16 - COMMON/PUNCC/IPUNCHeYCL( l9) lYCU{ $t Tl T6 HHY(501 ZONE i? COMMON/FR0MZN/NBDWLZ ( 1°0) ZONE l8 DIMENSION Nfl)#Y12.3.)#X3.hMC ? ZONE 19 IFCLSOLV,Ego0) GO TO 10 ZONE 20 IF2(LSOLV,EQ9,),ANDK$O;V,E9,4)? GO TO-IQ ZONE t2 lImNHH+MDL ZONE 22 00 12 J;1tI ZONE 23 JZi,DJLC3.JI ZONE 24 PW? DL(7#j)*ZiJZ) ZONE 25 C NOTE THAT IN PHASEl OF PROGRAM DWL(.J4) 15 PRESPECIFJED, ZONE . 6 C THIS RESTR1CTION TO BE LIFTNP MV§T Q MOPIFIEQ H9R4 ANP AT Pl ZONE 27 O BOTTOM AS WELL, ZONE 28 ZONg 29 ' a0 DO . );j,iJ. ZONE 30 30 NCI)0 ' ZONE 3a M I)zO ZONE 32 Xt 11 O, ZONE 33 ZONE 04 DO 2 J;3.MPL ZON E 35 35 IF(WL(8,fJ),GE,4)1G TQ ? ZONE 36 K1;DWLC3,J) ZONE 38 Y (Kl) -YK1)*L (wLJ Jt (P WL If #) . K?I ZONE ) fNCKl)4NCXKlVt3 ZONE 39 2 CONTINUE ZONE 40 40 DO 6 J;1eNHH 1ONE 41 K4;HHt@,Ji) ZONe 42 K5.H13. J) ZONE 43 RRiDT CK5) ZONE 44 XCK4)zX.K4)RR ZONE 45 45 MIK4)Mt(KO)+3 ZONE 46 IFCK4.LE.NZE) GO TO 6 ZONE 47 Y(K4)2Y(K4)$HH8u,4)#(HH(?7#,WPZCI4)1 ZONE 49 N(K4)wNCK4)+1 ZONe 49 O CONTINUE ZONE 53 O SU M1J)-= lt1vurODp ZONE 52 A;, ZONE * 5 ZONE 53 DO 9 l1,1NZONES ZONE 54 AZAYI ONE 55 5i93;SzBN t |) ZON4E 54 4 SU8RoUTINE ZONE RJOEE2O FTN V3,0 P38O oPT:1 07/28/76 O8,56,0, pAGE YAVzA/B Z0lNE 57 Do 7 191,NZONES ZONE 58 0 -Na4DWLZ(I)=N(i) ZONE 59 - IF(N().GTo)Go To 3 ZONE 60 60 Y(JI?YAV ZONE 61 0 Go TO 4 ZONE 62 3 Y( I)-Yt )/N ( I) ZONE 63 4 IFCM(I).GTi0)GO To 7 ZONE 6j 9 Mt )1 ZONE 65 ZONE 66 (- 7 CONTINUE ZONE 67 JIl;NHH+MDL ZONE 68 DO 5 J4z,ill ZONE 69 K:24DWL(3,J) ZONE 70 79 DWLt5,J); Y(tK2)/YAV ZONE 7i IFMLSOLVILT11) qO TO 5 ZONE 72 IF(K(SSOLV,LT94) GO TQ . ZONE 73 9 \ DWL(7,J)zX(K2)/MJ) ZONE 74 5 CONT1I1UE ZONE 75 iF(LSOLVtLT4,) GO TO t1 ZONE 76 , IF I KSO)LV,LT- 4) GO TO ZONE 77 WRITEg6,1.0) ZONE 78 WRITE(16,102) ZONg 7? . - DO 3 j-1 NZCNES ZONE 80 xRJoUI)JM( ZONE S. . , TEST FOR EQUAkITY QF FINA4 HRAPIA qHAR4*TUR8TjS ZONE ,2 X0$ XR z ZONE fF(XX GT,0iO$? 2PNCHw ZONE Ai. * < WIpY(/ Z---tE a 6 -Z -14 WRITE(6,10 ~,t%)Ni.R ZONE 101 FORMAT(1Hl,?QX VACTUAL ZQNP CHARACT9RISTIjS#/) ZONE 0 * FORMAT(02tHOO EHO4P$#,4XqWV N;ff025?PV!t WFHlT*9 ZONg 8a, 15X OREL. WLTH,*/) ZONE 8- M04 FOAT(4X,I2,2Ul5,2Fl5,4 ONE 0 * RETURN ZONE 9; E ND . ZONE - . ,4 . ,-0 S . . '.. 0 ' 'I . , ' o . . 5 I " ' ' Y re SUBROUTiI'E MAXIM RJOEE2U FTN V3,U-P3I8O OPTa1, 07/28/76, 08,56,00, PAG -SUHBROUTINE MAXIM(M) MAXIM 2 COM-MON /UF/ UC!55,41)#USTARC55)sM0L1 MAXIMI COMMON /PiRl2/klHH#MOL MAXIM 4 COMMON/LJWELNG/ DWL(11,90) MAXIM 5 * COMMloN /HOSHLP/ HIIC1O,55) MAXIM * COMMON I/SOLV/LSOLV#KSOLV MAXIM 7 COMMON /FROZEN/FHEEZE(55) MAXIM8 COMMON /MULT/ DWL,(1),DW(9) MAXIM9 Ki;;O MAXIM IQ4 4. 14 DO iOO I;;iNHH MAXIM ~ ON .FREEZE(I)20, MAX IM 12 . IF(M ,EQ. MDLI) J.,? MAXIM 13 MAXIM 14 * JK;MoL+I GOT0 3 MAXIM is 1GOo3MAXIM 17 K;M- MAXIM 08. . MAXIM 19 IFCJ.LEIMjQL) I. I TOWL14) MAXIM 2 I 000~5LVY 1 qojQ 008 MAXIM IFHH6i)GeL)O TO 0,98 MAXIM 2? MAXIM 23 * . T,L) FRE~ZE(I) ;1 MAXIM ~ IFM~,LT,Ia 0OT O MAXIM ? * MAXIM 26. M4AXIM 4 THIS BIT OF LOGIC 14 TRICVI~Y $9 hl§Tgh QARPFUI4 Y MAXIM 20 &wT WE ARE TRYgN~ TO DoAODw IUTIe NR Y THg TAME MAXIM 29 W.HAV PRO ESS D ALL POUTA THE5 A ST FEWW 0WLLINGSt , A US -0 MAXIM 30 30C HAVE SELECTED HIl3HER QUALIITY'DWLLINS~ LEAVIN2 TOO MANY LOE*UL MAXIM 3 CC TYDWELLINGS VAOANT, (Tl qEFT OF THIS IS TO' DRAG THE AVERA4q MAXIM 32 a LEVEL OF PRICES DOWN PELOW, THAT 1§ AwAY FROM THE UPPER RIGHT* MAXIM a IIAND CORNER WHERE PROFIT$ AR9 MAXIMIZED,) AT THE SAME TIME, WE MAXIM 34 C WAt TO ELIMINATE# IN PHASEto THE 5TEP DY STgP PRICE COOMPETITION~ MAXIM C A OQDWELLINQS FOR ATTRACTION OF HOUSEHOLDS, WHICH NOTORIOUSLY MAXIM 36 C SLOWS DOWN CONVE-RGENCE. WE TiiEREFORE AFREEZEA HOUSE~IOLD SELECT7 MAXIM 37 C IONS WHICH ARE LESS THAN THE CURRENT PWE1~LING INDEX jp AND RE, MAXIM4 38 *C STRICT SUCH PAIR-INGS TO NO MORE THAN Jo-1 HOUSEHOLDS,. THUS, BY TH& MAXIM 396 C TIMiE THE LAST (Mo'..) DWELLINQ IS PROCE5SED. THERE MUST 09 AT IEAT AM 40 40 . C ONE HOUSEHOLU AVAII.ABLE FOR CONSIDERATION (ASSUMING NHH',Gf, Ml)I,e MAXIJ4 41 C MAXIM 42 *IN ESSENCE. WHAT WE HAVE TRIEP TO DO IS TO PAIR HOU59HOLDS ,TO MAXIM 43 C OWE LLIMIGSBY DESCENDING ORDER OF INCOME AND INITIAL SERVICES, REW MAXIM 44. C SPECTIVELY! PURI'NG PHASEI OF TH~ ALGORITHM, MAXIM 45 C. MAXIM 46 c MAXIM 4? 889 CONTINUE MAXIM 44 IF(U(I,K) .GT, USTAP(I)) IQ.?O MAXIM 49 ~Q ~JJ;HH(6#J) MAXIM SO 9 I(6,I) 4 ;J jJ MAXIM HH(7,I);OW(6J)MAXIM 53 HH(8,I)QDWL(O,J- MAXIM ORHH (8 1);;DWL( 4 J) M4AX IM 5 O5 WL(9,J)4DWL(9tq)t4,Q MAXIM SUBROUTINr MAXIM RJOEE20 FTHI V3,D-P380 OpT;I 07/28/76 08,56,00, PAGE USTAR( I )zU( I .K) MAXIM 57 Go TO iOo MAXIM 58 20 IF(HHC6,I) ,EO, J) 30,1IU MA&IM 59 30 I)WL(9#J)QDWL(9,J)-1, 0 MAXIM 60 6U V-1z.E+6 MAXIM 61 LS TAR- MMLl HAXIM MX;l MAXIM 63 IF(LSOLV .EQ, 1) GO TO 777 MAXIM 64 MXzK MAXIM 65 777 CONT IUE IMAXIM 66 00 40 LzHXtMDQl1 MAXIM 67 IF(U(I,L) tGT, V) 32,40 MAXIM 68 3? 4;U(I,L) MAXIM 69 LS TA R;L MAXIM 70 40 CONTINUE MAXIM 71 IF(LSTAR ,EU, HlIL1) q4l42 MAXIM 7f 4 LsTARzML+. MAXIM 73 42 CONTINUE MAXIM 74 HH(6,I)yLSTAR MAXIM 75 HH(7pI;0DWl(6oL$TAR) MAXIM 76 HH(8t 1;DWWL(4#,lTAR) MAXIM 77 HH(9aI )DWL0lS TAR) MAXIM 78 PWL(9 LSTAR)gPW q9fL qTAR)tX0- MAXIM 7i USTARII)(;V MAXIM 80 t'00 CONTjr!E MAXIM 81 RETUF.JN MAXIM MAXI *. 0 .,. l7 uBROuTjNE SrIPSHT RJOEE20 FTN V3,0-p380 OpTaI 07/?8/76 08,56,00, pAGE 9 SUBROUTINE SNPbHT(J) SNPSHOT 2 DINENJSICN .VV(l5) SNPSHOT 3 COMI101 /FR0ZEN/FREEZL(55) SNPSHOT i CtiMMON / ITT/ ITER SNPSHOT 5 * COMr;0r4/HOSHLIV/ HH(10,55) SNPaHOT 6 COMMON/DWELUG/ DWL(11#90) SNPSHOT 7 COMIMOIJ /UFl U(55.41)eUSTAR(55)iMlJL1 SNPSHOT a COMMON /PRm2/NHH ,MPDL SNPSHOT 9 4 COMMON /KKNT/ KOUNT t5 ) KITER SNPSHOT '0 ,q COMMON'J/PARAM/ KFLG,X1oX2,X;5X4 SNPSHOT 1t COMMON /SOLV/LSOLV,IKSOLV SNPSHOT 12 PRINT J°,KITER, ITER#LSOLVYKSOLVjJfX4 SNPSHO T 1s 1P FORtIAT1/ ,2X,6IJKITER4f13,2Xu*H ITER92l3#2XseHLS5t0 ,I3q?X, SNP5HOT 14 J6HKSOLV-13,2X 6H OWL; 3,2YX,3HX4i;Fl2q8§ SNPSHOT 15 PRI1NT 36 SNPSHOT 16 3a FORMAT(1HO,1OX,OUTILITY MATRIX*) SNPSHOT 17 PRJr T 39. (Kil;,w MflDL1) SNPSHOT 1to 39 FORMAT(7HOHH/UWL1iXN2HUtt X.11ZI SSNPSHOT , P0 200 SNPSHOT 20 * tnO 20 1SHOMU; . .NPSNOT fl t ' ' VV(K) -UI K) SNPSHOT 22 IF(VVrK) ILE, rEt$ VY(I9ip W SNPSHOT 23 -O CONTINUE SNPSHOT 25 25 IF(VSTARUSTA )sTAR SNPSHOT 25 27 - F-SIRL0 Et5) VsTARi0 - SPSHOT RX1 ;, .-., - HH60HH(I,t1RFHEZE(Ill-0. SNPSHOtT 27 27 - PRNT 4pjH6VTR(VKii4" NSiT 2 4 FORMAC i3F4,1qF#2 6,F1 ,*. SNPSHOT 29 120CONTINUE SPHT 3 *§ : - *PRI1NT 5*,l -;' !. t; 1 'DL q SNPSHOT' 31. - - ;--:. PRINT So 31 -Q FORA t 1PR C, 6) .NSHOT 32 P pRINT 5 (DWL(4, ),KI -NPSHOT - 51 F0RMAT(14X#7H$ERVjCE§F$0¢.2,* SNPSHOT 34 PRINT 52#4DWL(9,K)hKioeMDLI - . SNPSH$OT 35 * * '. 52 FORMAT!?X 9OCCUPANTS,8Fj?j4. * tNPSHQT 3 *. PRII4T 53, C4WL(8iK)tK;1,MDL) SNPSHOT 37 53 FoRMlAT(15Xq6HSTATUSf1F.-4) SNPSHOT 3 . RETURN .. BNpSHOT 39 * END * SNPSHOT 40 * 2A ." - , - - *.1 . t .' I. *.~ ** - * . ','t *. .. SUBROUTINE ERRTB RJOJEE2Q FTN V3.O-P380 OPT;1 07/28/76 08,56,00, PAGIE0 SU8ROUiTINE ERRTEO ERRTO3 2 i *C THIS~ UBROUTINE PMOCCESAN ADDITIONAL TABLE ATl THE LAST ITERATION ERRTB 3 C SPECS PROVIDED UY SUE MARSHALL ERlRTa 4 C PROCRA11MER J, GUERON NOV 75 THE URBAN INSTITUTE ERRTB S C BESIDES SOME COMMON ARRAYS THE SU8WOUTINE EPECTS THE FOLLOHINQ PARAM ERRTB 6 C TFLG s0 NOT TO CALL THE SUER ENRTB 7 C ;;I To CALL IT AT Ti ERRTB a C :2 TO CALL IT AT T2 ERRTB 9> C ERMIAT(I#J) ja ZONE S(±.) ERRTB 1a C FOR T FL Gr-1 FOR TFLG ;? ERRTO .1* c I1 11wCOME I ;1 IINCOME ERRTB 12 C =2 PRC WHITE ;2 PRC WHITE ERUiTB 13 C &3 RELATIVE PRICE s3 MEANEXNDTN ERRTLI 14 C ;4 EXPENDITURE RATIO 54 UTILIZATION ERRTB 159 q5 EXPENDITUN~ RATIO ERRTQ 16 ERRTB 17 C EXPLN EX~IUE1111 FR I.OQICAI, BREAKDOWN OF EXPEt1DIT4Fg ERRTB3 18 C HATIO0 ER RTB 1 19 ERRTO 20 COMMiON /HOSHLO/ HH(IQ 551 ERRTO 210 * . COMMON/DwELNG/ Pw 1~90) ERRTB 22 COMMON /PTO X01jjM99RO 2 COMMON /SOLV/ LSOLVhKSOLV ERIRTO 2 COMMON /JPARM/TFLGQEXPLM,ERMAT(617)qRM9AN%7I ERRTO 2 20COmMON/NHZ.' HYP NONRRTEN8 27 COMMON /PRM2'/NHH4IMDL .. - 2 DIMENSION KDW(7)RAT17) ERRT8 8d2 P IMENSIOMKN(LABK(4),A7;)o TEi(#)§El)Rl)RP1lY1? 9RTO . 9 DIMENSION LANI7)eLD?Z(5h E RRTB 9, .. 40 REjAL KDWtKZN#K74 . ERTQ 31 PATA LABI/10HINCOME #1QHFRG WUH1-T H170 0- EiF RA~TIO EI RTO 32 *DATA. LABZ/1OHINCOME #IOHPRP WHIITE jjQH= ERR--3 l oHUTILIZATUNi1QHEXP RATIQ ERRITFi 34 * IF(TFLG,EQ,O,)GO TO ?99 ERRTIE 35 B *. .,NZON;NZON~Sw1 ERRTB .... C INI IALIZE C~ NEST.R ERRYB 7 DO 0O Izi,NZON ER ~ 3 KDWil)m0 ERRTS 39 K, SUBROUTINE CROWD RjOeE2O FTN V3,0-P380 OpTl. 07/28/76 08,56100i PAGE SUBROUT114E CROW4D CROWD 2 C CR&WD 3 C THIS MODUL- IS CALLED ONCE AT THE END OF A CYCLE OF ITERAJIONS CROWD 4 C ITIS PURPOSE IS TO AGGBREGATE SUBDWELLINGS TO THE MASTER DWELLING LEVE CROWD 5 5 C SPECS FROM RAy STRUyK( CROWD 6 a PROGRAHER J, GUERO CRHOWD 7 C THE URBAN INSTITUTE MAY 76 CROWD a C CROWD 9 C THE DWL(11.J) WILL REPRESENT A MASTER DWELI1NG ID FLAG CROWD 10 * C INPIJT PARAMETERS TO THlIS MODULU ARE NCROWD OF MASTER DWEWJN S CROWD 11 C QMST ARRAY THE 0 OF EACH MASTER DWEL. CROWD 12 COMMON/BETHAS/ BETHlAI(10)BETHIA2(1O)"ETHA3(10) CROWD 13 COMION/DWELNG/ DWL(11,90. CROWD 14 * COQMON/PRICE/ P2(10) PC(1O)?pRC(jO) CROWD 15 COMMON/CRWD/ NCROW#dUMST(1O) CROWD 16 COMMON/PRM2/NHH,MDL CROWD 17 DIMENSION Q CROWD 1t F NCROWD,EQ;O) GO To 9q9 CROWD 19 C INITIALIZE SUS To ZERCI CROWD 20 * *DO 20 Lil#NCROWD CROWD 21 QSM(L);O . CROWD 22 * PS. !L)90 CROWD 23 * NSD(L);0 CROWD 24 K KZN(0);Q CROWD 25 ;0Q CONTINUE CROWD 26 ULATES SUM OF Q QROWD .7 pO 5D M4 - CROWD 2 DWL119DWL(14.M) CROWD 29 - } NX7INTtU)WL$1) .CROWD 30 3U IF(NX.EQo) qQ TO 5Q; CROWD QSM(NX).Q$M(NX)tDWL(?Mt: -ROWD 32 X t- NSBlNX)zN$8NX)t4' .CROWD 33 50 CONTINUE CROWD 34 O C ALULATES SVM OF W*P CROWD 35 * 00 70 Mzl,MDL CROWD 34 DWL11NDWL 1,iM) CROWD 37 NX;INT(DWL1.. CROWP 34 '! IFINX,EOO) 00 TO 7 . CROWP 39- WEIGHT..DWL(2fM[(QSM(NX) CROWD 40 40 * PSM(NX)zPSM(NX)t4E Ij T.pw 5§MH - CROWD 41 KZN(NX);IDWL(3#M) CROWD 42 70 CONTINUE CROWD 4 * Cp CROWD 44 WRI TE6,809) CROWD 45 PO 10Q L-1,NCROWP CROWD 46 KZz9KZN(L) CROWD 47 RATIO0(PSM(L)-Pl(KZZf)/PC(K?Z) CROWD 48 QAw BETHA1(KZZ)*tETHNA(Z? .OETd,fi KZ1?*RATI° CROWD 49 QFINPQA*QMST(L) CROWD 50 WRITE(6#801) OFJNqPSM(1)Q5M(LNSB(Li,K Z CROWD 100 QO NTINUE CROWD $2 999 CONTINUE CROWD 53 o000 FORATHlo.,Ox,VCROWPING RERULTOWD /QX,*9fIN^4tO5p#P§M d cROW- -(X *QSUl1§.3X ,SUiDWLSf#,4XstZONEo// CROWD 55 . 801 F0RMAT(5X,3(1xfF9?2)f2(9XfI5)f CROWD - 6 - I; " t -,.~ jl *SUpROUTINE CROWD RJOEE2Q FTN V3,O-P38g OPTal O7/28/76 O8.56,O0, pAGp RE TUR HI ROWP 57 oE ND CROW D 586 , s 9 . . . ,> - 9,,, . . Ž S, . 6 I.0'''':8 7 . , . . . .. . . , * * , - - o, . - ..- : . 6 - ...> . ;.:; -t. .- . . ., . * . * ,0 §VBUROUTIN~E POLIC RJOEE20 FTN V3O.0P380 OPT;I 07/?~8/76 Op,5610Q, PAGE 1 ~UGgUTNEPOLICY(KFLG) PLC *C THI~ MOHULE USE TO BE S POr-iu PruGNAM PHL ICY 39 CCC NEWPOL I CcPOLICY 4 *bC IT QOMPARES THE HH MATRIX WITH A BAsE MATRIx BOTH ARJE SIORED IN POLICY C THE MATRIX X(20i50) THE BASE IN XtIs55) WHERE I1*1.iO POLICY 6 C THlE POLICy MA,TRIX IN x(5.S5) wHERE I;1l,2o POLICY 7 *C THE TASK IS A PATCH UP JOB H4NCE IS NOT ELEGANT SORRY# POLICY 8 C SpECS R, sTRUyK PROGRAMMER J, GUERON POLICY C THE URBAN INSTITUE JUNE 76 OIY 1 SCOMMON /GRPH/NGRPHitNHG(4)fIXHHI4O,4) POLICY it COMMON /UF/ UI)Ut1(55.ilh~U$TARC§5)#MDLl PQLICY t2 COMMON eBA$EHIH/ HHBS(j1,5S)#U8S(55) POLICY 1 COMiMON /TTLE/ TITLEP(4hsTJTEB(i), POLICY $ * COMMON/HOSHL,D/ HH(19,S5s POLICY 1 C OMMON/PRN2/ NHHIeMDL * DMENSION XO 05!)# POLICY $ VT)PP(,Ii1t55 .POLICY $a 9 Q F(KFLG,N9,1) GQ TO 199 POLICY 1 C.. THIS MODULE MODIFIfli HH%4fJ) A5§1§TAI1q E4101I~WTY VARIAPQ OII& 21 * C >~~~~~8ASEQ ON A LIST OF HH 'IDENTIFICATQ NIHRSI A4RAIPOCY 2 C A MAXIMU" OF 4 RUST~ODPI 9 PD IN t4POIC DO 19? I1;tgRPH POI.ICY 24 D'O 194 J410NQ POLICY 2 IDXH.,IXHH(J#I) PWY 2 * ~~190 CONTIOIPCYU2 49? CONTINUF POLICY . 2 193 ~1til OIY POiCY POLICY 62 * HH(4sI)zjXH4jIj? $9 ONTiNUE .POLICY 33 'C POLICY 34. ~99 CoNThINU PLC DO 220 U;,H OIY 36 P ,. o 200 iB14O1 .POLICY 3 1J;J+i ., POLICY *: EQ X I eI);HHBSIOflPOLICY 0 ?0Q CONTINUE POLICY 42 UC1IdBS.POLICY 42 U(2# I );UST R )POLICY 44 442?O CONTINUE .POLICY 44 C INI'lIALIZTNQ SUMS TO ZR OIY 4 CREFERENCINQ PAT A F49$ POLICY 46 NUT 6 POLICY 47 DO ~POLICY 4 0 K 116 POIC QC (J,K) 0,0POIY 5 S A (1eK) 0.0, POLICY 52 DO 6 1- ,5 .POLIY . 00 6 J 1 .POIC 1. .-POLICY' 54 SUBROUTINE POLICY HJOEE2Q FTtj V3,U-p380 OpT=1 07/28/76 08.56,00 PAGE . J 6 S (IPJ#K) 0,0 POLICY 56 RATIO 0.0 POLICY 57 0 DO 7 K 1'5 POLICY 58 UT (K) 0,0 POLICY 59 -) DO 7 L -l2 POLICY 60 7 SU (t,K) - ,a0 POLICY 61 DO 134 J - INHHI POLICY 62 IF (X(Wl;J),NEvXltJ)) 'U TO 909 POLICY 43 i34 CON T HMUr POLICY 54 S F T53F;*G Q9:-:!- ,2. EAhCH 1 19.0 3 ii TY QRPUP P000-y 6$ - - POLICY 66 K Mx144J) POLICY 07 DO 66 L 1.2 POLICY 68 66 SU (L,K) w SU(L,K)-t EXP( U({..J) ) POLICY 69* SU(#)POLICY 70 C K REFERS TO ELIGI814ITY GROUP POLICY 70 C TRANSFER4 poS cy 76 A (JoK) ;A (1,K) + X(15qJ) POLICY 72 C EXPENDITURES POLICY 73 CC(i#K) 4 CC1I,K + W ,00MOM POLICY 74 70 CC(2,K9 ;CC22,K t X(j7#J)#Xtl4) POLICY 75 0 A(2,K) -A(2JK) t X 7.J).Xf30fJ?) - POLICY 76 C PRICES AND QUANTITY POLICY 77 - Xc3tK) CE(3.K) + X(8,4 POLICY 78 c CC(4#K * CC(4.I() + X(J.0.. P.OL ICY 7P .C(5qK) 9 CC(5ftK t $(7o4) POLICY P0 CC(6,K) ZcC6,K? + XN1T , ICY at COUNTE. PQOICY 82 A%6eK) A(6,(K).t Ito POLICY 83 C . CALCULATINQ 44 s1FT, POLICY 83 M a X91J)i POLICY 85 N ; X(19,J# , POLICY 86 POLICY 87 I F (l, EO ,N) ? O TO SOOD POL ICY 88 S(MhuNS) SM,h,& t ) POLICY 90; 9a P0P CONTINUE P. POLICY 89 -DO 501 K r 1,4 - POWICY 93 IF (K E0O4) PO TO rP POLICY 92 A(1t5l iA(l5) t A(11jK POLICY - A(#5) a A(?5) + A%?fK) POLICY 94 DO 501 H 1ef? POLICY 95 50 SU(M,5) r SUMHP5) t SOJMIN6 POLICY 96 DO 502 K z 14, POLICY. 97. POLICY 94 DO 502 h ; 1,6 POLICY 99 7p? CC(M#5) e CC(M,5) t CC(H.I) . POLICY 100 C CALCUlATIN6 AVYFRA69 POLICY 101 D PO 600 j f #b POLICY 10? . 0 IF (A(6eK),LE,0,) GO TO 0OO POLICY 103 I K Q5) M ) $i §K M011 POLICY- 304 KF CIEge5) A(1#1K) A(K POK) / A*Ki * A..i? ? POLICY $05 A(?2K) A(2 PK) / A(6 #K) POLICY 106t A(3,K) 3- 4000,*(CC2,K) t CCQ(t#K)l / GP39) POLICY 107 A*#,K) 100,0*((CC(4.Kj CSs3iK)) / MI(3iM) pot?CY 308 . , o .c A(5'(K) s00,o*((CC(6,K) v I Ct05t(} POLICY $09 . .*,. .. ..' -. 19 UTM(K P S(2LICY I / SU(4F(? r XsO?IlOO POLICY 110 ; ;- . * 5UI3ROUTINE POLICY RJOEE20 FTN V3 O-P8oD OPTal 017/28/76 08.56.9,Qot ~ ~0D cONTINUE POLICY 11II *Do 601. K ;1.5 -PO'..Icy 112 601 PPC(K) 100,O.(( (,01*A(3#K) 1.)f (,O1*A(4,K) + 1,.) . 1, PoLicy 11a IF(A(lo2),EQ.0.) GO TO 862 POLICY 114 Ilq1 IF (AC6#2).GT,0.0) RATIO A(2,0) ACI#?) POLICY 1156 862 DO 6171< * 1,5 POLICY 116 IF CA(6tK),LE,0.) QO TO 617 POLICY 117 * ..PO 616 M P ,6 POLICY l18e 016 CC(M,K) aCC(M1I<) /A(61K) POLICY 119 0?0o7 CON~TINUE POLICY .190 * c OUTPUT POLICY 121 WRITE (NUTt105) TITLEB POLICY 122 a'.WRITE (NUT,106) TITLEP POLICY 1-24 I1(F 1GE-i WNTHII~ POLICY 124 DO0 875 K;INHPOLICY 1R5 WRITE(6,120) (X((IOu,sIoit0*),MIH(4IK) POLICY 126 O075 cONT 1NUE POLICY' .127 I 20 FORMATIXF 0F1 F POLICY 18 T(1X3*F742,X,F1O5 50j03r5QFO,.1, F,,POLICY *29 FE4iI/9-HlQ VP ORM~p py GOl;!RoI. CARP #/ OIY 130 10 5 FORMAT(1H4,00ASE VASIQEic * FORMATC1fIH0 POL joy 1§ ; #14A101/) POLICY 131 WRITE (NUT.107) POLICY 133 *07 FORMAT (1O!GROUP AV# V,IAQ YCHANP~ AV,CHAN9 AVpANPLY 134 IgE NUMBER _ARMKGO/1I TAX ~X~P~~PP~ 9VOATIT POL ICY 135 PRIcE IGRU RA11/ POLICY 1?6 CY~79 K NGOP RAW1/ POLICY 137 KPOLICY 138 *. IF (K.EQ,2) GO TO 709 POLIPY 139 WRT NT1O)4 AJ~~i~;~POLICY 140 oGO TO 7O5 POLICy 141 o799 WRITE (NUT,108) I.., 1 Ii iq60 RATIO POLICY 142 705 CONTINUE .POLICY 143 WRITE (NUT0751) POLICY 144 *1 75* FORiMAT (IHC.' IIA N TL!y9A.pIP HN.I POLICY 4 DO 710 K ; 1 POLICY 146 L0K POICY 147 * < . .7*0 WRIT (NIUT111B) Li UT(K)t PPP(K? POLICY 148 lip ORMAT (1H0,I8,F1O?.?Ol(.F792) POLICY 149 I04 ORMAT (HPOLICY 150 .WRITE NT6O POLICY 151 .* 630 . .FoRmAT (IHOq6Xfi BASE 9 POLICY 9 RASE P POLICY Pl~ POLICY 152 PJO 6j, K w 1~5 POLICY 153 L;K .POLIQY 154 63j WRIlE (NUT,050) Li (cC(M,i%), o F36 POLICY *5o 050 FQRMAT (lse).142XF,) OIY 5 9WRITINe THiE MIORATION MATRIX POLICY 157 DO 999 K ;;1.5 POLICY 158 L;K POLICY 5 16DWRITE(NUToIO9 L POLICY 16 1 F0RMAT (iN04//q' M1GRATIQN MOVEI1NT$ FOR QRQUPII,* I PLC 16* WRITE (NUT#11? POIC 164 41FORMAT C11H00//.20X,0N POLICY 1 2 3 4 5g) 'POLICY 16 PO 777 M II *15 OL * VJrROUTINE POLICY RJOEE20 FTN V3,U-P380 OPT;1 07/?8/76 08,56,00, PAGE 4 777 WRITE (NUT,112) ti, (S(M,pNtK)# N 1#5) POLICY 166 * 11 FORMAT (lHO,I3,F14.0,4F5,O) POLICY 167 999 CONTINUE POLICY 168 QO TO 990 POLICY 169E * 70J 909 WRITE CNUT,910) POLICY 170 910 FORMAT (lHO.' HOUSEHOLDS DONT .MATCH#' POLICY 171 990 CoNTINUE POLICY 172 RETURN POLICY 174 E NP POLICY 174 52 4. SAMPLE PROGRAM 4' .v ; ' ( * * .. ;. ~ ~ ,' * c;vX; :J *L' r, . 3 {. si,,i a L .8 i isoN -.i t t'.''b , &, FllH -d4 'o 4 fC ..'. U§ ' Sga +w #0 t > #v- tf W l j ,y,*,.,j,*..wX,' i'"!s ,,.* ''ftu'Ai'. -, I n:* i4r - ; - l,, -: . ,C':, ' ; ,', , Iu; . - . ,- 00. - d 1.r S > . .. 9.b .: A .x . ,.. S I0 I . .... S.i * . * t 4 X . o . . '60 3sd . I v y I,g - \ ' +8 .' i7Sv . W {) ? . . iflhfl & S0.04-p0 ~ w TilldN6 - . '. .F' s''.il-'rl ''1. - ,' ,"- . . b8t;* o6 ibo 'o j' t$t ts8a It+8 1it+gg8tJ';. 0.';, ... ar;, *6& ,idFd .dd iz d6 oo j* b ii e* ;.;.." a' o^ DI& io' ' 004'61 '00tib *+o Od di',4,' , 4i t' to. .F O a 10 0 d o 0. 'O ,,, 3 S 0 sodo§ to| so o eb' 0 o0 0 04 '°'O 0o ,4- *0" ^,.,, g..rlr.. ' odewo'i o' ' ' ' ' bo.' .\..t§! l6&O ~~ Oo 004&G b a oN| t+ztn o#X0 0043t90 0o^ ^.0 039 :. io' O@ Oo eao 'oi je+oni 00+3i4a ggrZ5' 0o+3aG'd ioo00 &3 ' i' ,6-1,,, Wi M- IN 0Ff0'00'tO40 i' Of3' '03t'° 'tO 3t' to.3° T09iU jU'^° .IU~ -'3q . .... U; * 'Okk Ai 'oo6 I6'o j6b Eo4!o 'tdlt' 'T0o+eso no.3U'0 0sLt Ul '&6 td *Xo' oo oo h01o DD6'0 'U+ti 0 @00+36 iO ' ±v1416,, r 0010'o eo 'o'0 'o0 *b'o '0'0 8003'0 *0o o'0 'ao t0+bo' F. Ft1§.-,j,! * o'ln'o 0n'0 0'O 'o'o 'o00 'o 0063'(i 'o006+'0 @0+5' oo+i35.'10 i~ EHl I' ., ,, .,, t looO096 00C O 0IO'0-ODO~0~tH~ 01' < F 1 ~~~I - .%;* "''4 o i i ! 0 '2 ,'' *- TO un @o 1 Z' 0 'TO+3V 00 o+399 10.ev8" , x * 4. W ¾ i :to o o I , n ; l-5! ', a A"Nt 'A 'lP t ljD|u|t swwEda t hr.E V " , at,* !,Yp t it . S; t, 4 .p .. - ' X. , ;: ',, ' ' ',i i' v'-E i ,, . , * ., , , ,. , . ,, . ., . .'; .; '' A -. ,, : ; l ' . , x - - it\,'- ;T43' X09" '7 3' ' ' ' ' O " '.+3 'O 'TO3' T '3 T 'T*V "',O4: ''' : X X, X :' 1 ; ,' ' S * . , '2' 0 5- -E. - . : l' ; i ; o ' e$I 0 - I" , .'''.''''' .'' ; s. '. ,'' :;' : 0:E@M '. .-.,.,. .., . , .' * S3 ; Sn' 'o' .o' -o' .oo on oo oo e d oo - il: It 6'Nndi ". * t . *, i W~~~~~ 4..X;n' f ' S:f 00l We ! * , ; , F * ' ; * * * 8 ,* . ' ;*' , ', I.§* 6 S ;0 W ; 0S t6W loo 8 '0 '04-t ' 00 190 0 #0*9L' 'WItOdO SA 00 ',4it ' S , !, r'0' ; . , e -t~ E- 9 \ '. , s : ' 4 -in' . \ 4 ', . :....; .; £ Q . . 0+3£L6 0 £+395tt 0° .0+t39 @ 0'0 ' 003 9 t 6;0 i s ," : tf+3zt&' *of8fsolo o #1!O+j09V!0 'CO+1T9UZ'd 'W19g0990 '20+39TWO' #C0+J60940.iO_ itR3 ,..w . .,f i, , ; :. . . C . ;.u27i' ,;' J ,1. L, ,'gv 1,!- r8 W 1 z <"w . . o C o . o . °*t 0 ! . 2 . ; i 3 i - JC.J. ;C *' ,. C4 C A4 ' . 'o o ' . t tN ^ w 4W . IL OCf.A_ 'IC' -;\.'A h tou3zN 0i i,'lXd / s - :I . 7i ) ;. 'o' 'oc .no o oo 'oo 'o.t z0360 2o1t-oi § a b C & . 4 ;. C .2i7 1.8 5 .1J.I'd ,AJ I lLid HJVJ S,dIIUd!J h10HYSU114 90N .o'adolAd.o 00's CI 00'? oo'0t O0 I 0000 D 0I 'T 00 t oo'-0 ' oo'p. 0 Oo'~ 00'?0 00,?8 oo'fr~ 00o 00 00 Igi 00' 00' 00? 0'. 0 9 I oo oa' oo'o O'o oo'Z oao ao ool' 0 60- 00'? OBI?9 flo'T9 A 00 o ao 0 o '0 n0l at 00' O0'n0'?1 00'Tt T 00's 00'1 0'6 60'96 I1 001, 00's 00'9?9 c0'9TT E? o 0 t'9 ol'9 *00'09 Oot q l * 00V0 9'~ 00 0* 00 . Li ? oQ19? 9L't9 Oaa O...t I GE, g96t 3o9'9? 00 0 809T9 *99t L-'91~ 00O 1,O gl9 gt z 'q99 9'6T&' 00 0. 90gA i T * .I a ' -6 *- L8 .* 0'a . - 6II : z I T 9 o'909 £T'65a oo'a OSot 0'oer 000*o f o p.oo't t'0lle i9VgL oo'o 19? O * O' Lt'965 oof'q 00'0 td'Tf .e 9 00'? qg'6g~ ciga O~ I aAt a 00'? E6ga W119t 000 ,. 6'LS Coll? G?'09 919'M 00,0 19'IOt I 2 * 00'9 W6'401 L9'U6 0 0O L9VA0 0' O9 6 T 26'?96t O0' 0 U 86t.~.. 00p jam * * i 6 t '9 ' 9 9? m '. a ~ 'r~ frt 'O * la * Io lot lop :fl02A j W W Ao 'o V -0 W' @. NM CD 40 s '31 CD NF M C. CJ I - Vl 5 020 oo 0D o I C O, C: 00 I CD C IC. C0 II C o A C 00 O IS 0'0 -, - .00 0 0 . * S 0 0 00 00 00 00 1 5 00 IiJ 00 * - 00 o00 * I ,0-0 00 0-0 - S00 -, ,. .00 -. ,0 5.,t.~ . . --, :i , i 0° ld 16 9f . ' e9 OT ') * OO' u U J ; - *-* A'tf o't ;£9j9 iZ 0't't- or *o 00 Iz 99f Oi otl OZ S LyZ OD, e iO t ';.,.-i ' ab l .C £Z4ZO' Z t Oig 'Zol QOi £ If t-S t; oo 60 C6 . . ilJ dJa^ ftn i f o J g1taf6 , ''''' * s o0l 0d e0'.ttl. -*.,.', ...; - . . t'69t o t . . SNd ,' to, '@0,06I190 goT -to'o dO t§& oO' Z t 4L g oodts ad¢ 0£& 6oo af; 0 00,, 0 'L 0 to999 s99 9O'69 o' *L4 000 O'~0 t5r'9Z 96'...... Lc'rO 00...t'9t oO't *oL . o 0000 6£6t 009£ a n o *o0lo 000 0 0'8g 9 9'6£t if! I£'Z0;t O09 55l9 00c0 0&L9 IZiL4t tTL'9t O0'* o - 99'9 00'O ££^96? L£'tg? ZXZ *tg'? o0't oO'9 * qo6 o0o t£ZZ 5L'iO! zn' £. LL'? 0O't' . f00' ..: oZ6 00' a9' t' L9'Z£t T98 rt 00'a no9g *o oo '6 'L9Z 00b'1 90'L 0OO9b OOt .0T , 0 L OL'6 000 l0o19?D t lI 00' ' 9OaZ 0 ,t noI'S L 6 00 0 lo0LA IV9t tfof IC, L 00901 . . LLXo 0 o ,A tot 05 'lt£ oo6~ i 90t ooa 19 . .0 15 0 @o'ot oo'o . £9'LZf - 9'e tt Z Oo'o -no sh'0? oo'o. @A'69fr 9L5z R'S ot?s O-d't . oo'AT, ; 0o 10 0 5 Zoo t . oltf - @960? 00' ooo '089 269LOL t£o'Z0e oo'tt 'oo'Z i'@; CXT t@0? 00 10 B '09a& To'99 O z' P oo go OL DI ,g , .', t 00'06, 96'4G W?L0 £9 I I a 00'0 . ? 0 T * I'il 1iiNvi ~ Ni 40J OD, 50 a 0 M6 halT~i~ *~~ -dO~f1d ItJ3VA . 4 0 .i6 ~ 4 eg glo oi 6I 0d0d N6 *o 0 od0d iaYo A * flOnOo ienr 00~~ ~ ~ J,0r 0:0 3Ao oa 09 snivis AlNVdrJnOo Ag sEmItr11"d . tol ,0 a ,16 6t6 t 9 c o S NOIlVS6l A OW .; b Iot os o nL re d 0 a, gtg lltt l ta I~SfO add NI - w w w w w w a IF -< -( 0 ow -zo m O't- liA.% ~- !, :`P In m X l -C1 1U= 2 -C -C,- -C -C -C -C I- -) .* 40,= C)C C mCDC In m mA% oo o o o oo0 - CC CC C C C C C C C C , ~ C ~ ~ C 0. I C .L C - -: CC ~ l C I Tn 0z m, CD cw 4A %vn % n z M Im ,u . -cl - . . 10 CD t-. . . 1 ozu V .; :4 C ck 3m-o CD c ew ::o CD CD CD ) ul C: C, C )I C -C Z * - - Cjc - >- - :0 2 -C >lt C D t> .UI ZC O 0 0 oOOOOO; -I IO ULT C 1 CAo C) o s D ve¢D O <> ¢ " = O = 1C io C {h * S> -r ` c It- C- w to - C) 3 I Z 2 Ol ol tDs .0* Utx 0 o ................. C° ° 4 '0 n - - - - - - Z . Ow C4b OVrt. CD Z ir C D -0 - A0 'Z VI o trUnuvO tD on#1 %AgtDEt % _-J.1 1 - -4 CtD U1 O- CD C0 m m +- m n p. -ft - C 3 C rot CAl - %r OCD r, CD -4 *0 I., O CD CD CD - tCY, O\) J0 0. N :tD C -C - * .o D ooD a, cwr C C) C 111 CJ %il CDCDCDCD CC O C) C:) C- ,4 4.6 "to -4.- T CDC o CD o C- c C C C C C C-- C)~ C:. -4 , 4CP , CD 0o -' AO -, -1 O 1 t0 ..r 4 lc -o -at co om OM o. AA o O J t~ E 4C l4n lo C,l -J * *. * --- 40 C C C 1 C CD CD - - - - - > Cw0o0s 0 Cm B C) Z N Ca U IN) Oc*. 'C - 4.'4.6 4-# 4.'C - - IC = 'C C0 CD~ O C C C, = 4,9 O e * 0 *** * ** ou- oe OT - 0 - L9 9 O 00 a 00'T -0'&h &S'99 00l '69 OTZ 00 o a Oa 0'0 It 0 0 0' i 96 9gt 90ot §N'60Z Ot'T 0a'£ o t LN'0r dtO' Md * NidM0Oij i hoHd AYNh1iiO o d ANOZ ^H * t~ t l3Md 1 AO QL 009b 9 60' Z SNVYW * o0o $)0 IoTo 006 5't OO'Ob hO'T 0 l00o0 £99 'tt 00'6£ o0 . E5'L 00'0 os'got To£Yt 00g£ nd' 0I ;LIL 00'0 .L'Idt 65' St 006tL . 000 9'2T G25't 00'96 04 §0'@ 00'0 09, §t 9e'6£t L5'QZ 0£' 005S no0's 0I o00 0 05'it 5Z'LSt FtW 65'Z oO'T OO'L * g*§ 00'0 £ 6 6'Tt ' OOt no0g LL'g OOO w8at 4t'L0t ti'99 OI ooz 00'6 @ l Colo 00'0 '9I L9M,£Z 6t' TT'z -00'£ od'ot I Et, 00'0 08'b5Z TT'a OQ'£ O'Tt *oZo avi' OO'1 oo'Zt C 00 0 0*0'9t£ .t'Ot£ 00'fI 0t'9 00'9T 0t Lfr' 00'0 g5'Z£ 051e?£ L'Z5 00o9 o0' T oZ Olt po.Tt 0- .800 00'T t 09l6 00'0 it 0bt t Ztb O' 00't 00'95 o0os9I Q616 00 §9LZfr L9'gtf 85t 5s OQ'L ooiL I Colo frP'64 L0'00 Tg fl IT 9o'ot Colo L0'lg0 00490 0L199 ?9Z5 009 oO'8t tz'ot oo'o 90&690 §L'99 W§Z9 6£'Z 00'6oQ6 9£&0t 00'0 614Z4 96'TLS LVE 00'6T 00nzI 40* t 00'0 'C9 t9'LOt 9W'tL W OO'Ot 'It9 Wot 00'0 8qdw WM * 90'69, t' OO'Z 00o'A I00 *t ' titf O0NYJ * ~ t~ ON~Il 06 iN1iM £(tL'91i9 rt 0 0'£ a'44 0t 0t t&6t -6Z6 Q5'10HZf0- AtcO'S Ot ' 0d INY -. 0 odtd IN dA e 9~ 0F 0J 10~ ho . szlosn ojn rA * f0l A0Ndl3 ~ 9~t1M * I~) 1~oj4ud INSYA sni/is AomvdnIon As SoNIllgMd I.. ddo bo o0 ooL4 -6 go 360 la Doon i0 .t sn d NVJW 0 NYJW ON! Wgd nsHn dA1 idAi As a1oH3Sn10H 3181b91i a6h 8 hJW 0NVgW gh1H N n dAl Cdbb't d.Ob6got ddAf't doooWi'i tno'£t bo6Zltgb 0 too ag t oOt 6.4 I tf 601o'lz 29 '6VZ Et-zi . 05t 99 t g 0000 5t A9 02 dt t 6ttglolit oc 5 'c10 t9 t agev' e g129^ i3i dXi NV4W 4 NOlW b NVt31 3N I 18d 6NISns "0 3NOt Ad Wdo00t ilSogi ilb101 oo 1llXo 0z !10lO)Nvi5A 00-' 0 00i6 S£'T t 0 ' n0C 000 tO'£t 9T£'T 00' 00 o IBr 00' 5 t S t oo i G oF Col06'0 Stl' 5£' 005 00'Lt 00I I GI OZ 5£ 'T 0rO2 Do 0jo 00 '~ at 'z 00 'IV n o I 9z 00'0 00F b' 0' no'5 IG e 00'- O Z'Tt OO't WO g 6 5' '2 oo0t ooT OO'T XTVZ8 Gait 91o9'OZ $1Ia I' OOT dole 00'- O' aot' og - £ 0O'T Z5'9Z1 OO' Z' 904Z OO't- O' - OO'- 03tQt OD,? 'TF r'5 f8 gO'Z 00' 00'- 95'0£ OO'Z ££'Z 75 '; OO't O& no::: O0T9 00'2 t - 95I e; 9GI 0aI od. - *&N hi .-ZI .V 6i I 0Oa o aIa 6 00a&, 6 0igo 0 0't 0O' 0'0 g6 4 00o I 0 a I 4 o't amid gi A~tJid hi o . ) o IOIm to o 'a W-0t za. o o' no"t0 - ~q 0'00O'6c 96 L 00'0 66,11 WToo'9gO0 it 0006 £'0 O' O 00'0 ogw 4 9 ~98 s 00'~t~ no L 00'0 aa L9?91 go'9a .il oo't 01 go 0060 4L'L0tl 9L10 o06i 00O 06 S 00 Do'tv 0016 o oo'11 0060 11'96c .t t 0olo Wa 00'ga . 80~9T - 0 0 - L'f 06'frtol Coll 0061t 00'0 OO19 * 1't 00o0z0o0'6 2O06 ao06' do Tt 00 * 9't 0' . 0 c0o 196tOt 90'EL w oiot 00ll 0009906L,. 6it U16'a-6 o at '"g 00O 66, L,i6 0g d0'lc oit 0'Ott alto 01g L96t' 00blN de'I't * ALIIIjA 0Nvd ONtIflh d DNIiP%1a nii-m- * 0odid 1NnVA 0 d3d INI)VA *10 OdDd INOYA 0 dnon HIAO snivis AoJNvdnolo 0 1 SDJNt1IIM0 * 1Jgril1w 9r4!sn0H p -i al m 1 n. . 0 1 tte C 009CLet 0 ol6661~0 roolst 0ge 00*0 * . .I .08o66o 0I L S! 000 ~ O9~O 0 0fr v - OO0' 1 f s90Do 9 ooWl~~ 000Zaoa I . 00OO0 00 O oOO 0O00 005t16 ozg;0 00: oo olo: 0A. S. v 0o0. 00f1'z 0OiJ'fr1 ~ 0'91 ooo0izf f a * 0 90'rfr OO OO'6t 06 009 0'0 (~.0O DO' 96fr T9 OO9oO' 0o' 000 0I'01 bO'Z ; oa.#6L oo': Wt68 0: O' oO0Zt. * 00" . 96'96 00' la'9t Z£iZ 00'£ 0no'T, 004^ Lob5 ooZ 12tL'6tz ta z* o' nof et :o 00 I t5oEz o 66 19£ ES TIz aO* no' ^ 6NI0MO it J08d ONVAJ360 6 d 1NffZ lbihn §S5Nt1113Md AG ti& 0o'o 69,09 V0'W 99'+ 00'O tY'L9 00'T5 00'6 0NOQ'O no't t 00*0 WL4 AT9'0T 5£' 00'62 n ' l00o0 0'S0GT Zf'99 ?6ZT 5£'t 00'- 00' £ 00is o'o- r'. 00L 96 O'o . Tot 66'Lt 00'9£ .-. °O5 00'0 05'8s5 9960 r08 5£' OO'£ 0019 fZ§ oolo 05*5tt C-6'IL5 tgiZZ SS' o I0v or OO'L 6£W 000 . £'6 L' RO'9Z S£'T UO£¢ d0'R Z9'O 00a0 £'oLS 9t 00'1 00'6 08 0 9.,8 O 00 000t 50W6 00m0 906t t0 tD0o Z6 00 t ZO? 9't t 'Z 0 0' 00'ZT W~6 00#0 WM A':01p ooT' W 0090 *$ 00'0 OO'0t£ 0L'9t. 00I't O6Z 00b9 tlO'£t 9966 00'0 0.Or 0g ' *5 9t''t 00 9 . no' LG6 oolo 5Z'OL£ OS'§L£ OO'Zf - 06'8 . 00'5- boDt 996 oao REIM 00 00'0' Fao o' L&A te48 96'9 tole OO'L no'et- * LoOt 0000 cz'L8L O40Z9 £'ZL 0040 no6t it'ot 0 0 o O6Of 064t TE WE 0'016 0f'6 co9ot 00 0 Z6iz 6'?/A OVIOL We' 00'tT oo * qj't 00'0 5Z 029 * 9'LO .L 9Z £fr' 0O'0 no'tz * io~t oo'o e6'$ot tsq'Z6t 0t9'1: OV£T fl0£s -O ott 00 0 toL6tt 90946t *A'ttg ZT'Z 00'Nt 6O' odm0H§,loH Ad O -bald l0 N A a nod tJ h5VA 0 i I.) n *1 dflO3n isnh? o d?1of s3A0 *snivis -AONVdf 0A SJt1IO 1rdnOjo AR tDNtW If T9 SN(JvHlt89tJ .1.. .~ .. *~ .,I * ..* .,, 6 * OOoO'O o d6 boio'0 0 o 'Ot l§h NV Ni3 0 tw OhtN8 ONIHn3 3A * , ~i ~ t 6 *NW4 *N~n %,dA* i A . .o. . . A. o.sn6 . . .. . :; 4dll aI0HjS60W AS d ,?~ > ,.~ .NI4e 3N.i. ~ * . '. 6Od8it fltt t£t ' tt * oooo't s£ tto ' O9tO OoOQ '*# t£0 't000'fr' £ O~9£ . OOTat '60'6t Z 9T6LAtt Zf8'Z98T 0Z9'Z98T. * OOOOo t 'g 2I tt?tOv tt86'89 zzza' otg Os5' z - 86OO t . gtjB 2XXz. b£F't £8'~ F LZ ..TB;9Z.t-+ *li Ht idd dX d NVoU 6 * .NiWHd DNt N ,N0Z * N0Z Ad a £^6 tVo? 1IubuJ i~iioi 0o ' 1Yi0i,M!N 00, 'O ljoj/pJ4v3yA. no'o LZ'Ot §£t woa bQ'Nf 00'0 Z6'ZT 5£I'T Qd' 00b .; '00'0 oO1 50S'T 00a9 COLO oo 0i0 5 WT 2T 00'5 60iGt 00'0 WUZ 5£'t 00'0 00'0 z'Lt G'JT 001 hi'9t 00'0 oO'§b GE'Z 00 nt *oo'o oo't'0 WI oo'~ o 00'- 2T'FZ 00'Z tL '5Z, 9LI7 00'T 0 o't .. *00D'- IZbO' OO'Z 10'T tO'- G0T X9'h! OO'Z 9'b 8@ OO't 0O'~ 001 00'9'9 3" O00 00 o0'- 9e'vb 00' 00'00 06'Z 9OQZ oo'0 Z'tT 0O'Z ZT'55 981 00-'t no0 - oo T,s no a 19Z - oolz Po 'L R 0 O' L~0 0 2 0 0I 00' 00, U6 on *4 . i' o'o. 109 o0 NVI 0O0 it £'Oto '6 00 gt @1'i00 k'rO 9I0 00'L 0 aoNV 090O' 0 2'96? 19'19 909 a ?001 fO 0g1i 00O 2'9,0r 9L09 I' T'e 00's 0'0 *.1,@1P6 000-aoOIa V 0'9f D0t 'Ilo IP n OLI0010 . 91 . £99b PO T ti 0 o 'L . no to 9600 00OD'o L?'91' *. 0(t9eltT 6'T 5' oo'e noo U it.4' 6?0?I ge C 00'0 O I Oa I 'iu 00' nO 6? re0 o&lgoINWA * o: 0lot NCW PSN01IV9Jj 0 I4 * fi 0 1?O0 ~ NS06' -~O i' ' ' 0 -.N A .otizf 96.;I Ot6 tg * o EST tl 18ol dlM9nO !n 3!G - ' NN' *: ..,n . . 0;:,^ oDbD 'o . ozO'U t too 6o ooto Oof' E9 £2^ 6P *5E .) ..A 164n : dkg WN3 b-v4 Nk3u ONIWHd ONmbhl .- 90W> .iou .00 W ldiOiMt a01oH 0 ji9 .N i YA' cold 00T91L6 00?6 1 - ;; . oa ostfl . ooVts4 Oo0W NV stsW oo1I4H£ NIi. ** -^dAi "t j-O' .. 6O ' ' - 2,' ., , 0000 .D , * ' o et, 'd o T - t 00ooo; -ol 66 0 0t O t 00 5 no 9t : . 00'0 WZbZ 5£'t 00'4 . 0oo,. ,,, . w00*0 WUb5 g 'T 0 0000 o:tt 0000a?a 009'~ W0 09 0 :e800'0 66 - 2 . v6 1 . 00 Q 1.17 31 )h;., ., 4 4 1.0000 2 9532 2 2 1,0000 6266 a 1.0000 '1510-; r .* , I. ' 7... ..,,* ' > . i. --. 0 * ; . 9- .....r ., ,,5 . : !; .7 - .. . ..... f. * ~f U- ., S , ,;.7. 0, .7 ; . 1 .> , 7 . - . . . .~ 7. ~,* * * .; .;' 1 ;'' X ' " ' . - '' ' ,' ' ' ' ' * " . S ' ' ' S ., ; ',' V'' i ' . ' ' . *' $* * ';i' ,';- *.%v gali ttoot ggdo oodt I 1 62 o D o 00lt 9 toto 00oo9t ST It - oo o'o 'i' o t¶f o Vo IOt 00o41: e 0' I 'f.tt 61 6c gAla a t oo 004 0u0 A ur. u0u z. .i 0 . 0'091 Q04oL* '4 t@' 00041w' 6$t 00' OD SO*JL 4 .'. A 9O k:''~ 60O- t fr9106't 'tf 00b,0 I t? ~ ~ 'd 0* -UOU ?~fL9 to 99f9'a 't Oao f~ ,L 9 ~ ~ . 5dZo'9 ~ 6k 69Tg O~ oo' L ' 9 As6 0'L1 ~f 00t 01 S 't o~~ 490t 010 7 69 Do 4 0e I nI%S L otLfL c 0~f too'~t 'it £?'1 T~ s 000,0 IT ~ L4~6~': OS? 91:fT ' frf1' co~o 'r -I t6t 9V9 LI~t 9'co O'o I t , T o 0000 it 't - I ' Aiin O OY'H 0oN 1fl9 SY (I 6H elM i0o 'NQYlot coo Io ~JItVAUYH o~N10 - bNv"$ 961d'1v0 - mititin Ni. j6NvI4 M W A .N *iA !DVSA . YH' A n of ~ o'A o0g? , 05'$ 8I'L 6'0L S - *0 2 1 o o* A . o ' t~ ad~~6 ~ t't ~ f~I L'?. d' 00s *t ooo'o It * 00t * 00f~ 00L9?9 U 0' 9 f0fr'L~ * 0'L C0L 2*'f go a~'9 000! a0 - ' - " I ? *E . -o' 9LI 0~' 000'0 ; c0t *e Eg~~9'Qf o I O9 igito'e It 000' I - It £Q rLh 16? Ole ~~'~ ' - 000 '09 9OZI0't *o 0'aO 9 9 6'90i gI 000'D0 ' 6'L Ole 002 Olt iAoIlOi oaugT St 00li 90t' AoI0t It p6w * -e la~ -L 0'7it Ltl ai -s T 0 s0 n C~ jj si~V ," " , ' * +' . ? - t .- * - 3,02 aN3- v l,'t, ,.9 ., ,I nU. ::o o :o io *~00 dng o O 'oAO Nd 1d - . tW -00 36Z64 9h3 lo t * . .. . : * . *; sI . -7 ,0 ;W TO '" T *R ' 49 9' 6. *0 f0 ' -0si I 'C 'Q ' '9 'I T 0 I AO 1 1'0 '0 0A lo ' a 1 ye 'S -'0 '0 'o '0 -*- - * ;;'S * ~T d(Th8D S°9 SIN3W"A^I4 ~0ti~11t,tW *6toi'O ^6~T !b69'09 tAld^8 * -l--" -- ~ 1 er i. ;.s .- ' 4 ,!w7, c -r a l . , . - .-,. . *;d 1.. ;. . 5 ;.v . . ,'O ' d ' .t ': . - Z ' '~~~. :I .0 ,~ - ' ' 0> lb 10.1 N-d O 16 1 dA6 g '.' ' N *o' , C A O . , .... - .g. ;.''* ,n ;-.S d d3 L l3 , - h s i t1I Nt . I - e -, *~ ; I.. . 171 - dnom! WAs §lAiSv3A,ov hIIeL V ;6 -o io Io do ;0 '0 .'. i 9a ot il ; 1NW Aw *jOLV bt ~Ie, -0 , o o - - 0o ad o 'o 6O t' 0* - . -- .t . -..r.: .% ' . V . ; ''"} I , .I - n ,' . 'F -r A" W; fS00C-t i ; ., . ., ,*,:.. . . . .. , :* ; * ; 8 , 1t r' j..! 2 'l' ,~ ' " ' -''0'1h' ,.,j f e . g' . . -- ... . . i. 7 .;IVN' . ' '"' ' . , gn'N | k, IN . s v -n I. .. ....