- 积分
- 188
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-11-26
- 最后登录
- 1970-1-1

|
发表于 2020-2-1 10:29:57
|
显示全部楼层
clc,clear;
load coast
ncid = netcdf.open('E:\假期\Jtest.cam.h3.0000-01-01-00000.nc','NC_NOWRITE'); %打开nc文件
ncdisp('E:\假期\Jtest.cam.h3.0000-01-01-00000.nc'); %在命令窗中显示nc文件的变量
U850Data = ncread('E:\假期\Jtest.cam.h3.0000-01-01-00000.nc','U850'); %读入变量U850
V850Data = ncread('E:\假期\Jtest.cam.h3.0000-01-01-00000.nc','V850'); %读入变量V850
PRECCData = ncread('E:\假期\Jtest.cam.h3.0000-01-01-00000.nc','PRECC');%读入变量PRECC
LonData = ncread('E:\假期\Jtest.cam.h3.0000-01-01-00000.nc','lon'); %读入变量lon
LatData = ncread('E:\假期\Jtest.cam.h3.0000-01-01-00000.nc','lat'); %读入变量lat
mean(PRECCData,3);
mPRECCData=mean(PRECCData,3)*2;
[x,y]= meshgrid(LonData,LatData);
contourf(x',y',mPRECCData);
hold on
colorbar
mean(U850Data,3);
mU850Data=mean(U850Data,3)*2;
mean(V850Data,3);
mV850Data=mean(V850Data,3)*2;
h=quiver(x',y',mU850Data,mV850Data)
hold on
set(h,'color','g','linewidth',1);
hold on
axesm('eqaconic','MapLatLimit',[-20,40],'MapLonLimit',[40,180])
framem on;gridm on; plabel on;
labelx=[160 161]; labely=[-16 -15];
labelu=[10 0]; labelv=[0 0]; % set 10m/s as a standard
scale_auto=Get_Autoscale(LonData,LatData,mV850Data,mU850Data);
scale_label= Get_Autoscale(labelx, labely, labelu, labelv);
scale_factor=scale_auto/scale_label;
Handle = quiverm(labely, labelx, labelv, labelu, 2.0*scale_factor);
set(Handle, 'Color', 'r', 'LineWidth', 2.5);
text(160,-15,'10m/s','fontsize',12);
set(gca,'LineWidth',1,'FontSize',10,'Ylim',[-20,40],'Xlim',[40,180],'Position',[0.142,0.09,0.75,0.84]...
,'XTick',[40:20:180],'XTicklabel',{'40^°E','60^°E','80^°E','100^°E','120^°E','140^°E','160^°E','180^°E'}...
,'YTick',[-20:10:40],'YTicklabel',{'-20^°S','-10^°S','0','10^°N','20^°N','30^°N','40^°N'}); %添加经纬度信息
hold off;
title('P(mm/day)','position',[110 -25]);
netcdf.close(ncid);
楼主您好,我这样画出来后,不显示箭头,请问可以如何修改呀~ |
|