Commit 9fdf399d authored by Toni's avatar Toni
Browse files

Fixing plot tail, so that dead time is not plotted

parent 55631374
......@@ -12,7 +12,7 @@ $offtext
* Set TRUE or FALSE
$set archive FALSE
* Set file name prefix for archivation
$set archive_prefix 2022-01-12_Testing
$set archive_prefix 2022-01-XX_Testing
$set archive_exec 0
$if "%archive%"=="TRUE" $set archive_exec 1
......@@ -72,7 +72,8 @@ Parameter
;
* Prepare priceElspot
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, ts_priceElspot;
execute_loaddc "%backbone_output_GDX%", mSettings, ts_priceElspot ;
tparam(t)=(ord(t)<=(min(mSettings('building','dataLength'),mSettings('building','t_end')) - (mSettings('building','t_horizon')-mSettings('building','t_jump'))));
%GDXparam2%(tparam,%GDXparam2set2%)=ts_priceElspot(tparam) +eps ;
$set Ylabel2_alt1 "Price of electricity (EUR/MWh) :"
*Prepare outside temperature
......@@ -97,7 +98,7 @@ $set TMPparam ts_influx_absolute
$set Ylabel "Track exogenous commodities (kWh) :"
$set Ylabel2 %Ylabel2_alt2%
Parameter %TMPparam%(t,node) "External power (weather) influencing building" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node)=sum((grid,node_building2node(node_building,node),f)$(sameas(grid,"building") and sameas(f,'f00')),w2kw*%GDXparam%(grid,node,f,t) + eps);
execute_unload 'plots/%GDXfile%', tparam=t, node, %TMPparam%, %GDXparam2set2alt2%, %GDXparam2alt2% ;
......@@ -109,7 +110,7 @@ $set TMPparam ts_influx_normalized
$set Ylabel "Track exogenous commodities per m2 (kWh) :"
$set Ylabel2 %Ylabel2_alt2%
Parameter %TMPparam%(t,node) "External power (weather) influencing building - normalized by bulding m2" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node)=sum((grid,node_building2node(node_building,node),f)$(sameas(grid,"building") and sameas(f,'f00')), w2kw*%GDXparam%(grid,node,f,t)/building_squares(node_building) + eps);
execute_unload 'plots/%GDXfile%', tparam=t, node, %TMPparam%, %GDXparam2set2alt2%, %GDXparam2alt2% ;
......@@ -124,7 +125,7 @@ $set TMPparam r_gen_heat_absolute
$set Ylabel Electricity usage (kW) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,unit) "Electric power consumption for heating and cooling of air" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),unit_heat_and_cool(unit))=sum((grid,node,f)$(sameas(grid,"elec") and sameas(node,'link_node') and sameas(f,'f00')),sign* w2kw* %GDXparam%(grid,node,unit,f,t) + eps);
execute_unload 'plots/%GDXfile%', tparam=t, unit, %TMPparam%, %GDXparam2set2%, %GDXparam2% ;
......@@ -136,7 +137,7 @@ $set TMPparam r_gen_heat_normalized
$set Ylabel Electricity usage per m2(kW) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,unit) "Electric power consumption for heating and cooling of air - normalized by bulding m2" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),unit_heat_and_cool(unit))=sum((grid,node,node_building2unit(node_building,unit),f)$(sameas(grid,"elec") and sameas(node,'link_node') and sameas(f,'f00')),sign* w2kw* %GDXparam%(grid,node,unit,f,t)/building_squares(node_building) + eps);
execute_unload 'plots/%GDXfile%', tparam=t, unit, %TMPparam%, %GDXparam2set2%, %GDXparam2% ;
......@@ -148,7 +149,7 @@ $set TMPparam r_gen_dhw_absolute
$set Ylabel Electricity usage (kW) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,unit) "Electric power consumption for domestic hot water" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),unit_DHW(unit))=sum((grid,node,f)$(sameas(grid,"elec") and sameas(node,'link_node') and sameas(f,'f00')),sign* w2kw* %GDXparam%(grid,node,unit,f,t) + eps);
execute_unload 'plots/%GDXfile%', tparam=t, unit, %TMPparam%, %GDXparam2set2%, %GDXparam2% ;
......@@ -160,7 +161,7 @@ $set TMPparam r_gen_dhw_normalized
$set Ylabel Electricity usage per m2(kW) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,unit) "Electric power consumption for domestic hot water - normalized by bulding m2" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),unit_DHW(unit))=sum((grid,node,node_building2unit(node_building,unit),f)$(sameas(grid,"elec") and sameas(node,'link_node') and sameas(f,'f00')),sign* w2kw* %GDXparam%(grid,node,unit,f,t)/building_squares(node_building) + eps);
execute_unload 'plots/%GDXfile%', tparam=t, unit, %TMPparam%, %GDXparam2set2%, %GDXparam2% ;
......@@ -175,7 +176,7 @@ $set TMPparam r_state_interior_air
$set Ylabel Temperture (Celcius) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,*) "Temperature at interior_air_and_furniture (Celcius)" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node_building_interior_air_and_furniture(node))=sum((grid,f)$(sameas(grid,"building") and sameas(f,'f00')),%GDXparam%(grid,node,f,t) - c2k + eps);
*%TMPparam%(tparam(t),'out_temp')=16;
......@@ -188,7 +189,7 @@ $set TMPparam r_state_internal_mass
$set Ylabel Temperture (Celcius) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,*) "Temperature at internal_mass: inside walls (Celcius)" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node_building_internal_mass(node))=sum((grid,f)$(sameas(grid,"building") and sameas(f,'f00')),%GDXparam%(grid,node,f,t) - c2k + eps);
*%TMPparam%(tparam(t),'out_temp')=16;
......@@ -201,7 +202,7 @@ $set TMPparam r_state_envelope_mass
$set Ylabel Temperture (Celcius) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,*) "Temperature at envelope_mass: outside walls (Celcius)" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node_building_envelope_mass(node))=sum((grid,f)$(sameas(grid,"building") and sameas(f,'f00')),%GDXparam%(grid,node,f,t) - c2k + eps);
*%TMPparam%(tparam(t),'out_temp')=16;
......@@ -214,7 +215,7 @@ $set TMPparam r_state_DHWT
$set Ylabel Temperture (Celcius) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,*) "Temperature at Domestic Hot Water Tank (Celcius)" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node_building_DHWT(node))=sum((grid,f)$(sameas(grid,"building") and sameas(f,'f00')),%GDXparam%(grid,node,f,t) - c2k + eps);
*%TMPparam%(tparam(t),'out_temp')=16;
......@@ -229,7 +230,7 @@ $set TMPparam r_state_DH1
$set Ylabel Temperture (Celcius) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,*) "Temperature at Domestic Hot Water Tank (Celcius)" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node)$(node_building2node("IDA_ESBO_DH1",node)) = sum((grid,f)$(sameas(grid,"building") and sameas(f,'f00')),%GDXparam%(grid,node,f,t) - c2k + eps);
*%TMPparam%(tparam(t),'out_temp')=16;
......@@ -242,7 +243,7 @@ $set TMPparam r_state_DH2
$set Ylabel Temperture (Celcius) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,*) "Temperature at Domestic Hot Water Tank (Celcius)" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node)$(node_building2node("IDA_ESBO_DH2",node)) = sum((grid,f)$(sameas(grid,"building") and sameas(f,'f00')),%GDXparam%(grid,node,f,t) - c2k + eps);
*%TMPparam%(tparam(t),'out_temp')=16;
......@@ -255,7 +256,7 @@ $set TMPparam r_state_AB
$set Ylabel Temperture (Celcius) :
$set Ylabel2 %Ylabel2_alt1%
Parameter %TMPparam%(t,*) "Temperature at Domestic Hot Water Tank (Celcius)" ;
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, %GDXparam%;
execute_loaddc "%backbone_output_GDX%", %GDXparam%;
* Below line is adjusted for each paramGDX and paramPlot pair
%TMPparam%(tparam(t),node)$(node_building2node("IDA_ESBO_AB",node)) = sum((grid,f)$(sameas(grid,"building") and sameas(f,'f00')),%GDXparam%(grid,node,f,t) - c2k + eps);
*%TMPparam%(tparam(t),'out_temp')=16;
......
......@@ -34,7 +34,7 @@ error_unit_details(node,unit,t) "Show which unit supplies energy to grid "
;
*Load result file : r_gen(grid, node, unit,f,t)
execute_loaddc "%backbone_output_GDX%", tparam=t_realized, r_gen;
execute_loaddc "%backbone_output_GDX%", r_gen;
summary("elec_ideal_heating_kWh_by_m2",node_building) = sum((node_building2node(node_building,node),unit_heat(unit),tparam(t)), w2kw * r_gen("building", node, unit,'f00',t)/ building_squares(node_building)+ eps) ;
summary("elec_ideal_cooling_kWh_by_m2",node_building) = sum((node_building2node(node_building,node),unit_cool(unit),tparam(t)), w2kw * r_gen("building", node, unit,'f00',t)/ building_squares(node_building)+ eps);
......
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