将不同格点数计算结果画在同一张图上的问题，谢谢大家！

-逆光年年-

最近一直被一个问题困扰，在用wrf模式计算时，有两组数值试验中设置的水平格点数数不相同（file1、2、3的水平格点数都是602， file4是1602），其他的设置都一样的，想把两组实验计算得出的累计降水量结果用grads画在同一张图上，但是可能是由于x坐标定义不同的问题，总是不能画在一张图上，格点数不同的那个画不出来，我又在gs中间（'set dfile 4'前面）另外定义了下坐标，可是还是不成功，不知道是哪里出了问题，以下是我的gs文件，麻烦大家帮我看看哪里不对~谢谢大家了！！！！{:soso_e163:}{:soso_e183:}
'reinit'
'clear'
'open d:/wrf/test/test102000.ctl'
'open d:/wrf/test/test202000.ctl'
'open d:/wrf/test/test252000.ctl'
'open d:/wrf/test/test14.32000.ctl'

'set parea 0 11 0.0 8.5'
'set string 1 bc'
'set strsiz 0.15'
'set parea off'
'set csmooth on'
'q hardware'
rec = sublin(result,2)
_cols = subwrd(rec,4)
if (_cols=256)
  'set rgb 30 0 80 0'
  'set rgb 31 0 110 0'
  'set rgb 32 0 150 0'
  'set rgb 33 0 200 0'
  'set rgb 34 0 250 0'
  'set rgb 35 150 255 0'
endif
'set grid off'
'set mpdraw off'

'set xlopts 1 6 0.18'
'set ylopts 1 6 0.18'
'set string 2 c 6 0'
'set strsiz 0.15'
'set display color white'
t=19
while(t<22)
'c'
'set t ' t
'set x 221 381'

'set y 1 1'
'set z 1 1'
'set xaxis -80 80 10'
'set vrange 0 270'

*'set x 221 381'
*'set parea 1.2 9.0 1 8'
*'set t ' t
*'set z 1'
*'set xaxis -80 80 10'
*'set yaxis 0 270 1'

'set dfile 1'
'set cthick 6'
'set grads off'
'set ccolor 3'
'set cstyle 1'
'set cmin 0.01'
'set grads off'
'd RAINC+RAINNC'
pull tt
tt=(t-1)

'set dfile 2'
'set cthick 6'
'set grads off'
'set ccolor 5'
'set cstyle 3'
'set cmin 0.01'
'set grads off'
'd RAINC+RAINNC'
pull tt
tt=(t-1)

'set dfile 3'
'set cthick 6'
'set grads off'
'set ccolor 6'
'set cstyle 4'
'set cmin 0.01'
'set grads off'
'd RAINC+RAINNC'
pull tt
tt=(t-1)

'set x 721 881'
'set y 1 1'
'set z 1 1'
'set xaxis -80 80 10'
'set vrange 0 270'

'set dfile 4'
'set cthick 6'
'set grads off'
'set ccolor 4'
'set cstyle 2'
'set cmin 0.01'
'set grads off'
'd RAINC+RAINNC'
pull tt
tt=(t-1)

'set string 1 c 6 0'
'set strsiz 0.18'
'draw xlab  Distance(km)'
'draw ylab Total precipitation (mm)'
*'draw title Radar Echo (dBZ)'
'draw string 1.8 0.2 Time='tt'(Hour)'
if(rc='0')break;endif;
'set display color white'
'enable print precipitation.met'
'print'
pull dummy
'clear'
t=t+1
endwhile
pull dummy
'disable print'

传说中的谁

貌似不同格距的数据不能直接计算吧，先插值成相同格距？

-逆光年年-

传说中的谁 发表于 2012-7-3 11:43

                               
登录/注册后可看大图

貌似不同格距的数据不能直接计算吧，先插值成相同格距？

要做差值啊~是在grads里做吗？之前没有做过~望赐教~谢谢

mofangbao

插值的方法论坛太多了，直接用GrADS或者借助第三方工具，都有帖子可寻，楼主发挥下动手搜索的能力吧。

-逆光年年-

mofangbao 发表于 2012-7-3 18:50

                               
登录/注册后可看大图

插值的方法论坛太多了，直接用GrADS或者借助第三方工具，都有帖子可寻，楼主发挥下动手搜索的能力吧。

谢谢你的提醒，主要是我以前从来没有做过差值，不知道开始从哪里入手，我已经查过这里的帖子了，在网上也查过了，大家的帖子都针对各自问题，在讨论细节方面的东西，我这个新手看不懂阿～～所以还是想哪位能针对我这个问题指教一下，谢谢大家了！！！

-逆光年年-

顶一下~

-逆光年年-

哦~没动静~再顶一下~

-逆光年年-

心碎了

覃小覃

好像不同格点能在同一张图画吧，我画模式结果和实况就是不同的格点，但是grads一般会默认为第一个打开的文件的格点，下面打开的不同格点数据画图就会出错，我用set dfile这个命令强制了一下，下面交叉画图就没有报格点错误了，听说不用这个set dfile  直接用完一个文件就close那个文件就可以，这个命令的话顺序画图可以，但是交叉画图还是上面的命令号。

andrewsoong

俺是用的NCL画的~~~飘过~~~