The design of cooling systems for gas turbine engine blades and vanes calls for efficient simulation programs. The main purpose of the described program is to determine the complete boundary condition at the coolant side to support a temperature calculation for the solid. For the simulation of convection and heat pick up of the coolant flow, pressure loss, and further effects to be found in a rotating frame, the cooling systems are represented by networks of nodes and flow elements. Within each flow element the fluid flow is modeled by a system of ordinary differential equations based on the one-dimensional conservation of mass, momentum, and energy. In this respect, the computer program differs from many other network computation programs. Concerning cooling configurations in rotating systems, the solution for a single flow element or the entire flow system is not guaranteed to be unique. This is due to rotational forces in combination with heat transfer and causes considerable computational difficulties, which can be overcome by a special path following method in which the angular velocity is selected as the parameter of homotopy. Results of the program are compared with measurements for three applications.