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

uperrua

8 b, [( A, l' t- l最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
4 ?  a' R5 Y3 N: f) G$ Z; U
clear;
; ^7 F0 i( A# z! A  f
aa=load('rho_m2.dat');  %振幅数据
) t0 U9 M- h2 W; T& z7 N
# ]% u% w. h2 G$ g1 d/ {  x( u7 fcc=load('thet_m2.dat');   %迟角数据

  {0 F- y" _2 G; P9 h9 ]; L4 M, @/ jbb=flipud(aa);

3 y& q# l0 _) o6 y0 Ycc=flipud(cc);
! }, ~6 ~2 e6 G$ b. s1 M+ E/ P
bb(find(bb<0.01))=NaN;%陆地设为NaN
6 Y9 F7 r1 J) Y
cc(find(cc<0.01))=NaN;

9 y, R) |6 |2 R7 km_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

. ^+ N& u- f# U* Jlon1=linspace(117.50,126.83,56);

lat1=linspace(34,41,43);

[lon1,lat1]=meshgrid(lon1,lat1);
3 F- r6 H' K: h6 ]) ~2 D
3 l' \& C  a! h[c,h]=m_contour(lon1,lat1,bb,15);

m_coast('patch',[0 0 0])
; q6 l1 M6 b4 `5 B- g; @
hold on
$ `# N1 @  {. e$ i) ^. X
& F* U. \7 k5 X3 A& s& ydd=cc;

dd(find(dd>340|dd<20))=NaN;

' D  Z- Z6 o. o8 P( R[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

9 t6 T3 X* O, K5 _7 v: I7 `# {- tclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
  J6 `/ \4 ?  S: Q/ h
hold on
1 U2 M0 _% |7 l! O) T; k% M# k0 m
a2=find(cc>340&cc<360);
; l6 a* A% s0 _" y! P3 M- H
cc(a2)=cc(a2)-360;
3 o1 X: q5 ?7 |/ j! r5 k/ e
[c,h]=m_contour(lon1,lat1,cc,[0 0]);

( ]* P- |% z( m8 d3 Pclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
  G% p# o4 i) Y2 H. a
m_gshhs_h('save','topodata');

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

m_grid('box','fancy')

xlabel('Longitude')
4 C+ G. _0 `2 |. H
& r1 y4 l9 Y* i' J3 `& }2 kylabel('Latitude')

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

matlab绘制同潮时线和等振幅线
$ Y6 W8 q* C5 e/ J


从图中我们可以看到渤黄海的4个无潮点。

CCxiaom

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