若想要提取从1951-2006年56年nc文件中的某些数据,一个一个处理非常麻烦,这里介绍种较为简易的方法。例如想提取6-8月的位势高度资料。
'reinit'
t5=1951
*作文件名循环
while(t5<=2006)
'set gxout fwrite'
'set fwrite D:/sichuan/hgt1/'%t5%'.dat'
'sdfopen e:/ncep1/hgt/hgt.'%t5%'.nc'
t3=t5-1950
*判断是否为闰年
if(t3=2|t3=6|t3=10|t3=14|t3=18|t3=22|t3=26|t3=30|t3=34|t3=38|t3=42|t3=46|t3=50|t3=54)
to=153
else
to=152
endif
t4=to+91
while(to<=t4)
'set t 'to
t1=1
while(t1<=12)
'set z 't1
'set lon 80 140'
'set lat 15 55'
'd hgt'
t1=t1+1
endwhile
to=to+1
endwhile
*这里必须先观点上述运行的文件,grads最多同时可以打开20个文件左右。
'reinit'
t5=t5+1
endwhile
'reinit'
这样可以提取你想要的年数据,然后你大可运用fortran对数据进行随心所欲的处理。
能否直接生成一个文件还正在探索中。