40       subroutine bardmc (y,      shgl,      shp,   argument
71      &          shgl(MAXTOP,nsd,maxsh,MAXQPT), shp(MAXTOP,maxsh,MAXQPT),    
250      &               fres, hfres, shgl(lcsyst,:,1:nshl,:),
269      &               strl(iel:inum,:), shgl(lcsyst,:,1:nshl,:),
672       subroutine hfilter (y, x, ien, hfres, shgl, shp, Qwtf)  argument
683      &          shgl(nsd,nshl,maxsh),       shg(npro,nshl,nsd),
713             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
714             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
715             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
716             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
717             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
718             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
719             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
720             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
721             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
826      &     hfres, shgl, shp, Qwtf)  argument
837      &          shgl(nsd,nshl,maxsh),       shg(npro,nshl,nsd),
869             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
870             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
871             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
872             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
873             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
874             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
875             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
876             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
877             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
945         shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
946      &              + shgl(2,n,intp) * dxidx(:,2,1)
947      &              + shgl(3,n,intp) * dxidx(:,3,1))
948         shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
949      &              + shgl(2,n,intp) * dxidx(:,2,2)
950      &              + shgl(3,n,intp) * dxidx(:,3,2))
951         shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
952      &              + shgl(2,n,intp) * dxidx(:,2,3)
953      &              + shgl(3,n,intp) * dxidx(:,3,3))
1060       subroutine disPtion (y, x, ien, shgl, shp, eps, xmudmi)   argument
1069      &          shgl(nsd,nshl,maxsh),       shp(nshl,maxsh),
1102             dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
1103             dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
1104             dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
1105             dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
1106             dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
1107             dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
1108             dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
1109             dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
1110             dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
1146         shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
1147      &              + shgl(2,n,intp) * dxidx(:,2,1)
1148      &              + shgl(3,n,intp) * dxidx(:,3,1))
1149         shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
1150      &              + shgl(2,n,intp) * dxidx(:,2,2)
1151      &              + shgl(3,n,intp) * dxidx(:,3,2))
1152         shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
1153      &              + shgl(2,n,intp) * dxidx(:,2,3)
1154      &              + shgl(3,n,intp) * dxidx(:,3,3))