Lcos cost inputs

ssc/cmod_battery.cpp

@@ -399,18 +399,18 @@ battstor::battstor(var_table& vt, bool setup_model, size_t nrec, double dt_hr, c
 
 
             // Battery bank replacement
-            if (vt.is_assigned("om_replacement_cost1"))
+            if (vt.is_assigned("om_batt_replacement_cost"))
             {
                 std::vector<ssc_number_t> replacement_cost(nyears);
-                ssc_number_t* parr = vt.as_array("om_replacement_cost1", &cnt);
+                ssc_number_t*  parr = vt.as_array("om_batt_replacement_cost", &cnt);
                 if (cnt == 1)
                 {
                     for (i = 0; i < nyears; i++)
                         replacement_cost[i] = parr[0] * (ssc_number_t)pow((double)(inflation_rate + 1), (double)i);
                 }
                 else if (cnt < nyears)
                 {
-                    throw exec_error("battery", "Invalid number for om_replacement_cost1, must be 1 or equal to analysis_period.");
+                    throw exec_error("battery", "Invalid number for om_batt_replacement_cost, must be 1 or equal to analysis_period.");
                 }
                 else {
                     for (i = 0; i < nyears; i++)

ssc/cmod_cashloan.cpp

@@ -26,8 +26,8 @@ OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 #include "lib_financial.h"
 #include "lib_util.h"
 #include "common_financial.h"
-#include <sstream>
 using namespace libfin;
+#include <sstream>
 
 static var_info vtab_cashloan[] = {
 /*   VARTYPE           DATATYPE          NAME                        LABEL                                  UNITS         META                      GROUP            REQUIRED_IF                 CONSTRAINTS                      UI_HINTS*/
@@ -37,11 +37,13 @@ static var_info vtab_cashloan[] = {
 
 	{ SSC_INPUT,        SSC_NUMBER,      "total_installed_cost",     "Total installed cost",               "$",            "",                      "System Costs",            "*",                      "MIN=0",                                         "" },
 	{ SSC_INPUT,        SSC_NUMBER,      "salvage_percentage",       "Salvage value percentage",           "%",            "",                      "Financial Parameters",      "?=0.0",                  "MIN=0,MAX=100",                 "" },
-
+    //{ SSC_INPUT,        SSC_NUMBER,     "batt_salvage_percentage",                     "Net pre-tax cash battery salvage value",	                               "%",	 "",					  "Financial Parameters",             "?=0",                     "MIN=0,MAX=100",      			"" },
+
 	{ SSC_INPUT,        SSC_ARRAY,       "annual_energy_value",      "Energy value",                       "$",            "",                      "System Output",      "*",                       "",                                         "" },
 	{ SSC_INPUT,        SSC_ARRAY,       "annual_themal_value",      "Energy value",                       "$",            "",                      "System Output",      "",                       "",                                         "" },
 	{ SSC_INPUT,        SSC_ARRAY,       "gen",                      "Power generated by renewable resource",                    "kW", "", "System Output", "*", "", "" },
-	{ SSC_INPUT,        SSC_ARRAY,       "degradation",              "Annual degradation", "%", "", "System Output", "*", "", "" },
+
+    { SSC_INPUT,        SSC_ARRAY,       "degradation",              "Annual degradation", "%", "", "System Output", "*", "", "" },
 	{ SSC_INPUT,        SSC_NUMBER,      "system_use_lifetime_output","Lifetime hourly system outputs", "0/1", "0=hourly first year,1=hourly lifetime", "Lifetime", "*", "INTEGER,MIN=0", "" },
 
 	/* financial outputs */
@@ -52,12 +54,15 @@ static var_info vtab_cashloan[] = {
 	{ SSC_OUTPUT,        SSC_NUMBER,     "payback",                  "Payback period",                            "years",        "",                      "Cash Flow",      "*",                       "",                                         "" },
 	// added 9/26/16 for Owen Zinaman Mexico
 	{ SSC_OUTPUT, SSC_NUMBER, "discounted_payback", "Discounted payback period", "years", "", "Cash Flow", "*", "", "" },
-	{ SSC_OUTPUT, SSC_NUMBER, "npv", "Net present value", "$", "", "Cash Flow", "*", "", "" },
+
+    { SSC_OUTPUT, SSC_NUMBER, "npv", "Net present value", "$", "", "Cash Flow", "*", "", "" },
 
 	{ SSC_OUTPUT,        SSC_NUMBER,     "present_value_oandm",                      "Present value of O&M expenses",				   "$",            "",                      "Financial Metrics",      "*",                       "",                                         "" },
 	{ SSC_OUTPUT,        SSC_NUMBER,     "present_value_oandm_nonfuel",              "Present value of non-fuel O&M expenses",				   "$",            "",                      "Financial Metrics",      "*",                       "",                                         "" },
 	{ SSC_OUTPUT,        SSC_NUMBER,     "present_value_fuel",                      "Present value of fuel expenses",				   "$",            "",                      "Financial Metrics",      "*",                       "",                                         "" },
 	{ SSC_OUTPUT,        SSC_NUMBER,     "present_value_insandproptax",             "Present value of insurance and property tax",				   "$",            "",                      "Financial Metrics",      "*",                       "",                                         "" },
+
+    { SSC_OUTPUT, SSC_NUMBER, "npv_energy_lcos_nom", "Present value of annual stored energy (nominal)", "kWh", "", "LCOE calculations", "", "", "" },
 
 	{ SSC_OUTPUT, SSC_NUMBER, "adjusted_installed_cost", "Net capital cost", "$", "", "Financial Metrics", "*", "", "" },
 	{ SSC_OUTPUT, SSC_NUMBER, "loan_amount", "Debt", "$", "", "Financial Metrics", "*", "", "" },
@@ -162,7 +167,9 @@ extern var_info
 	vtab_battery_replacement_cost[],
 	vtab_fuelcell_replacement_cost[],
 	vtab_tax_credits[],
-	vtab_payment_incentives[];
+	vtab_payment_incentives[],
+    vtab_lcos_inputs[];
+
 
 enum {
 	CF_degradation,
@@ -171,16 +178,6 @@ enum {
 	CF_thermal_value,
 	CF_value_added,
 
-	CF_om_fixed_expense,
-	CF_om_production_expense,
-	CF_om_capacity_expense,
-	CF_om_fixed1_expense,
-	CF_om_production1_expense,
-	CF_om_capacity1_expense,
-	CF_om_fixed2_expense,
-	CF_om_production2_expense,
-	CF_om_capacity2_expense,
-	CF_om_fuel_expense,
 
 	CF_om_opt_fuel_2_expense,
 	CF_om_opt_fuel_1_expense,
@@ -253,7 +250,29 @@ enum {
 
 	CF_nte,
 
-	CF_max };
+    CF_om_fixed_expense,
+    CF_om_production_expense,
+    CF_om_capacity_expense,
+    CF_om_fixed1_expense,
+    CF_om_production1_expense,
+    CF_om_capacity1_expense,
+    CF_om_fixed2_expense,
+    CF_om_production2_expense,
+    CF_om_capacity2_expense,
+    CF_om_fuel_expense,
+    CF_energy_charged_grid,
+    CF_energy_charged_pv,
+    CF_energy_discharged,
+    CF_charging_cost_pv,
+    CF_charging_cost_grid,
+    CF_om_cost_lcos,
+    CF_salvage_cost_lcos,
+    CF_investment_cost_lcos,
+    CF_annual_cost_lcos,
+    CF_util_escal_rate,
+
+    CF_max,
+};
 
 
 
@@ -262,6 +281,7 @@ class cm_cashloan : public compute_module
 {
 private:
 	util::matrix_t<double> cf;
+    util::matrix_t<double> cf_lcos;
 	double ibi_fed_amount;
 	double ibi_sta_amount;
 	double ibi_uti_amount;
@@ -290,6 +310,7 @@ class cm_cashloan : public compute_module
 		add_var_info(vtab_battery_replacement_cost);
 		add_var_info(vtab_fuelcell_replacement_cost);
 		add_var_info(vtab_cashloan);
+        add_var_info(vtab_lcos_inputs);
 	}
 
 	void exec( )
@@ -299,7 +320,6 @@ class cm_cashloan : public compute_module
 		bool is_commercial = (as_integer("market")==1);
 		bool is_mortgage = (as_integer("mortgage")==1);
 
-
 //		throw exec_error("cmod_cashloan", "mortgage = " + util::to_string(as_integer("mortgage")));
 //		if (is_commercial) log("commercial market"); else log("residential market");
 //		if (is_mortgage) log("mortgage loan"); else log("standard loan");
@@ -309,6 +329,7 @@ class cm_cashloan : public compute_module
 
 		// initialize cashflow matrix
 		cf.resize_fill( CF_max, nyears+1, 0.0 );
+        cf_lcos.resize_fill(CF_max, nyears + 1, 0.0);
 
 		// initialize energy and revenue
 		size_t count = 0;
@@ -476,17 +497,17 @@ class cm_cashloan : public compute_module
 
 		if (add_om_num_types > 0)
 		{
-			escal_or_annual(CF_om_fixed1_expense, nyears, "om_fixed1", inflation_rate, 1.0, false, as_double("om_fixed_escal")*0.01);
-			escal_or_annual(CF_om_production1_expense, nyears, "om_production1", inflation_rate, 0.001, false, as_double("om_production_escal")*0.01);
-			escal_or_annual(CF_om_capacity1_expense, nyears, "om_capacity1", inflation_rate, 1.0, false, as_double("om_capacity_escal")*0.01);
-			nameplate1 = as_number("om_capacity1_nameplate");
+            escal_or_annual(CF_om_fixed1_expense, nyears, "om_batt_fixed_cost", inflation_rate, 1.0, false, as_double("om_fixed_escal") * 0.01);
+            escal_or_annual(CF_om_production1_expense, nyears, "om_batt_variable_cost", inflation_rate, 0.001, false, as_double("om_production_escal") * 0.01);
+            escal_or_annual(CF_om_capacity1_expense, nyears, "om_batt_capacity_cost", inflation_rate, 1.0, false, as_double("om_capacity_escal") * 0.01);
+            nameplate1 = as_number("ui_batt_capacity");
 		}
 		if (add_om_num_types > 1)
 		{
-			escal_or_annual(CF_om_fixed2_expense, nyears, "om_fixed2", inflation_rate, 1.0, false, as_double("om_fixed_escal")*0.01);
-			escal_or_annual(CF_om_production2_expense, nyears, "om_production2", inflation_rate, 0.001, false, as_double("om_production_escal")*0.01);
-			escal_or_annual(CF_om_capacity2_expense, nyears, "om_capacity2", inflation_rate, 1.0, false, as_double("om_capacity_escal")*0.01);
-			nameplate2 = as_number("om_capacity2_nameplate");
+			escal_or_annual(CF_om_fixed2_expense, nyears, "om_fuelcell_fixed_cost", inflation_rate, 1.0, false, as_double("om_fixed_escal")*0.01);
+			escal_or_annual(CF_om_production2_expense, nyears, "om_fuelcell_variable_cost", inflation_rate, 0.001, false, as_double("om_production_escal")*0.01);
+			escal_or_annual(CF_om_capacity2_expense, nyears, "om_fuelcell_capacity_cost", inflation_rate, 1.0, false, as_double("om_capacity_escal")*0.01);
+			nameplate2 = as_number("ui_fuelcell_capacity");
 		}
 
         // battery cost - replacement from lifetime analysis
@@ -509,7 +530,7 @@ class cm_cashloan : public compute_module
             }
             double batt_cap = as_double("batt_computed_bank_capacity");
             // updated 10/17/15 per 10/14/15 meeting
-            escal_or_annual(CF_battery_replacement_cost_schedule, nyears, "om_replacement_cost1", inflation_rate, batt_cap, false, as_double("om_replacement_cost_escal") * 0.01);
+            escal_or_annual(CF_battery_replacement_cost_schedule, nyears, "om_batt_replacement_cost", inflation_rate, batt_cap, false, as_double("om_replacement_cost_escal") * 0.01);
 
             for (i = 0; i < nyears && i < (int)count; i++) {
                 // batt_rep and the cash flow sheets are 1 indexed, replacement_percent is zero indexed
@@ -526,7 +547,7 @@ class cm_cashloan : public compute_module
 				fuelcell_rep = as_array("fuelcell_replacement", &count); // replacements per year calculated
 			else // user specified
 				fuelcell_rep = as_array("fuelcell_replacement_schedule", &count); // replacements per year user-defined
-			escal_or_annual(CF_fuelcell_replacement_cost_schedule, nyears, "om_replacement_cost2", inflation_rate, nameplate2, false, as_double("om_replacement_cost_escal")*0.01);
+			escal_or_annual(CF_fuelcell_replacement_cost_schedule, nyears, "om_fuelcell_replacement_cost", inflation_rate, nameplate2, false, as_double("om_replacement_cost_escal")*0.01);
 
 			for ( i = 0; i < nyears && i < (int)count; i++) {
 				cf.at(CF_fuelcell_replacement_cost, i + 1) = fuelcell_rep[i] *
@@ -1002,7 +1023,20 @@ class cm_cashloan : public compute_module
 	assign("lcoptc_sta_nom", var_data((ssc_number_t) lcoptc_sta_nom));
 	assign("lcoptc_sta_real", var_data((ssc_number_t) lcoptc_sta_real));
 
-
+    ///////////////////////////////////////////////////////////////////////
+//LCOS Calculations
+    if (is_assigned("battery_total_cost_lcos") && as_double("battery_total_cost_lcos") != 0) {
+        for (int y = 0; y <= nyears; y++) {
+            cf_lcos.at(0, y) = cf.at(CF_battery_replacement_cost, y);
+            cf_lcos.at(1, y) = cf.at(CF_battery_replacement_cost_schedule, y);
+            cf_lcos.at(6, y) = cf.at(CF_om_fixed1_expense, y); //Fixed OM Battery cost
+            cf_lcos.at(7, y) = cf.at(CF_om_production1_expense, y); //Produciton OM Battery cost
+            cf_lcos.at(8, y) = cf.at(CF_om_capacity1_expense, y); //Capacity OM Battery Cost
+        }
+        int grid_charging_cost_version = 0;
+        lcos_calc(this, cf_lcos, nyears, nom_discount_rate, inflation_rate, lcoe_real, total_cost, real_discount_rate, grid_charging_cost_version, 0);
+    }
+    /////////////////////////////////////////////////////////////////////////////////////////
 
 	double wacc = 0.0;
 	wacc = (1.0 - debt_frac)*nom_discount_rate + debt_frac*loan_rate*(1.0 - cf.at(CF_effective_tax_frac, 1));