call int_bounds(nct-nnt-2,iphi_start,iphi_end)
iphi_start=iphi_start+nnt+2
iphi_end=iphi_end+nnt+2
+ call int_bounds(nres-3,itau_start,itau_end)
+ itau_start=itau_start+3
+ itau_end=itau_end+3
if (lprint) then
write (iout,*) 'Processor:',MyID,
& ' loc_start',loc_start,' loc_end',loc_end,
ithet_end=nres
iphi_start=nnt+3
iphi_end=nct
+ itau_start=4
+ itau_end=nres
#endif
return
end