Lines Matching refs:intp
84 do intp = 1, ngaussf
95 pres(:) = pres(:) + shp(n,intp) * yl(:,n,1)
96 u1(:) = u1(:) + shp(n,intp) * yl(:,n,2)
97 u2(:) = u2(:) + shp(n,intp) * yl(:,n,3)
98 u3(:) = u3(:) + shp(n,intp) * yl(:,n,4)
106 aci(:,1) = aci(:,1) + shp(n,intp) * acl(:,n,2)
107 aci(:,2) = aci(:,2) + shp(n,intp) * acl(:,n,3)
108 aci(:,3) = aci(:,3) + shp(n,intp) * acl(:,n,4)
120 dxdxi(:,1,1) = dxdxi(:,1,1) + xl(:,n,1) * shgl(1,n,intp)
121 dxdxi(:,1,2) = dxdxi(:,1,2) + xl(:,n,1) * shgl(2,n,intp)
122 dxdxi(:,1,3) = dxdxi(:,1,3) + xl(:,n,1) * shgl(3,n,intp)
123 dxdxi(:,2,1) = dxdxi(:,2,1) + xl(:,n,2) * shgl(1,n,intp)
124 dxdxi(:,2,2) = dxdxi(:,2,2) + xl(:,n,2) * shgl(2,n,intp)
125 dxdxi(:,2,3) = dxdxi(:,2,3) + xl(:,n,2) * shgl(3,n,intp)
126 dxdxi(:,3,1) = dxdxi(:,3,1) + xl(:,n,3) * shgl(1,n,intp)
127 dxdxi(:,3,2) = dxdxi(:,3,2) + xl(:,n,3) * shgl(2,n,intp)
128 dxdxi(:,3,3) = dxdxi(:,3,3) + xl(:,n,3) * shgl(3,n,intp)
159 wght=Qwtf(intp)
165 shg(:,n,1) = (shgl(1,n,intp) * dxidx(:,1,1)
166 & + shgl(2,n,intp) * dxidx(:,2,1)
167 & + shgl(3,n,intp) * dxidx(:,3,1))
168 shg(:,n,2) = (shgl(1,n,intp) * dxidx(:,1,2)
169 & + shgl(2,n,intp) * dxidx(:,2,2)
170 & + shgl(3,n,intp) * dxidx(:,3,2))
171 shg(:,n,3) = (shgl(1,n,intp) * dxidx(:,1,3)
172 & + shgl(2,n,intp) * dxidx(:,2,3)
173 & + shgl(3,n,intp) * dxidx(:,3,3))
218 xmudmif(:,intp) = xmudmif(:,intp) +
219 & cdelsq(ien(:,n)) * Snorm2(:)*shp(n,intp)
223 xmudmif(:,intp)=min(xmudmif(:,intp),1000.0*rmu(:)) !
225 xmudmif(:,intp)=max(xmudmif(:,intp), zero) ! don't let (xmudmi) < 0
237 shpfun(i,n) = shp(n,intp)
284 Temp = Temp + shp(n,intp) * yl(:,n,5)
293 xx(:,1) = xx(:,1) + shp(n,intp) * xl(:,n,1)
294 xx(:,2) = xx(:,2) + shp(n,intp) * xl(:,n,2)
295 xx(:,3) = xx(:,3) + shp(n,intp) * xl(:,n,3)
321 xx(:,1) = xx(:,1) + shp(n,intp) * xl(:,n,1)
322 xx(:,2) = xx(:,2) + shp(n,intp) * xl(:,n,2)
323 xx(:,3) = xx(:,3) + shp(n,intp) * xl(:,n,3)
437 dissi(:,2) = xmudmif(:,intp)*Snorm(:)*rho(:)*WdetJ(:)