| /phasta/phSolver/incompressible/ ! |
| H A D | e3dc.f | 469 subroutine e3dcSclr ( gradS, giju, gGradS, argument
503 dimension gradS(npro,nsd), gGradS(npro,nsd),
511 gGradS(:,1) = GijU(:,1) * gradS(:,1)
512 1 + GijU(:,4) * gradS(:,2)
513 2 + GijU(:,6) * gradS(:,3)
514 gGradS(:,2) = GijU(:,4) * gradS(:,1)
515 1 + GijU(:,2) * gradS(:,2)
516 2 + GijU(:,5) * gradS(:,3)
517 gGradS(:,3) = GijU(:,6) * gradS(:,1)
518 1 + GijU(:,5) * gradS(:,2)
[all …]
|
| H A D | e3source.f | 210 subroutine e3sourceSclr ( Sclr, Sdot, gradS, argument
240 & gradS(npro,nsd), dwl(npro,nenl), local
446 uMod(:,1) = u1(:) - tmp * gradS(:,1)
447 uMod(:,2) = u2(:) - tmp * gradS(:,2)
448 uMod(:,3) = u3(:) - tmp * gradS(:,3)
569 mytmp = sign_levelset / sqrt( gradS(:,1)*gradS(:,1)
570 & + gradS(:,2)*gradS(:,2)
571 & + gradS(:,3)*gradS(:,3))
573 uMod(:,1) = mytmp * gradS(:,1) ! These are for the LHS
574 uMod(:,2) = mytmp * gradS(:,2)
[all …]
|
| H A D | e3ivar.f | 316 & Sclr, Sdot, gradS, argument
329 & gradS(npro,nsd), shpfun(npro,nshl),
385 gradS = zero
387 gradS(:,1) = gradS(:,1) + shg(:,n,1) * yl(:,n,id)
388 gradS(:,2) = gradS(:,2) + shg(:,n,2) * yl(:,n,id)
389 gradS(:,3) = gradS(:,3) + shg(:,n,3) * yl(:,n,id)
492 call e3sourceSclr ( Sclr, Sdot, gradS, dwl,
506 rLS(:) = ( Sdot(:) + (u1*gradS(:,1) +
507 & u2*gradS(:,2) +
508 & u3*gradS(:,3)) )
|
| H A D | e3res.f | 269 & Sclr, Sdot, gradS, argument
278 & Sclr(npro), Sdot(npro), gradS(npro,nsd), local
302 rNa = rcp*(tmps*Sdot + uMod(:,1) * gradS(:,1)
303 & + uMod(:,2) * gradS(:,2)
304 & + uMod(:,3) * gradS(:,3) )
313 rGNa(:,1) = diffus * gradS(:,1) + uMod(:,1) * tmp
314 rGNa(:,2) = diffus * gradS(:,2) + uMod(:,2) * tmp
315 rGNa(:,3) = diffus * gradS(:,3) + uMod(:,3) * tmp
|
| H A D | e3.f | 187 real*8 gradS(npro,nsd), shg(npro,nshl,nsd), local
233 & Sclr, Sdot, gradS,
254 call e3dcSclr ( gradS, giju, gGradS,
263 & Sclr, Sdot, gradS,
273 & Sclr, Sdot, gradS,
|
| H A D | e3bvar.f | 810 & gradSl(npro,nsd), gradS(npro,nsd)
924 gradS(:,1) = dxidxb(:,1,1) * gradSl(:,1) +
929 gradS(:,2) = dxidxb(:,1,2) * gradSl(:,1) +
933 gradS(:,3) = dxidxb(:,1,3) * gradSl(:,1) +
941 flux = diffus * ( gradS(:,1) * bnorm(:,1)
942 & + gradS(:,2) * bnorm(:,2)
943 & + gradS(:,3) * bnorm(:,3) )
|
| H A D | e3lhs.f | 271 & Sclr, Sdot, gradS, argument
281 & Sclr(npro), Sdot(npro), gradS(npro,nsd), local
|
| /phasta/phSolver/incompressible/e3source/ ! |
| H A D | e3source.fimplicit | 210 subroutine e3sourceSclr ( Sclr, Sdot, gradS,
225 c gradS : spatial gradient of Sclr
240 & gradS(npro,nsd), dwl(npro,nenl),
422 uMod(:,1) = u1(:) - tmp * gradS(:,1)
423 uMod(:,2) = u2(:) - tmp * gradS(:,2)
424 uMod(:,3) = u3(:) - tmp * gradS(:,3)
545 mytmp = sign_levelset / sqrt( gradS(:,1)*gradS(:,1)
546 & + gradS(:,2)*gradS(:,2)
547 & + gradS(:,3)*gradS(:,3))
549 uMod(:,1) = mytmp * gradS(:,1) ! These are for the LHS
[all …]
|
| H A D | e3source.fRHSonly | 210 subroutine e3sourceSclr ( Sclr, Sdot, gradS,
225 c gradS : spatial gradient of Sclr
240 & gradS(npro,nsd), dwl(npro,nenl),
408 uMod(:,1) = u1(:) - tmp * gradS(:,1)
409 uMod(:,2) = u2(:) - tmp * gradS(:,2)
410 uMod(:,3) = u3(:) - tmp * gradS(:,3)
531 mytmp = sign_levelset / sqrt( gradS(:,1)*gradS(:,1)
532 & + gradS(:,2)*gradS(:,2)
533 & + gradS(:,3)*gradS(:,3))
535 uMod(:,1) = mytmp * gradS(:,1) ! These are for the LHS
[all …]
|