from caes import ICAES, plot_series
import matplotlib.pyplot as plt
# create system
inputs = ICAES.get_default_inputs()
inputs['depth'] = 5550.0
system = ICAES(inputs=inputs)
# system.debug_perf()
#
# # run single cycle and analyze
system.single_cycle()
results = system.analyze_performance()
results.to_csv('single_cycle_performance.csv')
system.data.to_csv('single_cycle_timeseries.csv')
print(results)
# # plot results
system.plot_overview()
system.plot_pressures()
system.plot_pressure_losses()
inputs['ML_cmp3'] = 2.0
inputs['ML_exp1'] = 2.0
inputs['ML_exp2'] = 2.5
inputs['ML_exp3'] = 3.0
inputs['delta_p_cmp12'] = 0.0
inputs['delta_p_cmp23'] = 0.02
inputs['delta_p_cmp34'] = 0.0
inputs['delta_p_exp12'] = 0.0
inputs['delta_p_exp23'] = 0.02
inputs['delta_p_exp34'] = 0.0
# Turn-off leakage
# inputs['loss_m_air'] = 0.0
# ===========================
# create system and run
# ===========================
system = ICAES(inputs=inputs)
system.single_cycle()
results = system.analyze_performance()
print(results)
results.to_csv(casename + 'single_cycle_performance.csv')
system.data.to_csv(casename + 'single_cycle_timeseries.csv')
# # plot results
system.plot_overview(casename=casename)
system.plot_pressures(casename=casename)
system.plot_pressure_losses(casename=casename)