Matlab绘制同潮时线和等振幅线

haidaowang

最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
- L) Y) U; d0 }2 ^1 cclear;

6 O1 v5 h% Q$ i9 |* Naa=load('rho_m2.dat');
: f7 b  X! s" M, n/ ^4 |
cc=load('thet_m2.dat');

. M. ]+ m. [( T2 i3 R% Cbb=flipud(aa);

! t) e) ^3 q7 s, h, Jcc=flipud(cc);

bb(find(bb<0.01))=NaN;%陆地设为NaN
" Z& _/ V+ V$ x0 V6 d8 [+ g* X" g9 N
cc(find(cc<0.01))=NaN;

m_proj('miller','lon',[117.50 126.83],'lat',[34 41]);
/ @8 u% @1 o# z& g8 m
* G  ]% S; D1 slon1=linspace(117.50,126.83,56);

lat1=linspace(34,41,43);

: P( I. }1 }! P! E5 P1 N  y) d& Z[lon1,lat1]=meshgrid(lon1,lat1);

" Y3 X8 g4 x+ k4 u6 W! `. R[c,h]=m_contour(lon1,lat1,bb,15);
: O* H0 n3 [; i3 S  o
& N4 }7 D7 W+ B. A! Y+ @m_coast('patch',[0 0 0])

1 e9 N! g& s. @% i9 L2 jhold on
) l7 \3 ]' q9 E1 ?' B4 N! @: u
dd=cc;

/ l: N' W0 {9 P. ldd(find(dd>340|dd<20))=NaN;

[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);
2 ~! Q6 s& w9 F, M- k
clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
' N6 o" ?& W8 s  C
hold on
( u0 N, ^# S! ^+ x4 z
a2=find(cc>340&cc<360);

( I5 U/ S* J- B, D2 U# }cc(a2)=cc(a2)-360;
; C  |; V: Z0 V, V; J. ?
( R. ~( P7 t2 G0 r# G2 X[c,h]=m_contour(lon1,lat1,cc,[0 0]);

clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
# K; T; U/ J7 i" M0 P
m_gshhs_h('save','topodata');
1 d6 }* f$ G$ W: O! ?2 Y8 f) O% @
m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);
' z- ]& l, j- g# b7 m9 R) u
m_grid('box','fancy')
9 e9 \* j+ ?9 |5 J2 J0 d  y: Z6 K
xlabel('Longitude')
% F* v0 e5 P. u- c6 M3 B8 {6 A
- Y' z; Z! e* d/ |$ U" i  d7 `ylabel('Latitude')

title('等振幅线与同潮时线')




从图中我们可以看到渤黄海的4个无潮点。
, ~" I& _" L4 {2 k& s
. _' w$ M  M+ M* `7 G

yin123

Matlab绘制同潮时线和等振幅线