Plot air temperature data on pressure levels (Hovmöller diagram).

10 次查看（过去 30 天）

显示更早的评论

khan 2018-6-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/403749-plot-air-temperature-data-on-pressure-levels-hovmoller-diagram

回答： kafhi bachrudin 2023-5-29

采纳的回答： jonas

amj_w.zip

Hi, i need to plot air temperature at certain pressure levels. For reference a netcdf sample file is attached.

I want the data to be plotted for a certain latitude suppose (55) and a certain Longitude range suppose (60-90), and pressure levels suppose (1-17). In this case the X-axis will denote the Longitude range while Y-axis will show the pressure levels.

Any help will be appreciated.

3 个评论
显示 1更早的评论隐藏 1更早的评论

khan 2018-6-4

Sir, i think what you did is right, can you please share the procedure so i can see the steps and use. This data is Air temperature anomalies with unit of Kelvin. Thank You for helping

jonas 2018-6-4

Done! Didn't want to post an answer before it was correct. Happy to help.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

jonas 2018-6-4

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/403749-plot-air-temperature-data-on-pressure-levels-hovmoller-diagram#answer_323049

编辑：jonas 2018-6-4

在 MATLAB Online 中打开

Here you go! Note that you specify the latitude and range of longitudes as the indices of their respective vectors, and not in units of degrees. You can easily try other types of plots (like controuf) or change the colormap.

ncid = netcdf.open('amj_w.nc','NOWRITE')
ncdisp('amj_w.nc');
lon = ncread('amj_w.nc','longitude');
lat = ncread('amj_w.nc','latitude');
p = ncread('amj_w.nc','level');
T = ncread('amj_w.nc','t');
lat_const=50; %the latitude as the index in the lat vector
range=1:100; %specify the lon range as index... 
T=permute(T,[3 1 2]);
surf(lon(range),p,T(:,range,lat_const),'edgecolor','none');
cb=colorbar;
cb.Label.String='Temperature / K';
xlabel('longitude / degrees east')
ylabel('pressure / millibars')
set(gca,'XLim',[lon(min(range)) lon(max(range))])
view([0 90])
box on
title(['Latitude: ',num2str(lat(lat_const)),' degrees'])
set(gca,'layer','top')
colormap(parula)