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,