Skip to content

Commit

Permalink
Implementation review comments on Vanderkelen et al, 2018b
Browse files Browse the repository at this point in the history
  • Loading branch information
@Ivanderkelen
Ivanderkelen committed Oct 17, 2018
1 parent f044967 commit 3853edc
Show file tree
Hide file tree
Showing 84 changed files with 810,717 additions and 170 deletions.
4 changes: 2 additions & 2 deletions .gitignore
View file
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@

biascorrection/PFT_lin

data
biascorrection/QUANT
notused
23 changes: 14 additions & 9 deletions WBmodel/WBmodel.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,25 +5,31 @@
% define the starting lakelevel
define_L0;

% define minimum and maximum lake level
Lbounds = define_minmax(flag_outscen);

% determine Pel
determ_Pel;

if flag_run == 1
L0 = lakelevel(1);

[L_obs,Qout_ac] = solveWB(P_wb_obs, E_wb_obs, Qin_wb_obs, Qout, A_lake, L0, ndays,flag_run,flag_outscen,lakelevel);
[L_obs,Qout_ac] = solveWB(P_wb_obs, E_wb_obs, Qin_wb_obs, Qout, A_lake, L0, ndays,flag_run,flag_outscen,Pel,k,Lbounds);

elseif flag_run == 2
L0 = lakelevel(1);

for i = 1:nRCMs

[L_ev(i,:),~] = solveWB(P_wb_ev(i,:),E_wb_ev(i,:),Qin_wb_ev(i,:),Qout, A_lake, L0, ndays,flag_run,flag_outscen,lakelevel);
[L_ev(i,:),~] = solveWB(P_wb_ev(i,:),E_wb_ev(i,:),Qin_wb_ev(i,:),Qout, A_lake, L0, ndays,flag_run,flag_outscen,Pel,k,Lbounds);

end

elseif flag_run == 3 % historical
L0 = lakelevel_all(1);
L0 = lakelevel_all(1097);
for i = 1:nm
fprintf('Calculating WB for model %d of %d \n',i,nm);
[L_hist(i,:),Qout_hist(i,:)] = solveWB(P_wb_hist(i,:),E_wb_hist(i,:),Qin_wb_hist(i,:),Qout, A_lake, L0, ndays,flag_run,flag_outscen,lakelevel);
[L_hist(i,:),Qout_hist(i,:)] = solveWB(P_wb_hist(i,:),E_wb_hist(i,:),Qin_wb_hist(i,:),Qout, A_lake, L0, ndays,flag_run,flag_outscen,Pel,k,Lbounds);
end


Expand All @@ -34,32 +40,31 @@

fprintf('Calculating WB for model %d of %d \n',i,nm);

[L_rcp26(i,:), Qout_rcp26(i,:)] = solveWB(P_wb_rcp26(i,:),E_wb_rcp26(i,:),Qin_wb_rcp26(i,:),Qout, A_lake, L0(1,i), ndays,flag_run,flag_outscen,lakelevel);
[L_rcp26(i,:), Qout_rcp26(i,:)] = solveWB(P_wb_rcp26(i,:),E_wb_rcp26(i,:),Qin_wb_rcp26(i,:),Qout, A_lake, L0(1,i), ndays,flag_run,flag_outscen,Pel,k,Lbounds);
end

elseif flag_run == 5 % RCP 4.5

L0 = L0_rcp45;

for i = 1:nm

fprintf('Calculating WB for model %d of %d \n',i,nm);

[L_rcp45(i,:), Qout_rcp45(i,:)] = solveWB(P_wb_rcp45(i,:),E_wb_rcp45(i,:),Qin_wb_rcp45(i,:),Qout, A_lake, L0(1,i), ndays,flag_run,flag_outscen,lakelevel);
[L_rcp45(i,:), Qout_rcp45(i,:)] = solveWB(P_wb_rcp45(i,:),E_wb_rcp45(i,:),Qin_wb_rcp45(i,:),Qout, A_lake, L0(1,i), ndays,flag_run,flag_outscen,Pel,k,Lbounds);
end



elseif flag_run == 6 % RCP 8.5

L0 = L0_rcp85;

for i = 1:nm

fprintf('Calculating WB for model %d of %d \n',i,nm);

[L_rcp85(i,:), Qout_rcp85(i,:)] = solveWB(P_wb_rcp85(i,:),E_wb_rcp85(i,:),Qin_wb_rcp85(i,:),Qout, A_lake, L0(1,i), ndays,flag_run,flag_outscen,lakelevel);
[L_rcp85(i,:), Qout_rcp85(i,:)] = solveWB(P_wb_rcp85(i,:),E_wb_rcp85(i,:),Qin_wb_rcp85(i,:),Qout, A_lake, L0(1,i), ndays,flag_run,flag_outscen,Pel,k,Lbounds);


end

end
Expand Down
2 changes: 1 addition & 1 deletion WBmodel/calcWBterms_obs.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -23,7 +23,7 @@
for t = 1:ndays
% Calculate terms of WB
P_mean(t) = nanmean(nanmean(P_lake(:,:,t))); % average all precip [m]
if (isnan(P_mean(t))) P_mean(t) = 0; end
if (isnan(P_mean(t))||P_mean(t) <0) P_mean(t) = 0; end
E_mean(t) = nanmean(nanmean(E_lake(:,:,t))); % average all evap [m]
end

Expand Down
11 changes: 6 additions & 5 deletions WBmodel/def_timebounds.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -21,8 +21,8 @@
time_end = [2008,12,31,23,0,0];
end
elseif flag_run == 3
time_begin = [1950, 1, 1, 0,0,0];
time_end = [2004,12,31,23,0,0];
time_begin = [1951, 1, 1, 0,0,0];
time_end = [2005,12,31,23,0,0];
elseif flag_run == 4 || flag_run == 5 || flag_run == 6
time_begin = [2006, 1, 1, 0,0,0];
time_end = [2100,12,31,23,0,0];
Expand All @@ -40,7 +40,8 @@
% Observations
time_begin_obs = [1993, 1, 1, 0,0,0]; % datum van eerste dahiti data beschikb
time_end_obs = [2014,12,31,23,0,0];



date_vec_obs= datevec(datenum(time_begin_obs ):1:datenum(time_end_obs ));
date_obs = date_vec_obs(:,1:3);

Expand All @@ -59,8 +60,8 @@
date_bc= date_vec_bc(:,1:3);

% Historical
time_begin_hist = [1950, 1, 1, 0,0,0];
time_end_hist = [2004,12,31,23,0,0];
time_begin_hist = [1951, 1, 1, 0,0,0];
time_end_hist = [2005,12,31,23,0,0];

date_vec_hist= datevec(datenum(time_begin_hist):1:datenum(time_end_hist));
date_hist = date_vec_hist(:,1:3);
Expand Down
29 changes: 29 additions & 0 deletions WBmodel/define_minmax.m
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,29 @@
% ------------------------------------------------------------------------
% Function to define the minimum and maximum lake level for the different
% outflow scenarios
% ------------------------------------------------------------------------

function [Lbounds] = define_minmax(flag_outscen)

load diff_abs_jinja
dam_height = diff_abs_jinja + 31;
safety = 7;
hmax = dam_height-safety;
% Agreed Curve
if flag_outscen == 4
Lmin = 10 + diff_abs_jinja; % Source: figure
Lmax = 13.5 + diff_abs_jinja;
elseif flag_outscen == 1 % constant outflow max
Lmin = diff_abs_jinja;
Lmax = hmax ;
elseif flag_outscen == 2 % constant outflow min
Lmin = diff_abs_jinja;
Lmax = hmax ;
elseif flag_outscen == 3 % constant lake level
Lmin = 0;
Lmax = hmax ;
end

Lbounds = [Lmin Lmax];

end
47 changes: 47 additions & 0 deletions WBmodel/determ_Pel.m
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,47 @@
% ------------------------------------------------------------------------
% Script to calculate Pel (kW/day) and determine k: efficiency factor (kN/m³)
% ------------------------------------------------------------------------

% load absolute base eigt dam
load diff_abs_jinja.mat

Pel_m_nalu = 180000; %[kW] (value from globalenergyobservatory.org)
Pel_m_kiira = 200000; %[kW] (value from globalenergyobservatory.org)
Pel_m_sum = Pel_m_nalu + Pel_m_kiira;

Qmax_nalu = 1200; % [m³/s] (turbines; Kizza et al., 2006 )
Qmax_kiira = 1100; % [m³/s] (turbines; Kizza et al., 2006 )
Qmax_mean = (Qmax_nalu + Qmax_kiira)/2;

% simplification: water head = relative height of the water at the dam.
h_max = 24; % maximum water head (m) (similar for nalu and kiira)

% calculate efficiency factor k [kN/m³]
k = Pel_m_sum./ (Qmax_mean * h_max);

% Pel elektricity generation corresponding to mean observed outflow of historical
% period
Qout_mean_ms = mean(outflow)/(24*60*60); % in m³/s
h_mean = mean(lakelevel_all) - diff_abs_jinja;
Pel_hist = mean(Qout_mean_ms*k*h_mean);

% Pel for period for 1955 (minimum outflow)
[isdate_1955, date_loc_date] = ismember(date_all(:,1),1955);
Qout_mean_1955_ms = 55*10^6./(24*60*60);
h_mean_1955 = lakelevel_all(isdate_1955') - diff_abs_jinja;
Pel_1955 = mean(Qout_mean_1955_ms*k*h_mean_1955);

% Pel for period for 1964 (maximum outflow)
[isdate_1964, date_loc_date] = ismember(date_all(:,1),1964);
Qout_mean_1964_ms = 138*10^6./(24*60*60);
h_mean_1964 = lakelevel_all(isdate_1964') - diff_abs_jinja;
Pel_1964 = mean(Qout_mean_1964_ms*k*h_mean_1964);


if flag_outscen == 1
Pel = Pel_hist;
elseif flag_outscen == 2
Pel = Pel_1964;
else
Pel = 0;
end
14 changes: 7 additions & 7 deletions WBmodel/initialisation.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -41,9 +41,9 @@

elseif flag_run == 3 % historical

RCM(:,1) = [{'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'HIRHAM5_'} ; {'CRCM5_'} ; {'CRCM5_'} ; {'RACMO22T_'}; {'RACMO22T_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'}];
GCM(:,1) = [{'CNRM-CM5'} ; {'EC-EARTH'} ; {'HadGEM2-ES'} ; {'MPI-ESM-LR'} ; {'EC-EARTH'} ; {'MPI-ESM-LR'} ; {'CanESM2'}; {'EC-EARTH'} ; {'HadGEM2-ES'}; {'CanESM2'}; {'CM5A-MR'}; {'CNRM-CM5'}; {'EC-EARTH'}; {'GFDL-ESM2M'}; {'HadGEM2-ES'}; {'MIROC5'}; {'MPI-ESM-LR'}; {'NorESM1-M'} ; {'CSIRO-Mk3-6-0'}; {'HadGEM2-ES'} ; {'MPI-ESM-LR'} ; {'EC-EARTH'} ; {'CM5A-LR'} ; {'GFDL-ESM2G'} ; {'MIROC5'} ];
RCM_text(:,1) = [{'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'HIRHAM5 '}; {'CRCM5 '} ; {'CRCM5 '} ; {'RACMO22T '}; {'RACMO22T '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '}];
RCM(:,1) = [{'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CRCM5_'} ; {'CRCM5_'} ; {'HIRHAM5_'} ; {'RACMO22T_'}; {'RACMO22T_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'} ; {'REMO2009_'}];
GCM(:,1) = [{'CNRM-CM5'} ; {'EC-EARTH'} ; {'HadGEM2-ES'} ; {'MPI-ESM-LR'} ; {'MPI-ESM-LR'} ; {'CanESM2'}; {'EC-EARTH'} ; {'EC-EARTH'} ; {'HadGEM2-ES'}; {'CanESM2'}; {'CM5A-MR'}; {'CNRM-CM5'}; {'EC-EARTH'}; {'GFDL-ESM2M'}; {'HadGEM2-ES'}; {'MIROC5'}; {'MPI-ESM-LR'}; {'NorESM1-M'} ; {'CSIRO-Mk3-6-0'}; {'HadGEM2-ES'} ; {'MPI-ESM-LR'} ; {'EC-EARTH'} ; {'CM5A-LR'} ; {'GFDL-ESM2G'} ; {'MIROC5'} ];
RCM_text(:,1) = [{'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CRCM5 '} ; {'CRCM5 '} ; {'HIRHAM5 '}; {'RACMO22T '}; {'RACMO22T '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '} ; {'REMO2009 '}];

elseif flag_run == 4 % rcp 2.6

Expand All @@ -53,9 +53,9 @@

elseif flag_run == 5 % rcp 45

RCM(:,1) = [{'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'HIRHAM5_'} ; {'CRCM5_'} ; {'CRCM5_'} ;{'RACMO22T_'} ; {'RACMO22T_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'REMO2009_'} ; {'REMO2009_'} ];
GCM(:,1) = [{'CNRM-CM5'} ; {'EC-EARTH'} ; {'HadGEM2-ES'} ; {'MPI-ESM-LR'} ; {'EC-EARTH'} ; {'MPI-ESM-LR'} ; {'CanESM2'} ;{'EC-EARTH'} ; {'HadGEM2-ES'}; {'CanESM2'}; {'CM5A-MR'}; {'CNRM-CM5'}; {'EC-EARTH'}; {'GFDL-ESM2M'}; {'HadGEM2-ES'}; {'MIROC5'}; {'MPI-ESM-LR'}; {'NorESM1-M'} ; {'CSIRO-Mk3-6-0'}; {'MPI-ESM-LR'} ; {'EC-EARTH'} ];
RCM_text(:,1) = [{'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'HIRHAM5 '} ; {'CRCM5 '} ; {'CRCM5 '} ;{'RACMO22T '} ; {'RACMO22T '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'REMO2009 '} ; {'REMO2009 '} ];
RCM(:,1) = [{'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CCLM4-8-17_'}; {'CRCM5_'} ; {'CRCM5_'} ;{'HIRHAM5_'} ;{'RACMO22T_'} ; {'RACMO22T_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'RCA4_'} ; {'REMO2009_'} ; {'REMO2009_'} ];
GCM(:,1) = [{'CNRM-CM5'} ; {'EC-EARTH'} ; {'HadGEM2-ES'} ; {'MPI-ESM-LR'} ; {'MPI-ESM-LR'} ; {'CanESM2'} ; {'EC-EARTH'};{'EC-EARTH'} ; {'HadGEM2-ES'}; {'CanESM2'}; {'CM5A-MR'}; {'CNRM-CM5'}; {'EC-EARTH'}; {'GFDL-ESM2M'}; {'HadGEM2-ES'}; {'MIROC5'}; {'MPI-ESM-LR'}; {'NorESM1-M'} ; {'CSIRO-Mk3-6-0'}; {'MPI-ESM-LR'} ; {'EC-EARTH'} ];
RCM_text(:,1) = [{'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CCLM4-8-17 '}; {'CRCM5 '} ; {'CRCM5 '} ;{'HIRHAM5 '} ; {'RACMO22T '} ; {'RACMO22T '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'RCA4 '} ; {'REMO2009 '} ; {'REMO2009 '} ];

elseif flag_run == 6 % rcp 85

Expand All @@ -68,7 +68,7 @@
RCM_all(:,1)= [{'CCLM4-8-17'}; {'CRCM5'} ; {'HIRHAM5'} ; {'RACMO22T'}; {'RCA4'} ; {'REMO2009'}];


% detemine number of models
% determine number of models
if flag_run == 1
nm = 1;
else
Expand Down
96 changes: 65 additions & 31 deletions WBmodel/load_WBterms.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -46,68 +46,102 @@

if flag_run == 2 % evaluation

E_wb_ev = csvread('QUANT\E_wb_ev_bc.csv',1)';
P_wb_ev = csvread('QUANT\P_wb_ev_bc.csv',1)';
Qin_wb_ev = csvread('QUANT\Qin_wb_ev_bc.csv',1)';
E_wb_ev = csvread('data\biascorrection\QUANT\E_wb_ev_bc.csv',1)';
P_wb_ev = csvread('data\biascorrection\QUANT\P_wb_ev_bc.csv',1)';
Qin_wb_ev = csvread('data\biascorrection\QUANT\Qin_wb_ev_bc.csv',1)';

elseif flag_run == 3 % historical

E_wb_hist = csvread('QUANT\E_wb_hist_bc.csv',1)';
P_wb_hist = csvread('QUANT\P_wb_hist_bc.csv',1)';
Qin_wb_hist = csvread('QUANT\Qin_wb_hist_bc.csv',1)';
E_wb_hist = csvread('data\biascorrection\QUANT\E_wb_hist_bc.csv',1)';
P_wb_hist = csvread('data\biascorrection\QUANT\P_wb_hist_bc.csv',1)';
Qin_wb_hist = csvread('data\biascorrection\QUANT\Qin_wb_hist_bc.csv',1)';

elseif flag_run == 4 % RCP 2.6

E_wb_rcp26 = csvread('QUANT\E_wb_rcp26_bc.csv',1)';
P_wb_rcp26 = csvread('QUANT\P_wb_rcp26_bc.csv',1)';
Qin_wb_rcp26 = csvread('QUANT\Qin_wb_rcp26_bc.csv',1)';
E_wb_rcp26 = csvread('data\biascorrection\QUANT\E_wb_rcp26_bc.csv',1)';
P_wb_rcp26 = csvread('data\biascorrection\QUANT\P_wb_rcp26_bc.csv',1)';
Qin_wb_rcp26 = csvread('data\biascorrection\QUANT\Qin_wb_rcp26_bc.csv',1)';

elseif flag_run == 5 % RCP 4.5

E_wb_rcp45 = csvread('QUANT\E_wb_rcp45_bc.csv',1)';
P_wb_rcp45 = csvread('QUANT\P_wb_rcp45_bc.csv',1)';
Qin_wb_rcp45 = csvread('QUANT\Qin_wb_rcp45_bc.csv',1)';
E_wb_rcp45 = csvread('data\biascorrection\QUANT\E_wb_rcp45_bc.csv',1)';
P_wb_rcp45 = csvread('data\biascorrection\QUANT\P_wb_rcp45_bc.csv',1)';
Qin_wb_rcp45 = csvread('data\biascorrection\QUANT\Qin_wb_rcp45_bc.csv',1)';

elseif flag_run == 6 % RCP 8.5

E_wb_rcp85 = csvread('QUANT\E_wb_rcp85_bc.csv',1)';
P_wb_rcp85 = csvread('QUANT\P_wb_rcp85_bc.csv',1)';
Qin_wb_rcp85 = csvread('QUANT\Qin_wb_rcp85_bc.csv',1)';
E_wb_rcp85 = csvread('data\biascorrection\QUANT\E_wb_rcp85_bc.csv',1)';
P_wb_rcp85 = csvread('data\biascorrection\QUANT\P_wb_rcp85_bc.csv',1)';
Qin_wb_rcp85 = csvread('data\biascorrection\QUANT\Qin_wb_rcp85_bc.csv',1)';

end

E_wb_hist = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\E_wb_hist_bc.csv',1)';
P_wb_hist = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\P_wb_hist_bc.csv',1)';
Qin_wb_hist = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\Qin_wb_hist_bc.csv',1)';

E_wb_rcp26 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\E_wb_rcp26_bc.csv',1)';
P_wb_rcp26 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\P_wb_rcp26_bc.csv',1)';
Qin_wb_rcp26= csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\Qin_wb_rcp26_bc.csv',1)';

E_wb_rcp45 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\E_wb_rcp45_bc.csv',1)';
P_wb_rcp45 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\P_wb_rcp45_bc.csv',1)';
Qin_wb_rcp45 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\Qin_wb_rcp45_bc.csv',1)';

E_wb_rcp85 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\E_wb_rcp85_bc.csv',1)';
P_wb_rcp85 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\P_wb_rcp85_bc.csv',1)';
Qin_wb_rcp85 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\QUANT\Qin_wb_rcp85_bc.csv',1)';


elseif flag_bc == 2 % LINEAR PARAMETRIC TRANSFORMATION

if flag_run == 2 % evaluation

E_wb_ev = csvread('PFT_lin\E_wb_ev_bc.csv',1)';
P_wb_ev = csvread('PFT_lin\P_wb_ev_bc.csv',1)';
Qin_wb_ev = csvread('PFT_lin\Qin_wb_ev_bc.csv',1)';
E_wb_ev = csvread('data\biascorrection\PFT_lin\E_wb_ev_bc.csv',1)';
P_wb_ev = csvread('data\biascorrection\PFT_lin\P_wb_ev_bc.csv',1)';
Qin_wb_ev = csvread('data\biascorrection\PFT_lin\Qin_wb_ev_bc.csv',1)';

elseif flag_run == 3 % historical

E_wb_hist = csvread('PFT_lin\E_wb_hist_bc.csv',1)';
P_wb_hist = csvread('PFT_lin\P_wb_hist_bc.csv',1)';
Qin_wb_hist = csvread('PFT_lin\Qin_wb_hist_bc.csv',1)';
E_wb_hist = csvread('data\biascorrection\PFT_lin\E_wb_hist_bc.csv',1)';
P_wb_hist = csvread('data\biascorrection\PFT_lin\P_wb_hist_bc.csv',1)';
Qin_wb_hist = csvread('data\biascorrection\PFT_lin\Qin_wb_hist_bc.csv',1)';

elseif flag_run == 4 % RCP 2.6

E_wb_rcp26 = csvread('PFT_lin\E_wb_rcp26_bc.csv',1)';
P_wb_rcp26 = csvread('PFT_lin\P_wb_rcp26_bc.csv',1)';
Qin_wb_rcp26 = csvread('PFT_lin\Qin_wb_rcp26_bc.csv',1)';
E_wb_rcp26 = csvread('data\biascorrection\PFT_lin\E_wb_rcp26_bc.csv',1)';
P_wb_rcp26 = csvread('data\biascorrection\PFT_lin\P_wb_rcp26_bc.csv',1)';
Qin_wb_rcp26 = csvread('data\biascorrection\PFT_lin\Qin_wb_rcp26_bc.csv',1)';

elseif flag_run == 5 % RCP 4.5

E_wb_rcp45 = csvread('PFT_lin\E_wb_rcp45_bc.csv',1)';
P_wb_rcp45 = csvread('PFT_lin\P_wb_rcp45_bc.csv',1)';
Qin_wb_rcp45 = csvread('PFT_lin\Qin_wb_rcp45_bc.csv',1)';
E_wb_rcp45 = csvread('data\biascorrection\PFT_lin\E_wb_rcp45_bc.csv',1)';
P_wb_rcp45 = csvread('data\biascorrection\PFT_lin\P_wb_rcp45_bc.csv',1)';
Qin_wb_rcp45 = csvread('data\biascorrection\PFT_lin\Qin_wb_rcp45_bc.csv',1)';

elseif flag_run == 6 % RCP 8.5

E_wb_rcp85 = csvread('PFT_lin\E_wb_rcp85_bc.csv',1)';
P_wb_rcp85 = csvread('PFT_lin\P_wb_rcp85_bc.csv',1)';
Qin_wb_rcp85 = csvread('PFT_lin\Qin_wb_rcp85_bc.csv',1)';
E_wb_rcp85 = csvread('data\biascorrection\PFT_lin\E_wb_rcp85_bc.csv',1)';
P_wb_rcp85 = csvread('data\biascorrection\PFT_lin\P_wb_rcp85_bc.csv',1)';
Qin_wb_rcp85 = csvread('data\biascorrection\PFT_lin\Qin_wb_rcp85_bc.csv',1)';

end
end


E_wb_hist = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\E_wb_hist_bc.csv',1)';
P_wb_hist = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\P_wb_hist_bc.csv',1)';
Qin_wb_hist = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\Qin_wb_hist_bc.csv',1)';

E_wb_rcp26 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\E_wb_rcp26_bc.csv',1)';
P_wb_rcp26 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\P_wb_rcp26_bc.csv',1)';
Qin_wb_rcp26 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\Qin_wb_rcp26_bc.csv',1)';

E_wb_rcp45i = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\E_wb_rcp45_bc.csv',1)';
P_wb_rcp45i = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\P_wb_rcp45_bc.csv',1)';
Qin_wb_rcp45i = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\Qin_wb_rcp45_bc.csv',1)';

E_wb_rcp85 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\E_wb_rcp85_bc.csv',1)';
P_wb_rcp85 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\P_wb_rcp85_bc.csv',1)';
Qin_wb_rcp85 = csvread('C:\Users\ivand\Documents\MATLAB\Thesis\Rscripts\data\PFT_lin\Qin_wb_rcp85_bc.csv',1)';

end
9 changes: 5 additions & 4 deletions WBmodel/load_data.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -32,15 +32,16 @@

% 3. Load longitude and latitude of CCLM grid
% ------------------------------------------------------------------------
% load lon_P.mat
% load lat_P.mat
load lon_grid.mat
load lat_grid.mat

nc = 0;
[lat_CCLM, lon_CCLM, depth_CCLM] = mf_load('lffd1996010100c.nc', 'DEPTH_LK', nc);



% 4. Load outflow
% ------------------------------------------------------------------------
load outflow.mat

% 5. Load DEM
% ------------------------------------------------------------------------
%load_DEM
Loading

0 comments on commit 3853edc

Please sign in to comment.