| 1 |
% |
|---|
| 2 |
% load in nod, ele, and bat files |
|---|
| 3 |
% |
|---|
| 4 |
[ele,x,y,z,bnd]=loadgrid('MESHH/g2ssubgom'); |
|---|
| 5 |
% |
|---|
| 6 |
% initiate subplot and plot desired bathymetric contours |
|---|
| 7 |
% |
|---|
| 8 |
subplot(2,2,1),plotbnd(x,y,bnd); |
|---|
| 9 |
axis('equal'); |
|---|
| 10 |
set(gca,'XLim',[100000 300000]); |
|---|
| 11 |
set(gca,'YLim',[-50000 150000]); |
|---|
| 12 |
set(gca,'Box','on'); |
|---|
| 13 |
set(gca,'XTick',[]); |
|---|
| 14 |
set(gca,'YTick',[]); |
|---|
| 15 |
cval=[50 100 200]; |
|---|
| 16 |
lcontour2(ele,x,y,z,cval); |
|---|
| 17 |
hold |
|---|
| 18 |
% set(gca,'visible','off'); |
|---|
| 19 |
% |
|---|
| 20 |
% declare title, load and plot desired drogs |
|---|
| 21 |
% |
|---|
| 22 |
title('JA/N95/fixed/22days') |
|---|
| 23 |
[gridname,ndrog,ndts,tsec,pth]=read_pth('JA.22h',4); |
|---|
| 24 |
npthm=0 |
|---|
| 25 |
for nd=1:ndrog |
|---|
| 26 |
for ndt=1:ndts |
|---|
| 27 |
npthm=npthm+1; |
|---|
| 28 |
npth=nd+(ndt-1)*ndrog; |
|---|
| 29 |
for i=1:4 |
|---|
| 30 |
pthm(npthm,i)=pth(npth,i); |
|---|
| 31 |
end |
|---|
| 32 |
end |
|---|
| 33 |
npthm=npthm+1; |
|---|
| 34 |
for i=1:4 |
|---|
| 35 |
pthm(npthm,i)=0/0; |
|---|
| 36 |
end |
|---|
| 37 |
end |
|---|
| 38 |
plot(pthm(:,1),pthm(:,2),'y-') |
|---|
| 39 |
plot(pth(1:ndrog,1),pth(1:ndrog,2),'b*') |
|---|
