下了ncep的1948年到2011年9月的位势高度场逐年1到12月平均资料，已经利用GrADS提取需要的逐年500hPa的位势高度场1到12月月平均资料(hgt.mon.mean.grd)，现在想求1948到2011年1到12月的分别的位势高度平均值，并写入新的文件(比如hgt.grd中)，请问fortran程序该怎么写，万分感谢！