grads 批量处理nc文件

若想要提取从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对数据进行随心所欲的处理。

 

能否直接生成一个文件还正在探索中。