source/cluster/wham/src-M-SAXS/rmsnat.f

   1       double precision function rmsnat(jcon)
   2       implicit real*8 (a-h,o-z)
   3       include 'DIMENSIONS'
   4       include 'sizesclu.dat'
   5       include 'DIMENSIONS.COMPAR'
   6       include 'COMMON.IOUNITS'
   7       include 'COMMON.CHAIN'
   8       include 'COMMON.INTERACT'
   9       include 'COMMON.VAR'
  10       include 'COMMON.CONTROL'
  11       integer ipermmin
  12       double precision przes(3),obrot(3,3)
  13       rmsnat=rmscalc(c(1,1),cref_pdb(1,1),przes,obrot,ipermmin)
  14       return
  15       end
  16 c-----------------------------------------------------------------------------
  17       double precision function gyrate()
  18       implicit real*8 (a-h,o-z)
  19       include 'DIMENSIONS'
  20       include 'COMMON.INTERACT'
  21       include 'COMMON.CHAIN'
  22       double precision cen(3),rg
  23
  24       do j=1,3
  25        cen(j)=0.0d0
  26       enddo
  27
  28       ii=0
  29       do i=nnt,nct
  30         if (itype(i).eq.ntyp1) cycle
  31         ii=ii+1
  32         do j=1,3
  33           cen(j)=cen(j)+c(j,i)
  34         enddo
  35       enddo
  36       do j=1,3
  37         cen(j)=cen(j)/dble(ii)
  38       enddo
  39       rg = 0.0d0
  40       do i = nnt, nct
  41         if (itype(i).eq.ntyp1) cycle
  42         do j=1,3
  43          rg = rg + (c(j,i)-cen(j))**2
  44         enddo
  45       end do
  46       gyrate = dsqrt(rg/dble(ii))
  47       return
  48       end