# fig-EETscaleup-A.py
# Calculates efficiency of EET based rewired carbon fixation
#
# Farshid Salimijazi and Buz Barstow
# Last updated by Buz Barstow on 2019-10-28
# ------------------------------------------------------------------------------------------------ #
from rewiredcarbon.scenario import ImportScenarioTable, CalculateScenarioEfficiencies, \
Plot_Efficiency_Scattergraphs, Export_Efficiency_Scattergraphs
from rewiredcarbon.utils import ensure_dir
scenarioTableFileName = 'input/fig-EETscaleup-A.csv'
outputFileDirname = 'output/fig-EETscaleup/fig-EETscaleup-A/'
ensure_dir(outputFileDirname)
outputFilePrefix = 'VariableResistivity'
scenarioDict = ImportScenarioTable(scenarioTableFileName)
efficienciesDict = CalculateScenarioEfficiencies(scenarioDict,
mode='scattergraph')
Plot_Efficiency_Scattergraphs(efficienciesDict, 'effTotalElectricalToFuel')
Export_Efficiency_Scattergraphs(outputFileDirname, outputFilePrefix, efficienciesDict, \
'effTotalElectricalToFuel', keysToPlot=None)
# Farshid Salimijazi and Buz Barstow # Last updated by Buz Barstow on 2019-10-25 # ------------------------------------------------------------------------------------------------ # from rewiredcarbon.scenario import ImportScenarioTable, CalculateScenarioEfficiencies, \ Plot_Efficiency_Bargraph, Generate_EfficienciesDict_Keys_Sorted_by_Efficiency, \ Export_Efficiency_Bargraph from rewiredcarbon.utils import ensure_dir scenarioTableFileName = 'input/fig-cbb_n2_to_amino_acids.csv' outputFilenameEff = 'output/fig-cbb-n2_to_amino_acids/fig-cbb_n2_to_amino_acids_eff.csv' outputFilenameFuelMassEff = 'output/fig-cbb-n2_to_amino_acids/fig-cbb_n2_to_amino_acids_fuel_mass.csv' ensure_dir(outputFilenameEff) ensure_dir(outputFilenameFuelMassEff) scenarioDict = ImportScenarioTable(scenarioTableFileName) efficienciesDict = CalculateScenarioEfficiencies(scenarioDict) # keysArray = \ # Generate_EfficienciesDict_Keys_Sorted_by_Efficiency(efficienciesDict, 'effTotalElectricalToFuel') keysArray = list(efficienciesDict.keys()) Plot_Efficiency_Bargraph(efficienciesDict, 'effTotalElectricalToFuel', \ 'effTotalElectricalToFuel_lowerError', 'effTotalElectricalToFuel_upperError', keysToPlot=keysArray)
# Last updated by Buz Barstow on 2019-10-25 # ------------------------------------------------------------------------------------------------ # from rewiredcarbon.scenario import ImportScenarioTable, CalculateScenarioEfficiencies, \ Plot_Efficiency_Bargraph, Generate_EfficienciesDict_Keys_Sorted_by_Efficiency, \ Export_Efficiency_Bargraph from rewiredcarbon.utils import ensure_dir scenarioTableFileName = 'input/fig-co2fixation.csv' outputFilename = 'output/fig-co2fixation/fig-co2fixation.csv' ensure_dir(outputFilename) scenarioDict = ImportScenarioTable(scenarioTableFileName) efficienciesDict = CalculateScenarioEfficiencies(scenarioDict) # keysArray = \ # Generate_EfficienciesDict_Keys_Sorted_by_Efficiency(efficienciesDict, 'effTotalElectricalToFuel') keysArray = list(efficienciesDict.keys()) Plot_Efficiency_Bargraph(efficienciesDict, 'effTotalElectricalToFuel', \ 'effTotalElectricalToFuel_lowerError', 'effTotalElectricalToFuel_upperError', keysToPlot=keysArray)
for targetEfficiencyRatio in targetEfficiencyRatios:
outputDict = \
Calculate_Density_and_Area_for_Target_to_Peak_Efficiency_Ratio_for_Array_of_Input_Powers(\
scenarioDict, efficienciesDict, efficienciesDictKeysForCalc, targetEfficiencyRatio, \
peakElectricalToFuelEfficiency, solarConstant=solarConstant)
collectedOutputDict[str(targetEfficiencyRatio)] = deepcopy(outputDict)
# Plot out the system footprint calculations and export them to a CSV file
keys = collectedOutputDict.keys()
vectorList_tankAreaRelativeToSolarPVArea = []
headerList_tankAreaRelativeToSolarPVArea = []
outputFilename_tankAreaRelativeToSolarPVArea \
= outputFileDirname + '/' + 'fig-h2scaleup-E-area_relative_to_PV_area.csv'
ensure_dir(outputFilename_tankAreaRelativeToSolarPVArea)
figure()
for key in keys:
outputDict = collectedOutputDict[key]
totalElectricalPowerArray = outputDict['totalElectricalPowerArray']
tankAreaRelativeToSolarPVAreaArray = outputDict[
'tankAreaRelativeToSolarPVAreaArray']
loglog(totalElectricalPowerArray,
tankAreaRelativeToSolarPVAreaArray,
label=key)
headerList_tankAreaRelativeToSolarPVArea.append('totalElectricalPower_' +
key)
