| /phasta/phSolver/compressible/ |
| H A D | getdiff.f | 46 & xx(npro)
54 xx=zero
56 xx(:)=xx(:) + shp(:,n) * xl(:,n,1)
90 where(xx(:).le. 0.42) !halfway btwn AIP and exit
92 elsewhere(xx(:).ge. 0.75) !2/3 of the way to the exit
96 & (55.65294821-55.65294821*fmax)*xx(:)*xx(:)*xx(:)
97 & +(-97.67092412+97.67092412*fmax)*xx(:)*xx(:)
98 & +(52.59203606-52.59203606*fmax)*xx(:)
112 where(xx(:).le.xstart)
114 elsewhere(xx(:).ge.xmidwy)
[all …]
|
| H A D | e3source.f | 34 & xl(npro,nenl,3), xx(npro,3)
80 xx=zero
82 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
83 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
84 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3)
95 rsq=xx(id,1)*xx(id,1) + xx(id,2)*xx(id,2)
132 radsqr=xx(id,2)**2+xx(id,1)**2
133 if(xx(id,3).lt. zinSponge) then ! map this into big outflow
137 xx(id,3)=(zinSponge-xx(id,3))*CoefRatioI2O
145 if((xx(id,3).gt.zoutSponge).or.(radsqr.gt.radsts)) then
[all …]
|
| /phasta/phSolver/incompressible/ |
| H A D | e3res.f | 353 & aci, u1, u2, u3, Temp, rho, xx, argument
359 & aci, xx local
392 call e3source(xx, src)
416 & -omega(2)*(omega(1)*xx(:,2)-omega(2)*xx(:,1))
417 & -omega(3)*(omega(1)*xx(:,3)-omega(3)*xx(:,1))
420 & -omega(1)*(omega(2)*xx(:,1)-omega(1)*xx(:,2))
421 & -omega(3)*(omega(2)*xx(:,3)-omega(3)*xx(:,2))
424 & -omega(1)*(omega(3)*xx(:,1)-omega(1)*xx(:,3))
425 & -omega(2)*(omega(3)*xx(:,2)-omega(2)*xx(:,3))
|
| H A D | e3ivar.f | 66 & src(npro,nsd), Temp(npro),xx(npro,nsd)
100 xx = zero
102 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
103 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
104 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3)
246 & aci, u1, u2, u3, Temp, rho, xx,
338 & umod(npro,nsd), Temp(npro),xx(npro,nsd),
416 xx = zero
418 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
419 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
[all …]
|
| H A D | e3source.f | 1 subroutine e3source(xx, src) argument
13 real*8 xx(npro,nsd), src(npro,nsd) local
26 r = sqrt( xx(iel,1)**2 + xx(iel,2)**2)
31 src(iel,1) = -xx(iel,2)/r * Stheta
32 src(iel,2) = xx(iel,1)/r * Stheta
170 x = xx(iel,1)
171 y = xx(iel,2)
172 z = xx(iel,3)
|
| H A D | advLES.f | 30 & xx(npro,nsd),
291 xx = zero
293 xx(:,1) = xx(:,1) + shp(n,intp) * xl(:,n,1)
294 xx(:,2) = xx(:,2) + shp(n,intp) * xl(:,n,2)
295 xx(:,3) = xx(:,3) + shp(n,intp) * xl(:,n,3)
298 call e3source(xx, src)
319 xx = zero
321 xx(:,1) = xx(:,1) + shp(n,intp) * xl(:,n,1)
322 xx(:,2) = xx(:,2) + shp(n,intp) * xl(:,n,2)
323 xx(:,3) = xx(:,3) + shp(n,intp) * xl(:,n,3)
[all …]
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 1 subroutine e3source(xx, src)
13 real*8 xx(npro,nsd), src(npro,nsd)
26 r = sqrt( xx(iel,1)**2 + xx(iel,2)**2)
31 src(iel,1) = -xx(iel,2)/r * Stheta
32 src(iel,2) = xx(iel,1)/r * Stheta
49 c$$$ x = xx(iel,1)
50 c$$$ y = xx(iel,2)
51 c$$$ z = xx(iel,3)
170 x = xx(iel,1)
171 y = xx(iel,2)
[all …]
|
| H A D | e3source.fRHSonly | 1 subroutine e3source(xx, src)
13 real*8 xx(npro,nsd), src(npro,nsd)
26 r = sqrt( xx(iel,1)**2 + xx(iel,2)**2)
31 src(iel,1) = -xx(iel,2)/r * Stheta
32 src(iel,2) = xx(iel,1)/r * Stheta
49 c$$$ x = xx(iel,1)
50 c$$$ y = xx(iel,2)
51 c$$$ z = xx(iel,3)
170 x = xx(iel,1)
171 y = xx(iel,2)
[all …]
|