Calculate derivative of Tauangle
#ifdef PARINTDER
- do i=iphi1_start,iphi1_end
+ do i=iphi1_start-1,iphi1_end
#else
- do i=4,nres
+ do i=3,nres
#endif
+cc dtauangle(j,intertyp,dervityp,residue number)
cc INTERTYP=1 SC...Ca...Ca..Ca
c the conventional case
sint=dsin(theta(i))
CCC third case SC...Ca...Ca...SC
#ifdef PARINTDER
- do i=iphi1_start,iphi1_end
+ do i=iphi1_start-1,iphi1_end
#else
- do i=4,nres
+ do i=3,nres
#endif
c the conventional case
sint=dsin(omicron(1,i))