20       subroutine cmass (shp, shgl, xl, em)  argument
36      &     xl(npro,nenl,nsd)
56          call e3metric( xl, shdrv,dxidx,shg,WdetJ)
73       subroutine cmassl (shp, shgl, shpf, shglf, xl, em, Qwtf)  argument
90      &     xl(npro,nenl,nsd)
120          call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf)
155          call e3metric( xl, shdrv,dxidx,shg,WdetJ)
196       subroutine cmasstl (shp, shgl, shpf, shglf, xl, em, Qwtf)  argument
213      &     xl(npro,nenl,nsd)
237          call e3metricf( xl, shdrv,dxidx,shg,WdetJ,Qwtf)
272          call e3metric( xl, shdrv,dxidx,shg,WdetJ)
316       subroutine e3metricf(  xl,      shgl,     dxidx,  argument
321       real*8     xl(npro,nenl,nsd),    shgl(npro,nsd,nshl),  local
333           dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(:,1,n)
334           dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(:,2,n)
335           dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(:,3,n)
336           dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(:,1,n)
337           dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(:,2,n)
338           dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(:,3,n)
339           dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(:,1,n)
340           dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(:,2,n)
341           dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(:,3,n)