the first working version of multichain homology
[unres.git] / source / wham / src-M / promienie.f
index 12a2e80..d04927d 100644 (file)
@@ -1,6 +1,7 @@
       subroutine promienie(*)
       implicit none
       include 'DIMENSIONS'
+      include 'DIMENSIONS.FREE'
       include 'COMMON.CONTROL'
       include 'COMMON.INTERACT'
       include 'COMMON.IOUNITS'
@@ -36,7 +37,7 @@
       enddo
       close (isidep1)
       do i=1,ntyp1
-        if (i.eq.10 .or. i.eq.21) then
+        if (i.eq.10 .or. i.eq.ntyp1) then
           dsc_inv(i)=0.0d0
         else
           dsc_inv(i)=1.0d0/dsc(i)