| /phasta/phSolver/compressible/ |
| H A D | e3visc.f | 3 & rmu, rlm, rlm2mu,
51 & rlm2mu(npro), con(npro),
69 stiff(:, 2, 2) = rlm2mu
72 stiff(:, 5, 2) = rlm2mu * u1
103 stiff(:, 8, 8) = rlm2mu
106 stiff(:,10, 8) = rlm2mu * u2
135 stiff(:,14,14) = rlm2mu
138 stiff(:,15,14) = rlm2mu * u3
230 compK(:, 1) = f1 * T * rlm2mu + f3 * T * rmu
234 compK(:, 3) = f1 * T * rmu + f3 * T * rlm2mu
[all …]
|
| H A D | e3ql.f | 38 & rlm(npro), rlm2mu(npro),
88 & rmu, rlm, rlm2mu, con, shape,
95 qdi(:,1) = rlm2mu * g1yi(:,2)
102 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3)
114 & + rlm2mu * g2yi(:,3)
119 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3)
132 & + rlm2mu * g3yi(:,4)
136 & + rlm2mu * u3 * g3yi(:,4)
|
| H A D | e3q.f | 39 & rlm(npro), rlm2mu(npro),
93 & rmu, rlm, rlm2mu, con, shape,
104 qdi(:,1) = rlm2mu * g1yi(:,2)
111 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3)
123 & + rlm2mu * g2yi(:,3)
128 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3)
141 & + rlm2mu * g3yi(:,4)
145 & + rlm2mu * u3 * g3yi(:,4)
|
| H A D | getdiff.f | 2 & rmu, rlm, rlm2mu, con, shp, argument
43 & rlm2mu(npro), con(npro),
263 rlm2mu = rlm + two * rmu
274 & rlm2mu, con, rho, Sclr) argument
307 & rlm2mu(npro), con(npro),
353 rlm2mu = rlm + two * rmu
|
| H A D | e3b.f | 85 & rlm2mu(npro), con(npro),
191 & rmu, rlm, rlm2mu, con, shape,
198 tau1n = bnorm(:,1) * (rlm2mu* g1yi(:,2) + rlm *g2yi(:,3)
203 & + bnorm(:,2) * (rlm * g1yi(:,2) + rlm2mu*g2yi(:,3)
209 & + rlm2mu*g3yi(:,4))
450 & rlm2mu(npro), con(npro),
574 & rlm, rlm2mu, con, rho, Sclr)
|
| H A D | e3bdg_nd.f | 6 & rmu, rlm2mu, con) argument
21 & rlm2mu(npro), Atau2(npro,nflow,nflow),
566 A0(:,2,2) = A0(:,2,2) + rlm2mu !rlm2mu
569 A0(:,5,2) = A0(:,5,2) + rlm2mu * u1
783 A0(:,3,3) = A0(:,3,3) + rlm2mu !rlm2mu
786 A0(:,5,3) = A0(:,5,3) + rlm2mu * u2
1001 A0(:,4,4) = A0(:,4,4) + rlm2mu !rlm2mu
1004 A0(:,5,4) = A0(:,5,4) + rlm2mu * u3
1342 A0(:,2,3) = A0(:,2,3) + rlm2mu - rmu
1343 A0(:,3,2) = A0(:,3,2) + rlm2mu - rmu
[all …]
|
| H A D | e3bdg.f | 6 & rmu, rlm2mu, con) argument
21 & rlm2mu(npro),
603 A0(:,2,2) = A0(:,2,2) + rlm2mu !rlm2mu
606 A0(:,5,2) = A0(:,5,2) + rlm2mu * u1
864 A0(:,3,3) = A0(:,3,3) + rlm2mu !rlm2mu
867 A0(:,5,3) = A0(:,5,3) + rlm2mu * u2
1126 A0(:,4,4) = A0(:,4,4) + rlm2mu !rlm2mu
1129 A0(:,5,4) = A0(:,5,4) + rlm2mu * u3
1610 A0(:,2,3) = A0(:,2,3) + rlm2mu - rmu
1611 A0(:,3,2) = A0(:,3,2) + rlm2mu - rmu
[all …]
|
| H A D | e3.f | 52 & shp(nshl,ngauss), rlm2mu(npro),
140 & rmu, rlm, rlm2mu,
179 & rmu, rlm, rlm2mu,
271 & rmu, rlm2mu, con)
280 & rmu, rlm2mu, con)
|
| H A D | e3ivar.f | 11 & rmu, rlm, rlm2mu, argument
81 & rlm2mu(npro), con(npro),
230 & rmu, rlm, rlm2mu, con, shape,
|
| /phasta/phSolver/incompressible/ |
| H A D | getfld.f | 1 subroutine getFld (T, cp, rmu, rlm, rlm2mu, argument
35 & rlm2mu(npro), con(npro)
68 rlm2mu = rlm + two * rmu
|