【请教】怎么将站点资料缺测补全啊？

mofangbao

孤鹜 发表于 2012-4-19 16:24

                               
登录/注册后可看大图

我写了两个小程序来分步输出，可能是读取、判断语句写得复杂了，程序运行了十几分钟才结束。不过，总算完 ...

你可以把程序贴出来，这样喜欢fortran的朋友可以改进改进

孤鹜

程序一
integer ten,twenty,thirty,forty,fifty,seventy,hundred,v,q
character(5) number,stid
character(11) c
character(30) file
character(6) lon,lat
character(10) zuowu
open (unit = 1,file="F:\data\turangshidu\stationlist.txt")
open (unit = 2,file="tttt.txt")
open (unit = 3,file="F:\data\turangshidu\AGME_AB2_CHN_TEN_STATION.txt")


do while(.not. eof(1))
        read(1,10)stid
        10 format(a5)
        do while(stid.ne.number)
                read(3,30)quzhanhao,number,didian,lon,lat
                30 format(a8,a5,a27,2a6)
        enddo
        file="F:\data\turangshidu\"//stid//".txt"
        open(unit = 20,file=file)
        do while(.not. eof(20))
                q=0
                read(20,200)year,mon,xun,zuowu,ten,twenty,thirty,forty,fifty,seventy,hundred,c,v
                200 format(i4,i4,i4,a10,i9,6i7,a11,i8)
                if(c.eq."  6(未灌溉)") then
                        q=1
                endif
        enddo
        if(q.eq.1)then
                ten=-9999
                twenty=-9999
                thirty=-9999
                forty=-9999
                fifty=-9999
                seventy=-9999
                hundred=-9999
                v=-9999
                do i=2007,2011
                        do j=1,12
                                do k=1,3
                                        write(2,20)stid,lat,lon,i,j,k,ten,twenty,thirty,forty,fifty,seventy,hundred,c,v
                                        20 format(a5,2a9,i6,i4,i4,i9,6i7,a11,i8)
                                enddo
                        enddo
                enddo
        endif
enddo
close(1)
close(2)
close(3)
close(20)
end

孤鹜

程序二
character(9) lat,lon,b,p
character(5) a,stid
character(11) c,n
integer year,mon,xun,d,e,f,ten,twenty,thirty,forty,fifty,seventy,hundred,g,h,i,j,k,l,m,v,o,q


open(2,file="tttt.txt")
open(3,file="end.txt")
do while(.not. eof(2))
        q=0
        read(2,10)a,b,p,d,e,f,g,h,i,j,k,l,m,n,o
        open(1,file="2.txt")
        do while(.not. eof(1))
                read(1,10)stid,lat,lon,year,mon,xun,ten,twenty,thirty,forty,fifty,seventy,hundred,c,v
                if (stid.eq.a.and.year.eq.d.and.mon.eq.e.and.xun.eq.f)then
                        write(3,10)stid,lat,lon,year,mon,xun,ten,twenty,thirty,forty,fifty,seventy,hundred,c,v
                        q=1
                endif
        enddo
        close(1)
        if(q.eq.0)then
                write(3,10)a,b,p,d,e,f,g,h,i,j,k,l,m,n,o
        endif
       
enddo

10 format(a5,2a9,i6,i4,i4,i9,6i7,a11,i8)
end

孤鹜

其中中间运算输出文件“tttt.txt”为

546243727