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

uperrua

; s8 Q  a# e7 J2 N! z最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考

! N$ h# O( [$ }# i. @" K3 rclear;

5 }2 G2 u3 Y. ^% Z: paa=load('rho_m2.dat');  %振幅数据

cc=load('thet_m2.dat');   %迟角数据
8 B+ J6 y) ^0 I  c2 V0 C1 ?* |
bb=flipud(aa);
6 Y. F0 P  H! j; i9 `
cc=flipud(cc);
& ]8 i5 y- ^- v  F
$ c$ J6 G+ C8 e- a0 cbb(find(bb<0.01))=NaN;%陆地设为NaN

cc(find(cc<0.01))=NaN;
/ T, J2 d+ \' D4 v9 j% R6 @: K$ A
; y  Q  r' z6 z. `5 ?' K! j3 T! |m_proj('miller','lon',[117.50 126.83],'lat',[34 41]);
3 x+ k6 K% M& ~, `0 D+ j: a* p
lon1=linspace(117.50,126.83,56);
. o/ I% W6 K# _  c! J
' e4 L( J' T% O2 N) }6 mlat1=linspace(34,41,43);

[lon1,lat1]=meshgrid(lon1,lat1);

$ I0 d3 I& ^4 p$ H- F+ G$ ][c,h]=m_contour(lon1,lat1,bb,15);
/ h/ a& g7 t& k) m0 @1 Q$ l. d
m_coast('patch',[0 0 0])

' F$ d. X+ [% [hold on

dd=cc;

dd(find(dd>340|dd<20))=NaN;
0 |6 n* I( M" V0 b- G7 K) c
- R( i1 Z9 w" _- E[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

( H" Z* u4 X% t$ J. aclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

1 Y/ I6 F! _3 L5 [2 R" `hold on
" T5 Z4 B1 K; o2 U% N; N5 r, G5 a
a2=find(cc>340&cc<360);
: t& q: \/ r' ]
4 V6 d8 l) _( H- z4 c) qcc(a2)=cc(a2)-360;

4 r0 {) d0 C8 R0 @, S[c,h]=m_contour(lon1,lat1,cc,[0 0]);
; v/ }9 n0 B! I; C2 T, T6 L
9 ]" T  G/ N2 k$ a+ o5 t9 ^1 w, jclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

! v5 U3 F: _( N* f) P2 cm_gshhs_h('save','topodata');
, g/ b7 H2 Y& k: s" L
& m" d: g1 ]* F1 |# I" Jm_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);
9 \: B$ M  v+ ]6 ?2 h
m_grid('box','fancy')
" e' N+ D9 T- c' F
: D' \3 R  w; Y3 O( c8 |# Vxlabel('Longitude')
# B( ?  b' C- v$ z/ ]
ylabel('Latitude')

* c; y7 {  k* a# Z+ Ltitle('等振幅线与同潮时线')

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


( E& m9 u0 R6 D6 L3 U. u* A
从图中我们可以看到渤黄海的4个无潮点。

CCxiaom

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