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