使用TMD软件读取TPXO潮汐数据

这里提供一个使用TMD软件读取渤海潮汐数据的程序，使用了TPXO为中国渤黄海单独制作的模型：

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振幅，可以很清楚看到在黄海三角洲和秦皇岛附近的无潮点

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 振幅')