& dcosomicron(j,1,2,i)-fac0*(dc_norm2(j,i-2+nres)-scalp*
& dc_norm(j,i-1+nres))/vbld(i-1+nres)
dtauangle(j,3,3,i)=-1/sing*dcostau(j,3,3,i)
& dcosomicron(j,1,2,i)-fac0*(dc_norm2(j,i-2+nres)-scalp*
& dc_norm(j,i-1+nres))/vbld(i-1+nres)
dtauangle(j,3,3,i)=-1/sing*dcostau(j,3,3,i)