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

[mw_shl_code=fortran,true]Program removezero

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

也许这个代码是你需要的

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

查看代码

打印?

01	Program removezero

02	Implicit none

03	Integer::i,j

04	Character(len=20)::name,name1,name2,name3(2),name4(1)

05	Real(kind=kind(0.0d0))::a,b,c

06	Integer::d,e,f,g,h,k,ierr

07	write(*,*) 'Input the name of xyz:'

08	Read(*,*) name

09	name1=trim(adjustl(name))//'.xyz'

10	name2=trim(adjustl(name))//'1.xyz'

11	Open(12,file=name1)

12	Open(13,file=name2)

13

14

Do i=1,1

15	read(12,*) name3(1),name3(2)

16	write(13,"(2x,a,2x,a)") trim(name3(1)),trim(name3(2))

17

End Do

18

19

Do

20	read(12,*,iostat=ierr) k,name4(1),a,b,c,d,e,f,g,h

21	If(ierr==0) then

22	If (f==0 .or. g==0 .or. h==0 ) cycle

23	write(13,"(2x,I4,2x,a2,2x,F12.6,f12.6,f12.6,2x,I4,2x,I4,2x,I4,2x,I4,2x,I4)")k,&

24	trim(name4(1)),a,b,c,d,e,f,g,h

25

End if

26	If(ierr/=0) exit

27

End Do

28

Read(*,*)

29

Close(12)

30

Close(13)

31

End