| /phasta/phSolver/compressible/ |
| H A D | e3mtrx.f | 4 & u1, u2, u3, 55 & u1(npro), u2(npro), 110 A0(:,3,1) = drdp * u2 114 A0(:,3,5) = drdT * u2 124 A0(:,5,3) = rho * u2 144 A1(:,3,1) = drdp * u1 * u2 145 A1(:,3,2) = rho * u2 148 A1(:,3,5) = drdT * u1 * u2 158 A1(:,5,3) = rho * u1 * u2 164 A2(:,1,1) = drdp * u2 [all …]
|
| H A D | e3visc.f | 4 & u1, u2, u3, 45 & u1(npro), u2(npro), 73 stiff(:, 5, 3) = rmu * u2 81 stiff(:, 5, 7) = rmu * u2 96 stiff(:,10, 2) = rlm * u2 106 stiff(:,10, 8) = rlm2mu * u2 115 stiff(:,10,14) = rlm * u2 129 stiff(:,15, 9) = rmu * u2 137 stiff(:,15,13) = rmu * u2 242 compK(:, 7) = f1 * T * rlm2mu * u1 + f2 * T * (rlm + rmu) * u2 [all …]
|
| H A D | e3tau.f | 2 & u1, u2, u3, 41 & u2(npro), u3(npro), 68 & + u2*u2*(gijd(:,1)*gijd(:,6)-gijd(:,4)*gijd(:,4)) 70 & + two *(u1*u2*(gijd(:,4)*gijd(:,5)-gijd(:,2)*gijd(:,6)) 72 & + u2*u3*(gijd(:,4)*gijd(:,2)-gijd(:,1)*gijd(:,5))) 83 & + u2*u2*gijd(:,3) 85 & +(u1*u2*gijd(:,2) 87 & + u2*u3*gijd(:,5))*two + 1.0e-15 !FIX FOR INVALID MESHES 163 & + (u1*(u1*gijd(:,1) + two*(u2*gijd(:,2)+u3*gijd(:,4))) 164 & + u2*(u2*gijd(:,3) + two*u3*gijd(:,5)) [all …]
|
| H A D | e3conv.f | 5 & u2, u3, rLyi, argument 56 & u2(npro), u3(npro), 78 ri(:, 3) = (- u1) * rho * u2 82 ri(:, 6) = (- u2) * rho 83 ri(:, 7) = (- u2) * rho * u1 84 ri(:, 8) = (- u2) * rho * u2 - pres 85 ri(:, 9) = (- u2) * rho * u3 86 ri(:,10) = (- u2) * rho * (ei + rk) - u2 * pres 90 ri(:,13) = (- u3) * rho * u2 359 & u2, u3, rLyti, argument [all …]
|
| H A D | e3ql.f | 36 & u1(npro), u2(npro), 78 & cp, u1, u2, 102 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3) 104 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 118 qdi(:,8) = rlm * u2 * g1yi(:,2) + rmu * u1 * g1yi(:,3) 119 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3) 121 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4) 134 & + rlm * u3 * g2yi(:,3) + rmu * u2 * g2yi(:,4) 135 & + rmu * u1 * g3yi(:,2) + rmu * u2 * g3yi(:,3)
|
| H A D | e3.f | 74 & u1(npro), u2(npro), 137 & u1, u2, u3, 152 & u1, u2, u3, 167 & u2, u3, rLyi, 180 & u1, u2, u3, 190 & rho, u1, u2, 204 & u1, u2, u3, 269 & u1, u2, u3, 278 & u1, u2, u3, 381 & u1(npro), u2(npro), [all …]
|
| H A D | e3q.f | 37 & u1(npro), u2(npro), 84 & cp, u1, u2, 111 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3) 113 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 127 qdi(:,8) = rlm * u2 * g1yi(:,2) + rmu * u1 * g1yi(:,3) 128 & + rmu * u1 * g2yi(:,2) + rlm2mu * u2 * g2yi(:,3) 130 & + rmu * u3 * g3yi(:,3) + rlm * u2 * g3yi(:,4) 143 & + rlm * u3 * g2yi(:,3) + rmu * u2 * g2yi(:,4) 144 & + rmu * u1 * g3yi(:,2) + rmu * u2 * g3yi(:,3)
|
| H A D | e3qvar.f | 4 & cp, u1, u2, u3) argument 47 & u1(npro), u2(npro), 62 u2 = zero 74 u2 = u2 + shp(:,n) * ycl(:,n,3)
|
| H A D | e3b.f | 71 & u1(npro), u2(npro), 132 & u1, u2, u3, 153 un = bnorm(:,1) * u1 + bnorm(:,2) * u2 + bnorm(:,3) * u3 169 F3 = rou * u2 + bnorm(:,2) * p 217 Fv5 = u1 * Fv2 + u2 * Fv3 + u3 * Fv4 437 dimension u1(npro), u2(npro), 486 & u2, u3, g1yti, 546 u2 = mytmp * g2yti 558 un = bnorm(:,1)*u1 + bnorm(:,2)*u2 + bnorm(:,3)*u3
|
| H A D | e3bvar.f | 4 & u1, u2, u3, rho, ei, argument 62 & u1(npro), u2(npro), 91 u2 = zero 100 u2 = u2 + shpb(:,nodlcl) * yl(:,nodlcl,3) 107 rk = pt5 * ( u1**2 + u2**2 + u3**2 ) 378 & u1, u2, u3, argument 422 & u1(npro), u2(npro), 447 u2 = zero 460 u2 = u2 + shpb(:,nodlcl) * ycl(:,nodlcl,3) 468 rk = pt5 * ( u1**2 + u2**2 + u3**2 )
|
| H A D | e3ivar.f | 9 & u1, u2, u3, argument 77 & u1(npro), u2(npro), 144 u2 = zero 153 u2 = u2 + shape(:,n) * ycl(:,n,3) 206 rk = pt5 * ( u1**2 + u2**2 + u3**2 ) 497 & rk, u1, u2, argument 567 & u1(npro), u2(npro), 590 u2 = zero 603 u2 = u2 + shape(:,n) * ycl( :,n,3) 637 rk = pt5 * ( u1**2 + u2**2 + u3**2 )
|
| H A D | MachControl.f90 | 104 real*8 :: u2, T local 110 u2 = vartsbuff(MC_iPP, 2, ivartsBuff)**2 + 114 MC_Mt = sqrt(u2/(1.4*Rgas*T))
|
| H A D | e3source.f | 2 & rho, u1, u2, 29 & u1(npro), u2(npro), 54 src(:,5) = u1*src(:,2) + u2*src(:,3) + u3*src(:,4) 188 src(:,5) = src(:,5) + (u1*sforce(:,1)+u2*sforce(:,2) 230 & u2, u3, xl, elDwl) 265 & u2 (npro), u3 (npro), 485 u2 = mytmp(:) * g2yti(:)
|
| H A D | e3ls.f | 4 & u1, u2, u3, argument 65 & u1(npro), u2(npro), 191 & u1, u2, u3, 209 & u1, u2, u3, 777 & u1, u2, u3, argument 834 & u1(npro), u2(npro), 892 & u1, u2, u3,
|
| /phasta/phSolver/incompressible/ |
| H A D | e3stab.f | 1 subroutine e3stab (rho, u1, u2, argument 31 & u2(npro), u3(npro), 79 4 + gijd(:,4) * u2 81 6 + u2 * ( gijd(:,4) * u1 82 7 + gijd(:,2) * u2 85 a + gijd(:,5) * u2 121 2 + two * (gijd(:,5)*gijd(:,6)-gijd(:,4)*gijd(:,3)) * u2 123 1 + u2 * ( (gijd(:,1)*gijd(:,3)-gijd(:,6)*gijd(:,6)) * u2 143 4 + gijd(:,4) * u2 145 6 + u2 * ( gijd(:,4) * u1 [all …]
|
| H A D | e3res.f | 1 subroutine e3Res ( u1, u2, u3, argument 41 dimension u1(npro), u2(npro), u3(npro), 140 rGNa(:,1,2) = tmp1 - u1(:)*u2(:)*rho(:) 142 rGNa(:,2,1) = tmp1 - u1(:)*u2(:)*rho(:) 143 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp - u2(:)*u2(:)*rho(:) 144 rGNa(:,2,3) = tmp2 - u3(:)*u2(:)*rho(:) 146 rGNa(:,3,2) = tmp2 - u3(:)*u2(:)*rho(:) 166 rGNa(:,1,2) = rGNa(:,1,2) + tmp1 * u2 169 rGNa(:,2,2) = rGNa(:,2,2) + tmp2 * u2 172 rGNa(:,3,2) = rGNa(:,3,2) + tmp3 * u2 [all …]
|
| H A D | e3ivar.f | 6 & u1, u2, u3, argument 63 & u1(npro), u2(npro), 83 u2 = zero 89 u2 = u2 + shpfun(:,n) * yl(:,n,3) 246 & aci, u1, u2, u3, Temp, rho, xx, 294 dkei = (u1(:)-dkei)**2 +u2(:)**2 ! u'^2+v'^2 318 & u1, u2, u3, argument 333 & u1(npro), u2(npro), 355 u2 = zero 362 u2 = u2 + shpfun(:,n) * yl(:,n,3) [all …]
|
| H A D | stats.f | 217 real*8 u1, u2, u3, r0, r1, r2, r3, r4, r5, t3, t4, t5 local 251 u2 = yAlpha(i,2) 258 r1 = res(2) + reg * u2 269 stsVelReg(i,2) = stsVelReg(i,2) + u2 273 r1 = res(2) * u2 + reg * u2 * u2 275 r3 = res(1) * u2 + res(2) * u1 + reg * u2 * u1 276 r4 = res(2) * u3 + res(3) * u2 + reg * u3 * u2 336 u2 = stsVec(i,2) 343 stsVel(i,2) = stsVel(i,2) + u2 347 stsVelSqr(i,2) = stsVelSqr(i,2) + u2*u2 [all …]
|
| H A D | e3.f | 56 & u2(npro), u3(npro), 103 & u1, u2, u3, 110 call e3stab (rho, u1, u2, 117 call e3Res ( u1, u2, u3, 129 call e3LHS ( u1, u2, u3, 190 & u1(npro), u2(npro), u3(npro), local 235 & u1, u2, u3,
|
| H A D | filters.f | 13 & u1(npro), u2(npro), 107 u2 = zero 117 u2 = u2 + shp(i,intp)*yl(:,i,3) 137 fresli(:,2) = WdetJ * u2 !G * u2 * WdetJ 141 fresli(:,5) = WdetJ * u2 * u2 ! G*u2*u2*WdetJ 143 fresli(:,7) = WdetJ * u1 * u2 ! G*u1*u2*WdetJ 145 fresli(:,9) = WdetJ * u2 * u3 ! G*u2*u3*WdetJ 186 & u1(npro), u2(npro), 412 & u1(npro), u2(npro), 507 u2 = zero [all …]
|
| H A D | getdiff.f | 216 u2=yl(e,n,3) 218 if((u1.eq.zero).and.(u2.eq.zero).and.(u3.eq.zero)) 302 u2=yl(e,n,3) 304 if((u1.eq.zero).and.(u2.eq.zero).and.(u3.eq.zero))
|
| H A D | e3source.f | 213 & u1, u2, u3, xl, argument 244 & u2(npro), u3(npro), local 447 uMod(:,2) = u2(:) - tmp * gradS(:,2) 507 uMod(:,2) = u2(:) - zero 578 u2 = umod(:,2)
|
| H A D | e3lhs.f | 1 subroutine e3LHS ( u1, u2, u3, argument 35 dimension u1(npro), u2(npro), u3(npro), 72 uBar(:,2) = tmp1 * u2(:) 116 & + u2 * shg(:,b,2) 125 t3(:,2) = t2(:,2) + tmp1 * u2 + tmp2 * rLui(:,2)
|
| /phasta/phSolver/common/ |
| H A D | asithf.f | 9 & u1(npro), u2(npro), 152 u2=zero 156 u2 = u2 + shp(i,intp)*ycl(:,i,3) 161 fresli(:,2) = fresli(:,22) * u2 !rho u2 * WdetJ 165 fresli(:,5) = fresli(:,2) * u2 !rho u2 u2 *WdetJ 167 fresli(:,7) = fresli(:,1) * u2 !rho u1 u2 *WdetJ
|
| H A D | bardmc.f | 680 & u1(npro), u2(npro), 762 u2=zero 766 u2 = u2 + shp(i,intp)*yl(:,i,3) 774 fresli(:,i,2) = rho * u2 * shp(i,intp) !rho * u2 * WdetJ 780 fresli(:,i,7) = rho * u2 * u2 * shp(i,intp) !rho * u2 * u2 * WdetJ 782 fresli(:,i,9) = rho * u1 * u2 * shp(i,intp) !rho * u1 * u2 * WdetJ 784 fresli(:,i,11)= rho * u2 * u3 * shp(i,intp) !rho * u2 * u3 * WdetJ 834 & u1(npro), u2(npro), 1000 u2=zero 1004 u2 = u2 + shp(i,intp)*yl(:,i,3) [all …]
|