X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fintcor.f;h=9195f8a5cc4b5ab51fabfd2a8cb0083d74e097bc;hb=bf8b3ab324569ca7ae7bb9045ad9358f6826a150;hp=a3cd5d0bbbb65be53f175a1ef35f311524a4a959;hpb=6e5f3f4866eb982f17fb6d5c2bb30bd7bf0585be;p=unres.git diff --git a/source/unres/src_MD-M/intcor.f b/source/unres/src_MD-M/intcor.f index a3cd5d0..9195f8a 100644 --- a/source/unres/src_MD-M/intcor.f +++ b/source/unres/src_MD-M/intcor.f @@ -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