我试图创建一个Hovmoller图matlab但我很新的编程和不知道该做什么。我试图这样做的涡度。以下是我目前的代码。任何输入或建议,将不胜感激。谢谢! ! !

路径(路径“/用户/ zoleak /文件/”);lev_name =‘700’; % mb

disp (strcat (lev_name“计算意味着字段…”,“m”))

file1 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.T.00.2015.1p5.45N45S.nc ';file2 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.U.00.2015.1p5.45N45S.nc ';file3 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.V.00.2015.1p5.45N45S.nc ';file4 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.GEO.00.2015.1p5.45N45S.nc ';file5 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.VORT.00.2015.1p5.45N45S.nc ';file6 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.PV.00.2015.1p5.45N45S.nc ';

%读文件lat = ncread (file1、“纬度”);朗= ncread (file1、经度);水平= ncread (file1、“水平”);u_data = ncread (file2, ' u ');v_data = ncread (file3 ' v ');vort_data = ncread (file5“vo”);

%的格式文件

x =经度(28:121);y = lat一句子);[nlons, nlats ntimes] =大小(vort_data);u_data = u_data (28:121 1:31 4:);v_data = v_data (28:121 1:31 4:);vort_data = vort_data (28:121 1:31 4:);

u_data =挤压(紧缩(u_data));v_data =挤压(紧缩(v_data));vort_data =挤压(紧缩(vort_data));

% %

clf;lat_avg =意味着(vort_data (5:15,:,:), 1);LONLIMS = (-60 - -139.5);LATLIMS = 30 [0];[lon_plot, lat_plot] = meshgrid (x, y);项目(“墨卡托”、“朗”,LONLIMS, lat, LATLIMS);% m_quiver (lon_plot lat_plot u_data (:,:)”, v_data(:,:)”,“线型”,“-”,“颜色”,[0 0 0],' AutoScaleFactor ', 2);抓住;contourf (lon_plot lat_plot lat_avg (:,:)”、“LevelList”[6 e-4:1e-5:6e-4],“LineColor”、“没有”)海岸(“补丁”,(1 1 1)linewi, 2);网格(linewi, 0.5,“tickdir”,“出”,“行”,“没有一个”);集(gca、字形大小,14);集(gcf,“颜色”,“白”); hcb=colorbar('eastoutside'); title(['5N-20N,700mb Vorticity']); set(hcb,'YTick',[-1e-4:1e-5:1e-4]) caxis([-.5e-4 .5e-4]) save_name1=strcat('/Users/zoleak/Documents/figures/HOVMOLLER/PAT_VORT'); print(save_name1,'dpdf'-r600);

我试图创建一个Hovmoller图MATLAB但我很新的编程和不知道该做什么。我想做涡度。以下是我目前的代码。任何输入或建议,将不胜感激。

路径(路径“/用户/ zoleak /文件/”);lev_name = ' 700 '; % mb disp (strcat (lev_name“计算意味着字段…”,“m”)) file1 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.T.00.2015.1p5.45N45S.nc ';file2 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.U.00.2015.1p5.45N45S.nc ';file3 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.V.00.2015.1p5.45N45S.nc ';file4 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.GEO.00.2015.1p5.45N45S.nc ';file5 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.VORT.00.2015.1p5.45N45S.nc ';file6 = ' /用户/ zoleak /文档/数据/帕特丽夏/ ERAi.PV.00.2015.1p5.45N45S.nc ';%读文件lat = ncread (file1、“纬度”);朗= ncread (file1、经度);水平= ncread (file1、“水平”); u_data=ncread(file2,'u'); v_data=ncread(file3,'v'); vort_data=ncread(file5,'vo');%的格式文件x =经度(28:121);y = lat一句子);[nlons, nlats ntimes] =大小(vort_data);u_data = u_data (28:121 1:31 4:);v_data = v_data (28:121 1:31 4:);vort_data = vort_data (28:121 1:31 4:);u_data =挤压(紧缩(u_data));v_data =挤压(紧缩(v_data));vort_data =挤压(紧缩(vort_data));% %clf;lat_avg =意味着(vort_data (5:15,:,:), 1);LONLIMS = (-60 - -139.5);LATLIMS = 30 [0];[lon_plot, lat_plot] = meshgrid (x, y);项目(“墨卡托”、“朗”,LONLIMS, lat, LATLIMS);% m_quiver (lon_plot lat_plot u_data (:,:)”, v_data(:,:)”,“线型”,“-”,“颜色”,[0 0 0],' AutoScaleFactor ', 2);抓住;contourf (lon_plot lat_plot lat_avg (:,:)”、“LevelList”[6 e-4:1e-5:6e-4],“LineColor”、“没有”)海岸(“补丁”,(1 1 1)linewi, 2);网格(linewi, 0.5,“tickdir”,“出”,“行”,“没有一个”);集(gca、字形大小,14);集(gcf,“颜色”,“白”); hcb=colorbar('eastoutside'); title(['5N-20N,700mb Vorticity']); set(hcb,'YTick',[-1e-4:1e-5:1e-4]) caxis([-.5e-4 .5e-4]) save_name1=strcat('/Users/zoleak/Documents/figures/HOVMOLLER/PAT_VORT'); print(save_name1,'dpdf'-r600);