xyz数据当中最后3列有0的部分不写入，但是其它部分正常写入

978142355

如题，我想让最后3列数据中，只要有数字0的就不输入到另一个文件中，比如正文第二行 1  C      6.079287   -0.155411   16.237832    82     2     3     7     0
输入到另一个文件时，变成1  C      6.079287   -0.155411   16.237832    82     2     3     7
我自己尝试编写程序，但是输出的结果变成了遇到后三列为0就不写入文件中了，请教各位如何实现我说的那个功能。

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

    Program removezero
    Implicit none
    Integer::i,j
    Character(len=20)::name,name1,name2,name3(2),name4(1)
    Real(kind=kind(0.0d0))::a,b,c
    Integer::d,e,f,g,h,k,ierr
    write(*,*) 'Input the name of xyz:'
    Read(*,*) name
    name1=trim(adjustl(name))//'.xyz'
    name2=trim(adjustl(name))//'1.xyz'
    Open(12,file=name1)
    Open(13,file=name2)

    Do i=1,1 
        read(12,*) name3(1),name3(2)
        write(13,"(2x,a,2x,a)") trim(name3(1)),trim(name3(2))
    End Do

    Do
        read(12,*,iostat=ierr) k,name4(1),a,b,c,d,e,f,g,h
        If(ierr==0) then
            If (f==0 .or. g==0 .or. h==0 ) cycle
            write(13,"(2x,I4,2x,a2,2x,F12.6,f12.6,f12.6,2x,I4,2x,I4,2x,I4,2x,I4,2x,I4)")k,&
            trim(name4(1)),a,b,c,d,e,f,g,h
        End if
        If(ierr/=0) exit
    End Do
    Read(*,*)
    Close(12)
    Close(13)
    End

fcode

也许这个代码是你需要的

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

Program removezero
  Implicit none
  Character(len=64):: name
  Real(kind=kind(0.0d0))::a(3)
  Integer:: z(5) , k , ierr
  write(*,*) 'Input the name of xyz:'
  Read(*,*) name
  Open(12,file=trim(adjustl(name))//'.xyz')
  Open(13,file=trim(adjustl(name))//'1.xyz')
  read(12,'(a64)') name
  write(13,"(2x,a)") trim(name)
  Do
    read(12,*,iostat=ierr) k,name,a,z
    If(ierr/=0) exit
    write(13,"(2x,i4,2x,a2,2x,3F12.6,5(2x,I4))") &
      k , trim(name) , a , pack(z,z/=0)
  End Do
  Read(*,*)
  Close(12)
  Close(13)
End Program removezero

978142355

fcode 发表于 2015-11-24 12:26
也许这个代码是你需要的

[mw_shl_code=fortran,true]Program removezero

老大，这确实是我所需要的，还是你厉害啊。我在你写的这个代码中，有一项我不太懂，那就是pack(z,z/=0)，虽然我也猜到了它的用途，不过它具体的用法和意思是什么？