Merge branch 'bartek' of mmka.chem.univ.gda.pl:unres into bartek
[unres.git] / source / unres / src_MD / 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