22
22
use fdf, only : fdf_string
23
23
use units, only : Ang
24
24
use m_mpi_utils, only : broadcast
25
use atomlist, only : xa, isa, cisa
25
use siesta_geom, only : xa, isa, cisa
26
26
use alloc, only : re_alloc
27
27
use sys, only : die
28
28
use files, only : slabel, label_length
89
89
call broadcast(xa(1:3,1:na))
91
91
! Construct references
93
call re_alloc(cisa,1,na,name="cisa",routine="read_struct")
94
95
write(cisa(ia), '("siesta:e",i3.3)') isa(ia)
97
98
end subroutine read_struct
98
!--------------------------------------------------------------------------
99
subroutine write_struct(cell, na, isa, iza, xa, moved )
99
!---------------------------------------------------------------------
100
subroutine write_struct(cell, na, isa, iza, xa, moved)
101
102
c real*8 cell(3,3) : Unit cell vectors
102
103
c integer na : Number of atoms
103
104
c integer isa(na) : Atomic species index
104
105
c integer iza(na) : Atomic numbers
105
106
c real*8 xa(3,na) : Atomic positions
106
c logical moved : True if structure is the "predicted"
107
c logical moved : True if structure is the "predicted"
107
108
c one after application of forces/stress.
109
110
integer, intent(in) :: na, isa(na), iza(na)