Fortran Coder

标题: 嵌套循环，不报错，但是大循环只执行一次 [打印本页]
作者: Zhuming    时间: 2020-5-27 16:25
标题: 嵌套循环，不报错，但是大循环只执行一次
为了对比两个txt，以2013.txt为准，删除2014.txt中和2013.txt不同的地方。
代码：
[Fortran] 纯文本查看 复制代码
program climate

implicit none



!声明变量

real::lon,lon1

real::lat,lat1

real(kind=8)::dataa,dataa1

integer::status1=0

integer::i,j



!读取

OPEN(11,file='2013.txt',status='old')

OPEN(12,file='2014.txt',status='old')

OPEN(13,file='new2014.txt',status='replace')



do j =1,4189,1

    READ(11,*,iostat=status1)lon,lat,dataa

    if(status1/=0)exit

    do i =1,4548,1

        READ(12,*,iostat=status1)lon1,lat1,dataa1

        if(lat/=lat1.OR.lon/=lon1) cycle

            write(*,*)lon1,lat1,dataa1

            write(13,*)lon1,lat1,dataa1

        if(status1/=0)exit

    end do

end do



CLOSE(11)

CLOSE(12)

CLOSE(13)



end program climate

错误信息如图        
   


climate.png (3.12 KB, 下载次数: 213)

climate.png
作者: li913    时间: 2020-5-27 22:13
[Fortran] 纯文本查看 复制代码
program climate

implicit none



!声明变量

real::lon,lon1

real::lat,lat1

real(kind=8)::dataa,dataa1

integer::status1=0

integer::i,j



!读取

OPEN(11,file='2013.txt',status='old')

OPEN(12,file='2014.txt',status='old')

OPEN(13,file='new2014.txt',status='replace')



do j =1,4189,1

    READ(11,*,iostat=status1)lon,lat,dataa

    if(status1/=0)exit

    do i =1,4548,1

        READ(12,*,iostat=status1)lon1,lat1,dataa1

        if(status1/=0)exit !************************

        if(abs(lat-lat1)+abs(lon-lon1)>1.0e-6) cycle !************************

        write(*,*)lon1,lat1,dataa1

        write(13,*)lon1,lat1,dataa1

    end do

    rewind(12) !************************

end do



CLOSE(11)

CLOSE(12)

CLOSE(13)



end program climate
作者: Zhuming    时间: 2020-5-28 11:15
li913 发表于 2020-5-27 22:13
[mw_shl_code=fortran,true]program climate
implicit none

学到了！！感谢



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