-
Notifications
You must be signed in to change notification settings - Fork 0
/
SPLM_Example4.m
36 lines (33 loc) · 1.25 KB
/
SPLM_Example4.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
%%--------------------------Example 4 :SPLM-------------------------------%
% Author: Benjamin Campforts Katholieke Universiteit Leuven
%
%ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo%
%References
%
% Campforts, B., and G. Govers (2015),Keeping the edge: A numerical method
% that avoids knickpoint smearing when solving the stream power law, J.
% Geophys. Res. Earth Surf., 120, doi:10.1002/2014JF003376.
%
% Campforts, B., Schwanghart W, and G. Govers (2015),TTLEM 1.0: A numerical
% package for accurate simulation of transient landscape evolution in
% MATLAB. GMD
%ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo%
%
%-------------------------------------------------------------------------%
% clearvars;
clc;close all force;
numM=4;
K=5e-6; m=.42; n=2;kappa=0;
parameters=[K m n kappa];
dx=100; x=1:dx:15E3;hackFactor=2;DA=x.^hackFactor;
spatial={dx x DA} ;
t_end=1E6;
timing=[t_end nan];
UScen=1; maxElevation=1000;
upliftData={UScen maxElevation};
iniSurf = x.*0;baseLevelDescent=0;
oriBed={iniSurf baseLevelDescent};
visibleFlag=0;
plotOut=2;
[z, dt]=SPLM(numM,parameters,spatial,timing,upliftData,oriBed,visibleFlag,plotOut);