WRF输出格点数据，如何用GRADS画单站点逐小时累积降水量随时间变化图图

shuihu88888

WRF输出格点数据，如何用GRADS画单站点逐小时累积降水量随时间变化图，如图
我的gs文件是：
'reinit'
'open E:\test.ctl'
'enable print E:\test.jpg'
'set lon 124.8'
'set lat 52.2'
'set t 25 73'
'set gxout contour'
'define p=rainc+rainnc'
it=26
while(it<=73)
'define tmp=p(t='it')-p(t='it-1')'
it=it+1
'd tmp'
endwhile
'print'
'printim E:\test.jpg white'
'disable print'

做个霸气的木头

shuihu88888 发表于 2012-8-24 12:06

                               
登录/注册后可看大图

你好，我的意思其实很简单，WRF模式模拟的一次降水，水平方向为经纬度，垂直方向为气压分层，我要画的是某 ...

水壶同学，你的目的是求相邻时刻的降水差值随时间的变化图，但是你用的是循环，每次循环求一个插值画出来，加 上你前面定义的时间是25 73，相当于每次循环画一条直线，就成了你的图。如果你想画差值随时间的变化图 应当把这个差值定义到一个变量里，然后一下子画出来。只要定义一个存储差值的变量就行了：
比如这样：
set t 1 10
define tmp=p(t+3)-p(t+2)
d  tmp
就是画的从3到12的10个小时的累计降水随时间变化

做个霸气的木头

如果你想解决这个问题，你得告诉大家：你的格点数据具体怎样的，ctl 的的代码应该贴出来，你遇到什么样的问题。

冰枫

wrf是移动网格内的数据吗？如果是移动的，那就不能这么直接减呀，经纬度不同会出错。我也没明白你到底是要问啥呢？

Forcast

我觉得楼主已经实现了这个功能吧

shuihu88888

冰枫 发表于 2012-8-23 22:30

                               
登录/注册后可看大图

wrf是移动网格内的数据吗？如果是移动的，那就不能这么直接减呀，经纬度不同会出错。我也没明白你到底是要问 ...

我的意思其实很简单的，就是如何画某一个站点的逐小时累积降水量，横坐标是时间，纵坐标是一小时累积一次的降水量

shuihu88888

做个霸气的木头 发表于 2012-8-23 21:56

                               
登录/注册后可看大图

如果你想解决这个问题，你得告诉大家：你的格点数据具体怎样的，ctl 的的代码应该贴出来，你遇到什么样的问 ...

你好，我的意思其实很简单，WRF模式模拟的一次降水，水平方向为经纬度，垂直方向为气压分层，我要画的是某一站点地表逐小时累积降水量随着时间的变化，横坐标为时间，纵坐标为每个小时的累积降水量，谢谢

shuihu88888

Forcast 发表于 2012-8-24 07:32

                               
登录/注册后可看大图

我觉得楼主已经实现了这个功能吧

你好，结果没有实现呢，画出来的结果是这样的，而不是我一开始贴出来的那样，不知道哪里出错了，求赐教

shuihu88888

做个霸气的木头 发表于 2012-8-24 15:20

                               
登录/注册后可看大图

水壶同学，你的目的是求相邻时刻的降水差值随时间的变化图，但是你用的是循环，每次循环求一个插值画出来 ...

谢谢阿木啦~~问题是解决了

ZOE

shuihu88888 发表于 2012-8-24 20:03

                               
登录/注册后可看大图

谢谢阿木啦~~问题是解决了

你好 想问你这个WRF的输出数据时逐小时的还是总的呢？
谢谢