| /phasta/phSolver/compressible/ |
| H A D | e3dc.f | 1 subroutine e3DC (g1yi, g2yi, g3yi, A0, raLS, argument 35 & g3yi(npro,nflow), A0(npro,5,5), 114 A0gyi(:,11) = A0(:,1,1)*g3yi(:,1) 115 & + A0(:,1,2)*g3yi(:,2) 116 & + A0(:,1,3)*g3yi(:,3) 117 & + A0(:,1,4)*g3yi(:,4) 118 & + A0(:,1,5)*g3yi(:,5) 119 A0gyi(:,12) = A0(:,2,1)*g3yi(:,1) 120 & + A0(:,2,2)*g3yi(:,2) 121 & + A0(:,2,3)*g3yi(:,3) [all …]
|
| H A D | e3qvar.f | 2 & xl, g1yi, g2yi, g3yi, argument 44 & g3yi(npro,nflow), shg(npro,nshl,nsd), 111 g3yi = zero 132 g3yi(:,1) = g3yi(:,1) + shg(:,n,3) * ycl(:,n,1) 133 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * ycl(:,n,2) 134 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * ycl(:,n,3) 135 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * ycl(:,n,4) 136 g3yi(:,5) = g3yi(:,5) + shg(:,n,3) * ycl(:,n,5)
|
| H A D | e3ql.f | 33 & g3yi(npro,nflow), shg(npro,nshl,nsd), 76 & g2yi, g3yi, shg, 97 & + rlm * g3yi(:,4) 101 & + rmu * g3yi(:,2) 105 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) 115 & + rlm * g3yi(:,4) 117 & + rmu * g3yi(:,3) 121 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4) 127 & + rmu * g3yi(:,2) 129 & + rmu * g3yi(:,3) [all …]
|
| H A D | e3q.f | 34 & g3yi(npro,nflow), shg(npro,nshl,nsd), 82 & g2yi, g3yi, shg, 106 & + rlm * g3yi(:,4) 110 & + rmu * g3yi(:,2) 114 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) 124 & + rlm * g3yi(:,4) 126 & + rmu * g3yi(:,3) 130 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4) 136 & + rmu * g3yi(:,2) 138 & + rmu * g3yi(:,3) [all …]
|
| H A D | e3ivar.f | 4 & g1yi, g2yi, g3yi, argument 69 & g2yi(npro,nflow), g3yi(npro,nflow), 249 g3yi = zero 310 g3yi(:,1) = g3yi(:,1) + shg(:,1,3) * yl(:,1,1) 315 g3yi(:,2) = g3yi(:,2) + shg(:,1,3) * yl(:,1,2) 320 g3yi(:,3) = g3yi(:,3) + shg(:,1,3) * yl(:,1,3) 325 g3yi(:,4) = g3yi(:,4) + shg(:,1,3) * yl(:,1,4) 330 g3yi(:,5) = g3yi(:,5) + shg(:,1,3) * yl(:,1,5) 349 g3yi(:,1) = g3yi(:,1) + shg(:,n,3) * yl(:,n,1) 350 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) [all …]
|
| H A D | e3conv.f | 1 subroutine e3conv (g1yi, g2yi, g3yi, argument 50 & g3yi(npro,nflow), 111 & + A3(:,1,1) * g3yi(:,1) 114 & + A3(:,1,4) * g3yi(:,4) 115 & + A3(:,1,5) * g3yi(:,5) 127 & + A3(:,2,1) * g3yi(:,1) 128 & + A3(:,2,2) * g3yi(:,2) 130 & + A3(:,2,4) * g3yi(:,4) 131 & + A3(:,2,5) * g3yi(:,5) 143 & + A3(:,3,1) * g3yi(:,1) [all …]
|
| H A D | vorticity.f | 102 real*8 g3yi(npro,nflow) local 127 g3yi = zero 133 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) ! du/dz 137 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3) ! dv/dz 141 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * yl(:,n,4) ! dw/dz 154 & (g2yi(:,4)-g3yi(:,3)) !dw/dy-dv/dz 158 & (g3yi(:,2)-g1yi(:,4)) !du/dz-dw/dx 167 & ( 0.5*(g2yi(:,4)-g3yi(:,3))**2 !dw/dy-dv/dz 168 & +0.5*(g3yi(:,2)-g1yi(:,4))**2 !du/dz-dw/dx 170 & -g1yi(:,2)**2-g2yi(:,3)**2-g3yi(:,4)**2 [all …]
|
| H A D | e3visc.f | 1 subroutine e3visc (g1yi, g2yi, g3yi, argument 43 & g3yi(npro,nflow), 280 & + rlm * g3yi(:,4) 286 & + rmu * g3yi(:,2) 291 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) 308 & + rlm * g3yi(:,4) 311 & + rmu * g3yi(:,3) 316 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4) 328 & + rmu * g3yi(:,2) 331 & + rmu * g3yi(:,3) [all …]
|
| H A D | e3b.f | 64 & g3yi(npro,nflow), WdetJb(npro), 130 & g2yi, g3yi, WdetJb, 199 & + rlm *g3yi(:,4)) 201 & + bnorm(:,3) * (rmu *(g3yi(:,2) + g1yi(:,4))) 204 & + rlm *g3yi(:,4)) 205 & + bnorm(:,3) * (rmu *(g3yi(:,3) + g2yi(:,4))) 206 tau3n = bnorm(:,1) * (rmu *(g3yi(:,2) + g1yi(:,4))) 207 & + bnorm(:,2) * (rmu *(g3yi(:,3) + g2yi(:,4))) 209 & + rlm2mu*g3yi(:,4)) 225 & bnorm(:,3) * g3yi(:,5) )
|
| H A D | e3.f | 64 & g3yi(npro,nflow), shg(npro,nshl,nsd), 132 & g1yi, g2yi, g3yi, 163 call e3conv (g1yi, g2yi, g3yi, 177 call e3visc (g1yi, g2yi, g3yi, 218 call e3dc (g1yi, g2yi, g3yi,
|
| H A D | e3bvar.f | 3 & g3yi, WdetJb, bnorm, pres, T, argument 59 & g2yi(npro,nflow), g3yi(npro,nflow), 262 g3yi(:,2) = dxidxb(:,1,3) * gl1yi(:,2) + 272 g3yi(:,3) = dxidxb(:,1,3) * gl1yi(:,3) + 282 g3yi(:,4) = dxidxb(:,1,3) * gl1yi(:,4) + 292 g3yi(:,5) = dxidxb(:,1,3) * gl1yi(:,5) +
|
| /phasta/phSolver/incompressible/ |
| H A D | e3res.f | 3 & g1yi, g2yi, g3yi, argument 43 & g1yi(npro,nflow), g2yi(npro,nflow), g3yi(npro,nflow), 102 tmp = -pres + tauC * (g1yi(:,2) + g2yi(:,3) + g3yi(:,4)) 104 tmp2 = rmu * ( g3yi(:,3) + g2yi(:,4) ) 105 tmp3 = rmu * ( g1yi(:,4) + g3yi(:,2) ) 115 & + ubar(:,3) * g3yi(:,2) 119 & + ubar(:,3) * g3yi(:,3) 123 & + ubar(:,3) * g3yi(:,4) 133 rGNa(:,3,3) = two * rmu * g3yi(:,4) + tmp 147 rGNa(:,3,3) = two * rmu * g3yi(:,4) + tmp - u3(:)*u3(:)*rho(:) [all …]
|
| H A D | e3ivar.f | 4 & g3yi, shg, dxidx, argument 59 & g2yi(npro,ndof), g3yi(npro,ndof), 151 g3yi = zero 163 g3yi(:,1) = g3yi(:,1) + shg(:,n,3) * yl(:,n,1) 164 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) 165 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3) 166 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * yl(:,n,4) 247 & g1yi, g2yi, g3yi, 344 & g3yi(npro,nflow), local 431 g3yi = zero [all …]
|
| H A D | e3qvar.f | 3 & g2yi, g3yi, shg, argument 40 & g3yi(npro,nflow), shg(npro,nshl,nsd), 98 g3yi = zero 130 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) 131 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3) 132 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * yl(:,n,4)
|
| H A D | e3qGradV.f | 31 & g3yi(npro,nflow), shg(npro,nshl,nsd), 56 & g2yi, g3yi, shg, 65 ql(:,i,3 ) = ql(:,i,3 )+ shape(:,i)*WdetJ*g3yi(:,2 ) ! du/dz 69 ql(:,i,6 ) = ql(:,i,6 )+ shape(:,i)*WdetJ*g3yi(:,3 ) ! dv/dz 73 ql(:,i,9 ) = ql(:,i,9 )+ shape(:,i)*WdetJ*g3yi(:,4 ) ! dw/dz
|
| H A D | e3q.f | 33 & g3yi(npro,nflow), shg(npro,nshl,nsd), 75 & g2yi, g3yi, shg, 93 qdi(:,7) = rmu * (g1yi(:,4) + g3yi(:,2)) 99 qdi(:,8) = rmu * (g2yi(:,4) + g3yi(:,3)) 103 qdi(:,3) = rmu * (g1yi(:,4) + g3yi(:,2)) 104 qdi(:,6)= rmu * (g2yi(:,4) + g3yi(:,3)) 105 qdi(:,9)= two * rmu * g3yi(:,4)
|
| H A D | e3ql.f | 33 & g3yi(npro,ndof), shg(npro,nshl,nsd), 63 & g2yi, g3yi, shg, 72 qdi(:,7) = rmu * (g1yi(:,4) + g3yi(:,2)) 78 qdi(:,8) = rmu * (g2yi(:,4) + g3yi(:,3)) 82 qdi(:,3) = rmu * (g1yi(:,4) + g3yi(:,2)) 83 qdi(:,6)= rmu * (g2yi(:,4) + g3yi(:,3)) 84 qdi(:,9)= two * rmu * g3yi(:,4)
|
| H A D | e3.f | 52 & g3yi(npro,ndof), shg(npro,nshl,nsd), 101 & g3yi, shg, dxidx, 119 & g1yi, g2yi, g3yi,
|
| H A D | e3bvar.f | 55 & g2yi(npro,ndof), g3yi(npro,ndof), 311 g3yi(:,2) = dxidxb(:,1,3) * gl1yi(:,2) + 321 g3yi(:,3) = dxidxb(:,1,3) * gl1yi(:,3) + 331 g3yi(:,4) = dxidxb(:,1,3) * gl1yi(:,4) + 351 & + bnorm(:,3) * (rmu * (g3yi(:,2) + g1yi(:,4))) 354 & + bnorm(:,3) * (rmu * (g3yi(:,3) + g2yi(:,4))) 355 tau3n = bnorm(:,1) * (rmu * (g3yi(:,2) + g1yi(:,4))) 356 & + bnorm(:,2) * (rmu * (g3yi(:,3) + g2yi(:,4))) 357 & + bnorm(:,3) * two * rmu * g3yi(:,4)
|
| H A D | advLES.f | 28 & g3yi(npro,ndof), divqi(npro,3), 182 g3yi = zero 194 g3yi(:,1) = g3yi(:,1) + shg(:,n,3) * yl(:,n,1) 195 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) 196 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3) 197 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * yl(:,n,4) 204 Sij(:,3) = g3yi(:,4) 206 Sij(:,5) = (g3yi(:,2)+g1yi(:,4))*pt5 207 Sij(:,6) = (g3yi(:,3)+g2yi(:,4))*pt5 360 & + u3 * g3yi(:,2) - src(:,1) ) * rho [all …]
|
| H A D | e3b.f | 62 & g3yi(npro,ndof), WdetJb(npro),
|