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

haidaowang

最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
clear;

1 V  o5 w2 d' f4 ]' Taa=load('rho_m2.dat');
+ ^" m) [4 p/ o1 P. r4 l$ ^
cc=load('thet_m2.dat');

$ M0 J. f+ S2 a; P0 E( lbb=flipud(aa);

$ Y% N% U. [  ^cc=flipud(cc);
1 U9 f7 ?; Z8 }, t% l
# S8 p8 a% h3 o/ d* ^bb(find(bb<0.01))=NaN;%陆地设为NaN
* {1 f: N' Q) o; R8 a/ V
cc(find(cc<0.01))=NaN;

m_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

lon1=linspace(117.50,126.83,56);

, d& H) i2 t# Xlat1=linspace(34,41,43);
$ f' c! W8 E# ~% o- ?( A7 U
* y3 h  e# s' W6 v( q2 b9 N[lon1,lat1]=meshgrid(lon1,lat1);

' L( u: i9 q2 |. K; ^) ][c,h]=m_contour(lon1,lat1,bb,15);

m_coast('patch',[0 0 0])
7 P' K* C$ F: N" X/ ~+ d& L
& a5 Y8 W. t  c3 thold on
! y; y' z9 M# c! i- |
' V' I" n3 g  ]& O' r* sdd=cc;

& x" ]0 O# S) `. U5 r. P5 y; Bdd(find(dd>340|dd<20))=NaN;

[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

hold on
- ]( i( k3 G% O% U# \, a
a2=find(cc>340&cc<360);

cc(a2)=cc(a2)-360;
( n+ C/ S$ U- A
[c,h]=m_contour(lon1,lat1,cc,[0 0]);
" z' t3 l1 B6 |5 L3 x' a: H( P6 M
clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

+ N* N. v7 F* M( wm_gshhs_h('save','topodata');
4 D# b7 C; Z9 X* u3 M& r
, k8 C% g/ I. j4 l$ ]m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);

0 L5 a; H3 ?. n3 sm_grid('box','fancy')
* f7 T7 s& z. D) ~
xlabel('Longitude')
2 t* G0 }( Y  L( y4 x5 }
ylabel('Latitude')

" L" S0 o. U3 }0 Ntitle('等振幅线与同潮时线')


" i1 X) ]( h3 F! `7 X
: X1 B( j+ m! E  _
从图中我们可以看到渤黄海的4个无潮点。

% J* h) y: M! j$ G( A7 L$ E+ E

yin123

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