Adjusted code to run stationary period, followed by nonstationary period (currently in tauP).