--- /dev/null
+cccccccccccccccccccccccccccccccccc
+ subroutine get_diff12(aarray,barray,diff)
+ implicit real*8 (a-h,o-z)
+ include 'DIMENSIONS'
+ include 'COMMON.CSA'
+ include 'COMMON.BANK'
+ include 'COMMON.CHAIN'
+ include 'COMMON.GEO'
+ dimension aarray(mxang,maxres,mxch),
+ & barray(mxang,maxres,mxch)
+
+ diff=0.d0
+ do k=1,numch
+ do j=2,nres-1
+c do i=1,4
+c do i=1,2
+ do i=1,ndiff
+ dif=rad2deg*dabs(aarray(i,j,k)-barray(i,j,k))
+ if(dif.gt.180.) dif=360.-dif
+ if (dif.gt.diffcut) diff=diff+dif
+ enddo
+ enddo
+ enddo
+
+ return
+ end
+ccccccccccccccccccccccccccccccccccccccccccccccccc