function [AddPath,ForwardModel,ValueParameters,rho_radial_obs2,ni]=InputParametersForward


% Subroutine to impose the input parameters of the software:
%
% The path of the folder in AddPath in which to put the output results
%
% PARAMETERS TO SET THE FORWARD MODEL
%
% The parameter "SourceType" to choose between the different geometric shapes: 
% sphere (1), rectangular prism (2), prolate or oblate ellipsoid (3), triaxial ellipsoid (4), cylinder (5)
%
% The parameter "SourceRotation" to choose (for the ellipsoids and the rectangular prism) 
% whether to fix the position of the source (SourceRotation=0) or or rotate the source in space (SourceRotation=1)
%
% The values of the parameters that define the source 







%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% PARAMETERS TO SET THE FORWARD MODEL %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

ni=0.25; % Poisson's ratio


%% Load the folder path
AddPath='';


%% Choose the deformation source (forward model) 
% SourceType=1 --> sphere
% SourceType=2 --> rectangular prism
% SourceType=3 --> prolate or oblate ellipsoid
% SourceType=4 --> triaxial ellipsoid
% SourceType=5 --> cylinder
SourceType=3; 


%% Choose the orientation of the source
% SourceRotation=1 --> the source rotate
% SourceRotation=0 --> the source does not rotate
if SourceType==2||SourceType==3||SourceType==4
    SourceRotation=1;
    ForwardModel=[SourceType, SourceRotation];
elseif SourceType==1||SourceType==5
    ForwardModel=[SourceType];
end



%% Insert the values of the source parameters

switch SourceType

    case 1  % sphere               
        Xc=496494; % [m]
        Yc=4250781; % [m]
        Zc=-688; % [m]
        DeltaV=101540; % [m]
        ValueParameters=[Xc,Yc,Zc,DeltaV];  

    case 2  % rectangular prism
        Xc=0; % [m]
        Yc=0; % [m]
        depth=-500;  % upper limit of the prism % [m] 
        l_a=1000; % prism width % [m]
        l_b=400; % prism depth % [m]
        l_c=200; % prism height % [m]
        eps_0=0.0028; % stress-free strain 
        ValueParameters=[Xc,Yc,depth,l_a,l_b,l_c,eps_0];
        if SourceRotation==1 % yes rotation
            alpha=140; % rotation angle around the z axis
            delta=40; % rotation angle around the y axis
            gamma=50; % rotation angle around the x axis
            ValueParameters=[ValueParameters,alpha,delta,gamma];
        end

    case 3 % prolate or oblate ellipsoid
        Xc=0; % [m]
        Yc=0; % [m]
        a=500; % length of the semi-axis along the x direction % [m]
        b=300; % length of the semi-axis along the y direction % [m]
        depth=-123; % upper limit % [m]
        eps_0=0.003; % stress-free strain 
        ValueParameters=[Xc,Yc,depth,a,b,eps_0];
        if SourceRotation==1 % yes rotation
            alpha=125; % rotation angle around the z axis
            delta=80; % rotation angle around the y axis
            gamma=5; % rotation angle around the x axis
            ValueParameters=[ValueParameters,alpha,delta,gamma];              
        end

    case 4 % triaxial ellipsoid
        Xc=0; % [m]
        Yc=0; % [m]
        a=500; % length of the semi-axis along the x direction % [m]
        b=100; % length of the semi-axis along the y direction % [m]
        c=500; % length of the semi-axis along the z direction % [m]
        depth=-200; % upper limit % [m]
        eps_0=0.003; % stress-free strain 
        ValueParameters=[Xc,Yc,depth,a,b,c,eps_0];
        if SourceRotation==1 % yes rotation
            alpha=12; % rotation angle around the z axis
            delta=12; % rotation angle around the y axis
            gamma=12; % rotation angle around the x axis
            ValueParameters=[ValueParameters,alpha,delta,gamma];               
        end

    case 5 % cylinder
        Xc=0; % [m]
        Yc=0; % [m]
        depth=-365; % upper limit % [m]
        R=105; % radius % [m]
        h=1000; % height % [m]
        eps_0=0.003; % stress-free strain
        ValueParameters=[Xc,Yc,depth,R,h,eps_0];        

end



% observation points over a ground surface profile along the x-axis 
rho_radial_obs2=[0:100:5000]+Xc;