ABSTRACT: A combination pumped storage reservoir system was simulated by modifying the WRE deep reservoir model. Each of the two reservoirs was described by a copy of the WRE model program, the two programs were converted into subroutines and were called upon alternately by a main program. Operationally, the contributing reservoir, i.e., the reservoir from which flow was discharged, was simulated for one execution interval (1 hour), followed by simulating the receiving reservoir for the same execution interval. The main program directed the discharge temperature of the contributing reservoir for each execution interval as input for simulating the receiving reservoir. The two subroutines were run in this interactive mode for a simulation period of one year.Two simulations, labeled “Pump” and “Net,” were effected and differed basically in the distribution of flow volumes exchanged between the reservoirs. In the ‘Pump’ simulation the total hourly flow volumes were distributed into appropriate horizontal layers of the lakes as determined by temperature‐density relationships, i.e., the pumped discharges were distributed into the upper lake and the generation discharges were distributed into the lower lake. In the ‘Net’ simulation only the net daily discharges (daily difference between pump‐back and generation at the upper dam), distributed uniformly over 24 hours, were mixed into the horizontal layers of the two lakes.Both simulations produced annual thermal regimes that were apparently within reason for the geographical area and the nature of the input data. However, neither accurately reflected a generalized conditions for the reservoirs because the Pump simulation reflected conditions in the forebay and tailrace of the Upper dam while the Net simulation reflected conditions of the remaining parts of the reservoirs.