function [saveasfilename]=v25_asymshocks_specs2(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_asymshocks2';
sim1=1000;sim2=200;
num=3;
% 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=0;
usecommonshocks=0;

%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=1.3;par.Pbmax=1.3;
par.sca=0.02;par.scb=0.02;

%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=[40 40];loop=250;%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=13;      %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=0; %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
    sigma1(:,:,i)=[0.06 0; 0 0.01*i];
end;clear i;
%par.kb=par.ka;
xaxisp=squeeze(sigma1(2,2,:));


%Alternative model configurations to simulate:
num=num;          %Number of model configurations
    %Adjust varrow dimension so that it is equal to num!

save([filename4 '1']);
par.thetaab=10;par.thetaba=10;
save([filename4 '2']);
par.thetaab=0.05;par.thetaba=0.05;
par.phiamax=10;par.phibmax=10;
save([filename4 '3']);