source/unres/check_bond.F90

   1       module check_bond_
   2
   3       use names
   4       use geometry_data
   5       use energy_data
   6       use geometry, only: chainbuild
   7       use energy, only: etotal
   8       implicit none
   9 !-----------------------------------------------------------------------------
  10       contains
  11 !-----------------------------------------------------------------------------
  12 ! check_bond.f
  13 !-----------------------------------------------------------------------------
  14       subroutine check_bond
  15 ! Subroutine is checking if the fitted function which describs sc_rot_pot
  16 ! is correct, printing, alpha,beta, energy, data - for some known theta.
  17 ! theta angle is read from the input file. Sc_rot_pot are printed
  18 ! for the second  residue in sequance.
  19 !       include 'DIMENSIONS'
  20 !       include 'COMMON.VAR'
  21 !       include 'COMMON.GEO'
  22 !       include 'COMMON.INTERACT'
  23 !       include 'COMMON.CHAIN'
  24        real(kind=8),dimension(0:n_ene) :: energia
  25 !el local variables
  26        integer :: it,i
  27
  28        it=itype(2)
  29        do i=1,101
  30          vbld(nres+2)=0.5d0+0.05d0*(i-1)
  31          call chainbuild
  32          call etotal(energia)
  33          write (2,*) vbld(nres+2),energia(17)
  34        enddo
  35       return
  36       end subroutine check_bond
  37
  38       end module check_bond_