예제 #1
0
    'frout': ['frstorage'],
    'lwrpu': ['lwrmixer'],
    'frpu': ['frmixer'],
    'mixerout': ['fr'],
    'lwrreprocessingwaste': ['lwrsink'],
    'frreprocessingwaste': ['frsink']
}

for commod, facility in commod_dict.items():
    agent_entry_dict[commod] = tester.get_agent_dict(output_file, facility)

all_dict['power'] = tester.supply_demand_dict_driving(output_file, demand_eq,
                                                      'power')

plotter.plot_demand_supply_agent(all_dict['power'], agent_entry_dict['power'],
                                 'power', '0-' + calc_method + '-power', True,
                                 True, False, 1)

front_commods = ['sourceout', 'enrichmentout']
mid_commods = ['mixerout']
back_commods = [
    'lwrstorageout', 'frstorageout', 'lwrout', 'frout', 'lwrreprocessingwaste',
    'frreprocessingwaste', 'frpu', 'lwrpu'
]

for commod in front_commods:
    all_dict[commod] = tester.supply_demand_dict_nondriving(
        output_file, commod, True)
    name = '0-' + calc_method + '-' + commod
    plotter.plot_demand_supply_agent(all_dict[commod],
                                     agent_entry_dict[commod], commod, name,
예제 #2
0
    'frout': ['frstorage'],
    'lwrtru': ['pumixerlwr'],
    'frtru': ['pumixerfr'],
    'lwrreprocessingwaste': ['lwrsink'],
    'frreprocessingwaste': ['frsink']
}

for commod, facility in commod_dict.items():
    agent_entry_dict[commod] = tester.get_agent_dict(output_file, facility)

# get supply deamnd dict
all_dict['power'] = tester.supply_demand_dict_driving(output_file, demand_eq,
                                                      'power')

plotter.plot_demand_supply_agent(all_dict['power'], agent_entry_dict['power'],
                                 'power', 'eg01-eg24-flatpower-d3ploy_power',
                                 True, True, False, 1)

front_commods = ['sourceout', 'enrichmentout']
back_commods = [
    'lwrstorageout', 'frstorageout', 'lwrout', 'frout', 'lwrreprocessingwaste',
    'frreprocessingwaste', 'frtru', 'lwrtru'
]

for commod in front_commods:
    all_dict[commod] = tester.supply_demand_dict_nondriving(
        output_file, commod, True)
    name = 'B2000-' + commod
    plotter.plot_demand_supply_agent(all_dict[commod],
                                     agent_entry_dict[commod], commod, name,
                                     True, True, False, 1)
        "reactor10",
        "newreactor",
    ],
)
agent_entry_dict["fuel"] = tester.get_agent_dict(
    output_file, ["source", "initialsource"]
)
agent_entry_dict["spent_fuel"] = tester.get_agent_dict(output_file, [
                                                       "sink"])
# plots demand, supply, calculated demand, calculated supply for the
# scenario for each calc method
plotter.plot_demand_supply_agent(
    all_dict_power,
    agent_entry_dict["power"],
    "power",
    name + " Power",
    True,
    False,
    False,
)
plotter.plot_demand_supply_agent(
    all_dict_fuel,
    agent_entry_dict["fuel"],
    "fuel",
    name + " Fuel",
    True,
    False,
    False)
plotter.plot_demand_supply_agent(
    all_dict_spent_fuel,
    agent_entry_dict["spent_fuel"],
    output_file = name + '.sqlite'
    with open(input_file, 'w') as f:
        json.dump(scenario_1_input[calc_method], f)
    s = subprocess.check_output(['cyclus', '-o', output_file, input_file],
                                universal_newlines=True,
                                env=ENV)

    # Initialize dicts
    all_dict_fuel = {}
    agent_entry_dict = {}

    all_dict_fuel = tester.supply_demand_dict_driving(output_file, demand_eq,
                                                      'fuel')
    agent_entry_dict['fuel'] = tester.get_agent_dict(output_file, ['source'])
    plotter.plot_demand_supply_agent(all_dict_fuel, agent_entry_dict['fuel'],
                                     'fuel', name + '_fuel', True, False,
                                     False)

    metric_dict = tester.metrics(all_dict_fuel, metric_dict, calc_method,
                                 'fuel', True)

    df = pd.DataFrame(metric_dict)
    df.to_csv('scenario_1_output.csv')

##########################################################################

######################################SCENARIO 2##########################
# scenario 2, source -> reactor (cycle time = 1, refuel time = 0) -> sink
scenario_2_input = {}
demand_eq = "1000*t"
    'frout': ['frstorage'],
    'lwrreprocessingoutpu': ['pumixerlwr'],
    'frreprocessingoutpu': ['pumixerfr'],
    'lwrreprocessingwaste': ['lwrsink'],
    'frreprocessingwaste': ['frsink']
}
for commod, facility in commod_dict.items():
    agent_entry_dict[commod] = tester.get_agent_dict(output_file, facility)
###########################

# get supply deamnd dict
# and plot
all_dict['power'] = tester.supply_demand_dict_driving(output_file, demand_eq,
                                                      'power')
plotter.plot_demand_supply_agent(all_dict['power'], agent_entry_dict['power'],
                                 'power', 'transitionscenario_1_input_power',
                                 True)

front_commods = ['sourceout', 'enrichmentout']
back_commods = [
    'lwrstorageout', 'frstorageout', 'lwrout', 'frout', 'lwrreprocessingoutpu',
    'frreprocessingoutpu', 'lwrreprocessingwaste', 'frreprocessingwaste'
]
for commod in front_commods:
    all_dict[commod] = tester.supply_demand_dict_nondriving(
        output_file, commod, True)
    name = 'transitionscenario_1_input_' + commod
    plotter.plot_demand_supply_agent(all_dict[commod],
                                     agent_entry_dict[commod], commod, name,
                                     True)
    metric_dict = tester.metrics(all_dict[commod], metric_dict, calc_method,
예제 #6
0
    agent_entry_dict['sourceoutput'] = tester.get_agent_dict(output_file, [
                                                             'source'])
    agent_entry_dict['reactor1output'] = tester.get_agent_dict(output_file, [
                                                               'separations1'])
    agent_entry_dict['separations1Pu'] = tester.get_agent_dict(output_file, [
                                                               'mixer'])
    agent_entry_dict['reactor2output'] = tester.get_agent_dict(output_file, [
                                                               'separations2'])

    # plots demand, supply, calculated demand, calculated supply for the
    # scenario for each calc method
    name1 = "scenario_7_input_" + calc_method + "_power"
    plotter.plot_demand_supply_agent(
        all_dict['power'],
        agent_entry_dict['power'],
        'power',
        name1,
        True,
        False,
        False)
    name2 = "scenario_7_input_" + calc_method + "_sourceoutput"
    plotter.plot_demand_supply_agent(
        all_dict['sourceoutput'],
        agent_entry_dict['sourceoutput'],
        'sourceoutput',
        name2,
        True,
        False,
        False)
    name3 = "scenario_7_input_" + calc_method + "_reactor1output"
    plotter.plot_demand_supply_agent(
        all_dict['reactor1output'],
예제 #7
0
    'frstorageout': ['frreprocessing'],
    'moxstorageout': ['moxreprocessing'],
    'frmixerout': ['frmixer'],
    'moxmixerout': ['moxmixer'],
    'lwrtru': ['frmixer', 'moxmixer'],
    'frtru': ['frmixer', 'moxmixer'],
    'moxtru': ['frmixer', 'moxmixer']
}

for commod, facility in commod_dict.items():
    agent_entry_dict[commod] = tester.get_agent_dict(output_file, facility)

all_dict['power'] = tester.supply_demand_dict_driving(output_file, demand_eq,
                                                      'power')
plotter.plot_demand_supply_agent(all_dict['power'], agent_entry_dict['power'],
                                 'power', name + '_power', True, True, False,
                                 1)

front_commods = ['sourceout', 'enrichmentout']
mid_commods = ['lwrtru', 'frtru', 'moxtru']
back_commods = ['lwrstorageout', 'frstorageout', 'moxstorageout']

for commod in front_commods:
    all_dict[commod] = tester.supply_demand_dict_nondriving(
        output_file, commod, True)
    name = 'flatpower-' + calc_method + '-' + commod
    plotter.plot_demand_supply_agent(all_dict[commod],
                                     agent_entry_dict[commod], commod, name,
                                     True, True, False, 1)
    metric_dict = tester.metrics(all_dict[commod], metric_dict, calc_method,
                                 commod, True)
    all_dict_power = tester.supply_demand_dict_driving(output_file, demand_eq,
                                                       'power')
    all_dict_fuel = tester.supply_demand_dict_nondriving(
        output_file, 'fuel', True)
    all_dict_spentfuel = tester.supply_demand_dict_nondriving(
        output_file, 'spentfuel', False)

    agent_entry_dict['power'] = tester.get_agent_dict(output_file, ['reactor'])
    agent_entry_dict['fuel'] = tester.get_agent_dict(output_file, ['source'])
    agent_entry_dict['spentfuel'] = tester.get_agent_dict(
        output_file, ['sink'])

    # plots demand, supply, calculated demand, calculated supply for the
    # scenario for each calc method
    plotter.plot_demand_supply_agent(all_dict_power, agent_entry_dict['power'],
                                     'power', name + "_power", True, False,
                                     False)
    plotter.plot_demand_supply_agent(all_dict_fuel, agent_entry_dict['fuel'],
                                     'fuel', name + "_fuel", True, False,
                                     False)
    plotter.plot_demand_supply_agent(all_dict_spentfuel,
                                     agent_entry_dict['spentfuel'],
                                     'spentfuel', name + "_spentfuel", False,
                                     False, False)

    metric_dict = tester.metrics(all_dict_power, metric_dict, calc_method,
                                 'power', True)
    metric_dict = tester.metrics(all_dict_fuel, metric_dict, calc_method,
                                 'fuel', True)
    metric_dict = tester.metrics(all_dict_spentfuel, metric_dict, calc_method,
                                 'spentfuel', False)
                                env=ENV)
    # Initialize dicts
    all_dict = {}
    agent_entry_dict = {}
    all_dict['power'] = tester.supply_demand_dict_driving(
        output_file, demand_eq, 'power')
    all_dict['sourceoutput'] = tester.supply_demand_dict_nondriving(
        output_file, 'sourceoutput', True)

    agent_entry_dict['power'] = tester.get_agent_dict(output_file,
                                                      ['reactor1', 'reactor2'])
    agent_entry_dict['sourceoutput'] = tester.get_agent_dict(
        output_file, ['source'])

    # plots demand, supply, calculated demand, calculated supply for the scenario for each calc method
    plotter.plot_demand_supply_agent(all_dict['power'],
                                     agent_entry_dict['power'], 'power', name,
                                     True)
    name2 = "scenario_7_input_" + calc_method + "_sourceoutput"
    plotter.plot_demand_supply_agent(all_dict['sourceoutput'],
                                     agent_entry_dict['sourceoutput'],
                                     'sourceoutput', name2, True)

    metric_dict = tester.metrics(all_dict['power'], metric_dict, calc_method,
                                 'power', True)
    metric_dict = tester.metrics(all_dict['sourceoutput'], metric_dict,
                                 calc_method, 'sourceoutput', True)

    df = pd.DataFrame(metric_dict)
    df.to_csv('scenario_7_output.csv')