524      &          nsons(nshg),                   iBC(nshg),
942      &                   nsons,  ifath,     x, cdelsq1)  argument
963      &          nsons(nshg),
1272      &                  nsons,       ifath,      x,  argument
1296      &          ifath(nshg),                   nsons(nshg),
1552             xavegr(:,m) = xavegr(:,m)/nsons(:)
1568             xaveg(:,m) = xaveg(:,m)/nsons(:)
1607      &                  nsons,       ifath,      x,  argument
1631      &          ifath(nshg),                   nsons(nshg),
1988             xavegr(:,m) = xavegr(:,m)/nsons(:)
2002             xaveg(:,m) = xaveg(:,m)/nsons(:)
2110      &                   nsons,  ifath,     x,   cdelsq1, fwr4,   argument
2132      &          nsons(nshg),           xfac(nshg,5),
2460      &                   nsons,  ifath,     x, fwr2, fwr3)  argument
2491      &          nsons(nshg),
2832          xnuder(:,1) = xnuder(:,1) / nsons(:)
2833          xnuder(:,2) = xnuder(:,2) / nsons(:)
2836          ynuder(:,m) = ynuder(:,m)/nsons(:)
2839          xmr(:,m) = xmr(:,m)/nsons(:)
2840          xlr(:,m) = xlr(:,m)/nsons(:)
2841          xl1r(:,m) = xl1r(:,m)/nsons(:)
2842          xl2r(:,m) = xl2r(:,m)/nsons(:)
2845          uir(:,1) = uir(:,1)/nsons(:)
2846          uir(:,2) = uir(:,2)/nsons(:)
2847          uir(:,3) = uir(:,3)/nsons(:)
2849          snormr(:) = snormr(:)/nsons(:)
2885          xnude(:,1) = xnude(:,1)/nsons(:)
2886          xnude(:,2) = xnude(:,2)/nsons(:)
2889          ynude(:,m) = ynude(:,m)/nsons(:)
2892          xm(:,m) = xm(:,m)/nsons(:)
2893          xl(:,m) = xl(:,m)/nsons(:)
2894          xl1(:,m) = xl1(:,m)/nsons(:)
2895          xl2(:,m) = xl2(:,m)/nsons(:)
2898          ui(:,1) = ui(:,1)/nsons(:)
2899          ui(:,2) = ui(:,2)/nsons(:)
2900          ui(:,3) = ui(:,3)/nsons(:)
2902          snorm(:) = snorm(:)/nsons(:)
3114          xnutf=xnuder(:,1)*xnuder(:,2)/nsons(:)
3126          xnuder(:,1)=xnuder(:,1)/nsons(:)
3127          xnuder(:,2)=xnuder(:,2)/nsons(:)
3143      &                   nsons,  ifath,     x,    fwr2, fwr3)  argument
3172      &          nsons(nshg),
3578          xnuder(:,1) = xnuder(:,1) / nsons(:)
3579          xnuder(:,2) = xnuder(:,2) / nsons(:)
3582          ynuder(:,m) = ynuder(:,m)/nsons(:)
3585          xmr(:,m) = xmr(:,m)/nsons(:)
3586          xlr(:,m) = xlr(:,m)/nsons(:)
3587          xl1r(:,m) = xl1r(:,m)/nsons(:)
3588          xl2r(:,m) = xl2r(:,m)/nsons(:)
3591          uir(:,1) = uir(:,1)/nsons(:)
3592          uir(:,2) = uir(:,2)/nsons(:)
3593          uir(:,3) = uir(:,3)/nsons(:)
3595          snormr(:) = snormr(:)/nsons(:)
3632          xnude(:,1) = xnude(:,1)/nsons(:)
3633          xnude(:,2) = xnude(:,2)/nsons(:)
3636          ynude(:,m) = ynude(:,m)/nsons(:)
3639          xm(:,m) = xm(:,m)/nsons(:)
3640          xl(:,m) = xl(:,m)/nsons(:)
3641          xl1(:,m) = xl1(:,m)/nsons(:)
3642          xl2(:,m) = xl2(:,m)/nsons(:)
3645          ui(:,1) = ui(:,1)/nsons(:)
3646          ui(:,2) = ui(:,2)/nsons(:)
3647          ui(:,3) = ui(:,3)/nsons(:)
3649          snorm(:) = snorm(:)/nsons(:)
3871          xnutf=xnuder(:,1)*xnuder(:,2)/nsons(:)
3883          xnuder(:,1)=xnuder(:,1)/nsons(:)
3884          xnuder(:,2)=xnuder(:,2)/nsons(:)