wrf输出数据提取物理量画出来的图与用wrf输出数据直接画的图错位

443753891

wrf输出数据提取物理量画出来的图与用wrf输出数据直接画的图错位
提取前画的正确图
提取数据画的图

不知为什么，请高手指导一下，以下贴出的提取和描述文件
wrf改动后的outd02.dat描述ctl:
dset F:\20090728\wrf\wrfout\outd02.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.2 MODEL
*pdef  186 192 lcc  30.567  117.229   93.500   96.500  60.00000  30.00000  121.00000  10000.000  10000.000
*xdef  507 linear  105.33969   0.04504504
*ydef  419 linear   20.93499   0.04504504
xdef  186 linear  105.33969   0.04504504
ydef  192 linear   20.93499   0.04504504
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   37 linear 18Z26JUL2009      60MN      
VARS  115
LU_INDEX       1  0  LAND USE CATEGORY (-)
U             19  0  x-wind component (m s-1)
V             19  0  y-wind component (m s-1)
W             19  0  z-wind component (m s-1)

提取gs文件：
'reinit'
'open  F:\20090728\wrf\wrfout\outd02.ctl'
'set fwrite  F:\20090728\wrf\rain\rain1.grd'
'set gxout fwrite'
'set t 1'
'set z 1'
'set x 1 186'
'set y 1 192'
'd RAINC.1(T=1)+RAINNC.1(T=1)'
i =1
while(i<37)
j=i+1

'set z 1'
'set x 1 186'
'set y 1 192'
'd RAINC.1(T='j')+RAINNC.1(T='j')-RAINC.1(T='i')-RAINNC.1(T='i')'
say 'u  'i'   'ok
i=i+1
endwhile

say 'u  'i'   '
'disable fwrite'
;

提取的的描述文件ctl
dset F:\20090728\wrf\rain\rain1.grd
undef 1.e30
title  OUTPUT FROM WRF V3.2 MODEL
pdef  186 192 lcc  30.567  117.229   93.500   96.500  60.00000  30.00000  121.00000  10000.000  10000.000
xdef  507 linear  105.33969   0.04504504
ydef  419 linear   20.93499   0.04504504
*xdef  186 linear  105.33969   0.04504504
*ydef  192 linear   20.93499   0.04504504
zdef   1 levels  1000.00000
tdef   37 linear 18Z27JUL2009      60MN      
VARS   1
RAIN1          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
ENDVARS

Forcast

建议你试试
*pdef  186 192 lcc  30.567  117.229   93.500   96.500  60.00000  30.00000  121.00000  *10000.000  10000.000
*xdef  507 linear  105.33969   0.04504504
*ydef  419 linear   20.93499   0.04504504
xdef  186 linear  105.33969   0.04504504
ydef  192 linear   20.93499   0.04504504

443753891

还是不行，pdef  186 192 lcc  30.567  117.229   93.500   96.500  60.00000  30.00000  121.00000  *10000.000  10000.000这行我在提取时打上星号了，不用这行

443753891

Forcast 发表于 2012-6-20 12:33

                               
登录/注册后可看大图

建议你试试
*pdef  186 192 lcc  30.567  117.229   93.500   96.500  60.00000  30.00000  121.00000  *1 ...

还是不行，
pdef  186 192 lcc  30.567  117.229   93.500   96.500  60.00000  30.00000  121.00000  *10000.000  10000.000这行我在提取时打上星号了，不用这行

catmiaow

我提取数据没设x和y的范围，直接全部提取，数据没有错，你试试。

443753891

catmiaow 发表于 2012-6-24 10:09

                               
登录/注册后可看大图

我提取数据没设x和y的范围，直接全部提取，数据没有错，你试试。

请问你ARWpost处理输出的CTL中xdef,ydef的参数修改没？我不设x和y的范围，不改动的话提取还是出错，画出的图更差，呈断续状。 如果改动的话还是老样子错位

catmiaow

443753891 发表于 2012-6-25 15:01

                               
登录/注册后可看大图

请问你ARWpost处理输出的CTL中xdef,ydef的参数修改没？我不设x和y的范围，不改动的话提取还是出错，画出的 ...

ARWpost处理我没修改，但是用grads提取的数据，编写ctl时改动了一些地方，把原来的ctl注释掉*options  byteswapped和*pdef  222 168 lcc  32.509  102.263  111.500   84.500  60.00000  30.00000  100.00000  10000.000  10000.000这两行。

443753891

catmiaow 发表于 2012-6-26 09:51

                               
登录/注册后可看大图

ARWpost处理我没修改，但是用grads提取的数据，编写ctl时改动了一些地方，把原来的ctl注释掉*options   ...

按照你的提取方法，画出的图不错位了，谢谢。

茉莉花

我用ARWpost处理wrf输出的数据，然后用grads显示，发现只有fields中list出的诊断变量（ fields = 'height,pressure,tk,tc'）可以正常显示，其他变量均不正常。这是什么原因？

443753891

茉莉花 发表于 2012-12-15 17:30

                               
登录/注册后可看大图

我用ARWpost处理wrf输出的数据，然后用grads显示，发现只有fields中list出的诊断变量（ fields = 'height,p ...

什么意思啊，会不会是其他的物理量你没输出吧？你最好贴上后处理生成的CTL文件。