| /phasta/phSolver/compressible/ |
| H A D | e3dc.f | 1 subroutine e3DC (g1yi, g2yi, g3yi, A0, raLS, argument 34 dimension g1yi(npro,nflow), g2yi(npro,nflow), 86 A0gyi( :,6) = A0(:,1,1)*g2yi(:,1) 87 & + A0(:,1,2)*g2yi(:,2) 88 & + A0(:,1,3)*g2yi(:,3) 89 & + A0(:,1,4)*g2yi(:,4) 90 & + A0(:,1,5)*g2yi(:,5) 91 A0gyi( :,7) = A0(:,2,1)*g2yi(:,1) 92 & + A0(:,2,2)*g2yi(:,2) 93 & + A0(:,2,3)*g2yi(:,3) [all …]
|
| H A D | e3qvar.f | 2 & xl, g1yi, g2yi, g3yi, argument 43 & g1yi(npro,nflow), g2yi(npro,nflow), 110 g2yi = zero 126 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * ycl(:,n,1) 127 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * ycl(:,n,2) 128 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * ycl(:,n,3) 129 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * ycl(:,n,4) 130 g2yi(:,5) = g2yi(:,5) + shg(:,n,2) * ycl(:,n,5)
|
| H A D | e3ql.f | 32 dimension g1yi(npro,nflow), g2yi(npro,nflow), 76 & g2yi, g3yi, shg, 96 & + rlm * g2yi(:,3) 99 & + rmu * g2yi(:,2) 104 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 112 & + rmu * g2yi(:,2) 114 & + rlm2mu * g2yi(:,3) 116 qdi(:,7) = rmu * g2yi(:,4) 119 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3) 120 & + rmu * u3 * g2yi(:,4) [all …]
|
| H A D | e3q.f | 33 dimension g1yi(npro,nflow), g2yi(npro,nflow), 82 & g2yi, g3yi, shg, 105 & + rlm * g2yi(:,3) 108 & + rmu * g2yi(:,2) 113 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 121 & + rmu * g2yi(:,2) 123 & + rlm2mu * g2yi(:,3) 125 qdi(:,7) = rmu * g2yi(:,4) 128 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3) 129 & + rmu * u3 * g2yi(:,4) [all …]
|
| H A D | e3ivar.f | 4 & g1yi, g2yi, g3yi, argument 69 & g2yi(npro,nflow), g3yi(npro,nflow), 248 g2yi = zero 285 g2yi(:,1) = g2yi(:,1) + shg(:,1,2) * yl(:,1,1) 290 g2yi(:,2) = g2yi(:,2) + shg(:,1,2) * yl(:,1,2) 295 g2yi(:,3) = g2yi(:,3) + shg(:,1,2) * yl(:,1,3) 300 g2yi(:,4) = g2yi(:,4) + shg(:,1,2) * yl(:,1,4) 305 g2yi(:,5) = g2yi(:,5) + shg(:,1,2) * yl(:,1,5) 343 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * yl(:,n,1) 344 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) [all …]
|
| H A D | e3conv.f | 1 subroutine e3conv (g1yi, g2yi, g3yi, argument 49 dimension g1yi(npro,nflow), g2yi(npro,nflow), 106 & + A2(:,1,1) * g2yi(:,1) 108 & + A2(:,1,3) * g2yi(:,3) 110 & + A2(:,1,5) * g2yi(:,5) 122 & + A2(:,2,1) * g2yi(:,1) 123 & + A2(:,2,2) * g2yi(:,2) 124 & + A2(:,2,3) * g2yi(:,3) 126 & + A2(:,2,5) * g2yi(:,5) 138 & + A2(:,3,1) * g2yi(:,1) [all …]
|
| H A D | vorticity.f | 101 real*8 g2yi(npro,nflow) local 126 g2yi = zero 132 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) ! du/dy 136 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) ! dv/dy 140 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4) ! dw/dy 154 & (g2yi(:,4)-g3yi(:,3)) !dw/dy-dv/dz 162 & (g1yi(:,3)-g2yi(:,2)) !dv/dx-du/dy 167 & ( 0.5*(g2yi(:,4)-g3yi(:,3))**2 !dw/dy-dv/dz 169 & +0.5*(g1yi(:,3)-g2yi(:,2))**2 !dv/dx-du/dy 170 & -g1yi(:,2)**2-g2yi(:,3)**2-g3yi(:,4)**2 [all …]
|
| H A D | e3visc.f | 1 subroutine e3visc (g1yi, g2yi, g3yi, argument 42 dimension g1yi(npro,nflow), g2yi(npro,nflow), 279 & + rlm * g2yi(:,3) 283 & + rmu * g2yi(:,2) 290 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 304 & + rmu * g2yi(:,2) 307 & + rlm2mu * g2yi(:,3) 310 rmi(:, 9) = rmu * g2yi(:,4) 314 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3) 315 & + rmu * u3 * g2yi(:,4) [all …]
|
| H A D | e3b.f | 63 dimension g1yi(npro,nflow), g2yi(npro,nflow), 130 & g2yi, g3yi, WdetJb, 198 tau1n = bnorm(:,1) * (rlm2mu* g1yi(:,2) + rlm *g2yi(:,3) 200 & + bnorm(:,2) * (rmu *(g2yi(:,2) + g1yi(:,3))) 202 tau2n = bnorm(:,1) * (rmu *(g2yi(:,2) + g1yi(:,3))) 203 & + bnorm(:,2) * (rlm * g1yi(:,2) + rlm2mu*g2yi(:,3) 205 & + bnorm(:,3) * (rmu *(g3yi(:,3) + g2yi(:,4))) 207 & + bnorm(:,2) * (rmu *(g3yi(:,3) + g2yi(:,4))) 208 & + bnorm(:,3) * (rlm * g1yi(:,2) + rlm *g2yi(:,3) 224 & bnorm(:,2) * g2yi(:,5) +
|
| H A D | e3.f | 63 dimension g1yi(npro,nflow), g2yi(npro,nflow), 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 | 2 & xlb, lnode, g1yi, g2yi, argument 59 & g2yi(npro,nflow), g3yi(npro,nflow), 259 g2yi(:,2) = dxidxb(:,1,2) * gl1yi(:,2) + 269 g2yi(:,3) = dxidxb(:,1,2) * gl1yi(:,3) + 279 g2yi(:,4) = dxidxb(:,1,2) * gl1yi(:,4) + 289 g2yi(:,5) = dxidxb(:,1,2) * 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)) 103 tmp1 = rmu * ( g2yi(:,2) + g1yi(:,3) ) 104 tmp2 = rmu * ( g3yi(:,3) + g2yi(:,4) ) 114 & + ubar(:,2) * g2yi(:,2) 118 & + ubar(:,2) * g2yi(:,3) 122 & + ubar(:,2) * g2yi(:,4) 129 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp 143 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp - u2(:)*u2(:)*rho(:) [all …]
|
| H A D | e3ivar.f | 3 & aci, g1yi, g2yi, argument 59 & g2yi(npro,ndof), g3yi(npro,ndof), 150 g2yi = zero 158 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * yl(:,n,1) 159 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) 160 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) 161 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4) 247 & g1yi, g2yi, g3yi, 343 & g1yi(npro,nflow), g2yi(npro,nflow), local 430 g2yi = zero [all …]
|
| H A D | e3qvar.f | 3 & g2yi, g3yi, shg, argument 39 & g1yi(npro,nflow), g2yi(npro,nflow), 97 g2yi = zero 126 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) 127 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) 128 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4)
|
| H A D | e3qGradV.f | 30 dimension g1yi(npro,nflow), g2yi(npro,nflow), 56 & g2yi, g3yi, shg, 64 ql(:,i,2 ) = ql(:,i,2 )+ shape(:,i)*WdetJ*g2yi(:,2 ) ! du/dy 68 ql(:,i,5 ) = ql(:,i,5 )+ shape(:,i)*WdetJ*g2yi(:,3 ) ! dv/dy 72 ql(:,i,8 ) = ql(:,i,8 )+ shape(:,i)*WdetJ*g2yi(:,4 ) ! dw/dy
|
| H A D | e3q.f | 32 dimension g1yi(npro,nflow), g2yi(npro,nflow), 75 & g2yi, g3yi, shg, 92 qdi(:,4) = rmu * (g1yi(:,3) + g2yi(:,2)) 97 qdi(:,2) = rmu * (g1yi(:,3) + g2yi(:,2)) 98 qdi(:,5) = two * rmu * g2yi(:,3) 99 qdi(:,8) = rmu * (g2yi(:,4) + g3yi(:,3)) 104 qdi(:,6)= rmu * (g2yi(:,4) + g3yi(:,3))
|
| H A D | e3ql.f | 32 dimension g1yi(npro,ndof), g2yi(npro,ndof), 63 & g2yi, g3yi, shg, 71 qdi(:,4) = rmu * (g1yi(:,3) + g2yi(:,2)) 76 qdi(:,2) = rmu * (g1yi(:,3) + g2yi(:,2)) 77 qdi(:,5) = two * rmu * g2yi(:,3) 78 qdi(:,8) = rmu * (g2yi(:,4) + g3yi(:,3)) 83 qdi(:,6)= rmu * (g2yi(:,4) + g3yi(:,3))
|
| H A D | e3.f | 51 dimension g1yi(npro,ndof), g2yi(npro,ndof), 100 & aci, g1yi, g2yi, 119 & g1yi, g2yi, g3yi,
|
| H A D | e3bvar.f | 55 & g2yi(npro,ndof), g3yi(npro,ndof), 308 g2yi(:,2) = dxidxb(:,1,2) * gl1yi(:,2) + 318 g2yi(:,3) = dxidxb(:,1,2) * gl1yi(:,3) + 328 g2yi(:,4) = dxidxb(:,1,2) * gl1yi(:,4) + 350 & + bnorm(:,2) * (rmu * (g2yi(:,2) + g1yi(:,3))) 352 tau2n = bnorm(:,1) * (rmu * (g2yi(:,2) + g1yi(:,3))) 353 & + bnorm(:,2) * two * rmu * g2yi(:,3) 354 & + bnorm(:,3) * (rmu * (g3yi(:,3) + g2yi(:,4))) 356 & + bnorm(:,2) * (rmu * (g3yi(:,3) + g2yi(:,4)))
|
| H A D | advLES.f | 27 & g1yi(npro,ndof), g2yi(npro,ndof), 181 g2yi = zero 189 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * yl(:,n,1) 190 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) 191 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) 192 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4) 203 Sij(:,2) = g2yi(:,3) 205 Sij(:,4) = (g2yi(:,2)+g1yi(:,3))*pt5 207 Sij(:,6) = (g3yi(:,3)+g2yi(:,4))*pt5 359 & + u2 * g2yi(:,2) [all …]
|
| H A D | e3b.f | 61 dimension g1yi(npro,ndof), g2yi(npro,ndof),
|