This project has involved development of a hydraulic fracture simulator which calculates fracture height as a function of distance from the wellbore in a situation in which a payzone is bounded by two zones in which the minimum in-situ stress is higher (the fracture is vertical). The fracture must be highly elongated (length/height ratio approximately greater than 4) and variations in elastic modulus across zones are ignored. First, we describe the leakoff and spurt loss calculations employed in the modeling. Second, we discuss a revised version of the vertically symmetric simulator (bounding zone stresses equal). The addition of non-Newtonian flow and leakoff (including spurt loss) is described in detail. An illustrative result is given. Third, we describe in detail the vertically asymmetric simulator (bounding zone stresses not equal). To illustrate the last results, we present design calculations for a 30,000 gallon fracture, which was the first stimulation in the Multi-Well Experiment. The 80 ft fracture interval in the Paludal zone has at its upper edge a 520 psi stress contrast, and at its lower edge a 1195 psi contrast. Computed fracture height growth above and below the perforated interval, bottomhole pressure, and width profiles in vertical sections are displayed. Comparison is made with diagnostic measurements of fracture length, height, and bottomhole pressure. The appropriate computer codes are included in this report. 21 references, 11 figures, 4 tables.