求助!!!Fortran读取grd文件……

Fortran菜鸟……求助一下大家!!!
现在有一个grd文件,是用grads转的nc资料,nc资料是1900-2009逐月的2米气温。转的GRD数据是1900-2009年1320个月的逐月全球平均2米气温资料。想要做出1900-2009年逐年平均的资料,要怎么写这个程序呀!!!不太会写循环呀,大家帮帮忙吧……非常感谢!!!
这是写的gs:
'sdfopen f:\nc\1.nc'
'set gxout fwrite'
'set fwrite f:\1gb.grd'
'set x 1'
'set y 1'
'set z 1 1'
'set t 1'
i=1
while(i<=1320)
'd t2m'
i=i+1
endwhile
'define var=tloop(aave(t2m,x=1,x=721,y=1,y=361))'
'd var'
'disable fwrite'

之后要怎么用Fortran做呀?``