clear all
clc

% Nominal swash zone
x0=(-100:1:0)';
h0=-0.05*x0;

% Nominal surf zone 
x1=(1:1:1000)';
h1=-0.15*(x1.^(0.5));

% Nominal shoreface
x2=(1001:1:10000)';
h2=-0.01888*(x2.^(0.8));


h=[h0;h1;h2];
x=[x0;x1;x2];
exp1_profile=[x h];

plot(x,h)

save exp1_profile.dat exp1_profile -ascii