function [koefAtonacijsko]=SP1996hor(epiFLT,Period,ecSoilCode)
%%%%sabeta pugliese 1996 koeficienti

%%%%%% Soil type%%%%%%%
%%%%%% Soil type%%%%%%%
if ecSoilCode==1
    %Rock
    SB=0;SC=0;SD=0;
elseif ecSoilCode==2
    %Stiff soil 
    SB=1;SC=0;SD=0;
elseif ecSoilCode==3
    %Soft soil
    SB=0;SC=1;SD=0;
elseif ecSoilCode==4
    SB=0;SC=0;SD=1;
else
    SB=0;SC=0;SD=1;
end


if epiFLT==1
    %epicentral distance
Tvec=[0.00 0.040 0.067 0.100 0.150 0.200 0.300 0.400 0.500 0.752 1.000 1.493 2.000 3.030 4.000];
b1vec=[-4.248 -3.712 -3.060 -2.791 -2.641 -2.759 -3.616 -4.780 -5.727 -7.067 -7.997 -9.243 -10.118 -11.339 -12.192];
b2vec=[0.836 0.760 0.700 0.700 0.714 0.744 0.868 1.025 1.151 1.312 1.409 1.520 1.582 1.646 1.669];
b3vec=[-1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000];
b4vec=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
h0vec=[5.0 4.7 6.3 6.2 5.9 5.7 5.4 5.2 5.0 4.7 4.4 4.0 3.6 3.0 2.6];
r0vec=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
bSBvec=[0.000 0.000 0.000 0.000 0.000 0.000 0.046 0.180 0.286 0.437 0.479 0.403 0.345 0.249 0.230];
bSCvec=[0.449 0.371 0.371 0.371 0.371 0.371 0.426 0.511 0.530 0.276 0.115 0.023 0.000 0.000 0.000];
bSDvec=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
sigmavec=[0.437 0.449 0.461 0.479 0.507 0.539 0.599 0.645 0.668 0.698 0.709 0.725 0.735 0.735 0.735];

b1=interp1(Tvec,b1vec,Period);
b2=interp1(Tvec,b2vec,Period);
b3=interp1(Tvec,b3vec,Period);
b4=interp1(Tvec,b4vec,Period);
h0=interp1(Tvec,h0vec,Period);
r0=interp1(Tvec,r0vec,Period);
bSB=interp1(Tvec,bSBvec,Period);
bSC=interp1(Tvec,bSCvec,Period);
bSD=interp1(Tvec,bSDvec,Period);
sigma=interp1(Tvec,sigmavec,Period);

else
    %fault distance
Tvec=[0.00 0.040 0.067 0.100 0.150 0.200 0.300 0.400 0.500 0.752 1.000 1.493 2.000 3.030 4.000];
b1vec=[-3.597 -2.994 -2.422 -2.079 -2.001 -2.174 -2.926 -4.311 -5.278 -6.722 -7.629 -8.927 -9.887 -11.155 -11.962];
b2vec=[0.705 0.629 0.587 0.564 0.587 0.629 0.725 0.921 1.048 1.220 1.312 1.428 1.497 1.554 1.577];
b3vec=[-1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000 -1.000];
b4vec=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
h0vec=[5.800 5.800 7.200 7.300 6.700 6.300 5.700 5.300 5.000 4.400 4.000 3.500 3.100 2.500 2.100];
r0vec=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
bSBvec=[0.000 0.000 0.000 0.000 0.000 0.000 0.150 0.265 0.359 0.534 0.557 0.484 0.424 0.348 0.299];
bSCvec=[0.389 0.299 0.299 0.299 0.299 0.299 0.380 0.484 0.507 0.276 0.115 0.023 0.000 0.000 0.000];
bSDvec=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
sigmavec=[0.398 0.419 0.419 0.437 0.461 0.495 0.564 0.608 0.624 0.652 0.668 0.679 0.684 0.684 0.684];

b1=interp1(Tvec,b1vec,Period);
b2=interp1(Tvec,b2vec,Period);
b3=interp1(Tvec,b3vec,Period);
b4=interp1(Tvec,b4vec,Period);
h0=interp1(Tvec,h0vec,Period);
r0=interp1(Tvec,r0vec,Period);
bSB=interp1(Tvec,bSBvec,Period);
bSC=interp1(Tvec,bSCvec,Period);
bSD=interp1(Tvec,bSDvec,Period);
sigma=interp1(Tvec,sigmavec,Period);
    
end


koefAtonacijsko.b1=b1;
koefAtonacijsko.b2=b2;
koefAtonacijsko.b3=b3;
koefAtonacijsko.b4=b4;
koefAtonacijsko.h0=h0;
koefAtonacijsko.r0=r0;
koefAtonacijsko.bSB=bSB;
koefAtonacijsko.bSC=bSC;
koefAtonacijsko.bSD=bSD;
koefAtonacijsko.sigma=sigma;

koefAtonacijsko.SB=SB;
koefAtonacijsko.SC=SC;
koefAtonacijsko.SD=SD;









