Home
last modified time | relevance | path

Searched refs:yl (Results 1 – 25 of 30) sorted by relevance

12

/phasta/phSolver/compressible/
H A De3juel.f1 subroutine e3juel (yl, acl, sls, A0, argument
25 dimension yl(npro,nshl,nflow), acl(npro,nshl,ndof),
55 rk = pt5 * (yl(:,inod,2)**2 + yl(:,inod,3)**2 + yl(:,inod,4)**2)
58 call getthm (yl(:,inod,1), yl(:,inod,5), sls,
64 yl(:,inod,1) = rho
65 yl(:,inod,2) = rho * yl(:,inod,2)
66 yl(:,inod,3) = rho * yl(:,inod,3)
67 yl(:,inod,4) = rho * yl(:,inod,4)
68 yl(:,inod,5) = rho * (ei + rk)
70 ub(:,:)=yl(:,1,:)+yl(:,2,:)
[all …]
H A De3ivar.f1 subroutine e3ivar (yl, ycl, acl, argument
63 dimension yl(npro,nshl,nflow), ycl(npro,nshl,ndof),
97 dui(:,1) = dui(:,1) + shape(:,n) * yl(:,n,1) ! p
98 dui(:,2) = dui(:,2) + shape(:,n) * yl(:,n,2) ! u1
99 dui(:,3) = dui(:,3) + shape(:,n) * yl(:,n,3) ! u2
100 dui(:,4) = dui(:,4) + shape(:,n) * yl(:,n,4) ! u3
101 dui(:,5) = dui(:,5) + shape(:,n) * yl(:,n,5) ! T
260 g1yi(:,1) = g1yi(:,1) + shg(:,1,1) * yl(:,1,1)
261 & + shg(:,2,1) * yl(:,2,1)
262 & + shg(:,3,1) * yl(:,3,1)
[all …]
H A Dvorticity.f89 real*8 yl(npro,nshl,ndof) local
108 call localy (y, yl, ien, ndof, 'gather ')
131 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2) ! du/dx
132 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2) ! du/dy
133 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2) ! du/dz
135 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3) ! dv/dx
136 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3) ! dv/dy
137 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3) ! dv/dz
139 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4) ! dw/dx
140 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4) ! dw/dy
[all …]
H A De3bvar.f1 subroutine e3bvar (yl, ycl, BCB, shpb, shglb, argument
53 dimension yl(npro,nshl,nflow), BCB(npro,nshlb,ndBCB),
98 pres = pres + shpb(:,nodlcl) * yl(:,nodlcl,1)
99 u1 = u1 + shpb(:,nodlcl) * yl(:,nodlcl,2)
100 u2 = u2 + shpb(:,nodlcl) * yl(:,nodlcl,3)
101 u3 = u3 + shpb(:,nodlcl) * yl(:,nodlcl,4)
102 T = T + shpb(:,nodlcl) * yl(:,nodlcl,5)
235 gl1yi(:,1) = gl1yi(:,1) + shglb(:,1,n) * yl(:,n,1)
236 gl1yi(:,2) = gl1yi(:,2) + shglb(:,1,n) * yl(:,n,2)
237 gl1yi(:,3) = gl1yi(:,3) + shglb(:,1,n) * yl(:,n,3)
[all …]
H A Dasbres.f22 dimension yl(npro,nshl,nflow), xlb(npro,nenl,nsd),
36 call localy(y, yl, ienb, nflow, 'gather ')
43 call e3b (yl, ycl, iBCB, BCB, shpb, shglb,
H A Dasires.f23 dimension yl(npro,nshl,nflow), ycl(npro,nshl,ndofl),
53 call localy(y, yl, ien, nflow, 'gather ')
78 call e3 (yl, ycl, acl, shp,
/phasta/phSolver/incompressible/
H A De3ivar.f1 subroutine e3ivar (yl, acl, shpfun, argument
55 dimension yl(npro,nshl,ndof), dwl(npro,nenl),
87 pres = pres + shpfun(:,n) * yl(:,n,1)
88 u1 = u1 + shpfun(:,n) * yl(:,n,2)
89 u2 = u2 + shpfun(:,n) * yl(:,n,3)
90 u3 = u3 + shpfun(:,n) * yl(:,n,4)
95 Temp = Temp + shpfun(:,n) * yl(:,n,5)
153 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * yl(:,n,1)
154 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
155 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3)
[all …]
H A Dgetdiff.f1 subroutine getDiff( dwl,yl, shape, xmudmi, xl, rmu, rho) argument
9 real*8 yl(npro,nshl,ndof), rmu(npro), xmudmi(npro,ngauss),
51 Sclr = Sclr + shape(:,n) * yl(:,n,isc)
89 call AddEddyViscSA(yl, shape, rmu)
92 call AddEddyViscKE(yl, dwl, shape, rho, sigmaInv, rmu)
129 subroutine getdiffsclr(shape, dwl, yl, diffus) argument
136 real*8 yl(npro,nshl,ndof), dwl(npro,nshl), shape(npro,nshl) local
143 call AddSAVar(yl, shape, diffus)
151 call AddEddyViscKE(yl, dwl, shape, rho, sigmaInv, diffus)
191 subroutine AddEddyViscSA(yl,shape,rmu)
[all …]
H A De3qvar.f1 subroutine e3qvar (yl, shgl, argument
37 dimension yl(npro,nshl,ndof),
122 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
123 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3)
124 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4)
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)
130 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2)
131 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3)
[all …]
H A De3source.f212 & yl, dxidx, rmu, argument
242 & yl(npro,nshl,ndof), dxidx(npro,nsd,nsd), local
293 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * yl(:,n,2)
294 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * yl(:,n,3)
295 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * yl(:,n,4)
297 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * yl(:,n,2)
298 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * yl(:,n,3)
299 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * yl(:,n,4)
301 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * yl(:,n,2)
302 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * yl(:,n,3)
[all …]
H A De3q.f1 subroutine e3q (yl, dwl, shp, shgl, argument
24 dimension yl(npro,nshl,ndof), dwl(npro,nenl),
73 call e3qvar (yl, shdrv,
87 call getdiff(dwl, yl, shape, xmudmi, xl, rmu, tmp)
168 g1yti(:) = g1yti(:) + shg(:,n,1) * yl(:,n,6)
169 g2yti(:) = g2yti(:) + shg(:,n,2) * yl(:,n,6)
170 g3yti(:) = g3yti(:) + shg(:,n,3) * yl(:,n,6)
206 subroutine e3qSclr (yl, dwl, shp, shgl, argument
219 dimension yl(npro,nshl,ndof), dwl(npro,nshl),
258 call e3qvarSclr (yl, shdrv,
[all …]
H A Dasiq.f32 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd),
53 call localy(y, yl, ien, ndof, 'gather ')
64 call e3q (yl, dwl, shp, shgl,
100 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd),
111 call localy(y, yl, ien, ndof, 'gather ')
122 call e3qSclr (yl, dwl, shp, shgl,
H A De3.f1 subroutine e3 (yl, acl, dwl, shp, argument
40 dimension yl(npro,nshl,ndof),
75 call e3ql (yl, dwl, shp, shgl,
94 call getdiff(dwl, yl, shpfun, xmudmi, xl, rmu, rho)
98 call e3ivar (yl, acl, shpfun,
163 subroutine e3Sclr (yl, acl, shp, argument
179 real*8 yl(npro,nshl,ndof), acl(npro,nshl,ndof), local
211 call e3qlSclr (yl, dwl, shp, shgl, xl, ql, sgn)
227 call getdiffsclr(shpfun,dwl,yl,diffus)
231 call e3ivarSclr(yl, acl, shpfun,
H A Dasbmfg.f27 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd),
42 call localy(y, yl, ienb, ndofl, 'gather ')
64 call e3b (ul, yl, acl, iBCB, BCB,
100 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd),
112 call localy(y, yl, ienb, ndofl, 'gather ')
122 call e3bSclr (yl, iBCB, BCB, shpb, shglb,
H A Dasigmr.f30 dimension yl(npro,nshl,ndofl), acl(npro,nshl,ndofl),
58 call localy(y, yl, ien, ndofl, 'gather ')
86 call e3 (yl, acl, dwl, shp,
111 call timeseries(yl,xl,ien,sgn)
147 real*8 yl(npro,nshl,ndofl), acl(npro,nshl,ndofl), local
164 call localy(y, yl, ien, ndofl, 'gather ')
180 call e3Sclr (yl, acl, shp,
H A DasiqGradV.f30 dimension yl(npro,nshl,ndof), xl(npro,nenl,nsd),
50 call localy(y, yl, ien, ndof, 'gather ')
57 call e3qGradV (yl, shp, shgl,
H A De3ql.f1 subroutine e3ql (yl, dwl, shp, shgl, argument
25 dimension yl(npro,nshl,ndof), dwl(npro,nshl),
61 call e3qvar (yl, shdrv,
66 call getdiff(dwl, yl, shape, xmudmi, xl,rmu, tmp)
177 subroutine e3qlSclr (yl, dwl, shp, shgl, argument
191 dimension yl(npro,nshl,ndof), dwl(npro,nshl),
224 call e3qvarSclr (yl, shdrv, xl,
229 call getdiffsclr(dwl,shape,yl, diffus)
H A De3bvar.f1 subroutine e3bvar (yl, acl, ul, argument
51 dimension yl(npro,nshl,ndof), rmu(npro),
128 pres = pres + shpb(:,nodlcl) * yl(:,nodlcl,1)
129 u1 = u1 + shpb(:,nodlcl) * yl(:,nodlcl,2)
130 u2 = u2 + shpb(:,nodlcl) * yl(:,nodlcl,3)
131 u3 = u3 + shpb(:,nodlcl) * yl(:,nodlcl,4)
287 gl1yi(:,1) = gl1yi(:,1) + shglb(:,1,n) * yl(:,n,1)
288 gl1yi(:,2) = gl1yi(:,2) + shglb(:,1,n) * yl(:,n,2)
289 gl1yi(:,3) = gl1yi(:,3) + shglb(:,1,n) * yl(:,n,3)
290 gl1yi(:,4) = gl1yi(:,4) + shglb(:,1,n) * yl(:,n,4)
[all …]
H A De3qGradV.f1 subroutine e3qGradV (yl, shp, shgl, argument
23 dimension yl(npro,nshl,ndof),
54 call e3qvar (yl, shdrv,
H A Dfilters.f11 dimension ien(npro,nshl), yl(npro,nshl,5),
25 call local (y, yl, ien, 5, 'gather ')
116 u1 = u1 + shp(i,intp)*yl(:,i,2)
117 u2 = u2 + shp(i,intp)*yl(:,i,3)
118 u3 = u3 + shp(i,intp)*yl(:,i,4)
120 S11 = S11 + shg(:,i,1)*yl(:,i,2)
121 S22 = S22 + shg(:,i,2)*yl(:,i,3)
122 S33 = S33 + shg(:,i,3)*yl(:,i,4)
124 S12 = S12 + shg(:,i,2)*yl(:,i,2)
125 & +shg(:,i,1)*yl(:,i,3)
[all …]
H A Dasbflx.f27 dimension yl(npro,nshl,ndofl), xlb(npro,nenl,nsd),
49 call localy(y, yl, ienb, ndofl, 'gather ')
62 call e3b (ul, yl, acl, iBCB, BCB,
/phasta/phSolver/common/
H A Dgetstrl.f7 dimension ien(npro,nshl), yl(npro,nshl,5),
15 call localy (y, yl, ien, 5, 'gather ')
20 yl (:,:,1) = yl(:,:,1) / (Rgas * yl(:,:,5))
22 yl (:,:,1) = one
78 fresli(:,22) = fresli(:,22)+shp(i,intp)*yl(:,i,1) ! density at qpt
103 fresli(:,j) = fresli(:,j)+shg(:,i,ig)*yl(:,i,iv)
111 fresli(:,13) = fresli(:,13)+shg(:,i,2)*yl(:,i,2)
112 & +shg(:,i,1)*yl(:,i,3)
113 fresli(:,14) = fresli(:,14)+shg(:,i,3)*yl(:,i,2)
114 & +shg(:,i,1)*yl(:,i,4)
[all …]
H A Dbardmc.f678 dimension ien(npro,nshl), yl(npro,nshl,nflow),
689 call local (y, yl, ien, 5, 'gather ')
694 yl (:,:,1) = yl(:,:,1) / (Rgas * yl(:,:,5)) !get density
696 yl (:,:,1) = one
756 rho = rho+shp(i,intp)*yl(:,i,1) !density at qpt
765 u1 = u1 + shp(i,intp)*yl(:,i,2)
766 u2 = u2 + shp(i,intp)*yl(:,i,3)
767 u3 = u3 + shp(i,intp)*yl(:,i,4)
832 dimension ien(npro,nshl), yl(npro,nshl,nflow),
844 call local (y, yl, ien, 5, 'gather ')
[all …]
/phasta/phSolver/incompressible/e3source/
H A De3source.fimplicit212 & yl, dxidx, rmu,
242 & yl(npro,nenl,ndof), dxidx(npro,nsd,nsd),
289 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * yl(:,n,2)
290 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * yl(:,n,3)
291 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * yl(:,n,4)
293 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * yl(:,n,2)
294 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * yl(:,n,3)
295 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * yl(:,n,4)
297 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * yl(:,n,2)
298 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * yl(:,n,3)
[all …]
H A De3source.fRHSonly212 & yl, dxidx, rmu,
242 & yl(npro,nenl,ndof), dxidx(npro,nsd,nsd),
289 gradV(:,1,1) = gradV(:,1,1) + shg(:,n,1) * yl(:,n,2)
290 gradV(:,2,1) = gradV(:,2,1) + shg(:,n,1) * yl(:,n,3)
291 gradV(:,3,1) = gradV(:,3,1) + shg(:,n,1) * yl(:,n,4)
293 gradV(:,1,2) = gradV(:,1,2) + shg(:,n,2) * yl(:,n,2)
294 gradV(:,2,2) = gradV(:,2,2) + shg(:,n,2) * yl(:,n,3)
295 gradV(:,3,2) = gradV(:,3,2) + shg(:,n,2) * yl(:,n,4)
297 gradV(:,1,3) = gradV(:,1,3) + shg(:,n,3) * yl(:,n,2)
298 gradV(:,2,3) = gradV(:,2,3) + shg(:,n,3) * yl(:,n,3)
[all …]

12