source/unres/src_MD/src/check_bond.f

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