| /phasta/phSolver/common/ |
| H A D | genbc.f | 105 nn=idxslpw(i)
106 call slpwBC(mpslpw(nn,1),mpslpw(nn,2),iBCg(nn),
107 & BCg(nn,:), BCtmpg(nn,:),
108 & iBC(nn), BC(nn,:),
109 & wlnorm(nn,:,:) )
129 nn = 0
132 nn = nn + 1
133 if(mod(nn,50).eq.1)
258 nn=ienb(j) ! global node number
259 if(firstvisit(nn)) then
[all …]
|
| H A D | setSuction_Duct3.f90 | 12 integer i, j, k, nn, nL local
48 nn = suctionNodeMap(i) ! map face node index to global node
49 xcoor = x(nn,1)
50 ycoor = x(nn,2)
54 wnormy = wnorm(nn, 2)
55 wnormz = wnorm(nn, 3)
62 BC(nn, 3) = nvel_su*wnorm(nn, 1) ! set xVel
63 BC(nn, 4) = nvel_su*wnorm(nn, 2) ! set yVel
64 BC(nn, 5) = nvel_su*wnorm(nn, 3) ! set zVel
66 BC(nn, 3) = nvel_sl*wnorm(nn, 1) ! set xVel
[all …]
|
| H A D | setBlowing_Duct3.f | 10 integer i, nn local
38 nn=jetinletf(i)
39 if(iTurbWall(nn).eq.0)then
42 xcoor=x(nn,1)
43 ycoor=x(nn,2)
44 zcoor=x(nn,3)
68 BC(nn,2) = Temp ! Temp
69 BC(nn,3) = veloMag*jetnx ! set x velocity
70 BC(nn,4) = veloMag*jetny ! set y velocity
71 BC(nn,5) = 0 ! set zVel=0
[all …]
|
| H A D | mod_wallData.f90 | 60 nn=ienb(j) !nn is the global index of suction node
62 aedg(:)=x(ienb(j+1),:)-x(nn,:)
63 bedg(:)=x(ienb(j-1),:)-x(nn,:)
65 aedg(:)=x(ienb(j+1),:)-x(nn,:)
66 bedg(:)=x(ienb(nenbl),:)-x(nn,:)
68 aedg(:)=x(ienb(1),:)-x(nn,:)
69 bedg(:)=x(ienb(j-1),:)-x(nn,:)
75 call addElemNormal(wNorm, eNorm, nn)
107 do nn = 1,nshg
108 wn(:) = wNorm(nn,:)
[all …]
|
| H A D | findslpw.f | 39 nn=ienb(j)
40 if(ifirstvisit(nn).eq.1)then
41 ifirstvisit(nn)=0
43 idxslpw(nslpwnd)=nn
45 if(.not.btest(mpslpw(nn,2),islpw))then
46 mpslpw(nn,1)=mpslpw(nn,1)+1
47 mpslpw(nn,2)=mpslpw(nn,2)+2**islpw
50 aedg(:)=x(ienb(j+1),:)-x(nn,:)
51 bedg(:)=x(ienb(j-1),:)-x(nn,:)
53 aedg(:)=x(ienb(j+1),:)-x(nn,:)
[all …]
|
| H A D | setBlowing_Duct.f | 14 integer i, nn local
27 nn=jetinletf(i)
29 if (iTurbWall(nn).eq.0)then
30 BC(nn,2) = qold(nn,5) ! set Temp
31 BC(nn,3) = 0 ! set xVel=0
32 BC(nn,4) = qold(nn,3) ! set and scale y velocity
33 BC(nn,5) = 0 ! set zVel=0
34 BC(nn,7) = qold(nn,6) ! set SA variable
|
| H A D | INIprofile.f | 37 nn=inlf(i)
38 if(y(nn,1).gt.u_max)u_max=y(nn,1)
109 do nn=1,nshg
110 xdel=(x(nn,1)-xinlet_all)/(xoutlet_all-xinlet_all)
112 y(nn,1)= xvel_ini + (inletVelX-xvel_ini)*xdel
113 y(nn,4)= outPres1 + (pres_ini-outPres1)*xdel
176 integer nn local
201 do nn=1,nshg
203 xcoor=x(nn,1)
204 ycoor=x(nn,2)
[all …]
|
| H A D | BCprofile2.f | 19 nn=BCfaceNode(i)
20 BC(nn,1)=outPres1
36 integer i,nn local
51 nn=BCfaceNode(i)
52 if(iTurbWall(nn).eq.0)then
54 xcoor=x(nn,1)
55 ycoor=x(nn,2)
56 zcoor=x(nn,3)
86 BC(nn,2) = Temp !Temp
87 BC(nn,3) = xVel ! set and scale x velocity
[all …]
|
| H A D | proces.f | 50 integer nn local
182 do nn = 1,nshg
183 … if(y(nn,6) .eq. 0) cycle !don't change wall boundary conditions, should be iTurbWall == 1
185 if(point2x(nn,1) .gt. EVrampXmax) then !downstream of the ramp
186 y(nn,6) = EVrampMax
187 elseif(point2x(nn,1) .gt. EVrampXmin) then !and x(:,1) <= EVrampXmax
191 y(nn,6) = c0 + (point2x(nn,1) - x1)*(
192 & c1 + (point2x(nn,1) - x1)*(
193 & (c2 - (x2 - x1)*c3) + (point2x(nn,1) - x1)*c3))
195 y(nn,6) = EVrampMin
|
| H A D | sfID2np.f | 39 …nn=ienb(j) !nn is a global node index. The array ienb is probably the global indeces of the nod…
40 … if(ifirstvisit(nn).eq.1)then !need this check because nodes are shared by other elements
41 ifirstvisit(nn)=0
43 …nmap(nsfnp)=nn !nmap is the mapping from nodal surface index to the global nodal numbering for o…
|
| H A D | genibc.f | 52 nn = 0
56 nn = nn + 1
57 if (mod(nn,50).eq.1) write(iecho,1000)ititle,(j,j=1,ndof)
|
| H A D | findTurbWall.f | 24 nn=ienb(j)
25 iTurbWall(nn)=1
|
| /phasta/phSolver/compressible/ |
| H A D | setBlowing_Duct2.f | 14 integer i, nn local
68 nn=jetinletf(i)
69 if(iTurbWall(nn).eq.0)then
82 BC(nn,2) = 317 ! Temp
83 BC(nn,4) = blowingVel ! set and scale y velocity
84 BC(nn,3) = 0
85 BC(nn,5) = 0
86 BC(nn,7) = 1.825e-5
|