include 'COMMON.GEO'
include 'COMMON.CHAIN'
include 'COMMON.VAR'
- include 'COMMON.REFSYS'
include 'COMMON.IOUNITS'
include 'COMMON.MCM'
- dimension x(3),e(3,3),rot(3,3),trans(3,3)
+ dimension x(3),e(3,3),e1(3),e2(3),e3(3),rot(3,3),trans(3,3)
error=.false.
nend=nstart+nbond
if (print_mc.gt.2) then
i2=nend
i3=nstart
i4=nstart+1
- call refsys(error)
+ call refsys(i2,i3,i4,e1,e2,e3,error)
C Return, if couldn't define the reference system.
if (error) return
C Compute the transformation matrix.