Fortran Coder

标题: 运行完出现fort.28文件 [打印本页]

作者: Zhuming    时间: 2020-6-20 18:52
标题: 运行完出现fort.28文件
说明:我是想将年数据,按1-12月提取出来,写入12个txt中。运行结果是:8月的数据和原始年数据一模一样。其他月份看上去正常。接着我将数据写入month8.txt的命令注释掉,运存出来,month8.txt依旧和原始数据文件一样。再接着,我把open()语句中month8.txt的文件通道号改成了33(原先是28),在运行时,month8.txt这次是没有数据的,但是多出来一个fort.28文件,大小和原数据文件一样大。因为原始文件比较大,传不了,我上传一个小一点的版本。
[Fortran] 纯文本查看 复制代码
program fkk
implicit none

    character(len=512)::cRead
!    character(len=20)::filename,DYfilename
    character(len=20)::wwlln_date,wwlln_time
    real::lat,lon
    character(len=5)::timingerror
    character(len=2)::station
    integer::status1=0
    integer::month,month1
    integer::i,j
   
    open(11,file='2014.txt',status='old',form='formatted')
    open(12,file='record.txt',status='replace')
    OPEN(21,file='month1.txt',status='replace',form='formatted')
    OPEN(22,file='month2.txt',status='replace',form='formatted')
    OPEN(23,file='month3.txt',status='replace',form='formatted')
    OPEN(24,file='month4.txt',status='replace',form='formatted')
    OPEN(25,file='month5.txt',status='replace',form='formatted')
    OPEN(26,file='month6.txt',status='replace',form='formatted')
    OPEN(27,file='month7.txt',status='replace',form='formatted')
    OPEN(33,file='month8.txt',status='replace',form='formatted')
    OPEN(29,file='month9.txt',status='replace',form='formatted')
    OPEN(30,file='month10.txt',status='replace',form='formatted')
    OPEN(31,file='month11.txt',status='replace',form='formatted')
    OPEN(32,file='month12.txt',status='replace',form='formatted')
   
        do while(.true.)
            READ(11,"(A512)",iostat=status1)cRead
            if (status1/=0) exit
            CALL parserRead(cRead)
!            write(*,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
!            READ(wwlln_date(1:4),"(I4)")year
!            write(*,*)year
            write(*,*) LEN_TRIM(wwlln_date)
            if(LEN_TRIM(wwlln_date)/=10) then    !这里是因为我的文件中大部分日期格式例如:2014/08/12,但有的月份是2014/8/10或者2014/8/1
                CALL specialread(wwlln_date)
                write(12,*)month1
            else
                READ(wwlln_date(6:7),"(I2)")month
                write(*,*)wwlln_date,month
            end if

            if(month==1) write(21,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==2) write(22,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==3) write(23,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==4) write(24,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==5) write(25,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==6) write(26,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==7) write(27,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==8) write(33,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==9) write(29,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==10) write(30,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==11) write(31,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            if(month==12) write(32,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            
            do i=21,32,1
                j=i-20
                if(month1==j) write(i,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            end do
            
        end do
    CLOSE(11)
    CLOSE(12)

        
        
        
contains
    !调用子程序
    subroutine parserRead(c)
        character(len=*)::c
        integer::j
        
        do j=1,LEN_TRIM(c)
            if (c(j:j)=="/") c(j:j)="|"
        end do
        READ(c,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
        do j=1,LEN_TRIM(wwlln_date)
            if (wwlln_date(j:j)=="|") wwlln_date(j:j)="/"
        end do
    end subroutine parserRead
   
    subroutine specialread(d)
        character(LEN=*)::d
        
        READ(d(6:6),"(I1)")month1

    end subroutine specialread
end program

QQ截图20200620184536.png (53.84 KB, 下载次数: 212)

QQ截图20200620184536.png

20140604.txt

914.83 KB, 下载次数: 0


作者: Zhuming    时间: 2020-6-20 18:58
另外改完文件号之后(month8.txt的,原先是28号,改成了33),运行结果month8.txt大小为0kb,就是没有数据写入。也就是说if(month==8) write(33,*)wwlln_date,wwlln_time,lat,lon,timingerror,station,这条命令无作用,是个摆设。
在线等,谢谢!
作者: Zhuming    时间: 2020-6-20 19:45
我笨的离谱哈哈哈....fort.28是因为我把文件号改了,但输出的还是28号文件。
另外,我尝试了把下面语句
            do i=21,32,1
                j=i-20
                if(month1==j) write(i,*)wwlln_date,wwlln_time,lat,lon,timingerror,station
            end do
这个语句放到else上边,然后将 if(month==1) write(21,*)wwlln_date,wwlln_time,lat,lon,timingerror,station....这些if都放到end if 前面,就没问题了。这个也可以用do 循环来写,一样的。
还有要加close,close也可以用do来写。我本来写了,后来为了测出到底哪里出问题,然后删了....
最后还是想请问下,为什么原先这样写,month8.txt会是原数据文件呢?





欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2