source/unres/src-HCD-5D/check_bond.f

   1        subroutine check_bond
   2        implicit none
   3        integer i,it
   4 C Subroutine is checking if the fitted function which describs sc_rot_pot
   5 C is correct, printing, alpha,beta, energy, data - for some known theta.
   6 C theta angle is read from the input file. Sc_rot_pot are printed
   7 C for the second  residue in sequance.
   8        include 'DIMENSIONS'
   9        include 'COMMON.VAR'
  10        include 'COMMON.GEO'
  11        include 'COMMON.INTERACT'
  12        include 'COMMON.CHAIN'
  13        double precision energia(0:n_ene)
  14        it=itype(2)
  15        do i=1,101
  16          vbld(nres+2)=0.5d0+0.05d0*(i-1)
  17          call chainbuild
  18          call etotal(energia)
  19          write (2,*) vbld(nres+2),energia(17)
  20        enddo
  21        return
  22        end