function [saveasfilename]=v25_defaultvalues_specs(filename4)

%General modelling specs (how to solve and simulate):
colorset='brgkcmy';usesolver='path';method='iter';%Solver & Method for solving & colorset for illustrating results.
par.simulmethod='interpolation';
shock_nodes=7;t_convergence=5;
saveasfilename='v25_defaultvalues';%Filename for the results
sim1=1000;sim2=200;%Simulation dimensions: sim1: How often to start from the steady state/how many series to simulate? sim2: How long to simulate each of these series which starts from the steady state?
num=13;%Number of scenarios
% 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 sigmas, new files need to be generated.
ESIM=ESIM(1:sim1,:,1:sim2);%Define shock variable ESIM (template is bigger in each dimension)
usecommonshocks=1;%Set to 1 if the same shocks are to be used by all different model specifications. Otherwise set to 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.4;par.Pbmax=1.4;
% par.sca=0.035;par.scb=0.035;
sca=0.035;scb=0.035;

%Model version parameters:
modelstructure_change=1;%Not used at the moment.
saveall=1;%What data should be saved? 1:All, 0: Only summary statistics are saved.
Mu=[1 1];filename(1).x='v25_standard.yaml';
sim_start=[1 1];%Starting point for simulations
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;%Define the number of state, response, expectation, and state+response variables

%Plot parameters:
%THIS WAY OF PLOTTING IS NOT USED ANY MORE AND THEREFORE THESE SETTINGS DO
%NOT HAVE ANY EFFECT
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=1;      %Number of xstates (x-axis)
xaxisp=[0.06];%How to plot states (x-axis values)
sigma1=[0.06 0; 0 0.06]; %Shocks 1
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


%Alternative model configurations to simulate:
%One spec file is generated for each scenario. That way different
%scenarios can be calculated.
num=num;          %Number of model configurations; Update above, this is just a reminder.
    %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']);%Scenario description: %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=sca;par.scb=scb;
save([filename4 '3']);%No trade
par.thetaab=0.05;par.thetaba=0.05;
par.sca=sca;par.scb=0.0;
save([filename4 '4']);%Trade but only one country (focus) has reserve
par.sca=0.0;par.scb=scb;
save([filename4 '5']);%Trade but only one country has reserve (the one out of focus)
par.sca=sca;par.scb=scb;
par.phiamax=2;par.phibmax=2;
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.02;par.scb=0.02;
par.Pamax=1.3;par.Pbmax=1.3;
save([filename4 '9']);%Tiny common reserve/two ind. reserves with lower trigger price
par.sca=0.0;par.scb=0.0;
par.Pamax=1.4;par.Pbmax=1.4;
par.ka=0.03; par.kb=0.03;
save([filename4 '10']);%Storage subsidiy 0.03 instead of 0.06
par.ka=0.00; par.kb=0.00;
save([filename4 '11']);%Storage subsidiy 0.00 instead of 0.06
par.sca=sca;par.scb=scb;
par.ka=0.03; par.kb=0.03;
save([filename4 '12']);%Common reserve/two ind. reserves + storage subsidy of 0.03 instead 0.06
par.sca=0.02;par.scb=0.02;
par.Pamax=1.3;par.Pbmax=1.3;
par.ka=0.04; par.kb=0.04;
save([filename4 '13']);%Common reserve/two ind. reserves + storage subsidy of 0.04 instead 0.06