rr0(i)=0.0D0
a0thet(i)=0.0D0
do j=1,2
- athet(j,i)=0.0D0
- bthet(j,i)=0.0D0
+ do k=-1,1
+ do kk=-1,1
+ athet(j,i,k,kk)=0.0D0
+ bthet(j,i,k,kk)=0.0D0
+ enddo
+ enddo
enddo
do j=0,3
polthet(j,i)=0.0D0
include 'COMMON.IOUNITS'
include 'COMMON.SETUP'
include 'COMMON.CONTROL'
+c write(2,*)"hpb_partition: nhpb=",nhpb
#ifdef MPI
call int_bounds(nhpb,link_start,link_end)
if (.not. out1file)
link_start=1
link_end=nhpb
#endif
+c write(2,*)"hpb_partition: link_start=",nhpb," link_end=",link_end
return
end