if __name__ == '__main__':
pLow = 10*10**3
pHigh = 100*10**5
results = []
int_p = []
for x in range(pLow+1, pHigh-1, 10*10**3):
int_p.append(x/1000)
results.append(cfwh(x))
thermal = [res.eta for res in results]
exergetic = [res.E for res in results]
idx = thermal.index(max(thermal))
print 'Max Thermal Efficiency of {} with an Intermediate pressure of {} kPa'.format(
max(thermal), int_p[thermal.index(max(thermal))])
print 'Max Exergetic Efficiency of {} with an Intermediate pressure of {} kPa'.format(
max(exergetic), int_p[exergetic.index(max(exergetic))])
print 'Turbine one: {}'.format(results[idx].turb_one.ef)
print 'Turbine two: {}'.format(results[idx].turb_two.ef)
print 'Pump: {}'.format(results[idx].pump_one.ef)
print 'Y {}'.format(results[idx].y)
print 'CFWH: {}'.format(results[idx].cfwh.ef)
plot.plotData('Thermal Efficiencies of a CFWH Cycle', 'Intermediate Pressure (kPa)',
'Thermal Efficiency', [int_p, thermal])
plot.plotData('Exergetic Efficiencies of a CFWH Cycle', 'Intermediate Pressure (kPa)',
'Exergetic Efficiency', [int_p, exergetic])
plot.plotComponent('CFWH Turbine One', [res.turb_one for res in results], int_p)
plot.plotComponent('CFWH Turbine Two', [res.turb_two for res in results], int_p)
plot.plotComponent('CFWH Pump', [res.pump_one for res in results], int_p)
plot.plotComponent('CFWH Feed Water', [res.cfwh for res in results], int_p)
int_p = []
for x in range(pLow+1000, pHigh-1000, 10000):
int_p.append(x/1000)
results.append(reheater(x))
thermal = [res.eta for res in results]
exergetic = [res.E for res in results]
idx = thermal.index(max(thermal))
print 'Max Thermal Efficiency of {} with an Intermediate pressure of {} kPa'.format(
max(thermal), int_p[thermal.index(max(thermal))])
print 'Max Exergetic Efficiency of {} with an Intermediate pressure of {} kPa'.format(
max(exergetic), int_p[exergetic.index(max(exergetic))])
print 'Turbine one: {}'.format(results[idx].turb_one.ef)
print 'Turbine two: {}'.format(results[idx].turb_two.ef)
print 'Pump: {}'.format(results[idx].pump.ef)
plot.plotData('Thermal Efficiencies of a Reheater Cycle', 'Intermediate Pressure (kPa)',
'Thermal Efficiency', [int_p, thermal])
plot.plotData('Exergetic Efficiencies of a Reheater Cycle', 'Intermediate Pressure (kPa)',
'Exergetic Efficiency', [int_p, exergetic])
plot.plotComponent('Reheater Turbine One', [res.turb_one for res in results], int_p)
plot.plotComponent('Reheater Turbine Two', [res.turb_two for res in results], int_p)
plot.plotComponent('Reheater Pump', [res.pump for res in results], int_p)
