include 'DIMENSIONS'
include 'COMMON.GEO'
include 'COMMON.CHAIN'
+ include 'COMMON.IOUNITS'
x12=c(1,i1)-c(1,i2)
x23=c(1,i3)-c(1,i2)
x34=c(1,i4)-c(1,i3)
if (dabs(scalar).gt.1.0D0)
&scalar=0.99999999999999D0*scalar/dabs(scalar)
angle=dacos(scalar)
-cd print '(2i4,10f7.3)',i2,i3,vx,vy,vz,wx,wy,wz,vnorm,wnorm,
-cd &scalar,angle
else
angle=pi
endif
scalar=tx*x23+ty*y23+tz*z23
if (scalar.lt.0.0D0) angle=-angle
beta=angle
+cd if ((vnorm.lt.0.01.or.wnorm.lt.0.01).and.angle.ne.pi
+cd & .and.i3.ne.maxres*2)
+cd & write(iout,'(a12,4i4,6f7.3,4f10.5)') 'beta warning',i1,i2,i3,i4,
+cd & vx,vy,vz,wx,wy,wz,vnorm,wnorm,scalar,angle*RAD2DEG
return
end
C