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

haidaowang

# n& A+ ~6 n. z, [最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
clear;

aa=load('rho_m2.dat');

% H6 S+ i, D7 }) Rcc=load('thet_m2.dat');
) l, a  b+ K' f" I2 I
bb=flipud(aa);
% Z( w0 }0 |) X1 _7 n2 y
3 g9 ?4 X/ ?( l9 @1 r+ y( }2 gcc=flipud(cc);

& \' S. K4 A& J" H: lbb(find(bb<0.01))=NaN;%陆地设为NaN
7 w0 R! G$ b) k. }
cc(find(cc<0.01))=NaN;
. p' w( }1 k& q  d  ?) J4 T
3 T) b# M" C' gm_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

lon1=linspace(117.50,126.83,56);

* x+ U: o. a5 y/ F  t. jlat1=linspace(34,41,43);
: K, A# B! G" q; U
[lon1,lat1]=meshgrid(lon1,lat1);
& |+ `. C$ G. @
8 e( `: H+ K& v1 k3 t, V[c,h]=m_contour(lon1,lat1,bb,15);

+ @3 e# |, D3 zm_coast('patch',[0 0 0])
/ u2 S; X. _, R  Y( _2 }
hold on
4 o" }; H* s2 x% {8 H7 a
dd=cc;
' o7 N3 @  `' R* k
dd(find(dd>340|dd<20))=NaN;

# l; z1 f) S# Y" M# g3 v) X; R' _[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

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

hold on

  O1 H8 d* l% y1 C) |% Da2=find(cc>340&cc<360);

& D" ^( C# ?4 L4 T% wcc(a2)=cc(a2)-360;

7 }% E3 V, G4 }% {/ x; v! i[c,h]=m_contour(lon1,lat1,cc,[0 0]);
/ |/ w% a& V) C
+ \8 b' ~( Y; x+ Y2 T$ e( n; qclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

. m1 O" q' {% Pm_gshhs_h('save','topodata');

m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);

m_grid('box','fancy')
8 U2 Y& O5 z1 }6 Z; Z/ H
xlabel('Longitude')

; x7 v7 P: u0 H) V$ A" ^: r- qylabel('Latitude')

/ B2 `* V0 p0 Ytitle('等振幅线与同潮时线')

5 d; M1 @# s0 O! Q


* Y0 u. u6 \( F从图中我们可以看到渤黄海的4个无潮点。
7 [+ \, l, A% ]) V: o

; m5 N" E" G! H; P0 T

yin123

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