|
|
@ -133,14 +133,14 @@ module io
|
|
|
|
do i = 1, ele_num
|
|
|
|
do i = 1, ele_num
|
|
|
|
do inod = 1, ng_node(lat_ele(i))
|
|
|
|
do inod = 1, ng_node(lat_ele(i))
|
|
|
|
do ibasis = 1, basisnum(lat_ele(i))
|
|
|
|
do ibasis = 1, basisnum(lat_ele(i))
|
|
|
|
write(11, '(i16, 3f23.15)') basis_type(ibasis,lat_ele(i)), r_node(:,ibasis,inod,i)
|
|
|
|
write(11, '(2i16, 3f23.15)') basis_type(ibasis,lat_ele(i)), 0, r_node(:,ibasis,inod,i)
|
|
|
|
end do
|
|
|
|
end do
|
|
|
|
end do
|
|
|
|
end do
|
|
|
|
end do
|
|
|
|
end do
|
|
|
|
|
|
|
|
|
|
|
|
!Write atom positions
|
|
|
|
!Write atom positions
|
|
|
|
do i = 1, atom_num
|
|
|
|
do i = 1, atom_num
|
|
|
|
write(11, '(i16, 3f23.15)') type_atom(i), r_atom(:,i)
|
|
|
|
write(11, '(2i16, 3f23.15)') type_atom(i), 1, r_atom(:,i)
|
|
|
|
end do
|
|
|
|
end do
|
|
|
|
|
|
|
|
|
|
|
|
!Finish writing
|
|
|
|
!Finish writing
|
|
|
|