c Obtaining the gamma derivatives from sine derivative
if (phi(i).gt.-pi4.and.phi(i).le.pi4.or.
& phi(i).gt.pi34.and.phi(i).le.pi.or.
c Obtaining the gamma derivatives from sine derivative
if (phi(i).gt.-pi4.and.phi(i).le.pi4.or.
& phi(i).gt.pi34.and.phi(i).le.pi.or.
call vecpr(dc_norm(1,i-1),dc_norm(1,i-2),vp1)
call vecpr(dc_norm(1,i-3),dc_norm(1,i-1),vp2)
call vecpr(dc_norm(1,i-3),dc_norm(1,i-2),vp3)
call vecpr(dc_norm(1,i-1),dc_norm(1,i-2),vp1)
call vecpr(dc_norm(1,i-3),dc_norm(1,i-1),vp2)
call vecpr(dc_norm(1,i-3),dc_norm(1,i-2),vp3)
scalp=scalar(dc_norm2(1,i-2+nres),dc_norm(1,i-1))
fac0=1.0d0/(sint1*sint)
fac1=cost*fac0
scalp=scalar(dc_norm2(1,i-2+nres),dc_norm(1,i-1))
fac0=1.0d0/(sint1*sint)
fac1=cost*fac0
dsintau(j,1,2,i)=
& -sing*(ctgt1*domicron(j,2,1,i-1)+ctgt*dtheta(j,1,i))
& -(fac0*vp2(j)+sing*dc_norm(j,i-2))*vbld_inv(i-1)
dsintau(j,1,2,i)=
& -sing*(ctgt1*domicron(j,2,1,i-1)+ctgt*dtheta(j,1,i))
& -(fac0*vp2(j)+sing*dc_norm(j,i-2))*vbld_inv(i-1)
-c write(iout,*) "dsintau", dsintau(j,1,2,i)
+c write(iout,*) "dsintau", dsintau(j,1,1,i),dsintau(j,1,2,i)
dtauangle(j,1,2,i)=cosg_inv*dsintau(j,1,2,i)
c Bug fixed 3/24/05 (AL)
dsintau(j,1,3,i)=-sing*ctgt*dtheta(j,2,i)
dtauangle(j,1,2,i)=cosg_inv*dsintau(j,1,2,i)
c Bug fixed 3/24/05 (AL)
dsintau(j,1,3,i)=-sing*ctgt*dtheta(j,2,i)
-c write(iout,*) "tu",i,k,(dtauangle(j,1,k,i),j=1,3)
+c write(iout,*) "tu",1,i,k,(dtauangle(j,1,k,i),j=1,3)