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

uperrua

, l; u9 E: g/ Q5 C) r7 `最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考

clear;

aa=load('rho_m2.dat');  %振幅数据
. S" r% t- b; D$ t$ p: w
cc=load('thet_m2.dat');   %迟角数据
- E7 `. V! y. q4 H$ E7 g/ N/ n
bb=flipud(aa);
# c% \3 N# K3 H- r' n' a
cc=flipud(cc);
& I0 D8 v. l$ O, g. g+ z9 O
, q# R  `$ l: _/ O' Bbb(find(bb<0.01))=NaN;%陆地设为NaN
2 ]4 ]- w' @7 B! k" q7 w2 d0 \+ l  I5 @
cc(find(cc<0.01))=NaN;

! D) d1 ^) V, _; i4 S/ Im_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

lon1=linspace(117.50,126.83,56);

lat1=linspace(34,41,43);

[lon1,lat1]=meshgrid(lon1,lat1);
% L: H4 E' j. D$ Y% O0 w
[c,h]=m_contour(lon1,lat1,bb,15);
+ b" ]9 T! u. u7 _# Z  W3 `
m_coast('patch',[0 0 0])

hold on

! i, S7 \, f2 j- H! F% ndd=cc;
/ x' g2 Q* n" `$ U4 _+ L
dd(find(dd>340|dd<20))=NaN;

[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);
5 s$ w& M; W5 M$ V  T7 W5 P7 \
3 R) o4 a" H/ U2 N' q! vclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
0 L4 T/ C+ [6 f5 Y. \/ S7 X
$ ]- K. B& W5 h. Phold on

2 G% x# H1 [' Ta2=find(cc>340&cc<360);
, a8 w+ z2 A% ?  {7 P/ S
cc(a2)=cc(a2)-360;
- A7 C) v1 _. ?' Y* K' y5 p# h
/ j* P/ T, N: [5 g+ G' s* u[c,h]=m_contour(lon1,lat1,cc,[0 0]);

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

m_gshhs_h('save','topodata');

m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);
6 d, ]  g$ j. T& I5 a$ U5 i5 v. k
m_grid('box','fancy')
: v+ u$ \' U  M' m- S8 R
xlabel('Longitude')

( m$ Y$ g$ f$ W# }0 _0 p$ h' sylabel('Latitude')

8 T0 L- }) t6 d! Wtitle('等振幅线与同潮时线')
/ U8 S/ I0 g( ?! ?) y1 E; W
matlab绘制同潮时线和等振幅线
3 T- R0 _- X7 ]5 a3 l9 F- U
% S+ f: R3 ]* @# ^' q8 x
( Q! q# v/ c- L. N% z( Z/ n; `
从图中我们可以看到渤黄海的4个无潮点。

CCxiaom

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