From 8779b5582b963757644d954a00fd4639e5934d98 Mon Sep 17 00:00:00 2001
From: irisR <zr92@cornell.edu>
Date: Sat, 17 Feb 2024 12:07:34 -0500
Subject: [PATCH] added two parameters including cost_perN and F_max to econ
 and simulation

---
 inputs/parameters.m       | 4 ++--
 simulation/modules/econ.m | 4 ++--
 simulation/simulation.m   | 4 ++--
 3 files changed, 6 insertions(+), 6 deletions(-)

diff --git a/inputs/parameters.m b/inputs/parameters.m
index 3acc2ee2..441f1b6b 100644
--- a/inputs/parameters.m
+++ b/inputs/parameters.m
@@ -50,6 +50,6 @@
             'LCOE_max', .5,...                  % maximum LCOE ($/kWh)
             'power_max', Inf,...                % maximum power (W)
             'eff_pto',   0.80,...               % PTO efficiency (-)
-            'eff_array', 0.95*0.98 );           % array availability and transmission efficiency (-)
-        
+            'eff_array', 0.95*0.98,...          % array availability and transmission efficiency (-)
+            'cost_perN', 0.1 );                 % cost per Newton   
 end
diff --git a/simulation/modules/econ.m b/simulation/modules/econ.m
index 578c7aa6..6a35aa99 100644
--- a/simulation/modules/econ.m
+++ b/simulation/modules/econ.m
@@ -1,5 +1,5 @@
 
-function [LCOE, capex, opex] = econ(m_m, M, cost_m, N_WEC, P_elec, FCR, efficiency)
+function [LCOE, capex, opex] = econ(m_m, M, cost_m, N_WEC, P_elec, FCR, cost_perN, F_max, efficiency)
 
 structural_cost = m_m.* cost_m;
 devicestructure =  N_WEC * structural_cost(M);
@@ -9,7 +9,7 @@
 development     = 4553000;
 infrastructure  = 990000;
 mooring         = N_WEC * 525000;
-pto             = N_WEC * 623000;
+pto             = N_WEC * (623000 + F_max * cost_perN);
 profitmargin    = 356000;
 installation    = 5909000;
 contingency     = 1590000;
diff --git a/simulation/simulation.m b/simulation/simulation.m
index f360e2a0..6f73b689 100644
--- a/simulation/simulation.m
+++ b/simulation/simulation.m
@@ -48,7 +48,7 @@
                                     in.M, in.h_s, in.T_s, in.rho_w, in.g, in.sigma_y, A_c, ...
                                     A_lat_sub, r_over_t, I, in.E);
 
-LCOE = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.eff_array);
+LCOE = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.cost_perN, in.F_max, in.eff_array);
 
 %% Assemble constraints g(x) >= 0
 g = zeros(1,18);
@@ -83,7 +83,7 @@
                                             in.t_fc, in.t_fb, in.t_sr, in.t_dt,...
                                             D_d, in.D_dt, in.theta_dt, in.T_s, in.h_d, ...
                                             in.M, in.rho_m, in.rho_w, in.m_scale);
-    [~,capex,opex] = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.eff_array);
+    [~,capex,opex] = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.cost_perN, in.F_max, in.eff_array);
     val.mass_f  = mass(1);
     val.mass_vc = mass(2);
     val.mass_rp = mass(3);