怎么实现循环读取呢

njzqxt

下载了NCEP的月平均气温资料 air.mon.mean.nc（1948年1月至2012年4月），想作历年5月平均气温的时间剖面图，请问各位XD，怎样读出每年5月的数据呢？万分感激！

'reinit'
'sdfopen D:\air.mon.mean.nc'
'set lat 20 30'
'set lon 97.5'
'set z 6 '
i=5
while(i<=21)
'set t 'i''
'set gxout contour'
'd air'
pull dummy
'print'
'c'
i=i+12
endwhile
;

mofangbao

不明白你的提问，麻烦再描述准确一些，你下面的gs和你的内容是什么关系？

做个霸气的木头

你是要读数据还是画图？

njzqxt

麻烦DX们了，我要作历年5月平均气温的时间剖面图。

做个霸气的木头

直接想画出来比较困难，我建议你利用循环将每年五月份的数据输出成grd文件，再这个ctl，然后就可以画时间剖面图了！

做个霸气的木头

直接想画出来比较困难，我建议你利用循环将每年五月份的数据输出成grd文件，再写个对应c的tl，然后就可以画时间剖面图了！

njzqxt

谢谢版主指点！
我先 利用循环将每年五月份的数据输出成grd文件，GS如下：
'reinit'
'set gxout fwrite'
'set fwrite D:\air.mon.mean.grd'
'sdfopen D:\air.mon.mean.nc'
'set z 6'     
t=5            
while(t<=761)  
'd air'
t=t+12         
endwhile
'disable fwrite'

-----------------------------------------------------------------------------
'set z 6'      *只取500hPa
t=5            *1948年5月（起始值）
while(t<=761)  *2011年5月（终止值）
t=t+12         *只取每年5月的值

运行后能生成GRD文件，请问我的脚本正确吗？

做个霸气的木头

njzqxt 发表于 2012-5-30 17:00

                               
登录/注册后可看大图

谢谢版主指点！
我先 利用循环将每年五月份的数据输出成grd文件，GS如下：
'reinit'

经纬度的设定要加上的吧！你要得到时间剖面，那么经纬度要有一个固定住~

njzqxt

经纬度的设定应该可以在输出图形前设定吧？
CTL文件，请帮忙看看有无问题，谢谢！
dset D:\air.mon.mean.grd
undef 32678
title NCEP/NCAR REANALYSIS PROJECT
xdef   180 linear   0  2
ydef  89 linear   -88  2
zdef   1 levels   0
tdef 1948 linear   may1948     1yr
vars   1
air 0 99 yearly Means of air Temperature
endvars

做个霸气的木头

njzqxt 发表于 2012-5-30 17:22

                               
登录/注册后可看大图

经纬度的设定应该可以在输出图形前设定吧？
CTL文件，请帮忙看看有无问题，谢谢！
dset D:\air.mon.mean. ...

当然可以啦！你的时间应该不对，时间序列应该是年的个数，如果为1948~2011，应该是tdef 64 linear may1948 1yr