fortran 文件讀寫

program main
	integer i,in,flag,n
	real*8 ::V(4,2000)
	character*88 fname
	
        call getarg(1,fname)
	open(7,file=trim(fname),status='OLD')
	i=0
	flag=0
	do while(flag==0) 
	i=i+1
	!read(7,IOSTAT=flag) n,V(1,i),V(2,i)
	read(7,*,IOSTAT=flag) V(:,i)
	enddo
	
	close(7)
	in=i-1
	open(7,file='rf_'//trim(fname))
      write(7,*)in
	do i=1,in
	write(7,*) V(1,i),V(2,i),0.0
	enddo
	do i=1,in
	write(7,*) V(3,i),V(4,i),0.0
	enddo
	close(7)
      
	print*,'Job done..'


endprogram