如何在Fortran 77中使用此代码?

时间:2017-04-23 20:23:44

标签: fortran header-files raster

这个问题Reading freely available UVR data using gfortran on mac OSX正是我想弄清楚的。我不在Mac上,所以我不知道他们是如何让它工作的。我使用的是g77 compiler

我不明白如何修改此脚本以使其适用于我的文件。该脚本可用here

    - Read header (size= pixel size *2byte):
character head*14400

read(10,rec=1) head
read(head,'(2i6,2f8.2,f8.4,2e12.5,a1,a8,a1,a40)') 
     & npixel,nline,lon_min,lat_max,reso,slope,offset,',',
     & para,',',outfile

 - Read data (e.g., fortran77)
parameter(nl=7200, ml=3601)
... open file by "unformatted", "recl=nl*2(byte)" (,"bytereclen")

integer*2 i2buf(nl,ml)
do m=1,ml
 read(10,rec=1+m) (i2buf(n,m), n=1,nl)
 do n=1,nl
   par=i2buf(n,m)*slope+offset
   write(6,*) 'PAR[Ein/m^2/day]=',par
 enddo
enddo

0 个答案:

没有答案