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

uperrua

最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
3 e# q' O2 v$ R/ F/ m
clear;
1 {1 o1 N. S! D' y% t4 l, T
( i( E* Y. w. B) X1 |% Eaa=load('rho_m2.dat');  %振幅数据

! D5 S6 m8 x$ X3 L. Ycc=load('thet_m2.dat');   %迟角数据

  h6 B" z0 Q# X( A# Gbb=flipud(aa);

( F6 s# ^8 C5 V$ }# C$ d7 j% Ncc=flipud(cc);

bb(find(bb<0.01))=NaN;%陆地设为NaN

' C) x* L+ j+ q- R' Ccc(find(cc<0.01))=NaN;

m_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

lon1=linspace(117.50,126.83,56);
4 V* U) M( ~9 P6 Z7 D4 f/ w# ^
! w0 j) u' J' X  F2 hlat1=linspace(34,41,43);

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

7 y9 r5 ~  X! A% ~. Y# R[c,h]=m_contour(lon1,lat1,bb,15);
( W  |! Y0 t" }+ G
m_coast('patch',[0 0 0])
3 p- m/ a* F& L3 w
hold on
7 a1 P  H+ O* h2 x/ T; y
6 R# D8 E$ n0 S0 l) `9 ~% Rdd=cc;

dd(find(dd>340|dd<20))=NaN;
" G1 b2 C9 @+ G6 }; s  b+ Y
- f( a3 `; _' ]0 G& Y; ]* l[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

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

7 q8 G' Y- {- \( R; ^( r) qhold on

a2=find(cc>340&cc<360);
- [/ T/ ?4 N+ r1 E* u. S& q
8 [4 x) a- [3 d$ M7 fcc(a2)=cc(a2)-360;
# [6 r0 }) F/ R  H. s
8 j- H) u( H/ }7 j* U6 _6 x[c,h]=m_contour(lon1,lat1,cc,[0 0]);
( K+ S: P, H- Y. a& w6 Y0 D
: z9 o5 R4 f# f7 Zclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

m_gshhs_h('save','topodata');

" @. V4 L0 k7 N1 [m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);

m_grid('box','fancy')
; [6 G8 n/ A7 \* I
2 b; w4 V1 K! g8 X; ]xlabel('Longitude')

' X# U; \& H9 U4 _ylabel('Latitude')

) N( M% ]1 l- m( Ztitle('等振幅线与同潮时线')
& M5 D6 B4 `: l* C$ w- r. H
matlab绘制同潮时线和等振幅线
/ S, f. A1 H# `4 \


6 M, O% E7 t3 B% F: f" o+ f从图中我们可以看到渤黄海的4个无潮点。

CCxiaom

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