| /phasta/phSolver/compressible/ |
| H A D | e3ivar.f | 2 & Sclr, shape, shgl, argument 82 & Sclr(npro) 112 Sclr = zero 115 Sclr = Sclr + shape(:,n) * ycl(:,n,isc) 120 call getthm (dui(:,1), dui(:,5), Sclr, 211 Sclr = zero 214 Sclr = Sclr + shape(:,n) * ycl(:,n,isc) 219 call getthm (pres, T, Sclr, 406 Sclr = zero 412 Sclr(i) = Sclr(i) + shape(i,n) * ycl(i,n,6) !scalar [all …]
|
| H A D | getthm.f | 1 subroutine getthm (pres, T, Sclr, rk, rho, argument 49 dimension pres(npro), Sclr(npro), 125 prop_blend(i) = 0.5*(one + Sclr(i)/epsilon_ls + 126 & (sin(pi*Sclr(i)/epsilon_ls))/pi ) 180 prop_blend(i) = 0.5*(one + Sclr(i)/epsilon_ls + 181 & (sin(pi*Sclr(i)/epsilon_ls))/pi ) 279 prop_blend(i) = 0.5*(one + Sclr(i)/epsilon_ls + 280 & (sin(pi*Sclr(i)/epsilon_ls))/pi )
|
| H A D | getdiff.f | 41 & rho(npro), Sclr(npro), 65 Sclr = zero 68 Sclr = Sclr + shp(:,n) * ycl(:,n,isc) 70 test_it = 0.5*(one + Sclr/epsilon_ls + 71 & (sin(pi*Sclr/epsilon_ls))/pi ) 274 & rlm2mu, con, rho, Sclr) 308 & rho(npro), Sclr(npro)
|
| H A D | e3qvar.f | 52 dimension dxdxi(npro,nsd,nsd), Sclr(npro) 86 Sclr = zero 89 Sclr = Sclr + shp(:,n) * ycl(:,n,isc) 95 call getthm (pres, T, Sclr,
|
| H A D | e3.f | 76 & Sclr(npro), dVdY(npro,15), 130 & Sclr, shape, shdrv, 245 call e3juel (yl, acl, Sclr, A0, WdetJ, rl, rml) 373 & Sclr(npro), srcp (npro) 432 & dxidx, Sclr, 445 & call e3sourceSclr (Sclr, rho, rmu, 470 & rho, u1, Sclr, 481 & rmu, Sclr, rho,
|
| H A D | e3conv.f | 358 & rho, u1, Sclr, 401 & g3yti(npro), Sclr(npro), 431 rti(:, 1) = rti(:,1) + (- u1) * rho * Sclr 432 rti(:, 2) = rti(:,2) + (- u2) * rho * Sclr 433 rti(:, 3) = rti(:,3) + (- u3) * rho * Sclr
|
| H A D | e3b.f | 443 dimension rho(npro), Sclr(npro), 489 & cp, rou, Sclr, 567 F = Sclr *rou 574 & rlm, rlm2mu, con, rho, Sclr) 589 …F = F + saSigmaInv*rho*((rmu/rho)+Sclr)*Sclrn!confirm the modificationc …
|
| H A D | e3visc.f | 44 & Sclr(npro), dxidx(npro,nsd,nsd), 362 & rmu, Sclr, rho, 396 & Sclr(npro), rho(npro), 416 stifft(:,1,1)=saSigmaInv*((rmu/rho)+max(zero,Sclr))
|
| H A D | e3bvar.f | 381 & Sclr, SclrF) 425 & rk(npro), Sclr(npro), 450 Sclr = zero 463 Sclr = Sclr + shpb(:,nodlcl) * ycl(:,nodlcl,id) 474 call getthm (pres, T, Sclr,
|
| H A D | e3source.f | 225 subroutine e3sourceSclr(Sclr, rho, rmu, argument 259 dimension Sclr (npro), ycl(npro,nshl,ndof), 297 sclrm=max(rnu/100.0,Sclr) ! sets a floor on SCLR
|
| H A D | e3bdg.f | 2600 & Sclr, 2610 & shp(npro,nshl), Sclr(npro), 2756 A0t(:) = A0t(:) +1.5*(rmu+Sclr) 2778 A0t(:) = A0t(:) +1.5*(rmu+Sclr) 2801 A0t(:) = A0t(:) +1.5*(rmu+Sclr)
|
| /phasta/phSolver/incompressible/ |
| H A D | e3ivar.f | 75 & Sclr(npro) 202 Sclr = zero 208 Sclr(i) = Sclr(i) + shpfun(i,n) * yl(i,n,6) !scalar 220 & + cos(pi*Sclr(i)/epsilon_ls)) * gyti(i,1) 222 & + cos(pi*Sclr(i)/epsilon_ls)) * gyti(i,2) 224 & + cos(pi*Sclr(i)/epsilon_ls)) * gyti(i,3) 316 & Sclr, Sdot, gradS, 328 & Sclr(npro), Sdot(npro), 357 Sclr = zero 364 Sclr = Sclr + shpfun(:,n) * yl(:,n,id) [all …]
|
| H A D | getdiff.f | 48 Sclr = zero 51 Sclr = Sclr + shape(:,n) * yl(:,n,isc) 57 prop_blend(i) = 0.5*(one + Sclr(i)/epsilon_ls + 58 & (sin(pi*Sclr(i)/epsilon_ls))/pi )
|
| H A D | e3.f | 188 & Sdot(npro), Sclr(npro), local 233 & Sclr, Sdot, gradS, 263 & Sclr, Sdot, gradS, 273 & Sclr, Sdot, gradS,
|
| H A D | e3source.f | 210 subroutine e3sourceSclr ( Sclr, Sdot, gradS, argument 239 real*8 Sclr(npro), Sdot(npro), local 324 SclrNN= max(Sclr(e),zero) 388 if(Sclr(e).gt.zero) ep=one
|
| H A D | e3res.f | 269 & Sclr, Sdot, gradS, 278 & Sclr(npro), Sdot(npro), gradS(npro,nsd), local
|
| H A D | e3lhs.f | 271 & Sclr, Sdot, gradS, 281 & Sclr(npro), Sdot(npro), gradS(npro,nsd), local
|
| /phasta/phSolver/common/ |
| H A D | solvecon.f | 184 dimension Sclr(npro), Sclrtmp(npro), 247 Sclr = Sclr + shape(:,i) * ycl(:,i,7) !d^kbar
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 210 subroutine e3sourceSclr ( Sclr, Sdot, gradS, 224 c Sclr : The scalar that is being transported in this pde. 225 c gradS : spatial gradient of Sclr 239 real*8 Sclr(npro), Sdot(npro), 320 SclrNN= max(Sclr(e),zero) 364 if(Sclr(e).gt.zero) ep=one 495 CAD If Sclr(:,1).gt.zero, result of sign_term function 1 496 CAD If Sclr(:,1).eq.zero, result of sign_term function 0 497 CAD If Sclr(:,1).lt.zero, result of sign_term function -1
|
| H A D | e3source.fRHSonly | 210 subroutine e3sourceSclr ( Sclr, Sdot, gradS, 224 c Sclr : The scalar that is being transported in this pde. 225 c gradS : spatial gradient of Sclr 239 real*8 Sclr(npro), Sdot(npro), 320 SclrNN= max(Sclr(e),zero) 350 if(Sclr(e).gt.zero) ep=one 481 CAD If Sclr(:,1).gt.zero, result of sign_term function 1 482 CAD If Sclr(:,1).eq.zero, result of sign_term function 0 483 CAD If Sclr(:,1).lt.zero, result of sign_term function -1
|