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

haidaowang

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

aa=load('rho_m2.dat');

1 h/ @4 k$ G  }2 J; c& \cc=load('thet_m2.dat');
: t3 r0 C7 O/ a  E
$ P+ m$ [4 j( B4 sbb=flipud(aa);
- E7 C- @, Z' u# H
cc=flipud(cc);

/ G; S$ ?+ t: L( H/ D+ |/ Obb(find(bb<0.01))=NaN;%陆地设为NaN
5 o4 R. M, F0 P! u
cc(find(cc<0.01))=NaN;

1 i2 M+ f* a6 y9 {1 Xm_proj('miller','lon',[117.50 126.83],'lat',[34 41]);
& S, Z7 @, H; C+ ^
6 E4 W9 J( c' P  D2 ?: _lon1=linspace(117.50,126.83,56);
) n+ {, \7 _4 m, i5 D* z- H0 I
  p& F; f; f  s$ {. Ilat1=linspace(34,41,43);
) [/ ]- F4 o' y: H. N) D1 @
- _. p1 s! l. L8 K6 n, J3 ?* u[lon1,lat1]=meshgrid(lon1,lat1);

[c,h]=m_contour(lon1,lat1,bb,15);
) K! p" `" C! d/ |; l4 }5 J
m_coast('patch',[0 0 0])
3 n' i. B* r9 c3 I* }  i
* |/ X. _9 N9 h' Q) e# |hold on

dd=cc;

, l, N$ Q; w: Rdd(find(dd>340|dd<20))=NaN;
" u5 d# }# h- T0 Z. f
8 _2 M, ?/ t. O# J6 ~! c[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);
3 V0 U% n/ p$ a) G& Y- B, Q  X
clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

hold on
( M* o: q! U- \8 O4 U' W+ d
, @! m. C( l9 z8 `a2=find(cc>340&cc<360);
$ C$ s, [+ k( S' d1 D
* \+ D  f# c  X4 b7 g5 S! A6 ecc(a2)=cc(a2)-360;

[c,h]=m_contour(lon1,lat1,cc,[0 0]);

clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
4 q- L4 \( _& T* g9 t: g6 u  i& \
m_gshhs_h('save','topodata');
7 E0 M3 ?+ Z/ K. V7 x; A
m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);
/ Y: c4 E( I' ~4 z% L9 x; R  p
m_grid('box','fancy')
6 }  P4 A* V3 P' u5 e
xlabel('Longitude')

ylabel('Latitude')
* b0 ]# ^7 i4 x1 b
title('等振幅线与同潮时线')

+ ~7 n" f: `# r# l3 H, ]$ n6 o9 Y
" q  T( Q9 M9 U, u& H' M6 @% _

8 A& k( y+ R/ l$ B/ l从图中我们可以看到渤黄海的4个无潮点。
& |% x3 ~) v* S; C; ]1 q

yin123

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