grads将wrf结果写成grd出错

alxnumb

      请求大神们帮助，我写了一个gs将模式转出的数据提取存放在grd里，结果所有时次的数据都一样，请问这是什么问题？非常感谢！！！

ctl文件dset Y:\2010\110d1.dat        
undef 1.e35
pdef  59  55 lcc   6.133   95.037    1.000    1.000  60.  30.  115.000  90000.  90000.
xdef  183 linear   77.9   0.40540546
ydef  102 linear    5.1   0.40540546
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        496 linear  0z01may2010  6hr
vars   46
U               19 0 U Compoment of wind                              
V               19 0 V Component of wind                              
W               19 0 W Component of wind                              
THETA           19 0 Theta                                             
TK              19 0 Temperature in K                                 
TC              19 0 Temperature in C                                 
P               19 0 Pressure (HPa)                                    
Z               19 0 Height (m)                                       
QVAPOR          19 0 Vapor                                             
QCLOUD          19 0 Cloud Water                                       
QRAIN           19 0 Rain Water                                       
QICE            19 0 -                                                
QSNOW           19 0 -                                                
QGRAUP          19 0 -                                                                  
RAINC           0  0 ACCUMULATED TOTAL CUMULUS PRECIPITATION           
RAINNC          0  0 ACCUMULATED TOTAL GRID SCALE PRECIPITATION        
                             
endvars
提取数据gs
'reinit'
'open Y:\110d1.ctl'
'set z 1'
'set t 1 80'
'set x 1 183'
'set y 1 102'
'set fwrite Y:\501.grd '
'set gxout fwrite'
'd rainnc'
'disable fwrite'
;
还想问一下应该用'set x 1 183'还是'set x 1 59'
183是经过插值后的格点数，不知道输出的话应该用哪个。感谢！

river

应该用循环吧
'reinit'
'open Y:\110d1.ctl'
'set fwrite Y:\501.grd '
'set gxout fwrite'
'set z 1'
'set x 1 183'
'set y 1 102'
tt=1
while(tt<=80)
'set t 'tt''
'd rainnc'
tt=tt+1
endwhile
'disable fwrite'
;

alxnumb

river 发表于 2012-12-14 08:27
应该用循环吧
'reinit'
'open Y:\110d1.ctl'

请问输出问输出文件x和y方向的格点数应该按照插值前的59 55，还是插值后的183 102呢

river

alxnumb 发表于 2012-12-14 08:49

                               
登录/注册后可看大图

请问输出问输出文件x和y方向的格点数应该按照插值前的59 55，还是插值后的183 102呢

你提取资料gs里的那个ctl是个怎么样的范围，就设置那个范围或者比那个范围小的范围就行啊。

alxnumb

river 发表于 2012-12-14 10:04
你提取资料gs里的那个ctl是个怎么样的范围，就设置那个范围或者比那个范围小的范围就行啊。

请问您知道降水的时间纬度图用gs应该怎么写么。我尝试了几次写出来都有问题。

alxnumb

river 发表于 2012-12-14 10:04
你提取资料gs里的那个ctl是个怎么样的范围，就设置那个范围或者比那个范围小的范围就行啊。

请问您知道降水的时间纬度图用gs应该怎么写么。我尝试了几次写出来都有问题。

river

alxnumb 发表于 2012-12-14 10:42

                               
登录/注册后可看大图

请问您知道降水的时间纬度图用gs应该怎么写么。我尝试了几次写出来都有问题。

那你可以把你的gs贴上来，让大家看看问题所在，一起交流。

alxnumb

river 发表于 2012-12-14 10:48

                               
登录/注册后可看大图

那你可以把你的gs贴上来，让大家看看问题所在，一起交流。

你好我发了一个纬度-时间图的帖子，里面有我写的gs请指教