Commit 0b7fd464 authored by @martinm's avatar @martinm
Browse files

Issues by Markus Millinger addressed

parent 09b03b88
......@@ -43,15 +43,15 @@ disp('Pre-calculations started');
%% Define Sets
load('SetList.mat');
time = 36; % t
tech = find(cellfun(@isempty,SetList.textdata.SetList(:,1))==1,1)-2; % i
techtype = find(cellfun(@isempty,SetList.textdata.SetList(:,9))==1,1)-2; % tt
modul = find(cellfun(@isempty,SetList.textdata.SetList(:,5))==1,1)-2; % m
time = 36; % t - time / number of years
tech = find(cellfun(@isempty,SetList.textdata.SetList(:,1))==1,1)-2; % i - number of technology concepts
techtype = find(cellfun(@isempty,SetList.textdata.SetList(:,9))==1,1)-2; % tt - number of technology types
modul = find(cellfun(@isempty,SetList.textdata.SetList(:,5))==1,1)-2; % m - number of technology modules
market = find(cellfun(@isempty,SetList.textdata.SetList(:,13))==1,1)-2; % j refers to the sub-sectors, described in the publications
cluster=3; %c
biotype = find(cellfun(@isempty,SetList.textdata.SetList(:,18))==1,1)-2; % bt
bioprod = find(cellfun(@isempty,SetList.textdata.SetList(:,22))==1,1)-2; % b
GHGtechtype = find(cellfun(@isempty,SetList.textdata.SetList(:,27))==1,1)-2; % gtt
biotype = find(cellfun(@isempty,SetList.textdata.SetList(:,18))==1,1)-2; % bt - biomass type
bioprod = find(cellfun(@isempty,SetList.textdata.SetList(:,22))==1,1)-2; % b - biomass product
GHGtechtype = find(cellfun(@isempty,SetList.textdata.SetList(:,27))==1,1)-2; % gtt - technology types defined for the calculation of the emission factors
% Reforming for GAMS communication
tims= strsplit(num2str(1:time));
......@@ -194,22 +194,22 @@ end
irgdx('idxdata')
% Renaming
v=vp;
vBio=vBiop; %#ok<*NASGU>
vGas=vGasp;
v3=v3p;
bu=bup;
bc=bcp;
ghgf=ghgfp;
ghgt=ghgtp;
ncap=ncapp;
ncap1=ncap1p;
ncap2=ncap2p;
next=nextp;
nprod=nprodp;
nxdec=nxdecp;
vBeh=vBehp;
tc=tcp;
v=vp; % Heat production [GJ]
vBio=vBiop; % Solid Biomass heat production [GJ]
vGas=vGasp; % gas_biogas_coal heat production [GJ]
v3=v3p; % Non-Biomass heat production [GJ]
bu=bup; % Actual converted biomass from biotype to bioprod [GJ]
bc=bcp; % Actual consumed biomass in the technology [GJ]
ghgf=ghgfp; % Actual feedstock GHG emissions [t]
ghgt=ghgtp; % Actual technology GHG emissions [t]
nprod=nprodp; % Number of HS producing heat
ncap=ncapp; % Number of existing HS per technology (Capacity)
ncap1=ncap1p; % Number of existing HS used for production
ncap2=ncap2p; % Overcapacity of existing HS
next=nextp; % Number of heating systems\modules extended
nxdec=nxdecp; % Number of HS of next that reach their lifetime
vBeh=vBehp; % Heat production in the clusters [GJ]
tc=tcp; % Total system costs
clear vp vBiop vGasp v3p bup bcp ncapp ncap1p ncap2p nextp nprodp nxdecp tcp ghgfp ghgtp vBehp
......
......@@ -72,8 +72,8 @@ $loadR BB
;
parameters
vc(t,i,j,b) Variable costs [:GJ]
inv(t,i,m,j) Investment cost per plant []
vc(t,i,j,b) Variable costs [EUR:GJ]
inv(t,i,m,j) Investment cost per plant [EUR]
pmBio(t,i,j) biomass share per technology [%]
pmGas(t,i,j) gas_coal share per technology [%]
pm3(t,i,j) non biomass_gas_coal share per technology [%]
......@@ -81,7 +81,7 @@ parameters
efGas(t,i,j) Conversion efficiency gas_biogas
efMethan(t,b) Conversion efficiency for the "biomethaneinspeiseanlage"
life(i,m,j) lifetime of heating system [a]
ba(t,bm) Available biomass [GJ and ha]
ba(t,bm) Available biomass from residues (bmwaste) [GJ] and available land for cultivation [ha]
bamaxw(t) Maximal allowed biomass usage from waste ba [%]
bamaxc(t) Maximal allowed biomass usage from cultivation ba [%]
yield(t,b) Yield of cultivation products [GJ:ha]
......@@ -94,9 +94,9 @@ parameters
nstart(i,j) Initial stock of HS
nsdec(t,i,m,j) Yearly decrease of initial stock of HS
culstart(b) Crop cultivation portfolio in the first 5 years
vcBeh(t,i,j,c) Intangible variable costs [:GJ]
invBeh(t,i,j,c) Intangible investment costs []
dBeh(t,j,c) Demand in the clusters
vcBeh(t,i,j,c) Intangible variable costs [EUR:GJ]
invBeh(t,i,j,c) Intangible investment costs [EUR]
dBeh(t,j,c) Demand in the clusters [GJ]
* definition of parameters for the export (equal the variables)
vp(t,i,j) Heat production [GJ]
......@@ -110,18 +110,18 @@ parameters
nprodp(t,i,j) Number of HS producing heat
ncapp(t,i,m,j) Number of existing HS per technology (Capacity)
ncap1p(t,i,m,j) Number of existing HS used for production
ncap2p(t,i,m,j) Overcapacity of existing HS
ncap2p(t,i,m,j) Overcapacity of existing HS [number of HS]
nextp(t,i,m,j) Number of heating systems\modules extended
nxdecp(t,i,m,j) Number of HS of next that reach their lifetime
vBehp(t,i,j,c) Heat production in the clusters [GJ]
tcp
tcp Total system costs [EUR]
returnStat(stat);
;
$loadR d, dcap, vc, inv, pmBio, pmGas, pm3, efBio, efGas, efMethan , life, ba, bamaxw, bamaxc, nstart, nsdec, yield, culstart, ghgr, ghgfeed, alloc, ghgmax, vcBeh, invBeh, dBeh
free variable
tc Total costs []
tc Total costs [EUR]
ghgtot GHG total emission [t]
;
......
This diff is collapsed.
......@@ -441,7 +441,7 @@ if fig6==1
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.FontSize=FT(4);
title('Biomasseverteilung nach Technologien','FontSize',FT(1));
title('Biomass distribution in the technologies','FontSize',FT(1));
ylabel(ylabbiomass,'FontSize',FT(2));
legend(h(end:-1:1),legtech(end:-1:1),'Location','bestoutside','FontSize',FT(3));
end
......
No preview for this file type
No preview for this file type
No preview for this file type
No preview for this file type
No preview for this file type
$gdxin matdata.gdx
$loadR t
$loadR i
$loadR m
$loadR j
$loadR c
$loadR bm
$loadR b
$loadR TB
$loadR MT
$loadR BB
$loadR d
$loadR dcap
$loadR vc
$loadR inv
$loadR pmBio
$loadR pmGas
$loadR pm3
$loadR efBio
$loadR efGas
$loadR efMethan
$loadR life
$loadR ba
$loadR bamaxw
$loadR bamaxc
$loadR nstart
$loadR nsdec
$loadR yield
$loadR culstart
$loadR ghgr
$loadR ghgfeed
$loadR alloc
$loadR ghgmax
$loadR dBeh
$loadR vcBeh
$loadR invBeh
$gdxin
% Copyright (c) 2004, Ben Barrowes
% All rights reserved.
%
% Redistribution and use in source and binary forms, with or without
% modification, are permitted provided that the following conditions are met:
%
% * Redistributions of source code must retain the above copyright notice, this
% list of conditions and the following disclaimer.
%
% * Redistributions in binary form must reproduce the above copyright notice,
% this list of conditions and the following disclaimer in the documentation
% and/or other materials provided with the distribution
% * Neither the name of Barrowes Consulting / Dartmouth College / US Army CRREL nor the names of its
% contributors may be used to endorse or promote products derived from this
% software without specific prior written permission.
% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
% DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE
% FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
% DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
% SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
% CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
% OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
% OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
function [ax,h]=suplabel(text,whichLabel,supAxes)
% PLaces text as a title, xlabel, or ylabel on a group of subplots.
% Returns a handle to the label and a handle to the axis.
......
Markdown is supported
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