Fortran Coder

标题: xyz数据当中最后3列有0的部分不写入，但是其它部分正常写入 [打印本页]

作者: 978142355 时间: 2015-11-24 10:57
标题: xyz数据当中最后3列有0的部分不写入，但是其它部分正常写入
如题，我想让最后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] 纯文本查看 复制代码

    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

AB22soakre.xyz

254.56 KB, 下载次数: 1

输入文件

作者: fcode 时间: 2015-11-24 12:26
也许这个代码是你需要的

[Fortran] 纯文本查看 复制代码

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 时间: 2015-11-24 14:30

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

[mw_shl_code=fortran,true]Program removezero

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

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)