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