1140 continue
1130 continue
- xmeana=xmeana/IHNT2(1)
- ymeana=ymeana/IHNT2(1)
- xmeanb=xmeanb/IHNT2(3)
- ymeanb=ymeanb/IHNT2(3)
- xmeanp=xmeanp/npart
- ymeanp=ymeanp/npart
- xm2=xm2/npart
- ym2=ym2/npart
- xym=xym/npart
+ IF (npart.GT.0) THEN
+ xmeana=xmeana/IHNT2(1)
+ ymeana=ymeana/IHNT2(1)
+ xmeanb=xmeanb/IHNT2(3)
+ ymeanb=ymeanb/IHNT2(3)
+ xmeanp=xmeanp/npart
+ ymeanp=ymeanp/npart
+ xm2=xm2/npart
+ ym2=ym2/npart
+ xym=xym/npart
- sx2=xm2-xmeanp*xmeanp
- sy2=ym2-ymeanp*ymeanp
- sxy=xym-xmeanp*ymeanp
+ sx2=xm2-xmeanp*xmeanp
+ sy2=ym2-ymeanp*ymeanp
+ sxy=xym-xmeanp*ymeanp
- delx=xmeanb-xmeana
- dely=ymeanb-ymeana
- dtmp=delx**2+dely**2
- bbtrue=sqrt(dtmp)
- dnumt=(sy2-sx2)*(delx**2-dely**2)-4D0*sxy*delx*dely
- ddent=(sy2+sx2)*bbtrue**2
- eccrp=dnumt/ddent
- dtmp=(sy2-sx2)*(sy2-sx2)+4D0*sxy*sxy
- eccpart=sqrt(dtmp)/(sx2+sy2)
- eccmc=(sy2-sx2)/(sy2+sx2)
- write(*,*),'HOUT: ',bb,' ',bbtrue,' ',ncolt,' ',npart,
- 1 ' ',eccrp,' ',eccpart
+ delx=xmeanb-xmeana
+ dely=ymeanb-ymeana
+ dtmp=delx**2+dely**2
+ bbtrue=sqrt(dtmp)
+ dnumt=(sy2-sx2)*(delx**2-dely**2)-4D0*sxy*delx*dely
+ ddent=(sy2+sx2)*bbtrue**2
+ eccrp=dnumt/ddent
+ dtmp=(sy2-sx2)*(sy2-sx2)+4D0*sxy*sxy
+ eccpart=sqrt(dtmp)/(sx2+sy2)
+ eccmc=(sy2-sx2)/(sy2+sx2)
+ write(*,*),'HOUT: ',bb,' ',bbtrue,' ',ncolt,' ',npart,
+ 1 ' ',eccrp,' ',eccpart
+ end if
end if
C ********total number interactions proj and targ has