grads画降雨日变化图的问题

carpenter3

６月降水资料，格式为.nc，tdef为２４０，每个3小时一个值。要怎么画降水日变化图，我自己写的脚本不对，求指教！
'reinit'
'sdfopen    I:＼tpr.06.nc'
'set gxout line'
'set grid off'
'set grads off'
'set x 1 '
'set y 1'
'set cthick 5'
'define tpra=aave(prcv,lon=50,lon=110,lat=20,lat=50)'
t=1
while(t<241)
'set t 't
'define tpraa=(tpra(t=t)+tpra(t=t+8))/2'
t=t+1
endwhile

'set xlabs 00||06||12||18||24'
'set ccolor 2'
'd tpraa'
;

njzqxt

目测可能出来的图是二维。可以固定经度为一个点，纬度设为一个范围，或者相反，这样能出三维的。

lu'an

时刻注意维数环境。
define tpra之前没有set t xx，默认set t 1，导致时间维只有一个时次。
define tpraa时在循环里面，每次都会覆盖前面的。看你的意思不是这样。

冬日情愫+

在定义tpra的时候应该设定时间维度 set t 1 240
求日变化 楼主的循环貌似有问题吧，还是说的不清楚
不知道是不是可以这么写
set  lon 50 110
set lat 20 50
set t 1 30
define tpra=ave(prcv,t+0,t=240,8)
然后根据需要的空间位置设定
set x 1
set y 1
d tloop(aave(tpra,lon=50,lon=110,lat=20,lat=50)

carpenter3

njzqxt 发表于 2013-8-23 09:13

                               
登录/注册后可看大图

目测可能出来的图是二维。可以固定经度为一个点，纬度设为一个范围，或者相反，这样能出三维的。

我就只想画二位的图，就是先做区域平均后，每个时间对应一个值，最后效果是横轴是00时刻、03时刻、06时刻~24时刻。纵轴是变量大小。这个要怎么实现？

carpenter3

lu'an 发表于 2013-8-23 10:55

                               
登录/注册后可看大图

时刻注意维数环境。
define tpra之前没有set t xx，默认set t 1，导致时间维只有一个时次。
define tpraa ...

同二楼回复

carpenter3

冬日情愫+ 发表于 2013-8-23 13:29

                               
登录/注册后可看大图

在定义tpra的时候应该设定时间维度 set t 1 240
求日变化 楼主的循环貌似有问题吧，还是说的不清楚
不知道 ...

这样算出来时逐日变化，横轴成了1号、2号....。我想先做区域平均，每个时刻对应一个值，然后同一时刻的再求平均，最后横轴是00时刻、03时刻、06时刻等，纵轴是变量大小。怎么画？

冬日情愫+

carpenter3 发表于 2013-8-23 15:07

                               
登录/注册后可看大图

这样算出来时逐日变化，横轴成了1号、2号....。我想先做区域平均，每个时刻对应一个值，然后同一时刻的再 ...

那只需要将tpra之前的时间定义成 set t 1 9  . 这样得到的tpra 就是同一时刻的平均。 先做区域平均和后做区域平均应该没有区别吧 。你也可以先做区域平均，不过这个时候定义tpra的时候应该是 set  t  1 240 .

lu'an

carpenter3 发表于 2013-8-23 15:02

                               
登录/注册后可看大图

同二楼回复

'reinit'
'sdfopen    I:＼tpr.06.nc'
'set grid off'
'set grads off'
'set x 1 '
'set y 1'
'set z 1'
'set t 1 240'
'define tpra=aave(prcv,lon=50,lon=110,lat=20,lat=50)'
******************************
'set gxout line'
'set cthick 5'
'set xlabs 00||06||12||18||24'
'set ccolor 2'
'd tpraa'
;
这样应该可行。