122 |
axis equal, view(2) |
axis equal, view(2) |
123 |
|
|
124 |
|
|
125 |
|
% Plot the bathymetry at 1/2 resolution |
126 |
|
k = 3; |
127 |
|
for k = 1:5 |
128 |
|
disp(sprintf(' k = %d',k)); |
129 |
|
ginfo(k).c_half = ginfo(k).cor(1:2:end,1:2:end,:); |
130 |
|
ginfo(k).b_half = ginfo(k).bathy(1:2:end,1:2:end); |
131 |
|
ginfo(k).b_half(find(ginfo(k).b_half >= 0.0)) = 1000; |
132 |
|
if k == 2 |
133 |
|
hold on |
134 |
|
end |
135 |
|
% figure(k) |
136 |
|
surf( ginfo(k).c_half(:,:,1), ... |
137 |
|
ginfo(k).c_half(:,:,2), ... |
138 |
|
ginfo(k).c_half(:,:,3), ginfo(k).b_half ) |
139 |
|
if k == 5 |
140 |
|
hold off |
141 |
|
end |
142 |
|
end |
143 |
|
axis equal, view(2) |
144 |
|
|
145 |
|
|
146 |
if do_print > 0 |
if do_print > 0 |
147 |
print -depsc llc_bathy.eps |
print -depsc llc_bathy.eps |
148 |
end |
end |