grads中nc中提取数据做风场合成分析

snowdrop917

本人grads入门菜鸟 第一次发帖 现在想做一下风场的合成分析，
1 从多年的uwnd.yyyy.nc和vwnd.yyyy.nc中选取有系统影响的时间段数据求合成分析，共53天的数据
2用'set fwrite G:\uwnd.dat'和用'set fwrite G:\vwnd.dat'将所需的数据放入dat中，并写了相应的uwnd.ctl和vwnd.ctl描述性文件，
3 最后想写个uv.gs画合成图
uv.gs如下
'reinit'
'open G:\uwnd.ctl'
'open G:\vwnd.ctl'
'set lat 10 30'
'set lon 90 130'
'set lev 500'
'set gxout barb'
'define u=ave(uwnd.1,t=1,t=53)'
'define v=ave(vwnd.2,t=1,t=53)'
'd u*2.5;v*2.5'
'printim G:\uv.png white'
'reinit'
;
出来的图如下：求解错误的原因？？？？？？？？ 望论坛各位帮助

uv

river

gs没什么问题么，你是觉得风速太大了是不是？是不是你的资料提取有问题呢。还有就是要看流场的形势你这个范围设置也不够大

ayzqs

可能是资料提取有问题，其他脚本贴出来让大家帮看看

随缘

为什么要放大2.5倍呢？

snowdrop917

第一次发帖就收到大家的回应 很是感动啊  真切的感受到了气象家园平台的力量，非常感谢各位
现在将我提取数据的方式放上来  各位看看问题出在哪？？？
1 uwnd1.gs
'reinit'
'set gxout fwrite'
'set fwrite G:\uwnd1.dat'
'sdfopen Y:\uwnd.1981.nc'
'set x 1 144'
'set y 1 73'
'set lev 500'
'set time 12z11jun1981'
'd uwnd'
'close 1'
*1
'sdfopen Y:\uwnd.1982.nc'
'set x 1 144'
'set y 1 73'
'set lev 500'
'set time 12z25jul1982'
'd uwnd'
'close 1'
*2
'sdfopen Y:\uwnd.1983.nc'
'set x 1 144'
'set y 1 73'
'set lev 500'
'set time 12z01nov1983'
'd uwnd'
'close 1'
*3
……
'sdfopen Y:\uwnd.2008.nc'
'set x 1 144'
'set y 1 73'
'set lev 500'
'set time 12z01nov2008'
'd uwnd'
'close 1'
*53
'disable fwrite'
'reinit'
;

2 uwnd1.ctl
DSET G:\uwnd1.dat
TITLE Upper Data
undef -9.99E33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1  levels 500
tdef 53 linear 12z11jun1982 6hr
vars 1
uwnd 1 99  
Endvars

3 uwnd2.gs
'reinit'
'set gxout fwrite'
'set fwrite G:\uwnd.dat'  
'open G:\uwnd1.ctl'
'set x 1 144'
'set y 1 73'
'set lev 500'
t=1
while(t<=53)
'set t 't
'd uwnd'
t=t+1
endwhile
'close 1'
'disable fwrite'
'reinit'
;
4 uwnd.ctl
DSET G:\uwnd.dat
TITLE Upper Data
undef -9.99E33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1  levels 500
tdef 53 linear 12z11jun1982 6hr
vars 1
uwnd 1 99  
Endvars

5 vwnd数据提取同uwnd

小草

对于500hPa上的风来说，确实是大了，看一下你原始数据文件中单位。看看是不是放大了10倍。

river

随缘 发表于 2012-12-10 21:25

                               
登录/注册后可看大图

为什么要放大2.5倍呢？

因为grads中一个风杆代表的是10节，如果你的风速是10m/s的话，在grads中才显示一个风杆。这个与我们平时一个风杆代表4m/s就不同了，所以要乘以2.5来使得grads中一个风杆代表4m/s

river

snowdrop917 发表于 2012-12-11 10:01

                               
登录/注册后可看大图

第一次发帖就收到大家的回应 很是感动啊  真切的感受到了气象家园平台的力量，非常感谢各位
现在将我提取数 ...

没看懂，怎么还提取了两遍

善良猪是我

学习学习了  这个我也有用  先留个脚印

xjliutao

NC的数据要怎么样处理呢