Discrepancy between PySAM and UI results

[color=#0033b3]import [/color]json [color=#0033b3]import [/color]pandas [color=#0033b3]as [/color]pd [color=#0033b3]import [/color]PySAM.Battery [color=#0033b3]as [/color]Battery [color=#0033b3]import [/color]PySAM.GenericSystem [color=#0033b3]as [/color]GenericSystem [color=#0033b3]import [/color]PySAM.Grid [color=#0033b3]as [/color]Grid [color=#0033b3]import [/color]PySAM.Singleowner [color=#0033b3]as [/color]SingleOwner [color=#0033b3]import [/color]PySAM.Utilityrate5 [color=#0033b3]as [/color]UtilityRate [color=#0033b3]import [/color]os [color=#0033b3]from [/color]numpy.testing [color=#0033b3]import [/color]assert_allclose [color=#0033b3]import [/color]numpy [color=#0033b3]as [/color]np [color=#0033b3]from [/color]pathlib [color=#0033b3]import [/color]Path [color=#0033b3]import [/color]matplotlib.pyplot [color=#0033b3]as [/color]plt [color=#8c8c8c]# directory where default inputs are stored [/color]working_dir = Path(r'<specify directory of saved files here!!!!!!!>') # open SAM UI result sam_ui_result = pd.read_csv(os.path.join(working_dir, [color=#067d17]"sam_ui_output_6b.csv"[/color])) [color=#8c8c8c]# Create a list of filenames. The order must be the same as the execution order [/color]module_names = [ [color=#067d17]"generic_system"[/color], [color=#067d17]"battery"[/color], [color=#067d17]"grid"[/color], [color=#067d17]"utilityrate5"[/color], [color=#067d17]"singleowner"[/color], ] file_names = [ [color=#067d17]"GenericBatterySingleOwner_6b_" [/color]+ x + [color=#067d17]".json" [/color][color=#0033b3]for [/color]x [color=#0033b3]in [/color]module_names ] paths = [os.path.join(working_dir, file_name) [color=#0033b3]for [/color]file_name [color=#0033b3]in [/color]file_names] [color=#8c8c8c]# read files [/color]module_inputs = [] [color=#0033b3]for [/color]f [color=#0033b3]in [/color]paths: [color=#0033b3]with [/color][color=#000080]open[/color](f, [color=#067d17]"r"[/color]) [color=#0033b3]as [/color]file: [color=#0033b3]try[/color]: module_inputs.append(json.load(file)) [color=#0033b3]except [/color][color=#000080]IOError [/color][color=#0033b3]as [/color]e: [color=#0033b3]raise [/color][color=#000080]IOError[/color]([color=#067d17]"Error when reading SAM json files" [/color]+ [color=#000080]str[/color](e)) [color=#8c8c8c]# initialize the primary module (Generic System) [/color]gs = GenericSystem.new() [color=#8c8c8c]# initialize the subsequent modules, linking them to the primary module [/color]batt = Battery.from_existing(gs) gr = Grid.from_existing(gs) ur = UtilityRate.from_existing(gs) so = SingleOwner.from_existing(gs) [color=#8c8c8c]# Create a list of modules. The execution order must be correct [/color]modules = [gs, batt, gr, ur, so] [color=#8c8c8c]# transfer input from module_inputs list to the modules [/color][color=#0033b3]for [/color]m, data [color=#0033b3]in [/color][color=#000080]zip[/color](modules, module_inputs): [color=#8c8c8c]# Loop through each key-value pair [/color][color=#8c8c8c] [/color][color=#0033b3]for [/color]k, v [color=#0033b3]in [/color]data.items(): [color=#0033b3]if [/color]k != [color=#067d17]"number_inputs"[/color]: k = k.replace( [color=#067d17]"adjust_"[/color], [color=#067d17]"" [/color][color=#067d17] [/color]) [color=#8c8c8c]# workaround for known bug, see https://github.com/NREL/pysam/issues/164 [/color][color=#8c8c8c] [/color][color=#0033b3]try[/color]: m.value(k, v) [color=#0033b3]except [/color][color=#000080]Exception [/color][color=#0033b3]as [/color]e: [color=#8c8c8c]# Catch any exception [/color][color=#8c8c8c] [/color][color=#0033b3]raise [/color][color=#000080]Exception[/color]( [color=#067d17]f"Error with key: [/color][color=#0037a6]{[/color]k[color=#0037a6]}[/color][color=#067d17], value: [/color][color=#0037a6]{[/color]v[color=#0037a6]}[/color][color=#067d17]. Original error: [/color][color=#0037a6]{[/color]e[color=#0037a6]}[/color][color=#067d17]" [/color][color=#067d17] [/color]) [color=#8c8c8c]# execute pySAM modules in sequence [/color][color=#0033b3]for [/color]m [color=#0033b3]in [/color]modules: m.execute() [color=#8c8c8c]## plotting [/color] sam_ui_result[[color=#067d17]'pysam_gen'[/color]] = modules[[color=#1750eb]1[/color]].SystemOutput.gen sam_ui_result[[color=#067d17]'pysam_ppa_price'[/color]] = modules[[color=#1750eb]1[/color]].Outputs.market_sell_rate_series_yr1 sam_ui_result[[color=#067d17]'pysam_soc'[/color]] = modules[[color=#1750eb]1[/color]].Outputs.batt_SOC datetimes_1_year = pd.date_range([color=#660099]start[/color]=[color=#067d17]'1990-1-1 0:30:00'[/color], [color=#660099]periods [/color]= [color=#1750eb]8760[/color], [color=#660099]freq [/color]= [color=#067d17]'H'[/color]) sam_ui_result[[color=#067d17]'timestamp'[/color]] = [color=#000080]list[/color](datetimes_1_year)*[color=#1750eb]25 [/color]year = [] [color=#0033b3]for [/color]i [color=#0033b3]in [/color]np.arange([color=#1750eb]25[/color]): year = year + [i]*[color=#1750eb]8760 [/color]sam_ui_result[[color=#067d17]'year'[/color]] = year sam_ui_result[[color=#067d17]'month'[/color]] = sam_ui_result[[color=#067d17]'timestamp'[/color]].dt.month sam_ui_result[[color=#067d17]'day'[/color]] = sam_ui_result[[color=#067d17]'timestamp'[/color]].dt.day sam_ui_result[[color=#067d17]'disc'[/color]] = [color=#000080]abs[/color]( sam_ui_result[[color=#067d17]"System power generated | (kW)"[/color]] - sam_ui_result[[color=#067d17]'pysam_gen'[/color]] ) max_disc = sam_ui_result.loc[sam_ui_result[[color=#067d17]'year'[/color]]==[color=#1750eb]0[/color],[color=#067d17]'disc'[/color]].idxmax() plot_year_wc, plot_month_wc, plot_day_wc = sam_ui_result.loc[max_disc, [[color=#067d17]'year'[/color],[color=#067d17]'month'[/color],[color=#067d17]'day'[/color]]] # the date of the worst-case discrepancy was May 1st # plot for days and months surrounding the worst-case date [color=#0033b3]for [/color]plot_month [color=#0033b3]in [/color][[color=#1750eb]4[/color], [color=#1750eb]5[/color]]: [color=#0033b3]for [/color]plot_day [color=#0033b3]in [/color][color=#000080]range[/color]([color=#1750eb]1[/color],[color=#1750eb]31[/color]): plot_df = sam_ui_result.loc[( (sam_ui_result[[color=#067d17]'year'[/color]] == plot_year) & (sam_ui_result[[color=#067d17]'month'[/color]] == plot_month) & (sam_ui_result[[color=#067d17]'day'[/color]] == plot_day) ), :].copy() fig, axs = plt.subplots([color=#1750eb]3[/color],[color=#1750eb]1[/color], [color=#660099]sharex[/color]=[color=#0033b3]True[/color], [color=#660099]figsize[/color]=([color=#1750eb]10[/color],[color=#1750eb]10[/color])) [color=#0033b3]for [/color]col [color=#0033b3]in [/color][[color=#067d17]'System power generated | (kW)'[/color], [color=#067d17]'Electricity to grid from battery AC | (kW)'[/color], [color=#067d17]'Electricity to battery from grid AC | (kW)'[/color], [color=#067d17]'pysam_gen'[/color]]: axs[[color=#1750eb]0[/color]].plot(plot_df[[color=#067d17]'timestamp'[/color]], plot_df[col], [color=#660099]label [/color]= col, [color=#660099]marker[/color]=[color=#067d17]'o'[/color]) axs[[color=#1750eb]1[/color]].plot(plot_df[[color=#067d17]'timestamp'[/color]], plot_df[[color=#067d17]'Power price for battery dispatch | ($/MWh)'[/color]], [color=#660099]marker[/color]=[color=#067d17]'o'[/color], [color=#660099]label[/color]=[color=#067d17]'UI'[/color]) axs[[color=#1750eb]1[/color]].plot(plot_df[[color=#067d17]'timestamp'[/color]], plot_df[[color=#067d17]'pysam_ppa_price'[/color]], [color=#660099]marker[/color]=[color=#067d17]'o'[/color], [color=#660099]label[/color]=[color=#067d17]'pysam'[/color]) axs[[color=#1750eb]1[/color]].set_ylabel([color=#067d17]'Power price for battery dispatch | ($/MWh)'[/color]) axs[[color=#1750eb]2[/color]].plot(plot_df[[color=#067d17]'timestamp'[/color]], plot_df[[color=#067d17]'Battery state of charge | (%)'[/color]], [color=#660099]marker[/color]=[color=#067d17]'o'[/color], [color=#660099]label[/color]=[color=#067d17]'UI'[/color]) axs[[color=#1750eb]2[/color]].plot(plot_df[[color=#067d17]'timestamp'[/color]], plot_df[[color=#067d17]'pysam_soc'[/color]], [color=#660099]marker[/color]=[color=#067d17]'o'[/color], [color=#660099]label[/color]=[color=#067d17]'pysam'[/color]) axs[[color=#1750eb]2[/color]].set_ylabel([color=#067d17]'Battery state of charge | (%)'[/color]) axs[[color=#1750eb]0[/color]].legend([color=#660099]prop[/color]={[color=#067d17]'size'[/color]: [color=#1750eb]8[/color]}) axs[[color=#1750eb]1[/color]].legend() axs[[color=#1750eb]2[/color]].legend() fig.autofmt_xdate() plt.title([color=#067d17]f'year[/color][color=#0037a6]{[/color]plot_year[color=#0037a6]}[/color][color=#067d17]_month[/color][color=#0037a6]{[/color]plot_month[color=#0037a6]}[/color][color=#067d17]_day[/color][color=#0037a6]{[/color]plot_day[color=#0037a6]}[/color][color=#067d17]'[/color]) plt.savefig(os.path.join(working_dir, [color=#067d17]f'troubleshoot_case6b_year[/color][color=#0037a6]{[/color]plot_year[color=#0037a6]}[/color][color=#067d17]_month[/color][color=#0037a6]{[/color]plot_month[color=#0037a6]}[/color][color=#067d17]_day[/color][color=#0037a6]{[/color]plot_day[color=#0037a6]}[/color][color=#067d17].jpg'[/color]), [color=#660099]dpi[/color]=[color=#1750eb]300[/color], [color=#660099]bbox_inches[/color]=[color=#067d17]'tight'[/color]) [color=#8c8c8c]# plt.show() [/color][color=#8c8c8c] [/color]plt.close() # the following raises an assertion error assert_allclose( modules[[color=#1750eb]1[/color]].SystemOutput.gen[:], sam_ui_result[[color=#067d17]"System power generated | (kW)"[/color]][:], [color=#660099]rtol[/color]=[color=#1750eb]1e-3[/color], [color=#660099]atol[/color]=[color=#1750eb]0.1[/color], )