Grads set fwrite问题

樱花阿喵

我有1948-2017年的1月500hPa位势高度场数据（hgt4817-500.grd），用set fwrite提取出1951-2012年1月的新grd文件(h5112.grd)，但用(hgt4817-500.grd)画1951年1月份和用（h5112.grd）画1951年1月份出来的图不一样，想知道原因是什么，求各位大神指教

Houckin

你的用来提取的gs文件里，可以看下吗

樱花阿喵

Houckin 发表于 2020-5-16 16:38
你的用来提取的gs文件里，可以看下吗

'reinit'
'open e:\dq\2\hgt4817-500.ctl'
'set fwrite e:\dq\3\h5112.grd'
'set gxout fwrite'
'set t 4 65'
'set lon 0 360'
'set lat -90 90'
'd hgt'
'disable fwrite'
;

Houckin

**这样试试？直接设set t 4 65，我不知道是不是这里有问题，我一般用下面的方式
'reinit'
'open e:\dq\2\hgt4817-500.ctl'
'set fwrite e:\dq\3\h5112.grd'
'set gxout fwrite'
**时间用循环
tt=4
while(tt<=65)
'set t  'tt
**一般坐标设成x,y不容易出问题
'set x 1  144'
'set y 1 73'
'd hgt'
tt=tt+1
endwhile
'disable fwrite'
;

樱花阿喵

谢谢你，问题已经解决啦

柚子

樱花阿喵 发表于 2020-5-22 10:25
谢谢你，问题已经解决啦

请问怎么解决的？用4楼的方法吗

樱花阿喵

用四楼方法也可以，我直接用的原始数据，只读需要用到的那一部分数据就可以了，希望可以帮到你。

yushaotian

Houckin 发表于 2020-5-16 21:20
**这样试试？直接设set t 4 65，我不知道是不是这里有问题，我一般用下面的方式
'reinit'
'open e:\dq\2\ ...

{:eb511:}{:eb511:}