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

uperrua

! A+ C9 @$ v- D5 m最近做海洋要素大作业模拟二维潮波需要绘制同潮时线和等振幅线，提供程序如下，仅供参考
9 L$ h. x7 ^  B6 m
) Q7 R' ?) A2 E0 ?clear;
' S% e# ]+ i, b# U, Y/ W
aa=load('rho_m2.dat');  %振幅数据

# f- P3 _' p, i) R0 G, ^cc=load('thet_m2.dat');   %迟角数据

bb=flipud(aa);

cc=flipud(cc);
  M* l: ^6 {" F. m
bb(find(bb<0.01))=NaN;%陆地设为NaN

7 i" G- Q' m( i1 Icc(find(cc<0.01))=NaN;

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

lon1=linspace(117.50,126.83,56);

- x4 Q+ U, V5 @" N$ V" Ulat1=linspace(34,41,43);

[lon1,lat1]=meshgrid(lon1,lat1);
$ f9 D% M! d5 S( S; w0 A6 Y0 p: X
[c,h]=m_contour(lon1,lat1,bb,15);
6 l' U! E3 A) @/ q! e. X7 E
m_coast('patch',[0 0 0])

, B. b1 {( d6 M* ~1 K/ j' P4 R3 Yhold on
4 J& D( S) ^6 |. i
dd=cc;
+ U: S: S  h) X7 R# G( O- [) h4 a
- T2 y! r& C" e1 w6 U; F  vdd(find(dd>340|dd<20))=NaN;

/ b6 n: O! M0 ]) h, D9 T$ C[c,h]=m_contour(lon1,lat1,dd,[30:30:330]);

- ?0 A& G7 u) x0 R8 y6 i' c' hclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);
0 Q# y1 }+ D8 n
hold on
2 |8 U9 i* q  R6 |# p
a2=find(cc>340&cc<360);

' |7 N1 e9 s! I2 X( u$ |cc(a2)=cc(a2)-360;

[c,h]=m_contour(lon1,lat1,cc,[0 0]);
$ F, k7 _/ i4 s& J0 k4 {
6 j/ k9 ~  [9 Z# Iclabel(c,h,'FontSize',10, 'Color','r','Rotation',0);

  h# C1 X4 q. J/ t& j, m* |5 H2 d( Zm_gshhs_h('save','topodata');

m_usercoast('topodata','patch',[0 0 0],'linewidth',1.2);
; m! `+ V2 R2 l
- c# @7 `1 _$ o4 dm_grid('box','fancy')
' I4 R7 J8 U6 T$ E. @) `8 {# m
xlabel('Longitude')

ylabel('Latitude')
. Q3 o. k8 b/ O- C+ ]) J
title('等振幅线与同潮时线')
7 d- h# Z! p8 L- c; @) i$ [4 F8 T6 z
. t6 p- c3 ^2 h1 ^$ Q& } matlab绘制同潮时线和等振幅线



从图中我们可以看到渤黄海的4个无潮点。

CCxiaom

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