+c do i=2,nres
+c vbld(i)=vbl
+c vbld_inv(i)=vblinv
+c enddo
+c if (itype(1).eq.ntyp1) then
+c vbld(2)=vbld(2)/2
+c vbld_inv(2)=vbld_inv(2)*2
+c endif
+c if (itype(nres).eq.ntyp1) then
+c vbld(nres)=vbld(nres)/2
+c vbld_inv(nres)=vbld_inv(nres)*2
+c endif
+c do i=2,nres-1
+c vbld(i+nres)=dsc(iabs(itype(i)))
+c vbld_inv(i+nres)=dsc_inv(iabs(itype(i)))