Program cooling tower calculation