973      &                  nsons,       ifath,      x,  argument
997      &          ifath(nshg),                   nsons(nshg),
1013       rnum(ifath(:)) = numNden(:,1)
1014       rden(ifath(:)) = numNden(:,2)
1178          xaveg(ifath(i),:) = xaveg(ifath(i),:) + xave(i,:)
1641      &                   nsons,  ifath,     x)  argument
1675      &          ifath(nshg),          iper(nshg),
1939          xnude(ifath(i),1) = xnude(ifath(i),1) + xnum(i)
1940          xnude(ifath(i),2) = xnude(ifath(i),2) + xden(i)
1942          ynude(ifath(i),1) = ynude(ifath(i),1) + xfac(i,1)
1943          ynude(ifath(i),2) = ynude(ifath(i),2) + xfac(i,2)
1944          ynude(ifath(i),3) = ynude(ifath(i),3) + xfac(i,3)
1945          ynude(ifath(i),4) = ynude(ifath(i),4) + xfac(i,4)
1946          ynude(ifath(i),5) = ynude(ifath(i),5) + xfac(i,5)
1948          xm(ifath(i),1) = xm(ifath(i),1) + xmij(i,1)
1949          xm(ifath(i),2) = xm(ifath(i),2) + xmij(i,2)
1950          xm(ifath(i),3) = xm(ifath(i),3) + xmij(i,3)
1951          xm(ifath(i),4) = xm(ifath(i),4) + xmij(i,4)
1952          xm(ifath(i),5) = xm(ifath(i),5) + xmij(i,5)
1953          xm(ifath(i),6) = xm(ifath(i),6) + xmij(i,6)
1955          xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
1956          xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
1957          xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
1958          xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
1959          xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
1960          xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)         
1962          xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
1963          xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
1964          xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
1965          xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
1966          xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
1967          xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)         
1969          xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1) 
1970          xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
1971          xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3) 
1972          xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
1973          xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
1974          xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
1976          ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
1977          ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
1978          ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
1980          snorm(ifath(i)) = snorm(ifath(i)) + strnrm(i)
2004          cdelsq2(:) = ynuder(ifath(:),6)  ! For comparison w/ cdelsq
2043             numNden(:,1) = whist*numNden(:,1)+wcur*xnuder(ifath(:),1)
2044             numNden(:,2) = whist*numNden(:,2)+wcur*xnuder(ifath(:),2)
2106             numNden(:,1) = whist*numNden(:,1)+wcur*xnude(ifath(:),1)
2107             numNden(:,2) = whist*numNden(:,2)+wcur*xnude(ifath(:),2)
2117           cdelsq2(:) = ynude(ifath(:),6)  ! For comparison w/ cdelsq
2288             xnude(ifath(i),2) = xnude(ifath(i),2) + strnrm(i)
2305             xnude(ifath(i),1) = xnude(ifath(i),1) + x(i,1)
2306             xnude(ifath(i),2) = xnude(ifath(i),2) + x(i,2)
2328      &                   nsons,  ifath,     x)  argument
2360      &          ifath(nshg),          iper(nshg),
2680          xnude(ifath(i),1) = xnude(ifath(i),1) + xnum(i)
2681          xnude(ifath(i),2) = xnude(ifath(i),2) + xden(i)
2683          ynude(ifath(i),1) = ynude(ifath(i),1) + xfac(i,1)
2684          ynude(ifath(i),2) = ynude(ifath(i),2) + xfac(i,2)
2685          ynude(ifath(i),3) = ynude(ifath(i),3) + xfac(i,3)
2686          ynude(ifath(i),4) = ynude(ifath(i),4) + xfac(i,4)
2687          ynude(ifath(i),5) = ynude(ifath(i),5) + xfac(i,5)
2689          xm(ifath(i),1) = xm(ifath(i),1) + xmij(i,1)
2690          xm(ifath(i),2) = xm(ifath(i),2) + xmij(i,2)
2691          xm(ifath(i),3) = xm(ifath(i),3) + xmij(i,3)
2692          xm(ifath(i),4) = xm(ifath(i),4) + xmij(i,4)
2693          xm(ifath(i),5) = xm(ifath(i),5) + xmij(i,5)
2694          xm(ifath(i),6) = xm(ifath(i),6) + xmij(i,6)
2696          xl(ifath(i),1) = xl(ifath(i),1) + xlij(i,1)
2697          xl(ifath(i),2) = xl(ifath(i),2) + xlij(i,2)
2698          xl(ifath(i),3) = xl(ifath(i),3) + xlij(i,3)
2699          xl(ifath(i),4) = xl(ifath(i),4) + xlij(i,4)
2700          xl(ifath(i),5) = xl(ifath(i),5) + xlij(i,5)
2701          xl(ifath(i),6) = xl(ifath(i),6) + xlij(i,6)         
2703          xl1(ifath(i),1) = xl1(ifath(i),1) + fres(i,4)
2704          xl1(ifath(i),2) = xl1(ifath(i),2) + fres(i,5)
2705          xl1(ifath(i),3) = xl1(ifath(i),3) + fres(i,6)
2706          xl1(ifath(i),4) = xl1(ifath(i),4) + fres(i,7)
2707          xl1(ifath(i),5) = xl1(ifath(i),5) + fres(i,8)
2708          xl1(ifath(i),6) = xl1(ifath(i),6) + fres(i,9)         
2710          xl2(ifath(i),1) = xl2(ifath(i),1) + fres(i,1)*fres(i,1) 
2711          xl2(ifath(i),2) = xl2(ifath(i),2) + fres(i,2)*fres(i,2)
2712          xl2(ifath(i),3) = xl2(ifath(i),3) + fres(i,3)*fres(i,3) 
2713          xl2(ifath(i),4) = xl2(ifath(i),4) + fres(i,1)*fres(i,2)
2714          xl2(ifath(i),5) = xl2(ifath(i),5) + fres(i,1)*fres(i,3)
2715          xl2(ifath(i),6) = xl2(ifath(i),6) + fres(i,2)*fres(i,3)
2717          ui(ifath(i),1) = ui(ifath(i),1) + fres(i,1)
2718          ui(ifath(i),2) = ui(ifath(i),2) + fres(i,2)
2719          ui(ifath(i),3) = ui(ifath(i),3) + fres(i,3)
2721          snorm(ifath(i)) = snorm(ifath(i)) + strnrm(i)
2745          cdelsq2(:) = ynuder(ifath(:),6)  ! For comparison w/ cdelsq
2785             numNden(:,1) = whist*numNden(:,1)+wcur*xnuder(ifath(:),1)
2786             numNden(:,2) = whist*numNden(:,2)+wcur*xnuder(ifath(:),2)
2847             numNden(:,1) = whist*numNden(:,1)+wcur*xnude(ifath(:),1)
2848             numNden(:,2) = whist*numNden(:,2)+wcur*xnude(ifath(:),2)
2858           cdelsq2(:) = ynude(ifath(:),6)  ! For comparison w/ cdelsq
3040             xnude(ifath(i),2) = xnude(ifath(i),2) + strnrm(i)
3057             xnude(ifath(i),1) = xnude(ifath(i),1) + x(i,1)
3058             xnude(ifath(i),2) = xnude(ifath(i),2) + x(i,2)