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

haidaowang

最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
8 |8 t0 H, G9 J( ~6 W5 _$ rclear;

aa=load('rho_m2.dat');
( w5 ?0 }" z( m: q* b0 [5 N
cc=load('thet_m2.dat');

5 b0 J  H- y0 z" E0 u9 h# Mbb=flipud(aa);

cc=flipud(cc);
8 P8 d# O8 c# F( I
bb(find(bb<0.01))=NaN;%陆地设为NaN

cc(find(cc<0.01))=NaN;
7 q' [) l# R/ J0 Z! u
m_proj('miller','lon',[117.50 126.83],'lat',[34 41]);

lon1=linspace(117.50,126.83,56);

/ J. d" v* B8 d, g, e7 hlat1=linspace(34,41,43);

[lon1,lat1]=meshgrid(lon1,lat1);
# d" _9 x4 e2 E
[c,h]=m_contour(lon1,lat1,bb,15);
( V4 x# Q) B* s9 z: G8 j' h
m_coast('patch',[0 0 0])
" w' \+ m3 A. h  A
hold on

( a# C5 o# Y2 ndd=cc;

dd(find(dd>340|dd<20))=NaN;
  d  y  C& |& y. E
% ?1 E4 x2 ^( w2 b( W5 c: @2 ]( T[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

- ^( l- j+ `) k% l4 |3 Z4 [clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

8 ^5 {! ~# t; O  n' Ehold on
. m* D# [# Q4 D
# w9 g+ T/ i- fa2=find(cc>340&cc<360);

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

" n% M, @- T. s/ {[c,h]=m_contour(lon1,lat1,cc,[0 0]);
3 }, S  ~+ y9 p; p5 E# Z0 J" C
clabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
8 E/ S* o+ a# U
m_gshhs_h('save','topodata');
( k% F2 j1 ]3 y
$ I0 o* K& Q7 d$ vm_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);

6 W: t3 a3 f% d5 l% E" Mm_grid('box','fancy')
3 v$ b0 t1 C( n
xlabel('Longitude')
3 J2 Q, f" e6 {3 M) J4 L" l
3 S3 f: s0 t$ p3 M" [8 _: wylabel('Latitude')

title('等振幅线与同潮时线')



" f) S* ^. w* P4 i" [
从图中我们可以看到渤黄海的4个无潮点。

& g) c* _; T3 N

yin123

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