-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