ERA 风场数据

夏夜

201401.nc (1.49 MB, 下载次数: 37)

2015-2-3 12:17 上传

点击文件名下载附件

感谢提供下载的详细介绍：http://bbs.06climate.com/forum.php?mod=viewthread&tid=13288&highlight=era

下载ERA风场发现uv两个变量计算得到的风速和从文件中直接读取的风速值有差异（如图）。这是什么原因？

下载2014年一月的风场数据，选择10 metre U wind component、10 metre V wind component和10 metre wind speed 这三个变量

附件是2014年1月的数据和得到图。

读取数据和绘图程序如下：

file='201401.nc';
lat=ncread(file,'latitude');
lon=ncread(file,'longitude');
u=ncread(file,'u10');
v=ncread(file,'v10');
velsi=ncread(file,'si10');

latmin=0;
latmax=65;
lonmin=115;
lonmax=285;
Nlat=find(lat<=latmax&lat>=latmin);
Nlon=find(lon<=lonmax&lon>=lonmin);
lat=lat(Nlat);
lon=lon(Nlon);
[Lon,Lat]=meshgrid(lon,lat);
u=u(Nlon,Nlat);
v=v(Nlon,Nlat);
velsi=velsi(Nlon,Nlat);
u=u';
v=v';
vel=sqrt(u.^2+v.^2);
velsi=velsi';

figure
subplot(211)
m_proj('Mercator','lon',[115 285],'lat',[0 65]);
m_coast('line','color',[0 0 0],'linewidth',1);
m_grid('box','fancy',...
'xtick',5,'ytick',7,'tickdir','out',...
'fontsize',10);
hold on
m_contourf(Lon,Lat,vel,[min(vel(:)):0.1:max(vel(:))],'linestyle','none')
colorbar
caxis([0 10])
title('wind speed calculated by u and v')
subplot(212)
m_proj('Mercator','lon',[115 285],'lat',[0 65]);
m_coast('line','color',[0 0 0],'linewidth',1);
m_grid('box','fancy',...
'xtick',5,'ytick',7,'tickdir','out',...
'fontsize',10);
hold on
m_contourf(Lon,Lat,velsi,[min(velsi(:)):0.1:max(velsi(:))],'linestyle','none')
colorbar
title('wind speed from ncfile')

超台十八

图呢? 差异有多大，算上下限看看

夏夜

超台十八 发表于 2015-2-3 13:50
图呢? 差异有多大，算上下限看看

上传时漏掉了。

超台十八

可能资料是 short integer, 要转到 float?

超台十八

我刚也算了，真的不一样。。。

可能 si 是每时次风速的平均？

这个真不知道

超台十八

超台十八 发表于 2015-2-3 16:14
我刚也算了，真的不一样。。。

可能 si 是每时次风速的平均？

看来是的。你看看下面那张风速值都比上面高的。

夏夜

超台十八 发表于 2015-2-3 16:16
看来是的。你看看下面那张风速值都比上面高的。

u v是什么的平均? 选择的是Monthly Means of Daily Means

夏夜

超台十八 发表于 2015-2-3 16:02
可能资料是 short integer, 要转到 float?

转为double也是这样的

MeteoInfo

你下载的应该是月平均数据，平均后的U/V分量和平均后的风速是不可能对应的（非线性关系），这是正常现象。这也是为什么网站会提供月平均的风速数据，如果是瞬时数据就不会有风速变量。

夏夜

MeteoInfo 发表于 2015-2-3 16:31
你下载的应该是月平均数据，平均后的U/V分量和平均后的风速是不可能对应的（非线性关系），这是正常现象。

多谢，是月平均的。