Skip to content

Commit

Permalink
Sensitivity analysis for SlipRatio coefficient
Browse files Browse the repository at this point in the history
Adding scenario to investigate sensitivity of repose slope flux method to changes in SlipRatio coefficient
  • Loading branch information
RegMeasures committed Mar 5, 2017
1 parent 66bdd25 commit d209710
Show file tree
Hide file tree
Showing 2 changed files with 12 additions and 10 deletions.
Binary file modified Inputs/Scenarios.xlsx
Binary file not shown.
22 changes: 12 additions & 10 deletions SensitivityAnalysis.m
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,7 @@

%% Load sensitivity scenarios table
% Read excel file
[~, ~, raw] = xlsread('Inputs\Scenarios.xlsx','Scenarios','A3:R49');
[~, ~, raw] = xlsread('Inputs\Scenarios.xlsx','Scenarios','A3:S50');
% Allocate imported array to column variable names
clear Scenarios
Scenarios.ID = raw(:,1);
Expand All @@ -25,15 +25,16 @@
Scenarios.BankFlux = cell2mat(raw(:,7));
Scenarios.StoredBE = cell2mat(raw(:,8));
Scenarios.UpwindBedload = cell2mat(raw(:,9));
Scenarios.Geometry = raw(:,10);
Scenarios.BankTestWL = cell2mat(raw(:,11));
Scenarios.lb = cell2mat(raw(:,12));
Scenarios.ub = cell2mat(raw(:,13));
Scenarios.dT = cell2mat(raw(:,14));
Scenarios.Vgeometry = raw(:,15);
Scenarios.Vradius = cell2mat(raw(:,16));
Scenarios.VBankTestWL = cell2mat(raw(:,17));
Scenarios.Comments = raw(:,18);
Scenarios.SlipRatio = cell2mat(raw(:,10));
Scenarios.Geometry = raw(:,11);
Scenarios.BankTestWL = cell2mat(raw(:,12));
Scenarios.lb = cell2mat(raw(:,13));
Scenarios.ub = cell2mat(raw(:,14));
Scenarios.dT = cell2mat(raw(:,15));
Scenarios.Vgeometry = raw(:,16);
Scenarios.Vradius = cell2mat(raw(:,17));
Scenarios.VBankTestWL = cell2mat(raw(:,18));
Scenarios.Comments = raw(:,19);

Scenarios = struct2table(Scenarios);
Scenarios.ID = strtrim(Scenarios.ID);
Expand Down Expand Up @@ -74,6 +75,7 @@
Inputs.Bank.Flux.Approach = Scenarios.BankFlux(ScenNo);
Inputs.Bank.Update.StoredBE = Scenarios.StoredBE(ScenNo);
Inputs.ST.UpwindBedload = Scenarios.UpwindBedload(ScenNo);
Inputs.Bank.Flux.SlipRatio = Scenarios.SlipRatio(ScenNo);
if exist(Scenarios.Geometry{ScenNo}, 'file')
Inputs.Hyd.InitialGeometry = csvread(Scenarios.Geometry{ScenNo});
else
Expand Down

0 comments on commit d209710

Please sign in to comment.