Merge branch 'devel' into AFM
[unres.git] / source / unres / src_MD-M / intcor.f
index a3cd5d0..9195f8a 100644 (file)
@@ -17,7 +17,11 @@ c
       z23=c(3,i3)-c(3,i2)
       vnorm=dsqrt(x12*x12+y12*y12+z12*z12)
       wnorm=dsqrt(x23*x23+y23*y23+z23*z23)
+      if ((vnorm.eq.0.0).or.(wnorm.eq.0.0)) then
+      scalar=1.0
+      else
       scalar=(x12*x23+y12*y23+z12*z23)/(vnorm*wnorm)
+      endif
       alpha=arcos(scalar)
       return
       end