'Returns the energy production [Wh] for the given release [m3/s]'

def energy_production(self, release):

    power = ((((1000 * 9.81) * self.head) * release) * self.efficiency)
    return (power * 3600)
