Skip to content

GitLab

Commit 96a7d692 authored by @martinm's avatar @martinm
Browse files

Added biomethane options from solid biomass for industrial processes. Changed... 

Added biomethane options from solid biomass for industrial processes. Changed the model starting year from 2015 to 2020. Added further climate targets (climate neutrality in 2045).
parent 57df7719
Expand all Hide whitespace changes
Inline Side-by-side
src/BioFeedCost.m
View file @ 96a7d692
......@@ -90,7 +90,7 @@ function[BP]=BioFeedCost(time,bioprod,priceDevFactor,feedstockInputData,BioResid
% Prizes for residues [/GJ]
% Prizes for residues [EUR/GJ]
if Sen==0
for b=1:9
for t=1:time
......@@ -151,8 +151,11 @@ function[BP]=BioFeedCost(time,bioprod,priceDevFactor,feedstockInputData,BioResid
BP(:,16)=BP(:,15)+7; % Miscanthus briquettes
BP(:,17)=BP(:,15)+5; % Miscanthus pellets
% Add invest costs for "Biomethaneinspeiseanlage" according Excel "BiomassData"
BP(:,[8:11 18:23])=BP(:,[8:11 18:23])+5.5;
% Add specific costs for the conversion of biogas to biomethane in /GJ
BP(:,[10:11 18:23])=BP(:,[10:11 18:23])+0.015/0.0036;
% Add specific costs for the feed-in of biomethane into the grid in /GJ
BP(:,[8:11 18:23])=BP(:,[8:11 18:23])+0.03/0.0036;
end
......
src/BiomassData.mat
View file @ 96a7d692
No preview for this file type
src/Biomass_Data.xlsx
View file @ 96a7d692
No preview for this file type
src/GHGTargetSetting_8095.m
View file @ 96a7d692
......@@ -61,6 +61,12 @@ ghgti80l = interp1(ghgin(1:7,25),ghgin(1:7,26),1990:1:2050,'linear');
% interpolate for Industry emissions 95% scenario
ghgti95l = interp1(ghgin(1:7,27),ghgin(1:7,28),1990:1:2050,'linear');
% interpolate for heat based emissions 100% scenario
ghgth100 = interp1(ghgin(1:7,29),ghgin(1:7,30),1990:1:2050,'linear');
% interpolate for heat based emissions 95% in 2045 scenario
ghgth95_2045 = interp1(ghgin(1:7,31),ghgin(1:7,32),1990:1:2050,'linear');
% Calculation of GHG reduction compared to prervious year
ghgrede80l=zeros(1,61);
ghgrede95l=zeros(1,61);
......@@ -70,6 +76,8 @@ ghgredh80l=zeros(1,61);
ghgredh95l=zeros(1,61);
ghgredo80l=zeros(1,61);
ghgredo95l=zeros(1,61);
ghgredh100=zeros(1,61);
ghgredh95_2045=zeros(1,61);
for t=2:61
ghgrede80l(t)=1-ghgte80l(t)/ghgte80l(t-1);
ghgrede95l(t)=1-ghgte95l(t)/ghgte95l(t-1);
......@@ -79,8 +87,13 @@ for t=2:61
ghgredh95l(t)=1-ghgth95l(t)/ghgth95l(t-1);
ghgredo80l(t)=1-ghgto80l(t)/ghgto80l(t-1);
ghgredo95l(t)=1-ghgto95l(t)/ghgto95l(t-1);
ghgredh100(t)=1-ghgth100(t)/ghgth100(t-1);
ghgredh95_2045(t)=1-ghgth95_2045(t)/ghgth95_2045(t-1);
end
% Replace Nan by zero
ghgredh100(isnan(ghgredh100))=0;
% Calculation of percentage compared to 1990
ghg90e80l=zeros(1,61);
ghg90e95l=zeros(1,61);
......@@ -103,10 +116,14 @@ end
% Save vector for model
ghg80=ghgredh80l(26:61);
save('ScenarioData.mat','ghg80','-append');
ghg95=ghgredh95l(26:61);
save('ScenarioData.mat','ghg95','-append');
ghg80=ghgredh80l(31:61);
save('ScenarioData.mat','ghg80','-append');
ghg95=ghgredh95l(31:61);
save('ScenarioData.mat','ghg95','-append');
ghg100=ghgredh100(31:61);
save('ScenarioData.mat','ghg100','-append');
ghg95_2045=ghgredh95_2045(31:61);
save('ScenarioData.mat','ghg95_2045','-append');
% Plot GHG Emission comparison
figure (1);
......
src/Main.m
View file @ 96a7d692
......@@ -43,7 +43,7 @@ disp('Pre-calculations started');
%% Define Sets
load('SetList.mat');
time = 36; % t - time / number of years
time = 31; % 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
......@@ -133,7 +133,7 @@ tb.type='set';
tb.form='full';
tb.val=bioprodtech';
%Definition of which biomass types go into which biomass products
% Definition of which biomass types go into which biomass products
biotypebioprod=SetList.data.BioType2BioProduct(3:end,3:end);
biotypebioprod(isnan(biotypebioprod))=0;
......
src/OptimizationModule.gms
View file @ 96a7d692
......@@ -26,7 +26,7 @@ sets
t Time
$loadR t
t2034(t) Time before 2035 /1*20/
t2034(t) Time before 2035 /1*15/
i Technologies
$loadR i
......@@ -93,7 +93,7 @@ parameters
dcap(t,j) Heat demand per house or HS [GJ]
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
culstart(b) Crop cultivation portfolio in the first 5 years [ha]
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]
......@@ -213,9 +213,9 @@ ba1fct Which residue biomass types can be used for which biomass produ
ba2fct Limitation of energy crops land potential to biomass products
ba3fct Which fossil biomass types can be used for which biomass products
bufct Produced amount biomass products = sum of consumed biomass per technologies
bufct Produced amount of biomass products = sum of consumed biomass over the technologies
bustart set portfolio of energy crops in starting years
bustart set portfolio of digestible energy crops in starting years
bustartKup set portfolio of KUP in starting years
bustartMis set portfolio of Miscanthus in the starting years
bumax max increase of energy crops
......@@ -242,18 +242,19 @@ totghgfct.. ghgtot=e=sum((t,i,j,b),ghgf(t,i,j,b))+sum((t,i,j),ghgt(
demandfct(t,j).. d(t,j)=e=sum((i),v(t,i,j));
dcapfct(t,j).. d(t,j)=e=sum(i,nprod(t,i,j))*dcap(t,j);
nfct(t+1,i,m,j) $ (ord(t)<36).. ncap(t+1,i,m,j)=e=ncap(t,i,m,j)+next(t+1,i,m,j)-ndec(t+1,i,m,j);
nfct(t+1,i,m,j) $ (ord(t)<31).. ncap(t+1,i,m,j)=e=ncap(t,i,m,j)+next(t+1,i,m,j)-ndec(t+1,i,m,j);
ncapfct(t,i,m,j).. ncap(t,i,m,j)=e=ncap1(t,i,m,j)+ncap2(t,i,m,j);
ncap2fct(t,i,m,j) $ (ord(t)>1).. ncap2(t,i,m,j)=e=ncap(t,i,m,j)-nprod(t,i,j);
ncap3fct(t+1,i,j) $ (ord(t)<36).. ncap2(t+1,i,"1",j)=g=ncap2(t,i,"1",j)-nsdec(t+1,i,"1",j)-nxdec(t+1,i,"1",j);
ncap3fct(t+1,i,j) $ (ord(t)<31).. ncap2(t+1,i,"1",j)=g=ncap2(t,i,"1",j)-nsdec(t+1,i,"1",j)-nxdec(t+1,i,"1",j);
ncap2ctrl1(t).. sum((i,j),ncap2(t,i,"1",j))=l=0.01*sum((i,j),nprod(t,i,j));
ndecfct(t,i,m,j) $ (ord(t)>1).. ndec(t,i,m,j)=e=nsdec(t,i,m,j)+nxdec(t,i,m,j);
nxdec1fct(t,i,m,j) $ (ord(t)+life(i,m,j)<37) .. nxdec(t+life(i,m,j),i,m,j)=e=next(t,i,m,j);
nxdec1fct(t,i,m,j) $ (ord(t)+life(i,m,j)<32) .. nxdec(t+life(i,m,j),i,m,j)=e=next(t,i,m,j);
nocfct(t,i,m,j).. nprod(t,i,j)=l=ncap1(t,i,m,j);
n1fct(t,i,m,j).. nprod(t,i,j)=e=ncap1(t,i,"1",j);
mbioprodfct(t,i,j).. nprod(t,i,j)*pmBio(t,i,j)*dcap(t,j)=e=sum(b,vBio(t,i,j,b));
mgasprodfct(t,i,j).. nprod(t,i,j)*pmGas(t,i,j)*dcap(t,j)=e=sum(b,vGas(t,i,j,b));
m3prodfct(t,i,j).. nprod(t,i,j)*pm3(t,i,j)*dcap(t,j)=e=v3(t,i,j);
......@@ -267,16 +268,17 @@ bamaxwastefct(t).. sum(bmwaste(bm)
ba1fct(t,bmwaste).. ba(t,bmwaste)=g=sum(b,bu(t,b,bmwaste));
ba2fct(t).. ba(t,"13")*bamaxc(t)=g=sum(bcult(b),bu(t,b,"13")/yield(t,b));
ba3fct(t).. ba(t,"14")=g=sum(b,bu(t,b,"13"));
ba3fct(t).. ba(t,"14")=g=sum(b,bu(t,b,"14"));
bufct(t,b).. sum(bm,bu(t,b,bm))=e=sum((i,j),bc(t,i,j,b));
bustart(bculst).. bu("1",bculst,"13")=e=yield("1",bculst)*culstart(bculst);
bustartKup.. sum((bkup,bm),bu("1",bkup,bm))=e=yield("1","12")*culstart("12");
bustartMis.. sum((bmis,bm),bu("1",bmis,bm))=e=yield("1","15")*culstart("15");
bumax(t+1,bculst) $ (ord(t)<36).. bu(t+1,bculst,"13")=l=2*bu(t,bculst,"13");
bumaxKup(t+1) $ (ord(t)<36).. sum(bkup,bu(t+1,bkup,"13"))=l=2*sum(bkup,bu(t,bkup,"13"));
bumaxMis(t+1) $ (ord(t)<36).. sum(bmis,bu(t+1,bmis,"13"))=l=2*sum(bmis,bu(t,bmis,"13"));
bumax(t+1,bculst) $ (ord(t)<31).. bu(t+1,bculst,"13")=l=2*bu(t,bculst,"13");
bumaxKup(t+1) $ (ord(t)<31).. sum(bkup,bu(t+1,bkup,"13"))=l=2*sum(bkup,bu(t,bkup,"13"));
bumaxMis(t+1) $ (ord(t)<31).. sum(bmis,bu(t+1,bmis,"13"))=l=2*sum(bmis,bu(t,bmis,"13"));
ghgffct(t,i,j,b).. ghgf(t,i,j,b)=e=alloc(i,j)*ghgfeed(b)*bc(t,i,j,b);
ghgtfct(t,i,j).. ghgt(t,i,j)=e=alloc(i,j)*ghgr(t,i,j)*v(t,i,j);
......
src/OptimizationModule.lst
View file @ 96a7d692
This diff is collapsed.
src/Plotting.m
View file @ 96a7d692
......@@ -184,8 +184,8 @@ if fig1==1
end
xlim([0 time+1])
ax=gca;
ax.XTick = [1 5:10:36];
ax.XTickLabel=[2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(titlefig1,'FontSize',FT(1));
ylabel(ylabheat,'FontSize',FT(2));
......@@ -213,8 +213,8 @@ if fig2==1
xlim([0 time+1])
ylim([0 max(sum(v(:,MT{j},j)/1000000,2))+max(sum(v(:,MT{j},j)/1000000,2))/20])
ax=gca;
ax.XTick = [1 6:10:36];
ax.XTickLabel= [2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(legmarket{j},'FontSize',FT(1));
l=legend(h(end:-1:1),legtech{MT{j}(end:-1:1)},'Location','southwest');
......@@ -240,8 +240,8 @@ if fig2==1
xlim([0 time+1])
ylim([0 max(sum(v(:,MT{j},j)/1000000,2))+max(sum(v(:,MT{j},j)/1000000,2))/20])
ax=gca;
ax.XTick = [1 6:10:36];
ax.XTickLabel= [2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(legmarket{j},'FontSize',FT(1));
l=legend(h(end:-1:1),legtech{MT{j}(end:-1:1)},'Location','southwest');
......@@ -275,8 +275,8 @@ if figBeh==1
xlim([0 time+1])
ax=gca;
if j==5
ax.XTick = [1 6:10:36];
ax.XTickLabel= [2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
else
ax.XTickLabel={[]};
......@@ -327,8 +327,8 @@ if fig3==1
hold on
xlim([0 time+1])
ax=gca;
ax.XTick = [1 6:10:36];
ax.XTickLabel=[2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(legtech{i},'FontSize',FT(1));
l=legend(h(end:-1:1),legbioprod{TB{i}(end:-1:1)},'Location','southwest');
......@@ -352,8 +352,8 @@ if fig4==1
xlim([0 time+1])
ylim([0 max(sum(ba.val(:,1:12)/1000000,2))+100])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title1fig3,'FontSize',FT(1));
ylabel(ylabbiomass,'FontSize',FT(2));
......@@ -369,8 +369,8 @@ if fig4==1
xlim([0 time+1])
ylim([0 max(sum(ba.val(:,1:12)/1000000,2))+100])
ax=gca;
ax.XTick = [1 6:10:36];
ax.XTickLabel=[2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title2fig3,'FontSize',FT(1));
ylabel(ylabbiomass,'FontSize',FT(2));
......@@ -385,8 +385,8 @@ if fig4==1
plot(squeeze(sum(squeeze(bu(:,:,13))./yield.val(:,:),2,'omitnan')/(1000000)),'--','color','black')
xlim([0 time+1])
ax=gca;
ax.XTick = [1 6:10:36];
ax.XTickLabel=[2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title3fig3,'FontSize',FT(1));
ylabel(ylabarea,'FontSize',FT(2));
......@@ -404,8 +404,8 @@ if fig5==1
h=area(squeeze(sum(sum(bc(:,:,:,1:23),2),4))/1000000);
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title1fig31,'FontSize',FT(1));
ylabel(ylabbiomass,'FontSize',FT(2));
......@@ -419,8 +419,8 @@ if fig5==1
end
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title2fig31,'FontSize',FT(1));
ylabel(ylabbiomass,'FontSize',FT(2));
......@@ -438,8 +438,8 @@ if fig6==1
end
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title('Biomass distribution in the technologies','FontSize',FT(1));
ylabel(ylabbiomass,'FontSize',FT(2));
......@@ -489,8 +489,8 @@ if fig7==1
plot(costyearly/1000000000,'color','k','LineWidth',2)
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title1fig4,'FontSize',FT(1));
xlabel(xlabtime,'FontSize',FT(2));
......@@ -511,8 +511,8 @@ if fig7==1
end
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
hold on
h=bar(Ypos/1000000000,'stack');
......@@ -536,8 +536,8 @@ if fig8==1
end
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
hold on;
line(1:time,ghgmax.val(:)/1000000,'color','black');
......@@ -558,8 +558,8 @@ if fig9==1
subplot(2,1,1)
plot(100*sum(sum(ncap2(:,:,1,:),2),4)./sum(sum(nprod(:,:,:),2),3))
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title('Overcapacity in %/a','FontSize',FT(1));
ylabel('%','FontSize',FT(2));
......@@ -571,8 +571,8 @@ if fig9==1
h(i).FaceColor=coltech{i};
end
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title('Number of heating systems','FontSize',FT(1));
ylabel('Number of model heating systems in Mio','FontSize',FT(2));
......@@ -592,8 +592,8 @@ if fig10==1
hold on
xlim([0 time+1])
ax=gca;
ax.XTick = [1 6:10:36];
ax.XTickLabel=[2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(legmarket{j},'FontSize',8);
ylabel('ncap2','FontSize',FT(2));
......@@ -619,8 +619,8 @@ if fig11==1
hold on
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(['Overcapacity (ncap2) - ' legmarket{j}],'FontSize',FT(1));
ylabel('Number heating systems','FontSize',FT(2));
......@@ -637,8 +637,8 @@ if fig11==1
hold on
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(['Used capacity (ncap1) - ' legmarket{j}],'FontSize',FT(1));
ylabel('Number heating systems','FontSize',FT(2));
......@@ -655,8 +655,8 @@ if fig11==1
hold on
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(['Decommisioned from the newly invested (nxdec) - ' legmarket{j}],'FontSize',FT(1));
ylabel('Number heating systems','FontSize',FT(2));
......@@ -673,8 +673,8 @@ if fig11==1
hold on
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(['Newly invested in (next) - ' legmarket{j}],'FontSize',FT(1));
ylabel('Number heating systems','FontSize',FT(2));
......@@ -694,8 +694,8 @@ if fig11==1
plot(sum(nxdec(:,:,1,j),2))
plot(1,sum(nstart.val(:,j)),'*')
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(['All n`s - ' legmarket{j}],'FontSize',FT(1));
ylabel('Number heating systems','FontSize',FT(2));
......@@ -713,8 +713,8 @@ if fig12==1
subplot(2,1,1)
bar(squeeze(sum(ncap(:,CapTech,TM{CapTech},:),4)),'DisplayName','ncap')
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(['Capacity of plant modules in ' legtech{CapTech}],'FontSize',FT(1));
ylabel('Number of plants','FontSize',FT(2));
......@@ -725,8 +725,8 @@ if fig12==1
subplot(2,1,2);
bar(squeeze(sum(next(:,CapTech,TM{CapTech},:),4)),'DisplayName','next')
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(['Investments in ' legtech{CapTech}],'FontSize',FT(1));
ylabel('Number of plants','FontSize',FT(2));
......@@ -745,8 +745,8 @@ if fig13==1
plot(GasPrice*0.36)
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title1fig9,'FontSize',FT(1));
ylabel('Cent/kWh','FontSize',FT(2));
......@@ -758,8 +758,8 @@ if fig13==1
plot(PowerPrice(:,[1 9 15])*0.36)
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title2fig9,'FontSize',FT(1));
ylabel('Cent/kWh','FontSize',FT(2));
......@@ -771,8 +771,8 @@ if fig13==1
h1=plot(COCert);
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title3fig9,'FontSize',FT(1));
ylabel(ylabcost2,'FontSize',FT(2));
......@@ -783,8 +783,8 @@ if fig13==1
plot(yield.val(:,[10 11 12 15 18:23]))
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title(title4fig9,'FontSize',FT(1));
ylabel(ylabyield,'FontSize',FT(2));
......@@ -794,53 +794,43 @@ end
%% Plot for investment and variable costs over time / markets/ technologies / biomassproducts
%% Plot for Gestehungskosten over time / markets/ technologies / biomassproducts
if fig14==1
for j=nMarketOp
figure (14)
set(gcf,'position',figsize)
k=1;
for i=[MT{j} 1000]
if length(MT{j})+1<=4
for i=MT{j}
if length(MT{j})<=4
subplot(2,2,k);
elseif length(MT{j})+1>4 && length(MT{j})+1<=6
elseif length(MT{j})>4 && length(MT{j})<=6
subplot(2,3,k);
elseif length(MT{j})+1>6 && length(MT{j})+1<=9
elseif length(MT{j})>6 && length(MT{j})<=9
subplot(3,3,k);
elseif length(MT{j})+1>9 && length(MT{j})+1<=12
elseif length(MT{j})>9 && length(MT{j})<=12
subplot(3,4,k);
elseif length(MT{j})>12 && length(MT{j})<=16
subplot(4,4,k);
end
k=k+1;
% Investment costs
if i==1000
plot(sum(inv.val(:,MT{j},:,j),3)/1000)
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.FontSize=FT(4);
title('Investment costs including subsidies (if selected)','FontSize',FT(1));
ylabel('Costs [t]','FontSize',FT(2));
l=legend(legtech{MT{j}},'Location','northeast');
% Gestehungskosten = variable costs + invest costs
if isempty(TB{i})==1
plot(3.6*squeeze(vc.val(:,i,j,1)++sum(inv.val(:,i,:,j),3)./dcap.val(:,j)))
else
plot(3.6*squeeze(vc.val(:,i,j,TB{i})+sum(inv.val(:,i,:,j),3)./dcap.val(:,j)))
end
xlim([0 time+1])
ax=gca;
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.FontSize=FT(4);
title([legtech{i}],'FontSize',FT(1));
ylabel('Costs [/MWh]','FontSize',FT(2));
if length(TB{i})>1
l=legend(legbioprod{TB{i}},'Location','northwest');
l.FontSize=FT(3);
else % variable costs
if isempty(TB{i})==1
plot(squeeze(vc.val(:,i,j,1)))
else
plot(squeeze(vc.val(:,i,j,TB{i})))
end
xlim([0 time+1])
ax=gca;
ax.XTick = 1:5:36;
ax.XTickLabel=2015:5:2050;
ax.FontSize=FT(4);
title([legtech{i}],'FontSize',FT(1));
ylabel('Costs [/GJ]','FontSize',FT(2));
if length(TB{i})>1
l=legend(legbioprod{TB{i}},'Location','northwest');
l.FontSize=FT(3);
end
end
end
end
......@@ -858,8 +848,8 @@ if fig15==1
hold on
xlim([0 time+1])
ax=gca;
ax.XTick = [1 6:10:36];
ax.XTickLabel=[2015 2020:10:2050];
ax.XTick = 1:10:31;
ax.XTickLabel=2020:10:2050;
ax.YTick = [5 10:10:100];
ax.FontSize=FT(4);
ylabel('Costs (/GJ)','FontSize',FT(2));
......
src/Results/Temp.mat
View file @ 96a7d692
No preview for this file type
src/ScenarioData.mat
View file @ 96a7d692
No preview for this file type
src/SetList.mat
View file @ 96a7d692
No preview for this file type
src/SetList.xlsx
View file @ 96a7d692
No preview for this file type
src/SetParameter.m
View file @ 96a7d692
This diff is collapsed.
src/Sets.mat
View file @ 96a7d692
No preview for this file type