GrADS画wrfout降水问题

weiguangjian

首先感谢各位看帖人！
我用GrADS画wrfout逐时降水，出现如下问题。贴出*.ctl, *.gs关键语句以及错误信息，请各位大虾指正！
********************************************
.ctl
********
dset ^case24_d03_lv19.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.3.1 MODEL
pdef  300 300 lcc  35.541  119.712  150.500  150.500  60.00000  30.00000  120.00000   1000.000   1000.000
xdef  784 linear  117.93935   0.00450450
ydef  633 linear   34.10793   0.00450450
zdef   19 levels
1000.00000
950.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
450.00000
400.00000
350.00000
300.00000
250.00000
200.00000
150.00000
100.00000
tdef   25 linear 00Z21SEP2012      60MN
VARS   35
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
********************************************************************
*.gs
*****
'open ./case24_d03_lv19.ctl'
'set lon 119.98'
'set lat 35.88'
'set z 1'
'set t 1 25'
'define p=rainc+rainnc'
'define p2=p(t+0)-p(t-1)'
'd p2'
**********************************************************************
错误信息
**********
No hardcopy metafile open
All files closed; all defined objects released
All GrADS attributes have been reinitialized
Notice: Implied interpolation for file ./case24_d03_lv19.ctl
Interpolation will be performed on any data displayed from this file
Data Request Error: Invalid grid coordinates
  World coordinates convert to non-integer grid coordinates
    Variable = rainc  Dimension = 0
  Error ocurred at column 1
DEFINE error:  Invalid expression.
Data Request Error: Invalid grid coordinates
  World coordinates convert to non-integer grid coordinates
    Variable = p  Dimension = 0
  Error ocurred at column 1
DEFINE error:  Invalid expression.
Syntax Error:  Invalid Operand
  'p2' not a variable or function name
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = p2
No hardcopy metafile open
All files closed; all defined objects released
All GrADS attributes have been reinitialized
*****
错误提示好像是坐标问题和定义变量出错，以前也这么做过，都能出图。整糊涂了！
******************************************************

weiguangjian

没人给意见呢{:soso_e150:}{:soso_e135:}{:soso_e118:}

天心云

没指定T，p2无法计算。

river

把'define p=rainc+rainnc'去掉，直接写成下面这样不就行了么
'define p2=rainnc(t+1)-rainnc+rainc(t+1)-rainc'

weiguangjian

river 发表于 2013-3-25 17:11

                               
登录/注册后可看大图

把'define p=rainc+rainnc'去掉，直接写成下面这样不就行了么
'define p2=rainnc(t+1)-rainnc+rainc(t+1)- ...

多谢指点！但是这样还是出不了图，错误提示如下，

Data Request Error: Invalid grid coordinates
  World coordinates convert to non-integer grid coordinates
    Variable = rainnc  Dimension = 0
  Error ocurred at column 1
DEFINE error:  Invalid expression.
Data Request Error: Invalid grid coordinates
  World coordinates convert to non-integer grid coordinates
    Variable = p  Dimension = 0
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = p
No hardcopy metafile open

weiguangjian

天心云 发表于 2013-3-25 16:21

                               
登录/注册后可看大图

没指定T，p2无法计算。

多谢指点！那是否可以写成下面的循环
i=1
while(i<25)
'set t 'i
'define p1=RAINC+RAINNC'
'set t 'i+1
'define p2=RAINC+RAINNC'
'define p=p2-p1'
i=i+1
endwile
'set t 2 25'
'd p'
我照这样试了试还是不行

river

weiguangjian 发表于 2013-3-25 18:20

                               
登录/注册后可看大图

多谢指点！但是这样还是出不了图，错误提示如下，

Data Request Error: Invalid grid coordinates

我不是说把'define p=rainc+rainnc'去掉了么，怎会还有P的错误啊

river

weiguangjian 发表于 2013-3-25 18:20

                               
登录/注册后可看大图

多谢指点！但是这样还是出不了图，错误提示如下，

Data Request Error: Invalid grid coordinates

用循环吧
tt=1
while(tt<=25)
'set t 'tt''
'define p2=rainnc(t+1)-rainnc+rainc(t+1)-rainc'
'd p2'
tt=tt+1
endwhile
但是这样求出来的不止一张图了，你想得到的就是一张图是吧？楼主应该把你想要计算的是什么说清楚啊

weiguangjian

river 发表于 2013-3-25 18:57

                               
登录/注册后可看大图

用循环吧
tt=1
while(tt

我画的是某一站点的逐时降水，经纬度都是固定的，就是需要得到一张曲线图，不知道下面这种写法对不对，试了试，还是有错：'set lon 119.98'
'set lat 35.88'
'set z 1'
i=1
while(i<25)
'set t 'i
'define p1=RAINC+RAINNC'
'set t 'i+1
'define p2=RAINC+RAINNC'
'define p=p2-p1'
i=i+1
endwile
'set t 2 25'
'd p'

weiguangjian

river 发表于 2013-3-25 18:40

                               
登录/注册后可看大图

我不是说把'define p=rainc+rainnc'去掉了么，怎会还有P的错误啊

是去掉了，我只是写了p这个变量名