matlab 读取ROMSnc画海表面温度分布图
为什么画？主要来自粉丝询问：
更多内容查看原文链接请点击查看更多内容
请点击查看原文，会持续分享相关内容
他/她知道如何画图？
但不知道表层，可能对模式不熟悉。我也不熟悉。我给出两种方法：第一个读取坐标系sigma或者其它。看depth分布，这里这里就是-1到0的比例：那么-1是海底了。因此，你画图使用end即可pcolor(lon,lat,nanmean(temp(:,:,end,:),4))
第二种方法，海表面温度，不是第一就是最后，两个都画出来，表层温度大就解决了问题：
后台很多粉丝留言。持续更新，敬请~欢迎和感谢更多粉丝留言和推荐本公众号
代码：加载颜色包，和使用jet版本
.rtcContent { padding: 30px; } .lineNode {font-size: 12pt; font-family: “Times New Roman”, Menlo, Monaco, Consolas, “Courier New”, monospace; font-style: normal; font-weight: normal; }clear;clc;close all;ncdisp(“roms_his.nc”)file = “roms_his.nc”;temp = double(ncread(file,‘temp’));lon = double(ncread(file,‘lon_rho’));lat = double(ncread(file,‘lat_rho’));depth = double(ncread(file,‘s_rho’));time = double(ncread(file,‘ocean_time’))./86400+datenum(1858,11,17);time_end = datevec(time);%% 颜色包load(‘morecolor1_500.mat’)close allfigure(‘Position’,[50 50 1200 850],‘color’,‘w’)pcolor(lon,lat,nanmean(temp(:,:,end,:),4))shading interpcolormap(morecolor11_500)colorbarcaxis([25 30])title(‘roms海表面温度分布图’)export_fig(‘海表面温度分布图.jpg’,’-r600’)close allfigure(‘Position’,[50 50 1200 850],‘color’,‘w’)pcolor(lon,lat,nanmean(temp(:,:,1,:),4))shading interpcolormap(morecolor11_500)colorbarcaxis([25 30])title(‘roms海底温度分布图’)export_fig(‘海底温度分布图.jpg’,’-r600’)close allfigure(‘Position’,[50 50 1200 850],‘color’,‘w’)pcolor(lon,lat,nanmean(temp(:,:,end,:),4))shading interpcolormap(‘jet’)colorbarcaxis([25 30])title(‘roms海表面温度分布图’)export_fig(‘海表面温度分布图1.jpg’,’-r600’)close allfigure(‘Position’,[50 50 1200 850],‘color’,‘w’)pcolor(lon,lat,nanmean(temp(:,:,1,:),4))shading interpcolormap(‘jet’)colorbarcaxis([25 30])title(‘roms海底温度分布图’)export_fig(‘海底温度分布图1.jpg’,’-r600’)