求助：fwrite出的grd文件画出来的图和直接画出来的不一致

一一一一

目标是从多个文件中提取满足某条件的数据，求平均值画图，先是直接画了一次，然后想做一下两个条件的差值，因为原文件比较多，决定分别把满足两个条件的数据写出来，再d差值，但是在写出来的那一步遇到了问题，画出来的图和直接画得不一样，如下：
两个图分别如下，左边是先画的图，右边是写出来的数据画的。

直接画的gs如下：
'reinit'

'sdfopen f:\data\hgt\hgt.interim.1982.nc'
……省略其中17个
'sdfopen f:\data\hgt\hgt.interim.2010.nc'

'enable print d:\grads\test.gmf'

'set lat 10 70'
'set lon 65 140'
'set lev 850'
'set grid off'
'set ccolor 1'

'set t 152'
'define a1=hgt'
……省略中间的37个
'set t 10414'
'define a39=hgt.19'

'define h=(a1+a2+a3+a4+a5+a6+a7+a8+a9+a10+a11+a12+a13+a14+a15+a16+a17+a18+a19+a20+a21+a22+a23+a24+a25+a26+a27+a28+a29+a30+a31+a32+a33+a34+a35+a36+a37+a38+a39)/3900'

'd h'
'print'
'disable print'

；

然后写grd的gs如下：
'reinit'

'sdfopen f:\data\hgt\hgt.interim.1982.nc'
省略如上
'sdfopen f:\data\hgt\hgt.interim.2010.nc'

'set gxout fwrite'
'set fwrite d:\data\hgt_n.grd'

'set lat 10 70'
'set lon 65 140'
'set lev 850'

'set t 152'
'define a1=hgt'
省略如上
'set t 10414'
'define a39=hgt.19'

'define h=(a1+a2+a3+a4+a5+a6+a7+a8+a9+a10+a11+a12+a13+a14+a15+a16+a17+a18+a19+a20+a21+a22+a23+a24+a25+a26+a27+a28+a29+a30+a31+a32+a33+a34+a35+a36+a37+a38+a39)/3900'

'd h'
'disable fwrite'
;


ctl文件如下：
dset d:\data\hgt_n.grd
title hgt
undef -9.99e08
xdef 51 linear 65 1.5
ydef 41 linear 10 1.5
zdef 1 levels 850
tdef 1 linear 00z01may1979 1440mn
vars 1
h 0 99   h
endvars


然后图是直接画的：

nolejkwq

应该是数据描述文件的格点和数据文件的没有对应好的原因吧，好好检查一下，重新写ctl文件应该能解决吧

男紫汗

据说要用x，y，而非lat，lon

一一一一

男紫汗 发表于 2018-5-7 19:19
据说要用x，y，而非lat，lon

谢谢！！！换了xy之后解决了！！！

LemonTree007

不明白楼主的tdef后面的1440mn什么意思，月份的话应该用mo