add to delete

eros_hochrhein.m

+%__________________________________________________________________________
+%                               PREPARE GRIDS
+%--------------------------------------------------------------------------
+%
+%+++++++++++++++++ Example data (1% inclined, flat surface) +++++++++++++++
+addpath('.\mfiles');
+% ALT (elevation model) 
+dem=GRIDobj('./Topo/channel_dx8_S0.01.tif');
+
+% RAIN (sources (>0) and sinks (-1))
+rain = dem*0;
+rain.Z(95:117)=-1;
+rain.Z([35712;35713;35714;35715;35716;35717;35718;35719;35720;35721;35722;35723;35724])=0.441;
+rain.Z([15476;15622;15768;15914;16060;16206;16352;16498;16644;16790;16936;17082;17228])=0.559;
+GRIDobj2grd(rain,['./Topo/',dem.name,'.rain']);
+
+%SED (sediment thickness in meters)
+sed = dem*0;
+sed = sed+10;
+GRIDobj2grd(sed,['./Topo/',dem.name,'.sed']);
+
+LEM.dem = dem;
+LEM.rain = rain;
+LEM.sed = sed;
+
+GRIDobj2grd(dem,['./Topo/',dem.name,'.alt']);
+GRIDobj2grd(rain,['./Topo/',dem.name,'.rain']);
+GRIDobj2grd(sed,['./Topo/',dem.name,'.sed']);
+
+%+++++++++++++++++++++++++++ Real topography ++++++++++++++++++++++++++++++
+% % ALT (elevation model)
+% dem = GRIDobj(['./Topo/','SantaMariaBasin_outlet_utm20s.tif']);
+% GRIDobj2grd(dem,['./Topo/',dem.name,'.alt']);
+% 
+% % RAIN (sources and sinks)
+% rain = dem*0;
+% rain.Z(38:43,1)=1;
+% rain.Z(1,44:60)=-1;
+% GRIDobj2grd(rain,['./Topo/',dem.name,'.rain']);
+% 
+% % SED (sediment thickness)
+% sed = dem*0;
+% sed = sed+10;
+% GRIDobj2grd(sed,['./Topo/',dem.name,'.sed']);
+% 
+% LEM.dem = dem;
+% LEM.rain = rain;
+% LEM.sed = sed;
+%__________________________________________________________________________
+%%                           DEFINE INPUT PARAMETERS
+%--------------------------------------------------------------------------
+
+LEM.experiment = 'single_thread';                % Project name
+LEM.ErosPath = 'C:\\Users\\mey\\SynologyDrive\\erosmatlabinterface';    % Path to Eros.exe
+LEM.outfolder = 'playingaround';              % folder to store results in
+
+LEM.ErosPath = 'C:\\Projects\\EROS\\Hochrhein';     % Path to bin-folder, where .exe resides
+LEM.outfolder = 'Results';                          % folder to store results in
+
+LEM.start = 0;                       % start time (negative time)
+LEM.stop = 2000;                      % [s] length of model run
+LEM.draw = 10; 						% output interval
+LEM.init = 2000;                    
+
+LEM.step = 5;
+LEM.stepmin = 0.3;
+LEM.stepmax = 0.3;
+
+
+LEM.TU = 100;                           % unknown parameter
+LEM.floodos = 'stationary:pow';
+
+%--------------------------------------------------------------------------
+% EROSION/DEPOSITION
+%--------------------------------------------------------------------------
+LEM.erosion_model = 'shear_stress';     % (stream_power, shear_stress, shear_mpm)
+LEM.deposition_model = 'constant';      % need to know whether there are other options!
+
+% ALLUVIAL
+LEM.fluvial_stress_exponent = 1.5;      % exponent in sediment flux eq. (MPM): qs = E(tau-tau_c)^a
+LEM.fluvial_erodability = 0.4;       % [kg-1.5 m-3.5 s-2] E in MPM equation
+LEM.fluvial_sediment_threshold = 4.0;   % [Pa] critical shear stress (tau_c) in MPM equation
+LEM.deposition_length = 2;              % [m] xi in vertical erosion term: edot = qs/xi
+
+% Lateral erosion/deposition
+LEM.fluvial_lateral_erosion_coefficient = 0.05;             % dimensionless coefficient (Eq. 17 in Davy, Croissant, Lague (2017))
+LEM.fluvial_lateral_deposition_coefficient = 0.5;
+
+LEM.outbend_erosion_coefficient                   = 1.0;
+LEM.inbend_erosion_coefficient                    = 1.0;
+
+% BEDROCK
+LEM.fluvial_basement_erodability = 0.2;
+LEM.fluvial_basement_threshold = 4;
+
+% Lateral erosion/deposition
+LEM.lateral_erosion = 0.05;             % dimensionless coefficient (Eq. 17 in Davy, Croissant, Lague (2017))
+LEM.lateral_deposition = 0.5;
+
+LEM.outbend_erosion_coefficient                   = 1.000000;
+LEM.inbend_erosion_coefficient                    = 1.00000;
+%--------------------------------------------------------------------------
+% FLOW MODEL
+%--------------------------------------------------------------------------
+LEM.diffusion_coefficient = 0.5;
+LEM.poisson_coefficient = 4;
+LEM.flood_model = 1;
+LEM.flow_model = 'manning';
+LEM.minimum_flow_coefficient = 0.0;
+LEM.friction_coefficient = 0.025;       % 
+LEM.water_depth_limit = -1;
+LEM.flow_depth_constant = 1.0;
+LEM.flow_only = 0;
+LEM.fictious_area = 1;
+
+
+%--------------------------------------------------------------------------
+% OUTPUTS TO WRITE
+%--------------------------------------------------------------------------
+LEM.stress = 1;
+LEM.water = 1;
+LEM.discharge = 1;
+LEM.downward = 1;
+LEM.slope = 1;
+LEM.qs = 1;
+LEM.capacity = 1;
+LEM.sediment = 1;
+LEM.rainfall = 0;
+
+LEM.str_write = '';
+LEM.str_nowrite = '';
+
+writeErosInputs(LEM);
+
+%% run model
+system([LEM.experiment,'.bat'])
+