Commit 10/1/12

[unres.git] / source / unres / src_MIN / readpdb.F

diff --git a/source/unres/src_MIN/readpdb.F b/source/unres/src_MIN/readpdb.F
index dacd12d..1d79870 100644 (file)
--- a/source/unres/src_MIN/readpdb.F
+++ b/source/unres/src_MIN/readpdb.F
@@ -258,6 +258,9 @@ c      endif
           enddo
           iti=itype(i)
           di=dist(i,nres+i)
+C 9/29/12 Adam: Correction for zero SC-SC bond length 
+          if (itype(i).ne.10 .and. itype(i).ne.21. and. di.eq.0.0d0) 
+     &     di=dsc(itype(i))
           vbld(i+nres)=di
           if (itype(i).ne.10) then
             vbld_inv(i+nres)=1.0d0/di