为什么不同时间的位势高度场图相同？？

sth.special.

我要画的是1966—2005年 清明那天40年的平均位势高度场，使用的是NC逐日再分析资料
用grads画图时发现 每年清明那天的位势高度场图完全一样，因此40平均位势高度场图也与每年的图相同。想知道是怎么回事
将nc数据写成grd的gs如下
y.1=365;y.2=365;y.3=366;y.4=365;y.5=365;y.6=365;y.7=366;y.8=365;y.9=365;y.10=365;y.11=366;y.12=365;y.13=365;y.14=365;y.15=366;y.16=365;y.17=365;y.18=365;y.19=366;y.20=365;y.21=365;y.22=365;y.23=366;y.24=365;y.25=365;y.26=365;y.27=366;y.28=365;y.29=365;y.30=365;y.31=366;y.32=365;y.33=365;y.34=365;y.35=366;y.36=365;y.37=365;y.38=365;y.39=366;y.40=365;
r.1=5;r.2=5;r.3=5;r.4=5;r.5=5;r.6=5;r.7=5;r.8=5;r.9=5;r.10=5;r.11=4;r.12=5;r.13=5;r.14=5;r.15=4;r.16=5;r.17=5;r.18
=5;r.19=4;r.20=5;r.21=5;r.22=5;r.23=4;r.24=5;r.25=5;r.26=5;r.27=4;r.28=5;r.29=5;r.30=5;r.31=4;r.32=5;r.33=5;r.34=5;r.35=4;r.36=5;r.37=5;r.38=5;r.39=4;r.40=5;      
i=1
'set fwrite d:\lulu\xg\hgt.grd'
while(i<=40)
  'sdfopen d:\lulu\hgt\hgt.'1965+i'.nc'
*****************'写成GRD格式'*******************
  d.i=y.i-275+r.i
*say d.i y.i r.i
  'set gxout fwrite'
  'set lev 500'   
  'set lon 70 180'
  'set lat 0 55'
  'set t 'd.i''  
* say d.i
  'd hgt'
  i=i+1
endwhile
'disable fwrite'
;
该grd文件的ctl
dset d:\lulu\xg\hgt.grd
undef -999.0
title sample hgt data
xdef 45 linear 70 2.5
ydef 23 linear 0 2.5
zdef 1 levels 500
tdef 40 linear jan1966 1yr
vars 1
hgt 1 99 hgt data
endvars

传说中的谁

你试试看时间变量有木有改变嘛

sth.special.

传说中的谁 发表于 2012-5-10 16:10

                               
登录/注册后可看大图

你试试看时间变量有木有改变嘛

使用过 say语句显示的时间变量是变的 可是图出的都一样

言深深

sth.special. 发表于 2012-5-10 16:12

                               
登录/注册后可看大图

使用过 say语句显示的时间变量是变的 可是图出的都一样

不太懂grads，但是对此你是否可以查看一下你的数据是否一样，如果数据也是完全一样的，恐怕需要考察是不是数据搞错了的

sth.special.

言深深 发表于 2012-5-10 17:31

                               
登录/注册后可看大图

不太懂grads，但是对此你是否可以查看一下你的数据是否一样，如果数据也是完全一样的，恐怕需要考察是不是 ...

好的，我试试看能不能把数据输出来看看

pqman3

程序总共打开了40个文件，但是却没有切换文件号。这样后面39次循环和第1次循环一样，读的都是第1个文件，即1966年的数据。

Forcast

可以确定的是你出错了，错在哪里仔细找。

冰枫

我以前也碰到过，后来是把提取后的grd放到循环内，即一年产生一个grd，再用fortran写到一起的，这样就不会一样了。笨方法

sth.special.

冰枫 发表于 2012-5-11 08:30

                               
登录/注册后可看大图

我以前也碰到过，后来是把提取后的grd放到循环内，即一年产生一个grd，再用fortran写到一起的，这样就不会一 ...

好的 我试试，能出正确的结果就是好的 哈哈