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