40       lfathe = 0
75       subroutine getnu (ien, strl, xmudmi, cdelsq, lfathe)  argument
81      &           lfathe(npro,ngauss)         
86          xmudmi(iel,intp) = cdelsq( lfathe(iel,intp) )*
208       subroutine getylist (ylist, lfathe, shp, x, ien)  argument
212       dimension lfathe(npro,maxsh),    shp(nshl,maxsh), 
240             lfathe(nel,intp) = nlist
247                   lfathe(nel,intp) = imatch
254             lfathe(nel,intp) = nlist
448       subroutine smagcoeff (fres, xden, xnum, ien, shp, lfathe)  argument
459      &          lfathe(npro,maxsh)
545          xden( lfathe(nel,intp) ) = xden( lfathe(nel,intp) ) +
548          xnum( lfathe(nel,intp) ) = xnum( lfathe(nel,intp) ) +