do j=1,3
dcostheta(j,1,i)=-(dc_norm(j,i-1)+cost*dc_norm(j,i-2))/
& vbld(i-1)
do j=1,3
dcostheta(j,1,i)=-(dc_norm(j,i-1)+cost*dc_norm(j,i-2))/
& vbld(i-1)
do j=1,3
CC Calculate derivative over first omicron (Cai-2,Cai-1,SCi-1)
dcosomicron(j,1,1,i)=-(dc_norm(j,i-1+nres)+
do j=1,3
CC Calculate derivative over first omicron (Cai-2,Cai-1,SCi-1)
dcosomicron(j,1,1,i)=-(dc_norm(j,i-1+nres)+
dcostau(j,1,1,i)=fac1*dcosomicron(j,2,2,i-1)+fac3*
& dcosomicron(j,2,2,i-1)-fac0*(dc_norm(j,i-1)-scalp*
& (dc_norm2(j,i-2+nres)))/vbld(i-2+nres)
dcostau(j,1,1,i)=fac1*dcosomicron(j,2,2,i-1)+fac3*
& dcosomicron(j,2,2,i-1)-fac0*(dc_norm(j,i-1)-scalp*
& (dc_norm2(j,i-2+nres)))/vbld(i-2+nres)
- dtauangle(j,1,1,i)=-1/sing*dcostau(j,1,1,i)
+ dtauangle(j,1,1,i)=-1.0d0/sing*dcostau(j,1,1,i)
dcostau(j,1,2,i)=fac1*dcosomicron(j,2,1,i-1)+fac2*
& dcostheta(j,1,i)+fac3*dcosomicron(j,2,1,i-1)+fac4*
& dcostheta(j,1,i)
dcostau(j,1,2,i)=fac1*dcosomicron(j,2,1,i-1)+fac2*
& dcostheta(j,1,i)+fac3*dcosomicron(j,2,1,i-1)+fac4*
& dcostheta(j,1,i)
- dtauangle(j,1,2,i)=-1/sing*dcostau(j,1,2,i)
+ dtauangle(j,1,2,i)=-1.0d0/sing*dcostau(j,1,2,i)
dcostau(j,1,3,i)=fac2*dcostheta(j,2,i)+fac4*
& dcostheta(j,2,i)-fac0*(-dc_norm(j,i-2+nres)-scalp*
& dc_norm(j,i-1))/vbld(i)
dcostau(j,1,3,i)=fac2*dcostheta(j,2,i)+fac4*
& dcostheta(j,2,i)-fac0*(-dc_norm(j,i-2+nres)-scalp*
& dc_norm(j,i-1))/vbld(i)
- dtauangle(j,1,3,i)=-1/sing*dcostau(j,1,3,i)
+ dtauangle(j,1,3,i)=-1.0d0/sing*dcostau(j,1,3,i)
dcostau(j,2,1,i)=fac1*dcostheta(j,1,i-1)+fac3*
& dcostheta(j,1,i-1)-fac0*(dc_norm(j,i-1+nres)-scalp*
& dc_norm(j,i-3))/vbld(i-2)
dcostau(j,2,1,i)=fac1*dcostheta(j,1,i-1)+fac3*
& dcostheta(j,1,i-1)-fac0*(dc_norm(j,i-1+nres)-scalp*
& dc_norm(j,i-3))/vbld(i-2)
- dtauangle(j,2,1,i)=-1/sing*dcostau(j,2,1,i)
+ dtauangle(j,2,1,i)=-1.0d0/sing*dcostau(j,2,1,i)
dcostau(j,2,2,i)=fac1*dcostheta(j,2,i-1)+fac2*
& dcosomicron(j,1,1,i)+fac3*dcostheta(j,2,i-1)+fac4*
& dcosomicron(j,1,1,i)
dcostau(j,2,2,i)=fac1*dcostheta(j,2,i-1)+fac2*
& dcosomicron(j,1,1,i)+fac3*dcostheta(j,2,i-1)+fac4*
& dcosomicron(j,1,1,i)
- dtauangle(j,2,2,i)=-1/sing*dcostau(j,2,2,i)
+ dtauangle(j,2,2,i)=-1.0d0/sing*dcostau(j,2,2,i)
dcostau(j,2,3,i)=fac2*dcosomicron(j,1,2,i)+fac4*
& dcosomicron(j,1,2,i)-fac0*(dc_norm(j,i-3)-scalp*
& dc_norm(j,i-1+nres))/vbld(i-1+nres)
dcostau(j,2,3,i)=fac2*dcosomicron(j,1,2,i)+fac4*
& dcosomicron(j,1,2,i)-fac0*(dc_norm(j,i-3)-scalp*
& dc_norm(j,i-1+nres))/vbld(i-1+nres)
- dtauangle(j,2,3,i)=-1/sing*dcostau(j,2,3,i)
+ dtauangle(j,2,3,i)=-1.0d0/sing*dcostau(j,2,3,i)
c write(iout,*) i,j,"else", dtauangle(j,2,3,i)
enddo
endif
c write(iout,*) i,j,"else", dtauangle(j,2,3,i)
enddo
endif
dcostau(j,3,1,i)=fac1*dcosomicron(j,2,2,i-1)+fac3*
& dcosomicron(j,2,2,i-1)-fac0*(dc_norm(j,i-1+nres)-scalp*
& dc_norm2(j,i-2+nres))/vbld(i-2+nres)
dcostau(j,3,1,i)=fac1*dcosomicron(j,2,2,i-1)+fac3*
& dcosomicron(j,2,2,i-1)-fac0*(dc_norm(j,i-1+nres)-scalp*
& dc_norm2(j,i-2+nres))/vbld(i-2+nres)
- dtauangle(j,3,1,i)=-1/sing*dcostau(j,3,1,i)
+ dtauangle(j,3,1,i)=-1.0d0/sing*dcostau(j,3,1,i)
dcostau(j,3,2,i)=fac1*dcosomicron(j,2,1,i-1)+fac2*
& dcosomicron(j,1,1,i)+fac3*dcosomicron(j,2,1,i-1)+fac4*
& dcosomicron(j,1,1,i)
dcostau(j,3,2,i)=fac1*dcosomicron(j,2,1,i-1)+fac2*
& dcosomicron(j,1,1,i)+fac3*dcosomicron(j,2,1,i-1)+fac4*
& dcosomicron(j,1,1,i)
- dtauangle(j,3,2,i)=-1/sing*dcostau(j,3,2,i)
+ dtauangle(j,3,2,i)=-1.0d0/sing*dcostau(j,3,2,i)
dcostau(j,3,3,i)=fac2*dcosomicron(j,1,2,i)+fac4*
& dcosomicron(j,1,2,i)-fac0*(dc_norm2(j,i-2+nres)-scalp*
& dc_norm(j,i-1+nres))/vbld(i-1+nres)
dcostau(j,3,3,i)=fac2*dcosomicron(j,1,2,i)+fac4*
& 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)
+ dtauangle(j,3,3,i)=-1.0d0/sing*dcostau(j,3,3,i)