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

uperrua

3 ?7 l! h  Q" S% ]最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
% L, F* s  E  p0 H
clear;
$ \& ]9 P2 \" s2 n2 I( g4 U  `
8 }/ V) I1 _0 `: B& C2 L4 A# Xaa=load('rho_m2.dat');  %振幅数据

* Y8 E1 H  r& f% Z: M7 ]6 \cc=load('thet_m2.dat');   %迟角数据

, m8 s- }8 Y+ K' q, T/ V1 R; f( O4 Nbb=flipud(aa);
' l; c1 l# }  ~# }
1 d! \& H8 R6 V# o  I  r- R3 g" Ucc=flipud(cc);

; R3 R9 {; j2 P3 H. e: g) D9 s; Y6 tbb(find(bb<0.01))=NaN;%陆地设为NaN

cc(find(cc<0.01))=NaN;
  A: @/ r1 M4 g$ r
( \6 y; f, ]+ Vm_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

lon1=linspace(117.50,126.83,56);

4 I2 S) g* `$ {; blat1=linspace(34,41,43);

[lon1,lat1]=meshgrid(lon1,lat1);
" E$ R$ d. ]0 |
[c,h]=m_contour(lon1,lat1,bb,15);
" m0 S( H: `5 X* {
" o0 [& I7 r3 I/ R( Lm_coast('patch',[0 0 0])

hold on
+ D2 M: ?$ ]* Q7 F' ~: C7 u
& ~$ K$ X0 w7 W( X! U  N: ^dd=cc;
6 @6 f: j( }) f; y2 N
' [/ _* J. ~& ~, I6 Add(find(dd>340|dd<20))=NaN;
6 h! M! b0 ~$ Q- `
0 m+ y* X3 T& t[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);
- m9 Q9 k$ T6 c7 k3 q: Z1 E& a, S
clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
  \) n4 I8 v! u1 N2 E. m
hold on

6 F: ?: f$ [0 t" w! ?3 j" G4 t. ^a2=find(cc>340&cc<360);

cc(a2)=cc(a2)-360;

' N. h( ^4 t5 F& _: e/ d2 Z[c,h]=m_contour(lon1,lat1,cc,[0 0]);

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

# ^+ N: v8 s& u( H8 Z. f/ |m_gshhs_h('save','topodata');
$ q0 U/ q& p! J1 q; z
m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);

m_grid('box','fancy')
+ c5 ~) l6 R9 p4 z, j3 r
' @3 t% k) L: {; }/ @3 v4 p4 _xlabel('Longitude')

( m  B( g* p) R3 T) |5 ~ylabel('Latitude')

title('等振幅线与同潮时线')
$ B) s- q+ @! U, C- Y
matlab绘制同潮时线和等振幅线



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

CCxiaom

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