| /phasta/phSolver/incompressible/ |
| H A D | e3res.f | 4 & rLui, rmu, rho, argument
44 & rLui(npro,nsd), rmu(npro), rho(npro),
93 rNa(:,1) = rNa(:,1) + (omega(2)-omega(3)) * tauM * rLui(:,1)
94 rNa(:,2) = rNa(:,2) + (omega(3)-omega(1)) * tauM * rLui(:,2)
95 rNa(:,3) = rNa(:,3) + (omega(1)-omega(2)) * tauM * rLui(:,3)
161 tmp1 = tauM * rLui(:,1)
162 tmp2 = tauM * rLui(:,2)
163 tmp3 = tauM * rLui(:,3)
193 & * ( rLui(:,1) * g1yi(:,2)
194 & + rLui(:,2) * g2yi(:,2)
[all …]
|
| H A D | e3stab.f | 2 & u3, dxidx, rLui, argument
33 & rLui(npro,nsd),
299 tauBar = rLui(:,1) * ( gijd(:,1) * rLui(:,1)
300 & + gijd(:,4) * rLui(:,2)
301 & + gijd(:,6) * rLui(:,3) )
302 & + rLui(:,2) * ( gijd(:,4) * rLui(:,1)
303 & + gijd(:,2) * rLui(:,2)
304 & + gijd(:,5) * rLui(:,3) )
305 & + rLui(:,3) * ( gijd(:,6) * rLui(:,1)
306 & + gijd(:,5) * rLui(:,2)
[all …]
|
| H A D | e3ivar.f | 7 & ql, rLui, src, argument
65 & ql(npro,nshl,idflx), rLui(npro,nsd),
248 & rLui, src, divqi)
265 & tmp1(:)*rLui(:,1)*rLui(:,1)
267 & tmp1(:)*rLui(:,2)*rLui(:,2)
269 & tmp1(:)*rLui(:,3)*rLui(:,3)
342 real*8 rLui(npro,nsd), aci(npro,nsd), local
460 & rLui, src, divqi)
466 call e3uBar(rho, src, dxidx, rLui, rmu, uBar)
|
| H A D | e3.f | 57 & rLui(npro,nsd), uBar(npro,nsd),
104 & ql, rLui, src,
111 & u3, dxidx, rLui,
120 & rLui, rmu, rho,
131 & rLui, rmu,
|
| H A D | e3lhs.f | 3 & rLui, rmu, argument
37 & rLui(npro,3), rmu(npro),
120 tmp2 = tauBar * ( rLui(:,1) * shg(:,b,1)
121 & + rLui(:,2) * shg(:,b,2)
122 & + rLui(:,3) * shg(:,b,3) )
124 t3(:,1) = t2(:,1) + tmp1 * u1 + tmp2 * rLui(:,1)
125 t3(:,2) = t2(:,2) + tmp1 * u2 + tmp2 * rLui(:,2)
126 t3(:,3) = t2(:,3) + tmp1 * u3 + tmp2 * rLui(:,3)
|
| H A D | advLES.f | 32 & rLui(npro,3),
358 rLui(:,1) =(aci(:,1) + u1 * g1yi(:,2)
363 rLui(:,2) =(aci(:,2) + u1 * g1yi(:,3)
368 rLui(:,3) =(aci(:,3) + u1 * g1yi(:,4)
375 rLui(:,1)=rlui(:,1)+u1*divu
376 rLui(:,2)=rlui(:,2)+u2*divu
377 rLui(:,3)=rlui(:,3)+u3*divu
384 & u3, dxidx, rLui,
391 stressli(:,1) = u1(:)*rLui(:,1)
392 stressli(:,2) = u1(:)*rLui(:,2)
[all …]
|