多个文件的降水站点资料转格点的fortran求解

Elric

是这样的，在论坛找了好久之后呢，会转单时次单要素的资料，但是我需要批量处理很多.000文件，所以看看了批量处理的方法，写了程序运行试了一下，发现不大对好像，所以发上来求助一个各位大佬~
program test
    implicit none
    integer i
    character(40)::filename='F:\bishe\1998-2016\1998\98010108.000'
    character*8 stid
    real lon,lat,var,rain,time
    integer nlev,flag
    do i=101,131
    write(filename(28:30),'(i3)')i
    time=0.0
    nlev=1
    flag=1
open(1,file=filename,status='old')
open(2,file='f:\bishe\sta.grd',form='binary')
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
read(1,*)
10  read(1,*,end=100)stid,lon,lat,var,rain
    write(2)stid,lon,lat,var,rain
    goto 10
100 continue
    close(1)
    nlev=0
    print*,stid,lon,lat,var,rain
    write(2)stid,lat,lon,time,nlev,flag
    close(2)
    end do
    end

Elric

跪求各位大佬指个明路

Elric

QAQ各位大大，看一看瞧一瞧啊

Elric

求助啊老铁们

Elric

各位大佬指个明路