1c========================================================== 2c Initialize the jet mdot when setting the jet mdot 3c dynamically based on the contraction mdot, accompanied 4c with setBlowing_Duct2.f 5c========================================================== 6 7 subroutine setBlowing_Duct(BC,iTurbWall) 8 9 use blowingDuct 10 use readarrays ! for qold 11 include "common.h" 12 include "mpif.h" 13 include "auxmpi.h" 14 integer i, nn 15 real*8 BC(nshg,ndofBC) 16 real*8 yVel, Temp 17 integer isfID 18 integer iTurbWall(nshg) 19 20 allocate(jetinletf(nshg)) 21 isfID=isetBlowing_Duct 22 call sfID2np(isfID,njetinlet,jetinletf) 23c... above, find jet inlet nodes 24 25 if(njetinlet .gt. 0)then 26 do i=1,njetinlet 27 nn=jetinletf(i) 28c... from IC (p u v w T nu) 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 35 endif 36 enddo 37 endif 38c... above, take BC from IC, before setting blowing vel dynamically in the time step loop 39 40 return 41 end 42 43