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

uperrua

最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考

clear;

aa=load('rho_m2.dat');  %振幅数据

0 D3 b, M, N2 P5 r5 {1 q) u% M, lcc=load('thet_m2.dat');   %迟角数据
! Y0 r+ l6 v3 e  h) B2 L, ^% }
bb=flipud(aa);
7 P3 {/ H) i! j- l) ^
cc=flipud(cc);

. m3 A# ~' r) D" S8 vbb(find(bb<0.01))=NaN;%陆地设为NaN
; l1 U2 i0 \0 t, [* E
. Y* z" {9 s0 _& o9 ?, hcc(find(cc<0.01))=NaN;

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

lon1=linspace(117.50,126.83,56);

+ Y" D" A/ a- B, t; c. T3 Jlat1=linspace(34,41,43);

" g$ S6 l4 J8 ~& H7 \( a[lon1,lat1]=meshgrid(lon1,lat1);
$ _  `' V# i8 i$ y1 o( X
$ K2 @4 k% R: x[c,h]=m_contour(lon1,lat1,bb,15);

m_coast('patch',[0 0 0])

hold on
3 K1 S6 P7 O9 j4 b
dd=cc;
# v( f/ c% }* V4 z# m
dd(find(dd>340|dd<20))=NaN;

[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

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

/ K  V8 \0 M1 u/ ], `+ y; fhold on
0 ?. j- ]: }: {1 ?
a2=find(cc>340&cc<360);

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

[c,h]=m_contour(lon1,lat1,cc,[0 0]);

' }/ ]# N! n) }( w3 Y& S4 n7 N8 A9 r* Uclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

m_gshhs_h('save','topodata');
6 O2 s0 p8 [" P; w& O# i
m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);

m_grid('box','fancy')

5 G4 N( R, H9 Y+ k8 A# C, J4 |6 Q7 gxlabel('Longitude')

; Q5 z9 g1 D2 w! G# N/ {* [/ x4 wylabel('Latitude')

title('等振幅线与同潮时线')
' ~3 n9 [6 u  M; ~0 A* B7 g
& a' n: ?3 D/ a: I3 Z4 D* {0 m matlab绘制同潮时线和等振幅线
) `  K$ K$ L' S
; P8 I- [' W& X7 W& m( T8 Y: S7 A
8 O0 B' [3 M0 s: N. w# q
* |% O, O1 |% j# I4 k, i6 Y/ G从图中我们可以看到渤黄海的4个无潮点。

CCxiaom

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