Commit 55631374 authored by Toni's avatar Toni
Browse files

Added p_price_el_select to select flat or elspot price

parent 46788de9
......@@ -207,9 +207,10 @@ Parameters
* --- Other time dependent parameters -----------------------------------------
Parameters
p_stepLength(mType, f, t) "Length of an interval in hours"
p_stepLength(mType, f, t) "Length of an interval in hours"
p_stepLengthNoReset(mType, f, t) "Length of an interval in hours - includes also lengths of previously realized intervals"
p_s_discountFactor(s) "Discount factor for samples when using a multi-year horizon"
p_s_discountFactor(s) "Discount factor for samples when using a multi-year horizon"
p_price_el_select "Use flat price 1 EUR/MWh (set 0). Use elspot price (set 1)"
;
* --- Energy modelling for buildings -----------------------------------------
......
......@@ -104,6 +104,7 @@ $loaddc unit_DHW
$loaddc unit_heat_and_cool
$loaddc node_building2unit
$loaddc building_squares c2k
$loaddc p_price_el_select
$gdxin
......
......@@ -20,14 +20,17 @@ $offtext
* --- Objective Function Definition -------------------------------------------
* =============================================================================
Parameter
ts_priceElspotNode(node,t);
ts_priceElspotNode(node,t)=1E6;
ts_priceElspotNode(node,t) /#node.#t 1E6/;
*node_building_DHWT node_building_envelope_mass node_building_interior_air_and_furniture node_building_internal_mass
*ts_priceElspotNode(node_building_DHWT,t)=ts_priceElspot(t);
*ts_priceElspotNode(node_building_interior_air_and_furniture,t)=ts_priceElspot(t);
ts_priceElspotNode("link_node",t)=ts_priceElspot(t);
if(0=p_price_el_select,
ts_priceElspotNode("link_node",t)=1$ts_priceElspot(t);
elseif 1=p_price_el_select,
ts_priceElspotNode("link_node",t)=ts_priceElspot(t);
else
abort "Error, invalid value for p_price_el_select, set 0 or 1";
);
q_obj ..
+ v_obj * 1e6 // Objective function valued in MEUR instead of EUR (or whatever monetary unit the data is in)
......@@ -48,9 +51,7 @@ q_obj ..
+ sum(inc_dec,
+ sum(gn(grid, node),
+ vq_gen(inc_dec, grid, node, s, f, t)
* * ts_priceElspotNode(node,t) // Line A: Use line A or B by commenting out the other one. Line A using Elspot price
* 1 // Line B: Use line A or B by commenting out the other one. Line B: Flat price
* ts_priceElspotNode(node,t)
* original *( PENALTY_BALANCE(grid, node)${not p_gnBoundaryPropertiesForStates(grid, node, 'balancePenalty', 'useTimeSeries')}
* original + ts_node_(grid, node, 'balancePenalty', s, f, t)${p_gnBoundaryPropertiesForStates(grid, node, 'balancePenalty', 'useTimeSeries')}
* original )
......
......@@ -125,6 +125,9 @@ $gdxIn input\elspot_prices_2013.gdx
$load tsIso elspotIsoBB ts_priceElspot=elspotBB
$gdxin
*Setting
Parameter p_price_el_select "Use flat price 1 EUR/MWh (set 0). Use elspot price (set 1)" /1/;
* Write input with domains at execution time to take in consideration modifications
execute_unloaddi "input\inputDataAdjusted1.gdx";
execute_unloaddi "input\inputDataAdjusted1Squuezed.gdx" node_building2node, node_building2unit, ts_priceElspot ts_influx effLevelGroupUnit p_gnBoundaryPropertiesForStates p_gnu_io p_gn p_unit node p_gnn unit unitUnittype grid p_s_discountFactor unittype
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment