nc转grd，经过画图验证与nc直接画图相差很大，求解

紫儿

我想把nc的文件转为grd，下列是我对应的gs文件，以及用来画图验证我的grd文件是否正确的ctl文件，但是验证结果和直接用nc资料话的相差很大，该检查的地方我都检查了没发现问题，请各位帮忙看下问题在哪里，谢谢
以下为nc自带的ctl文件

以下是转换grd的gs文件
'reinit'
'sdfopen d:\reason\q.nc'
'set fwrite d:\result\q.grd'
'set gxout fwrite'
'set lev 800'
num=1
while(num<=420)
'set t 'num
'set x 101 113'
'set y 32 44'
'd q'
num=num+1
endwhile
'disable fwrite'（其中经纬，高度都是我选取的需要的范围）
;
以下是画图用的ctl文件
dset d:\result\q.grd
title humidity
undef -9.99e+33
xdef 13 linear 100 1
ydef 13 linear 31 1
zdef 1 levels 800
tdef 420 linear 00Z01JAN1979 1mo
vars 1
q=>q 1 t,z,y,x  specific humidity
endvars
以下是用nc画图的gs文件
'reinit'
'sdfopen d:\reason\q.nc'
'set lon 100 112'
'set lat 31 43'
'set lev 800'
'define qn=ave(q,t=1,t=420)'
*'set mpdset cnworld'
'd qn'
;

麦田_smile

'set lev 800'
num=1
while(num<=420)
'set t 'num
'set x 101 113'
'set y 32 44'
'd q'
改为
'set x 101 113'
'set y 32 44'
'set lev 800'
num=1
while(num<=420)
'set t 'num
'd q'
试试看

river

图呢，原图和你自己的都贴上来。

mofangbao

图以及绘对比图的gs

紫儿

river 发表于 2015-4-30 12:34
图呢，原图和你自己的都贴上来。

用ctl画出的图

用nc画出的图

左边是ctl画出的图，右边是nc画出的图，都是求420个时次的平均
以下是用nc画图的gs文件
'reinit'
'sdfopen d:\reason\q.nc'
'set lon 100 112'
'set lat 31 43'
'set lev 800'
'define qn=ave(q,t=1,t=420)'
*'set mpdset cnworld'
'd qn'
;

紫儿

麦田_smile 发表于 2015-4-30 12:31
'set lev 800'
num=1
while(num

试过了，还是不行

lqouc

紫儿 发表于 2015-4-30 13:01
试过了，还是不行

输出的时候为什么是xy而不是lonlat？

紫儿

lqouc 发表于 2015-4-30 13:50
输出的时候为什么是xy而不是lonlat？

lon lat是用nc画图时设置用的，用我自编的ctl文件，因为ctl都已设置好了维度信息，直接就open ctl 然后d ave(q,t=1,t=420)画图了

lqouc

紫儿 发表于 2015-4-30 13:57
lon lat是用nc画图时设置用的，用我自编的ctl文件，因为ctl都已设置好了维度信息，直接就open ctl 然后d  ...

拜托好好看下我说的话啦，我是指你第一段gs。输出成grd是你用的是x和y，整个帖子你一直都没给出nc文件的ctlinfo，谁知道你设置的'set x 101 113' 'set y 32 44'是不是就能代表，lon 101 113和lat 32 44.

紫儿

lqouc 发表于 2015-4-30 14:12
拜托好好看下我说的话啦，我是指你第一段gs。输出成grd是你用的是x和y，整个帖子你一直都没给出nc文件的c ...

不好意思，我理解错了，图片为nc的ctl文件，分辨率是1
我是想取北纬31-43，东经100-112,经过计算设置的'set x 101 113'
'set y 32 44'