Skip to content
Snippets Groups Projects
Commit 35a1ad48 authored by Dr. rer. nat. Jürgen Mey's avatar Dr. rer. nat. Jürgen Mey :eye:
Browse files

added more templates

parent 295e1c6c
Branches baseline
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
mfiles/eros_baseline.m 0 → 100644
+ 117
− 0
View file @ 35a1ad48
%--------------------------------------------------------------------------
% PREPARE GRIDS
%--------------------------------------------------------------------------
%
% addpath('.\mfiles')
% ALT (elevation model)
clear
dem=GRIDobj('./Topo/HochRhein_LAKE_1000m.tif');
% RAIN (sources (>0) and sinks (-1))
rain = GRIDobj('.\Topo\HochRhein_MAP_1000m.tif');
% WATER
water = GRIDobj('.\Topo\HochRhein_WATER+LAKE_1000m.tif');
% UPLIFT
% uplift = GRIDobj('.\Topo\uplift_m_per_s_baselevel_Basel.tif');
uplift = GRIDobj('.\Topo\uplift_m_per_s_nagra_baselevel_Basel.tif');
uplift = resample(uplift,dem);
% SED (sediment thickness in meters)
sed = GRIDobj('.\Topo\Hochrhein_SEDLAKE_1000m.tif');
% sed = dem*0;
% sed.Z(~isnan(sed.Z))=10;
LEM.dem = dem;
LEM.rain = rain;
LEM.sed = sed;
LEM.uplift = uplift;
LEM.water = water;
GRIDobj2grd(dem,['./Topo/',dem.name,'.alt']);
GRIDobj2grd(rain,['./Topo/',dem.name,'.rain']);
GRIDobj2grd(sed,['./Topo/',dem.name,'.sed']);
GRIDobj2grd(uplift,['./Topo/',dem.name,'.uplift']);
GRIDobj2grd(water,['./Topo/',dem.name,'.water']);
%--------------------------------------------------------------------------
%% DEFINE INPUT PARAMETERS
%--------------------------------------------------------------------------
LEM.experiment = 'baseline_test'; % Project name
LEM.ErosPath = 'D:\\USER\\mey'; % Path to .exe
LEM.outfolder = 'baseline_test\\sediments'; % folder to store results in
% LEM.inflow = 1060; % [m3s-1]water inflow at source cells
LEM.rainfall = 7.6e-8*10; % Sets the precipitation rate per unit surface when multiplied by the rainfall map
LEM.initial_sediment_stock = 0; % [%] volumetric sediment inflow at source cells)
LEM.inertia = 0; % refers to inertia term in shallow water equation
LEM.end = 3e8; LEM.end_option = 'time'; % length of model run
LEM.draw = 3e4; LEM.draw_option = 'time'; % output interval
LEM.step = 1e3; LEM.step_option = 'volume';
LEM.stepmin = 1e2;
LEM.stepmax = 1e4;
LEM.initbegin = 1e+3; % initialization time (-)
LEM.initend = 1e+3;
LEM.initstep = 2;
LEM.TU_coefficient = 1; % sets the proportion of rain pixels that make up 1 TU
LEM.flow_model = 'stationary:pow';
LEM.erosion_multiply = 1000;%15654; % multiplying factor for erosion rates. Equivalent to consider an "erosion time" larger than the hydrodynamic time
LEM.uplift_multiplier = 425516;%15654;
LEM.limiter = 1e-1;
LEM.continue_run = -1; % continue previous simulation from the specified stage (-1 for new run)
%--------------------------------------------------------------------------
% EROSION/DEPOSITION
%--------------------------------------------------------------------------
LEM.erosion_model = 'MPM'; % (stream_power, shear_stress, shear_mpm)
LEM.deposition_model = 'constant'; % need to know whether there are other options!
% ALLUVIAL
LEM.deposition_length = 1000; % [m] xi in vertical erosion term: edot = qs/xi
% LATERAL EROSION/DEPOSITION
LEM.fluvial_lateral_erosion_coefficient = 1e-4; % dimensionless coefficient (Eq. 17 in Davy, Croissant, Lague (2017))
LEM.fluvial_lateral_deposition_coefficient = 0.5;
LEM.lateral_erosion_model = 1;
LEM.lateral_deposition_model = 'constant';
% BEDROCK
LEM.poisson_coefficient = 5;
LEM.diffusion_coefficient = 4;
LEM.sediment_grain = 0.025;
LEM.basement_grain = 0.25;
%--------------------------------------------------------------------------
% FLOW MODEL
%--------------------------------------------------------------------------
LEM.friction_model = 'manning';
LEM.friction_coefficient = 0.025; %
LEM.flow_boundary = 'free';
%--------------------------------------------------------------------------
% OUTPUTS TO WRITE
%--------------------------------------------------------------------------
LEM.stress = 1;
LEM.waters = 1;
LEM.discharge = 1;
LEM.downward = 0;
LEM.slope = 1;
LEM.qs = 1;
LEM.capacity = 1;
LEM.sediment = 1;
LEM.flux =1;
LEM.stock =1;
LEM.str_write = '';
LEM.str_nowrite = '';
writeErosInputs(LEM);
%% run model
system([LEM.experiment,'.bat'])
\ No newline at end of file
mfiles/eros_hillslope_test2.m 0 → 100644
+ 138
− 0
View file @ 35a1ad48
%--------------------------------------------------------------------------
% PREPARE GRIDS
%--------------------------------------------------------------------------
%
% addpath('.\mfiles')
% ALT (elevation model)
clear
dem=GRIDobj('./Topo/Alps_basin_3.tif');
% RAIN (sources (>0) and sinks (-1))
rain = GRIDobj('.\Topo\Alps_basin_3_rain.tif'); % MAP after WorldClim2 (mm/yr)
% rain = resample(rain,dem);
% rain = rain/1000; % convert to m/yr
% rain = rain/3600/24/365.25/dem.cellsize.^2; % convert to m/s
% rain.Z(1:end,1)=-1;
% rain.Z(1,1:end)=-1;
% rain.Z(end,1:end)=-1;
% rain.Z(1:end,end)=-1;
% WATER
water = GRIDobj('.\Topo\Alps_basin_3_water.tif');
% UPLIFT
% uplift = GRIDobj('.\Topo\uplift_m_per_s_baselevel_Basel.tif');
uplift = GRIDobj('.\Topo\uplift_m_per_s_nagra_baselevel_Basel.tif');
uplift = resample(uplift,dem);
% SED (sediment thickness in meters)
sed = GRIDobj('.\Topo\mqu_140715g_utm32n.tif');
sed = resample(sed,dem);
sed.Z(isnan(sed.Z))=0;
LEM.dem = dem;
LEM.rain = rain;
% LEM.sed = sed;
% LEM.uplift = uplift;
LEM.water = water;
GRIDobj2grd(dem,['./Topo/',dem.name,'.alt']);
GRIDobj2grd(rain,['./Topo/',dem.name,'.rain']);
% GRIDobj2grd(sed,['./Topo/',dem.name,'.sed']);
% GRIDobj2grd(uplift,['./Topo/',dem.name,'.uplift']);
GRIDobj2grd(water,['./Topo/',dem.name,'.water']);
%--------------------------------------------------------------------------
%% DEFINE INPUT PARAMETERS
%--------------------------------------------------------------------------
LEM.experiment = 'hillslope_test2'; % Project name
LEM.ErosPath = 'D:\\USER\\mey'; % Path to .exe
LEM.outfolder = 'hillslope_test\\rgqs'; % folder to store results in
% LEM.inflow = 1060; % [m3s-1]water inflow at source cells
LEM.rainfall = 2.8e-5*2; % Sets the precipitation rate per unit surface when multiplied by the rainfall map
LEM.initial_sediment_stock = 0; % % The total "stock" of sediment at the precipiton landing is: input_sediment_concentration*cs_map[i]*Precipiton_volume
LEM.inertia = 0; % refers to inertia term in shallow water equation
LEM.begin = 0; LEM.begin_option = 'time'; % start time
LEM.end = 9e+5; LEM.end_option = 'time'; % length of model run
LEM.draw = 3000; LEM.draw_option = 'time'; % output interval
LEM.step = 100; LEM.step_option = 'volume';
LEM.stepmin = 1e1;
LEM.stepmax = 1e4;
LEM.initbegin = 1e+1; % initialization time (-)
LEM.initend = 1e+1;
LEM.initstep = 2;
LEM.TU_coefficient = '1:dir'; % sets the proportion of rain pixels that make up 1 TU
LEM.flow_model = 'stationary:pow';
LEM.erosion_multiply = 100; % multiplying factor for erosion rates. Equivalent to consider an "erosion time" larger than the hydrodynamic time
LEM.uplift_multiplier = 1/3600/24/365.25/1000*700000;
LEM.limiter = 1e-1;
LEM.continue_run = -1;
%--------------------------------------------------------------------------
% EROSION/DEPOSITION
%--------------------------------------------------------------------------
LEM.erosion_model = 'MPM'; % (stream_power, shear_stress, shear_mpm)
LEM.deposition_model = 'constant'; % need to know whether there are other options!
LEM.eros_version = 'eros7.3.112';
LEM.stress_model = 'rgqs';
% ALLUVIAL
LEM.fluvial_stress_exponent = 1.5; % exponent in sediment flux eq. (MPM): qs = E(tau-tau_c)^a
LEM.fluvial_erodability = 2.6e-8; % [kg-1.5 m-3.5 s-2] E in MPM equation
LEM.fluvial_sediment_threshold = 0.05; % [Pa] critical shear stress (tau_c) in MPM equation
LEM.deposition_length = 30; % [m] xi in vertical erosion term: edot = qs/xi
% LATERAL EROSION/DEPOSITION
LEM.fluvial_lateral_erosion_coefficient = 1e-4; % dimensionless coefficient (Eq. 17 in Davy, Croissant, Lague (2017))
LEM.fluvial_lateral_deposition_coefficient = 0.5;
LEM.lateral_erosion_model = 1;
LEM.lateral_deposition_model = 'constant';
% BEDROCK
LEM.fluvial_basement_erodability = 0.1;
LEM.fluvial_basement_threshold = 0.5;
LEM.outbend_erosion_coefficient = 1.000000;
LEM.inbend_erosion_coefficient = 1.00000;
LEM.poisson_coefficient = 5;
LEM.diffusion_coefficient = 4;
LEM.sediment_grain = 0.0025;
LEM.basement_grain = 0.025;
%--------------------------------------------------------------------------
% FLOW MODEL
%--------------------------------------------------------------------------
LEM.friction_model = 'manning';
LEM.friction_coefficient = 0.025; %
LEM.flow_boundary = '2';
%--------------------------------------------------------------------------
% OUTPUTS TO WRITE
%--------------------------------------------------------------------------
LEM.stress = 1;
LEM.waters = 1;
LEM.discharge = 1;
LEM.downward = 0;
LEM.slope = 1;
LEM.qs = 1;
LEM.capacity = 1;
LEM.sediment = 1;
LEM.flux =1;
LEM.stock =1;
LEM.str_write = '';
LEM.str_nowrite = '';
writeErosInputs(LEM);
%% run model
system([LEM.experiment,'.bat'])
\ No newline at end of file
mfiles/eros_hochrhein.m 0 → 100644
+ 168
− 0
View file @ 35a1ad48
%--------------------------------------------------------------------------
% PREPARE GRIDS
%--------------------------------------------------------------------------
%
% addpath('.\mfiles')
% ALT (elevation model)
clear
dem=GRIDobj('./Topo/cop30DEM_utm32n_hochrhein_carved.tif');
% RAIN (sources (>0) and sinks (-1))
rain = GRIDobj('.\Topo\map_wc2.tif'); % MAP after WorldClim2 (mm/yr)
rain = resample(rain,dem);
rain = rain/1000; % convert to m/yr
rain = rain/3600/24/365.25/dem.cellsize.^2; % convert to m^3/s
% INFLOWS
inflow_rhine = 346*3; % (m^3/s)
inflow_rhine_y = [273:279]; % y-location (column) of inlet
inflow_rhine_x = ones(7,1)'*3266; % x-location (row) of inlet
inflow_rhine = inflow_rhine/length(inflow_rhine_x)/dem.cellsize.^2; % divide by number of inflow cells and by cellsize^2
rain.Z(inflow_rhine_y,inflow_rhine_x(1)) = ones(length(inflow_rhine_x),1)*inflow_rhine;
inflow_aare = 349*3; % (m^3/s)
inflow_aare_x = [1461:1462]; % y-location (column) of inlet
inflow_aare_y = ones(2,1)'*998; % x-location (row) of inlet
inflow_aare = inflow_aare/length(inflow_aare_x)/dem.cellsize.^2; % divide by number of inflow cells and by cellsize^2
rain.Z(inflow_aare_y,inflow_aare_x(1)) = ones(length(inflow_aare_x),1)*inflow_aare;
inflow_reuss = 140*3; % (m^3/s)
inflow_reuss_x = [1680:1682]; % y-location (column) of inlet
inflow_reuss_y = ones(3,1)'*998; % x-location (row) of inlet
inflow_reuss = inflow_reuss/length(inflow_reuss_x)/dem.cellsize.^2; % divide by number of inflow cells and by cellsize^2
rain.Z(inflow_reuss_y,inflow_reuss_x(1)) = ones(length(inflow_reuss_x),1)*inflow_reuss;
inflow_limmat = 114*3; % (m^3/s)
inflow_limmat_x = [1851:1852]; % y-location (column) of inlet
inflow_limmat_y = ones(2,1)'*998; % x-location (row) of inlet
inflow_limmat = inflow_limmat/length(inflow_limmat_x)/dem.cellsize.^2; % divide by number of inflow cells and by cellsize^2
rain.Z(inflow_limmat_y,inflow_limmat_x(1)) = ones(length(inflow_limmat_x),1)*inflow_limmat;
rain.Z(1:end,1)=-1;
rain.Z(1,1:end)=-1;
% rain.Z(end,1:end)=-1;
% rivers = GRIDobj('.\Topo\rivers2raster.tif');
% INITIAL SEDIMENT CONCENTRATION
cs = rain;
cs.Z(cs.Z<inflow_rhine)=0;
% WATER
% water = GRIDobj('.\Topo\HochRhein_WATER+LAKE_1000m.tif');
% UPLIFT
% uplift = GRIDobj('.\Topo\uplift_m_per_s_baselevel_Basel.tif');
uplift = GRIDobj('.\Topo\uplift_m_per_s_nagra_baselevel_Basel.tif');
uplift = resample(uplift,dem);
% SED (sediment thickness in meters)
sed = GRIDobj('.\Topo\mqu_140715g_utm32n.tif');
sed = resample(sed,dem);
sed.Z(isnan(sed.Z))=0;
LEM.dem = dem;
LEM.rain = rain;
LEM.sed = sed;
LEM.uplift = uplift;
% LEM.water = water;
LEM.cs = cs;
GRIDobj2grd(dem,['./Topo/',dem.name,'.alt']);
GRIDobj2grd(rain,['./Topo/',dem.name,'.rain']);
GRIDobj2grd(sed,['./Topo/',dem.name,'.sed']);
GRIDobj2grd(uplift,['./Topo/',dem.name,'.uplift']);
% GRIDobj2grd(water,['./Topo/',dem.name,'.water']);
GRIDobj2grd(cs,['./Topo/',dem.name,'.cs']);
%--------------------------------------------------------------------------
%% DEFINE INPUT PARAMETERS
%--------------------------------------------------------------------------
LEM.experiment = 'hochrhein'; % Project name
LEM.ErosPath = 'D:\\USER\\mey'; % Path to .exe
LEM.outfolder = 'hochrhein\\4rivers'; % folder to store results in
LEM.eros_version = 'eros7.3.112';
% LEM.inflow = 1060; % [m3s-1]water inflow at source cells
LEM.rainfall = 1; % Sets the precipitation rate per unit surface when multiplied by the rainfall map
LEM.initial_sediment_stock = 0.001; % % The total "stock" of sediment at the precipiton landing is: input_sediment_concentration*cs_map[i]*Precipiton_volume
LEM.inertia = 0; % refers to inertia term in shallow water equation
LEM.begin = 0; LEM.begin_option = 'time'; % start time
LEM.end = 90e5; LEM.end_option = 'time'; % length of model run
LEM.draw = 1000; LEM.draw_option = 'time'; % output interval
LEM.step = 0.05e2; LEM.step_option = 'volume';
LEM.stepmin = 0.5e1;
LEM.stepmax = 1e4;
LEM.initbegin = 1e+1; % initialization time (-)
LEM.initend = 1e+1;
LEM.initstep = 2;
LEM.TU_coefficient = 0.01; % sets the proportion of rain pixels that make up 1 TU
LEM.flow_model = 'stationary:pow';
LEM.erosion_multiply = 100; % multiplying factor for erosion rates. Equivalent to consider an "erosion time" larger than the hydrodynamic time
LEM.uplift_multiplier = 700000;
LEM.limiter = 1e-1;
LEM.continue_run = -1;
%--------------------------------------------------------------------------
% EROSION/DEPOSITION
%--------------------------------------------------------------------------
LEM.erosion_model = 'MPM'; % (stream_power, shear_stress, shear_mpm)
LEM.deposition_model = 'constant'; % need to know whether there are other options!
LEM.stress_model = 'rgqs:dir';
% ALLUVIAL
LEM.fluvial_stress_exponent = 1.5; % exponent in sediment flux eq. (MPM): qs = E(tau-tau_c)^a
LEM.fluvial_erodability = 2.6e-8; % [kg-1.5 m-3.5 s-2] E in MPM equation
LEM.fluvial_sediment_threshold = 0.05; % [Pa] critical shear stress (tau_c) in MPM equation
LEM.deposition_length = 30; % [m] xi in vertical erosion term: edot = qs/xi
% LATERAL EROSION/DEPOSITION
LEM.fluvial_lateral_erosion_coefficient = 1e-4; % dimensionless coefficient (Eq. 17 in Davy, Croissant, Lague (2017))
LEM.fluvial_lateral_deposition_coefficient = 0.5;
LEM.lateral_erosion_model = 1;
LEM.lateral_deposition_model = 'constant';
% BEDROCK
LEM.fluvial_basement_erodability = 0.1;
LEM.fluvial_basement_threshold = 0.5;
LEM.outbend_erosion_coefficient = 1.000000;
LEM.inbend_erosion_coefficient = 1.00000;
LEM.poisson_coefficient = 5;
LEM.diffusion_coefficient = 4;
LEM.sediment_grain = 0.0025;
LEM.basement_grain = 0.025;
%--------------------------------------------------------------------------
% FLOW MODEL
%--------------------------------------------------------------------------
LEM.friction_model = 'manning';
LEM.friction_coefficient = 0.025; %
LEM.flow_boundary = 'free';
%--------------------------------------------------------------------------
% OUTPUTS TO WRITE
%--------------------------------------------------------------------------
LEM.stress = 1;
LEM.waters = 1;
LEM.discharge = 1;
LEM.downward = 0;
LEM.slope = 1;
LEM.qs = 1;
LEM.capacity = 1;
LEM.sediment = 1;
LEM.flux =1;
LEM.stock =1;
LEM.str_write = '';
LEM.str_nowrite = '';
writeErosInputs(LEM);
%% run model
system([LEM.experiment,'.bat'])
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment