973      &                  nsons,       ifath,      x,  argument
997      &          ifath(nshg),                   nsons(nshg),
1185             xavegr(:,m) = xavegr(:,m)/nsons(:)
1199             xaveg(:,m) = xaveg(:,m)/nsons(:)
1641      &                   nsons,  ifath,     x)  argument
1672      &          nsons(nshg),
2017          xnuder(:,1) = xnuder(:,1) / nsons(:)
2018          xnuder(:,2) = xnuder(:,2) / nsons(:)
2021          ynuder(:,m) = ynuder(:,m)/nsons(:)
2024          xmr(:,m) = xmr(:,m)/nsons(:)
2025          xlr(:,m) = xlr(:,m)/nsons(:)
2026          xl1r(:,m) = xl1r(:,m)/nsons(:)
2027          xl2r(:,m) = xl2r(:,m)/nsons(:)
2030          uir(:,1) = uir(:,1)/nsons(:)
2031          uir(:,2) = uir(:,2)/nsons(:)
2032          uir(:,3) = uir(:,3)/nsons(:)
2034          snormr(:) = snormr(:)/nsons(:)
2070          xnude(:,1) = xnude(:,1)/nsons(:)
2071          xnude(:,2) = xnude(:,2)/nsons(:)
2074          ynude(:,m) = ynude(:,m)/nsons(:)
2077          xm(:,m) = xm(:,m)/nsons(:)
2078          xl(:,m) = xl(:,m)/nsons(:)
2079          xl1(:,m) = xl1(:,m)/nsons(:)
2080          xl2(:,m) = xl2(:,m)/nsons(:)
2083          ui(:,1) = ui(:,1)/nsons(:)
2084          ui(:,2) = ui(:,2)/nsons(:)
2085          ui(:,3) = ui(:,3)/nsons(:)
2087          snorm(:) = snorm(:)/nsons(:)
2299          xnutf=xnuder(:,1)*xnuder(:,2)/nsons(:)
2311          xnuder(:,1)=xnuder(:,1)/nsons(:)
2312          xnuder(:,2)=xnuder(:,2)/nsons(:)
2328      &                   nsons,  ifath,     x)  argument
2357      &          nsons(nshg),
2758          xnuder(:,1) = xnuder(:,1) / nsons(:)
2759          xnuder(:,2) = xnuder(:,2) / nsons(:)
2762          ynuder(:,m) = ynuder(:,m)/nsons(:)
2765          xmr(:,m) = xmr(:,m)/nsons(:)
2766          xlr(:,m) = xlr(:,m)/nsons(:)
2767          xl1r(:,m) = xl1r(:,m)/nsons(:)
2768          xl2r(:,m) = xl2r(:,m)/nsons(:)
2771          uir(:,1) = uir(:,1)/nsons(:)
2772          uir(:,2) = uir(:,2)/nsons(:)
2773          uir(:,3) = uir(:,3)/nsons(:)
2775          snormr(:) = snormr(:)/nsons(:)
2812          xnude(:,1) = xnude(:,1)/nsons(:)
2813          xnude(:,2) = xnude(:,2)/nsons(:)
2816          ynude(:,m) = ynude(:,m)/nsons(:)
2819          xm(:,m) = xm(:,m)/nsons(:)
2820          xl(:,m) = xl(:,m)/nsons(:)
2821          xl1(:,m) = xl1(:,m)/nsons(:)
2822          xl2(:,m) = xl2(:,m)/nsons(:)
2825          ui(:,1) = ui(:,1)/nsons(:)
2826          ui(:,2) = ui(:,2)/nsons(:)
2827          ui(:,3) = ui(:,3)/nsons(:)
2829          snorm(:) = snorm(:)/nsons(:)
3051          xnutf=xnuder(:,1)*xnuder(:,2)/nsons(:)
3063          xnuder(:,1)=xnuder(:,1)/nsons(:)
3064          xnuder(:,2)=xnuder(:,2)/nsons(:)