%script to run diffuse
%This script uplifts a hill in the middle of the profile (the uplift profile is edfined below in the Harb vector which is total vertical uplift for the model run for each node)


%You should change these below:
dx = 1;   %usually m

Hinitial = [1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1];

k = 1   %usually m^2/kyr
finaltime = 50;   %usually kyr

pt = 100 ; %every how many time steps you want an output profile drawn

%make sure that Harb and Hinitial have the same length (number of elements or there will be a problem.  Or just comment this line out with a precending %:
Harb = [0 0 0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0];



%don't mess with this:
diffuse(dx, Hinitial, k, finaltime, pt, Harb)