The UFZ services GitLab and Mattermost will be unavailable on Monday, July 4 from 06:00 AM to 08:00 AM due to maintenance work.

Commit 390cceaa authored by Markus Millinger's avatar Markus Millinger
Browse files

Updated plotting of merit order and added dummy colours

parent 1d734eb8
......@@ -27,9 +27,9 @@ close all;
%% Data Import
% Conversion options, biomass crops and residues and other data
techData = ...
importdata('../data/BENOPT_inputData_200409.xls'); %Save to variable - comment out if the .mat-file is up to date
save('../data/techData.mat','techData'); %Save to .mat-file - comment out if the .mat-file is up to date
% techData = ...
% importdata('../data/BENOPT_inputData_200409.xls'); %Save to variable - comment out if the .mat-file is up to date
% save('../data/techData.mat','techData'); %Save to .mat-file - comment out if the .mat-file is up to date
techData = ...
importdata('../data/techData.mat'); %reads data from .mat file
......
......@@ -497,7 +497,7 @@ end
%% Merit order GHG/cost
% j=0;
% figure()
for scenario=1:6%:6%[1]
for scenario=1:1%:6%[1]
figure()
yearVector=[2030,2050];
plotNum=size(yearVector,2);
......@@ -506,10 +506,10 @@ for scenario=1:6%:6%[1]
% j=j+1;
% subplot(1,6,j)
subplot(1,plotNum,i)
paretoIter = 2;
paretoIter = 1;
year = yearVector(i);%2050;
timePoint = year-2019;
techIndex = g.fuelTypeDef(:,3)==0 & g.fuelTypeDef(:,4)==0 & g.techsectors(:,end)==0 & g.fuelTypeDef(:,7)==0; %Exclude LNG, CH4 and intermediate H2 AND EVs
techIndex = g.fuelTypeDef(:,3)==0 & g.fuelTypeDef(:,4)==0 & g.techsectors(:,end)==0;% & g.fuelTypeDef(:,7)==0; %Exclude LNG, CH4 and intermediate H2 AND EVs
prod = paretoVar(scenario,paretoIter).ghgAbateTech(timePoint,techIndex);
prod(prod<0) = s.referenceGHGemission.*paretoVar(scenario,paretoIter).prd2(timePoint,prod<0) - prod(prod<0); %Abatement of intermediate CH4 is without GHG reference and thus negative, this has to be corrected here
prod = prod*10^-3;
......@@ -540,7 +540,7 @@ for scenario=1:6%:6%[1]
yLabel = [''];
end
fileName = ['../fig/meritOrderGHG_S' num2str(scenario) '_' num2str(year)];
yMax = 700;
yMax = 2000;
xMax = sum(prod);
% if i==2
% tit = ['Merit order of fuel production, Scenario ' num2str(scenario)];
......
......@@ -57,5 +57,11 @@ col = [2,175,198;...
50,50,50
71,94,251;...
250,250,20;...
50,50,50
71,94,251;...
250,250,20;...
50,50,50
71,94,251;...
250,250,20;...
50,50,50]/255;
end
\ No newline at end of file
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