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

uperrua

最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
% o, l& A: \& I3 z: A5 l( Y
; C9 E" h& }% c2 `clear;
, R6 i* H/ M5 g  h9 Z
# c" B" A9 z0 T% f6 kaa=load('rho_m2.dat');  %振幅数据

cc=load('thet_m2.dat');   %迟角数据

bb=flipud(aa);

  e% [: }/ ]6 V- m7 J# x  xcc=flipud(cc);

# Z# C+ |! p5 Xbb(find(bb<0.01))=NaN;%陆地设为NaN
8 w( W+ i* p  O
cc(find(cc<0.01))=NaN;
6 G7 P+ l3 y/ {; H" C# D  Y
m_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

lon1=linspace(117.50,126.83,56);

$ r) b! C2 }5 w% w5 t& A: nlat1=linspace(34,41,43);
- [0 K6 g3 W8 w4 t9 W( Q3 _) w
[lon1,lat1]=meshgrid(lon1,lat1);
7 d7 Y( k/ a. a9 q% I
[c,h]=m_contour(lon1,lat1,bb,15);
# N7 G( _1 f$ @% n6 j
: b0 r, o6 k$ G/ E, x6 Jm_coast('patch',[0 0 0])
: E, u4 K$ W& v" I
hold on

dd=cc;
5 U* i; n/ k" E4 A# G
dd(find(dd>340|dd<20))=NaN;
1 F' z  j2 m( }& S( s% e* x& V/ ^
! ]; P3 s2 g  {3 @0 I/ x[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);
. b* K# O9 A/ P3 }9 r0 X" T
clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

hold on

, W; t; W9 e# U& n# y: l' Z% fa2=find(cc>340&cc<360);
: P; y( \" K& j& M: G1 t# G- u
) L% T; i# l/ ~( \6 X; n  c3 l& Scc(a2)=cc(a2)-360;

6 {& r( U- C0 u4 Z) P8 S' [( d[c,h]=m_contour(lon1,lat1,cc,[0 0]);
* v. G4 V* M" h
clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
6 t) i; ]5 a, X" U4 U
* n: }& J4 Q0 l' ]2 F" l: ?m_gshhs_h('save','topodata');
5 D8 B4 h+ X' |0 L7 e
m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);
/ C+ I4 c9 D( F) w1 G
m_grid('box','fancy')
5 H: {. M8 E0 \8 E( m$ Z# p
xlabel('Longitude')

ylabel('Latitude')

title('等振幅线与同潮时线')
( u* ~, o* W- |* w, P
4 {' F% ?6 k9 {# O matlab绘制同潮时线和等振幅线
9 ~$ ~* l* T( K) a3 S
- h1 T6 l* q1 e* L* A# M
* C9 U: l! c9 j8 E6 {+ `: G$ i
, D: p$ c7 A; v! c( F1 H' }从图中我们可以看到渤黄海的4个无潮点。

CCxiaom

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