subroutine writepdb(pdbfilename,numberofatoms,kindofatom, &
                    resnum,coord)
  use defs; implicit none
  character(len=*),intent(in) :: pdbfilename
  integer(i4b),intent(in) :: numberofatoms
  character(len=4),dimension(:,:),intent(in) :: kindofatom
  integer(i4b),dimension(:),intent(in) :: resnum
  real(dp),dimension(:,:),intent(in) :: coord ! in angstroms
  integer(i4b) :: i
  character(len=6), parameter :: cipher='ATOM'
  open(27,file=pdbfilename)
  do i = 1,numberofatoms 
!  write(27,'(a4,2x,i5,1x,a4,a1,a4,1x,i4,4x,3(f8.3))') &  ! old
   write(27,'(a4,2x,i5,1x,a4,a4,1x,i5,4x,3(f8.3))') &
    cipher,i,kindofatom(:,i),resnum(i),coord(:,i) 
  end do
  close(27)
end subroutine writepdb