call int_bounds(nres-2,ithet_start,ithet_end)
ithet_start=ithet_start+2
ithet_end=ithet_end+2
+ write (iout,*) "ithet_start",ithet_start," ithet_end",ithet_end
call int_bounds(nct-nnt-2,iturn3_start,iturn3_end)
iturn3_start=iturn3_start+nnt
iphi_start=iturn3_start+2
include 'COMMON.CONTROL'
include 'COMMON.MD'
include 'COMMON.INTERACT'
- write(iout,*)"homology_partition: lim_odl=",lim_odl,
- & " lim_dih",lim_dih
+cd write(iout,*)"homology_partition: lim_odl=",lim_odl,
+cd & " lim_dih",lim_dih
#ifdef MPI
- write (iout,*) "MPI"
+ if (me.eq.king .or. .not. out1file) write (iout,*) "MPI"
call int_bounds(lim_odl,link_start_homo,link_end_homo)
call int_bounds(lim_dih-nnt+1,idihconstr_start_homo,
& idihconstr_end_homo)