implicit real*8 (a-h,o-z)
include 'DIMENSIONS'
include 'DIMENSIONS.ZSCOPT'
+ include 'DIMENSIONS.FREE'
#ifdef MPI
include 'mpif.h'
#endif
& " lim_dih",lim_dih
#ifdef MPL
call int_bounds(lim_odl,link_start_homo,link_end_homo)
- call int_bounds(lim_dih-nnt+1,idihconstr_start_homo,
+ call int_bounds(lim_dih,idihconstr_start_homo,
& idihconstr_end_homo)
- idihconstr_start_homo=idihconstr_start_homo+nnt-1
- idihconstr_end_homo=idihconstr_end_homo+nnt-1
+ idihconstr_start_homo=idihconstr_start_homo+nnt-1+3
+ idihconstr_end_homo=idihconstr_end_homo+nnt-1+3
if (me.eq.king .or. .not. out1file)
& write (iout,*) 'Processor',fg_rank,' CG group',kolor,
& ' absolute rank',MyRank,
#else
link_start_homo=1
link_end_homo=lim_odl
- idihconstr_start_homo=nnt
- idihconstr_end_homo=lim_dih
+ idihconstr_start_homo=nnt+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,