请问如果根据fortran程序写出GrADS所需的ctl文件呢？谢谢！

physzw

                               

数据是0.1o×0.1o的原数据提取的，覆盖范围为（120-123oE，36-38oN），是基于2012年排放清量，按减排控制计划，预测的2015年排放量。排放清单包括工业，电力、交通和居民源。

数据存储格式为无格式码，读取的fortran程序如下：

!howto compile:: pgf90 -byteswapio liu_read.f90 -o liu.out

!execution::./liu.out

!

parameter(nx=31,ny=21,np=4,nt=12)

character*80fname

realsx(nx,ny,np,nt)  !nx=e-w grids; ny=s-n grids; nt=12 months

                     !np=4::industry,energy, transportation, residential

fname='BC_monthly.grd'

open(11,file=trim(fname),status='old',&

form='unformatted',access='direct',recl=nx*ny*4)

kk=0

doit=1,nt

doip=1,np

   kk=kk+1

   read(11,rec=kk)((sx(i,j,ip,it),i=1,nx),j=1,ny)

enddo

enddo

close(11)

!writeout

end