求助用grads提取数据时发生数据漂移

bright

利用grads计算1000hPa上位温并存储，但画图之后发现数据发生了漂移，找了很久也没找到错误在哪里，恳请各位帮忙看看。下面是gs:
*存储1000hPa上的位温
'reinit'
'open E:\5mydata\2ctl\fnl_20190401_00_00.ctl'
'set gxout fwrite'
'set fwrite e:/5mydata/6data/1000.grd'
*--------------------------
ii=1
while(ii<=21)
'set t 'ii''
'set z 1'
'set x 1 360'
'set y 1 181'
'd TMPprs*pow(1000./lev,0.286)'
ii=ii+1
endwhile
*'set t 1 21'
*'set z 1'
*'set x 1 360'
*'set y 1 181'
*--------------------------------
*'d TMPprs*pow(1000./lev,0.286)'
*--------------------------------
'disable fwrite'
数据配套的ctl 是
dset E:/5mydata/6data/1000.grd
title qse
undef 9.999E+20
options sequential
xdef 360 linear 0.0 1.0
ydef 181 linear -90.0 1.0
zdef 1 levels 1000
tdef 21 linear 00Z01apr2019 6hr
vars 1
qse 0 99 potential temperature
ENDVARS

bright

漂移的图如下

苏拉苏拉

x方向错位，感觉可能是经度起始位置不对造成的，原始数据ctl长什么样

278803532

你画图的时候用的是经纬度，那么写数据的时候最好也用经纬度，要是用X和Y，那你要算好第一个格点的经纬度。然后对应的CTL是不是第一个格点。
'set x 1 360'
'set y 1 181'
就是这两行，设置成经纬度就可以了哈。

bright

非常感谢各位回复，原始数据ctl如下：


dset e:\5mydata\6data\4ncep\fnl_2019%m2%d2_%h2_00.grib2
index e:\5mydata\6data\4ncep\fnl_20190401_00_00.grib2.idx
undef 9.999E+20
title e:\5mydata\6data\4ncep\fnl_20190401_00_00.grib2
* produced by g2ctl v0.1.1
* command line options: e:\5mydata\6data\4ncep\fnl_20190401_00_00.grib2
* griddef=1:0:(360 x 181):grid_template=0:winds(N/S): lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 1.000000 lon 0.000000 to 359.000000 by 1.000000 #points=65160:winds(N/S)

dtype grib2
options template
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 21 linear 00Z01apr2019 6hr
* PROFILE hPa
zdef 31 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000 700 500 300 200 100
options pascals
vars 122

bright

改成经纬度，图完全就不对了