jgrad_start(i)=i+1
jgrad_end(i)=nres
enddo
+! THIS SHOULD BE FIXED
+ itmp=0
+ do i=1,4
+ itmp=itmp+nres_molec(i)
+ enddo
+ call int_bounds(nres_molec(5),icatb_start,icatb_end)
+ icatb_start=icatb_start+itmp
+ icatb_end=icatb_end+itmp
+
+
+
if (lprint) then
write (*,*) 'Processor:',fg_rank,' CG group',kolor,&
' absolute rank',myrank,&