Methods are presented to generate isothermal and nonisothermal pressure-transient data for injection and falloff periods affected by a propagating or closing fracture. Both changing and stationary fracture properties can be handled along with stationary and moving temperature and saturation profiles (which are based on Buckley-Leverett-type computations). The approach is based on general analytical solutions for the Laplace-transformed pressure in multizone elliptic systems under constant-rate production. It extends earlier results on flow in single- and two-zone elliptic systems. Rate distributions along fractures and rate transients when the wellbore pressure is kept constant can also be obtained.