- #1
Physicslad78
- 47
- 0
Hello All
I am trying to compile a program and I guess it compiles fine but there is a problem in execution. I get the following error:
End of file
apparent state: unit 5 (ngh.in)
last format: list io
lately reading direct formatted external IO
Aborted
I am not sure what the error is actually (I am sure it is in open related to file ngh.in). I will post the program here and I appreciate it if someone can give me tips. Thanks
program find_neighbor
c
implicit double precision (a-h,o-z)
parameter(natm=12,nlmx=1000)
dimension aa(3,3),xa(3,natm)
c
open(5,file='ngh.in')
open(6,file='ngh.list')
c
read (5,*) ishell
write (6,'(/a,i1)') ' ishell = ',ishell
read (5,*) rc,alat
write (6,'(/a,2(3f8.3,2x))') ' rc and alat = ',rc,alat
rc=rc/alat
c
write (6,'(/a)') ' real space lattice vectors:'
do i=1,3
read (5,*) (aa(j,i),j=1,3)
write (6,'(3f8.3)') (aa(j,i),j=1,3)
end do
c
read (5,*) nat
write (6,'(/a,i2)') ' nat =',nat
c
write (6,'(/a)') ' atom coordinates:'
do iat=1,nat
read (5,*) (xa(j,iat),j=1,3)
write (6,'(i2,3f10.5)') iat,(xa(j,iat),j=1,3)
end do
c
call neigh(ishell,nat,aa,xa,alat,rc)
call sort(n,ra,rb,iatm,lv,wksp,iwksp)
call indexx(n,arrin,indx)
call shell(nn,ns,arn,axn,it)
c
end
I am trying to compile a program and I guess it compiles fine but there is a problem in execution. I get the following error:
End of file
apparent state: unit 5 (ngh.in)
last format: list io
lately reading direct formatted external IO
Aborted
I am not sure what the error is actually (I am sure it is in open related to file ngh.in). I will post the program here and I appreciate it if someone can give me tips. Thanks
program find_neighbor
c
implicit double precision (a-h,o-z)
parameter(natm=12,nlmx=1000)
dimension aa(3,3),xa(3,natm)
c
open(5,file='ngh.in')
open(6,file='ngh.list')
c
read (5,*) ishell
write (6,'(/a,i1)') ' ishell = ',ishell
read (5,*) rc,alat
write (6,'(/a,2(3f8.3,2x))') ' rc and alat = ',rc,alat
rc=rc/alat
c
write (6,'(/a)') ' real space lattice vectors:'
do i=1,3
read (5,*) (aa(j,i),j=1,3)
write (6,'(3f8.3)') (aa(j,i),j=1,3)
end do
c
read (5,*) nat
write (6,'(/a,i2)') ' nat =',nat
c
write (6,'(/a)') ' atom coordinates:'
do iat=1,nat
read (5,*) (xa(j,iat),j=1,3)
write (6,'(i2,3f10.5)') iat,(xa(j,iat),j=1,3)
end do
c
call neigh(ishell,nat,aa,xa,alat,rc)
call sort(n,ra,rb,iatm,lv,wksp,iwksp)
call indexx(n,arrin,indx)
call shell(nn,ns,arn,axn,it)
c
end