一天四个时次求日平均，不知道两种方法的差异在哪里？

jolincai

对一天四个时次的NC资料求日平均，两天共731天2924个时次。以下两种方法提取资料，为什么画图的结果会不一样呢？不懂错的方法错在哪里？（和原始NC资料画图相同的方法算对）。求大神告诉错误在哪里。PS 两种方法画同一天的图 等值线相同，但是等值线所处的经纬度就不一样了，相差一定的经度（60个经度左右）和一定的纬度（30个纬度左右），只有set fwrite这一步不同，其他的CTL 和画图GS都是一样的，不知道为什么会出现这样的差别。
第一种：（和NC资料画图一样的）
'reinit'
'sdfopen E:\data\erai\temp\1995_19962mt.nc'
'set fwrite E:\data\erai\temp\1995_19962mtdaily.grd'
'set gxout fwrite'
'set x 1 144'
'set y 1 73'
i=1
while (i<=2921)
'd ave(t2m,t='i',t='i+3')'
i=i+4
endwhile
*'disable fwrite'
*'d ave(t2m,t=97,t=100)'
'disable fwrite'
;
第二种：（错误的）
'reinit'
'sdfopen E:\data\erai\temp\1995_19962mt.nc'
'set fwrite E:\data\erai\temp\1995_19962mtdaily.grd'
'set gxout fwrite'
'set lon 0 360'
'set lat -90 90'
i=1
while (i<=2921)
'd ave(t2m,t='i',t='i+3')'
i=i+4
endwhile
*'disable fwrite'
*'d ave(t2m,t=97,t=100)'
'disable fwrite'
;

这两种方法分别求出日平均后，分别画第25天的图（任何相同的一天都可以），即set t 25。 然后和原始NC资料 d ave(t2m,t=97,t=100)直接计算第25天的日平均画图，第一种方法图就一样；第二种方法图就不一样，即等值线相同，但相差一定的经纬度。实在想不通为什么，求解。
顺便问一下  ctl里 xdef linear 0 2.5和xdef linear -180 2.5有什么区别？

十里庐州啼

我按第一个gs做日平均，但是出图只有第一张跟原始一样，之后就不一样了，楼主做出来的是全部一样的么？

jolincai

十里庐州啼 发表于 2017-6-21 19:41
我按第一个gs做日平均，但是出图只有第一张跟原始一样，之后就不一样了，楼主做出来的是全部一样的么？

必须用格点

南南偶卖糕Oo

想问问楼主，t2m 是什么呀