X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M-newcorr%2Fcheck_bond.f;fp=source%2Funres%2Fsrc_MD-M-newcorr%2Fcheck_bond.f;h=c8a4ad1d78ba71aaf48e3878f7269f26ce3688f3;hb=7308760ff07636ef6b1ee28d8c3a67a23c14b34b;hp=0000000000000000000000000000000000000000;hpb=9a54ab407f6d0d9d564d52763b3e2136450b9ffc;p=unres.git

diff --git a/source/unres/src_MD-M-newcorr/check_bond.f b/source/unres/src_MD-M-newcorr/check_bond.f
new file mode 100644
index 0000000..c8a4ad1
--- /dev/null
+++ b/source/unres/src_MD-M-newcorr/check_bond.f
@@ -0,0 +1,20 @@
+       subroutine check_bond
+C Subroutine is checking if the fitted function which describs sc_rot_pot
+C is correct, printing, alpha,beta, energy, data - for some known theta. 
+C theta angle is read from the input file. Sc_rot_pot are printed 
+C for the second  residue in sequance.
+       include 'DIMENSIONS'
+       include 'COMMON.VAR'
+       include 'COMMON.GEO'
+       include 'COMMON.INTERACT'
+       include 'COMMON.CHAIN'
+       double precision energia(0:n_ene)
+       it=itype(2)
+       do i=1,101
+         vbld(nres+2)=0.5d0+0.05d0*(i-1)
+         call chainbuild
+         call etotal(energia)
+         write (2,*) vbld(nres+2),energia(17)
+       enddo
+       return
+       end