這里提供一個(gè)使用TMD軟件讀取渤海潮汐數(shù)據(jù)的程序，使用了TPXO為中國(guó)渤黃海單獨(dú)制作的模型：

x=117.5:0.1:122.5;

y=37:0.1:41;

[X,Y]=meshgrid(x,y);

[amp,Gph,Depth,conList]=extract_HC('DATA/Model_YS',Y,X,'z');

%讀取振幅amp,遲角Gph,水深Depth,分潮名在conList里

figure()

m_proj('Mercator','lon',[116.5 124],'lat',[36 41.5]);

m_grid('box','fancy')

set(gcf,'color',[1 1 1]);

m_coast('patch',[0.5 0.5 0.5],'linewidth',1);

set(gca,'Fontsize',13)

xlabel('Longitude','Fontsize',15)

ylabel('Latitude','Fontsize',15)

hold on;

[c,h]=m_contourf(X,Y,Depth);

ch=colorbar;set(get(ch,'title'),'string','[meter]');%?

title('水深')

%讀取M2振幅，可以很清楚看到在黃海三角洲和秦皇島附近的無(wú)潮點(diǎn)

figure()

m_proj('Mercator','lon',[116.5 124],'lat',[36 41.5]);

m_grid('box','fancy')

set(gcf,'color',[1 1 1]);

m_coast('patch',[0.5 0.5 0.5],'linewidth',1);

set(gca,'Fontsize',13)

xlabel('Longitude','Fontsize',15)

ylabel('Latitude','Fontsize',15)

hold on;

[c,h]=m_contourf(X,Y,squeeze(amp(1,:,:)));

ch=colorbar;set(get(ch,'title'),'string','[meter]');%?

title('M2 振幅')