function [saveasfilename]=v25_privstorsubsidies_specs(filename4)
%Different private storage subsidies

%General modelling specs (how to solve and simulate):
colorset='brgkcmy';usesolver='path';method='iter';
par.simulmethod='solve';
shock_nodes=7;t_convergence=5;
saveasfilename='v25_privstorsubsidies';
sim1=1000;sim2=200;
num=1;
% ESIM=0;
load('results/v25_speedtest_useESIM_0.06.mat','ESIM'); %THESE SHOCKS ARE FOR sigma1=sigma1=[0.06 0; 0 0.06]; For other sigma, new file needs to be generated.
ESIM=ESIM(1:sim1,:,1:sim2);
usecommonshocks=1;

%Specific model parameters:
%par.ka=0.06; par.kb=0.06;
par.ra=0.03; par.rb=0.03;
par.elastDa=-0.2;par.elastDb=-0.2;
par.elastSa=0.2;par.elastSb=0.2;
par.thetaab=0.05;par.thetaba=0.05;
par.gammaa=1;par.gammab=1;
par.phiamax=0;par.phibmax=0;
par.Xabmax=1E5;par.Xbamax=1E5;
par.Pamax=9.4;par.Pbmax=9.4;
par.sca=0.0;par.scb=0.0;

%Model version parameters:
modelstructure_change=1;saveall=1;
Mu=[1 1];filename(1).x='v25_standard.yaml';
sigma1=[0.06 0; 0 0.06]; %Shocks 1
sim_start=[1 1];
approx_order=[50 50];loop=300;%50 50 300
state_min=[0.6 0.6];state_max=[1.7 1.7];
xvarnum=12;svarnum=2;evarnum=4;varnum=xvarnum+svarnum;

%Plot parameters:
varrow(1:num)=9; %dim of varrow must equal num   %Response/State variable number to plot (index to choose row)
linesfor=[1 1 1 1 1];%Which vars to plot:mean,mean-sd,mean+sd,min,max (0:no; 1:yes)
line(1).type='-';line(2).type='--';line(3).type=':'; %Define line type to use for plotting
x_label='SD of shocks';y_label='Price';

%Parameters to vary:
xstates=20;      %Number of xstates (x-axis)
% xaxisp=[0.06 0.055 0.05 0.045 0.04 0.035 0.03 0.025 0.02 0.015 0.01 0.005 0];%How to plot states (y-axis values)
synchron=1; %If the same parameter is changed in both countries, then set synchron to 1 (i.e. k_a and k_b are changed in the same way for different scenarios), otherwise set synchron to zero
for i=1:xstates
    par.ka(i,1)=(i-1)*0.005;
end;clear i;
par.kb=par.ka;
xaxisp=par.ka';

% par.ka(1,1)=0.06; par.kb(1,1)=0.06;
% par.ka(2,1)=0.055; par.kb(2,1)=0.055;
% par.ka(3,1)=0.05; par.kb(3,1)=0.05;
% par.ka(4,1)=0.045; par.kb(4,1)=0.045;
% par.ka(5,1)=0.04; par.kb(5,1)=0.04;
% par.ka(6,1)=0.035; par.kb(6,1)=0.035;
% par.ka(7,1)=0.03; par.kb(7,1)=0.03;
% par.ka(8,1)=0.025; par.kb(8,1)=0.025;
% par.ka(9,1)=0.02; par.kb(9,1)=0.02;
% par.ka(10,1)=0.015; par.kb(10,1)=0.015;
% par.ka(11,1)=0.01; par.kb(11,1)=0.01;
% par.ka(12,1)=0.005; par.kb(12,1)=0.005;
% par.ka(13,1)=0.00; par.kb(13,1)=0.00;

% par.ka(8,1)=-0.01; par.kb(8,1)=-0.01;
% par.ka(9,1)=-0.02; par.kb(9,1)=-0.02;
% par.ka(10,1)=-0.03; par.kb(10,1)=-0.03;
% par.ka(11,1)=-0.04; par.kb(11,1)=-0.04;
% par.ka(12,1)=-0.05; par.kb(12,1)=-0.05;
% % par.ka(4,1)=-0.05; par.kb(12,1)=-0.05;

%Alternative model configurations to simulate:
num=num;          %Number of model configurations
    %Adjust varrow dimension so that it is equal to num!
% par.thetaab=1E5;par.thetaba=1E5;
% par.sca=0.0;par.scb=0.0;
save([filename4 '1']);%No trade, no reserve
% par.thetaab=0.05;par.thetaba=0.05;
% save([filename4 '2']);%No reserve
% par.thetaab=1E5;par.thetaba=1E5;
% par.sca=0.1;par.scb=0.1;
% save([filename4 '3']);%No trade
% par.thetaab=0.05;par.thetaba=0.05;
% par.sca=0.1;par.scb=0.0;
% save([filename4 '4']);%Trade but only one country (focus) has reserve
% par.sca=0.0;par.scb=0.1;
% save([filename4 '5']);%Trade but only one country has reserve (the one out of focus)
% par.sca=0.1;par.scb=0.1;
% par.phiamax=1E5;par.phibmax=1E5;
% save([filename4 '6']);%Trade only if no country uses reserve
% par.phiamax=0;par.phibmax=0;
% par.Xabmax=par.sca;par.Xbamax=par.scb;
% save([filename4 '7']);%Trade only up to capacity of reserve
% par.Xabmax=1E5;par.Xbamax=1E5;
% save([filename4 '8']);%Common reserve/two ind. reserves
% par.sca=0.0;par.scb=0.0;
% par.ka=0.04; par.kb=0.04;
% save([filename4 '9']);%Storage subsidiy 0.04 instead of 0.06
% par.ka=0.02; par.kb=0.02;
% save([filename4 '10']);%Storage subsidiy 0.02 instead of 0.06
% par.sca=0.1;par.scb=0.1;
% par.ka=0.04; par.kb=0.04;
% save([filename4 '11']);%Common reserve/two ind. reserves + storage subsidy of 0.04 instead 0.06