以2013.txt为基准文件,剔除2014.txt中多余的部分。
代码:
program climate
implicit none
!声明变量
real::lon,lon1
real::lat,lat1
real(kind=4)::dataa,dataa1
integer::status1=0
!读取
OPEN(11,file='2013.txt',status='old')
OPEN(12,file='2014.txt',status='old')
OPEN(13,file='new2014.txt',status='replace')
do while(.true.)
READ(11,*,iostat=status1)lon,lat,dataa
if(status1/=0)exit
!write(*,*)lon,lat,dataa
do while(.true.)
READ(12,*,iostat=status1)lon1,lat1,dataa1
if(status1/=0)exit
!write(*,*)lon1,lat1,dataa1
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