forked from amforte/Topographic-Analysis-Kit
-
Notifications
You must be signed in to change notification settings - Fork 0
/
PlotIndividualBasins.m
76 lines (63 loc) · 2.07 KB
/
PlotIndividualBasins.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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
function PlotIndividualBasins(location_of_data_files)
% Function takes outputs from 'ProcessRiverBasins' function and makes and saves plots for each basin with stream profiles, chi-z, and slope area
%
%Inputs:
% location_of_data_files - full path of folder which contains the mat files from 'ProcessRiverBasins'
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Function Written by Adam M. Forte - Last Revised Fall 2015 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
cd(location_of_data_files);
%% Build File List
% Get Basin Numbers
AllFullFiles=dir('*_Data.mat');
num_basins=numel(AllFullFiles);
basin_nums=zeros(num_basins,1);
for jj=1:num_basins
fileName=AllFullFiles(jj,1).name;
basin_nums(jj)=sscanf(fileName,'%*6s %i'); %%%
end
FileCell=cell(num_basins,1);
for kk=1:num_basins
basin_num=basin_nums(kk);
SearchAllString=['*_' num2str(basin_num) '_Data.mat'];
SearchSubString=['*_' num2str(basin_num) '_DataSubset*.mat'];
if numel(dir(SearchSubString))>0
Files=dir(SearchSubString);
else
Files=dir(SearchAllString);
end
FileCell{kk}=Files;
end
fileList=vertcat(FileCell{:});
num_files=numel(fileList);
for ii=1:num_files
load(fileList(ii,1).name,'DEMc','Chic','Sc','SAc','RiverMouth','drainage_area');
f1=figure(1);
set(f1,'Units','inches','Position',[1.0 1.5 10 10],'renderer','painters','PaperSize',[10 10],'PaperPositionMode','auto');
clf
subplot(3,1,1);
hold on
plotdz(Sc,DEMc,'dunit','km','smooth',true);
xlabel('Distance (km)');
ylabel('Elevation (m)');
hold off
subplot(3,1,2);
hold on
plot(Chic.chi,Chic.elev,'-k');
xlabel('Chi');
ylabel('Elevation (m)');
hold off
a1=subplot(3,1,3);
hold on
scatter(SAc.a,SAc.g,'o','MarkerFaceColor','b','MarkerEdgeColor','k');
set(a1,'XScale','log','YScale','log');
xlabel('Log Area');
ylabel('Log Slope');
hold off
suptitle(['Basin Number: ' num2str(RiverMouth(:,3)) ' - Drainage Area: ' num2str(drainage_area)]);
fileName=['BasinPlot_' num2str(RiverMouth(:,3)) '.pdf'];
print(f1,'-dpdf',fileName);
close all
end
end