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

%General modelling specs (how to solve and simulate):
colorset='brgkcmy';usesolver='path';method='iter';
par.simulmethod='interpolation';
shock_nodes=7;t_convergence=5;
saveasfilename='v25_reservevariation';
sim1=600;sim2=200;
num=20;
% 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.005;par.scb=0.005;

%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=21;      %Number of xstates (x-axis)
% yaxisp=[1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65];%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

par.Pamax(1:xstates,1)=1.075+(0.025:0.025:0.025*xstates);
par.Pbmax=par.Pamax;
xaxisp=par.Pamax';


%Alternative model configurations to simulate:
num=num;          %Number of model configurations
    %Adjust varrow dimension so that it is equal to num!
for i99=1:num%Set to 18 and adjust factor 2 for higher precision
    par.sca=i99*0.005;par.scb=par.sca;
    sim_start=[1+par.sca 1+par.scb];%[1.0004+0.014+par.sca 1.0004+0.014+par.scb];
    save([filename4 num2str(i99)]);
end;clear i99;
    
% par.sca=0.005;par.scb=0.005;
% save([filename4 '1']);
% par.sca=0.01;par.scb=0.01;
% save([filename4 '2']);
% par.sca=0.015;par.scb=0.015;
% save([filename4 '3']);
% par.sca=0.02;par.scb=0.02;
% save([filename4 '4']);
% par.sca=0.025;par.scb=0.025;
% save([filename4 '5']);
% par.sca=0.03;par.scb=0.03;
% save([filename4 '6']);
% par.sca=0.035;par.scb=0.035;
% save([filename4 '7']);
% par.sca=0.04;par.scb=0.04;
% save([filename4 '8']);
% par.sca=0.045;par.scb=0.045;
% save([filename4 '9']);
% par.sca=0.05;par.scb=0.05;
% save([filename4 '10']);
% par.sca=0.055;par.scb=0.055;
% save([filename4 '11']);
% par.sca=0.06;par.scb=0.06;
% save([filename4 '12']);