link_start_homo=1
link_end_homo=lim_odl
idihconstr_start_homo=nnt+3
- idihconstr_end_homo=lim_dih+3
+ idihconstr_end_homo=lim_dih+nnt-1+3
write (iout,*)
& ' lim_odl',lim_odl,' link_start=',link_start_homo,
& ' link_end',link_end_homo,' lim_dih',lim_dih,
link_start_homo=1
link_end_homo=lim_odl
idihconstr_start_homo=nnt+3
- idihconstr_end_homo=lim_dih+3
+ idihconstr_end_homo=lim_dih+nnt-1+3
write (iout,*)
& ' lim_odl',lim_odl,' link_start=',link_start_homo,
& ' link_end',link_end_homo,' lim_dih',lim_dih,
link_start_homo=1
link_end_homo=lim_odl
idihconstr_start_homo=nnt+3
- idihconstr_end_homo=lim_dih+3
+ idihconstr_end_homo=lim_dih+nnt-1+3
write (iout,*)
& ' lim_odl',lim_odl,' link_start=',link_start_homo,
& ' link_end',link_end_homo,' lim_dih',lim_dih,
link_start_homo=1
link_end_homo=lim_odl
idihconstr_start_homo=nnt+3
- idihconstr_end_homo=lim_dih+3
+ idihconstr_end_homo=lim_dih+nnt-1+3
write (iout,*)
& ' lim_odl',lim_odl,' link_start=',link_start_homo,
& ' link_end',link_end_homo,' lim_dih',lim_dih,
link_start_homo=1
link_end_homo=lim_odl
idihconstr_start_homo=nnt+3
- idihconstr_end_homo=lim_dih+3
+ idihconstr_end_homo=lim_dih+nnt-1+3
write (iout,*)
& ' lim_odl',lim_odl,' link_start=',link_start_homo,
& ' link_end',link_end_homo,' lim_dih',lim_dih,
& " the value read in: ",energia(0),eini," point",
& iii+1,indstart(me1)+iii," T",
& 1.0d0/(1.987D-3*beta_h(ib,ipar))
+#define DEBUG
#ifdef DEBUG
call enerprint(energia(0),fT)
#endif
+#undef DEBUG
#ifdef DEBUG
write (iout,'(4f10.5,2i5)') 0.0,energia(0),0.0,
& 1.0d0/(beta_h(ib,ipar)*1.987D-3),
link_start_homo=1
link_end_homo=lim_odl
idihconstr_start_homo=nnt+3
- idihconstr_end_homo=lim_dih+3
+ idihconstr_end_homo=lim_dih+nnt-1+3
write (iout,*)
& ' lim_odl',lim_odl,' link_start=',link_start_homo,
& ' link_end',link_end_homo,' lim_dih',lim_dih,