用gs以后和原始数据直接处理效果不一样

Jeanela

大家帮我看看是什么问题 个人看不出来了  下面是我的gs和ctl
'reinit'
'open e:\dg3\dg3.ctl'
'set gxout fwrite'
'set fwrite e:\dg1\x1.dat'
'set lon 0 358.125'
'set lat -88.542 88.542'
t1=6
while(t1<=360)
z1=1
  while(z1<=28)
    'set z 'z1
    'define a=ave((CNVHRsig+LWHRsig+SHAHRsig+SWHRsig)*24*3600,t='t1',t='t1+2')'   
    'd a'
    z1=z1+1
  endwhile
t1=t1+12
endwhile
'disable fwrite'


dset  e:\dg1\x1.dat
title  latent heat
undef 9.999E+20
xdef 192 linear 0.000000 1.875
ydef 94 levels
-88.542 -86.653 -84.753 -82.851 -80.947 -79.043 -77.139 -75.235 -73.331 -71.426
-69.522 -67.617 -65.713 -63.808 -61.903 -59.999 -58.094 -56.189 -54.285 -52.380
-50.475 -48.571 -46.666 -44.761 -42.856 -40.952 -39.047 -37.142 -35.238 -33.333
-31.428 -29.523 -27.619 -25.714 -23.809 -21.904 -20.000 -18.095 -16.190 -14.286
-12.381 -10.476  -8.571  -6.667  -4.762  -2.857  -0.952   0.952   2.857   4.762
  6.667   8.571  10.476  12.381  14.286  16.190  18.095  20.000  21.904  23.809
25.714  27.619  29.523  31.428  33.333  35.238  37.142  39.047  40.952  42.856
44.761  46.666  48.571  50.475  52.380  54.285  56.189  58.094  59.999  61.903
63.808  65.713  67.617  69.522  71.426  73.331  75.235  77.139  79.043  80.947
82.851  84.753  86.653  88.542
zdef 28 levels
0.9950 0.9820 0.9640 0.9420 0.9155 0.8835 0.8454 0.8009 0.7499 0.6934 0.6323 0.5678 0.5012 0.4352 0.3717 0.3122 0.2582 0.2102 0.1682 0.1323 0.1023 0.778 0.578 0.413 0.278 0.174 0.093 0.025
tdef 30 linear 00Z06jun1979 1yr
vars 1
a  28  99  latent heat
endvars

斥鷃

我补充一下，没用ave处理前的原始数据的CTL中有若干变量，有些变量有28层，但部分变量只有21层甚至更少。

Jeanela

斥鷃 发表于 2014-4-15 20:41
我补充一下，没用ave处理前的原始数据的CTL中有若干变量，有些变量有28层，但部分变量只有21层甚至更少。

第一层数据和原始数据处理的就不重合啊  gs哪里错了  Ave取夏季的值要怎样取  不对？

守望之月/;-D

小菜鸟看不出来问题~~
有没有试过如果只做一层平均的话会有问题么？