1*59599516SKenneth E. Jansen subroutine asbwmod (y, ac, x, BC, iBC, 2*59599516SKenneth E. Jansen & iper, ilwork, ifath, velbar) 3*59599516SKenneth E. Jansenc 4*59599516SKenneth E. Jansenc---------------------------------------------------------------------- 5*59599516SKenneth E. Jansenc 6*59599516SKenneth E. Jansenc This routine assembles traction BCs for a modeled wall 7*59599516SKenneth E. Jansenc 8*59599516SKenneth E. Jansenc---------------------------------------------------------------------- 9*59599516SKenneth E. Jansenc 10*59599516SKenneth E. Jansen use pointer_data 11*59599516SKenneth E. Jansen include "common.h" 12*59599516SKenneth E. Jansenc 13*59599516SKenneth E. Jansen dimension y(nshg,ndof), x(numnp, nsd), 14*59599516SKenneth E. Jansen & BC(nshg,ndofBC), iBC(nshg), 15*59599516SKenneth E. Jansen & iper(nshg), ilwork(nlwork), 16*59599516SKenneth E. Jansen & ifath(numnp), velbar(nfath,nflow), 17*59599516SKenneth E. Jansen & ac(nshg,ndof) 18*59599516SKenneth E. Jansenc 19*59599516SKenneth E. Jansenc.... compute and assemble the residuals corresponding to the 20*59599516SKenneth E. Jansenc boundary integral 21*59599516SKenneth E. Jansenc 22*59599516SKenneth E. Jansen call settauw (y, x, 23*59599516SKenneth E. Jansen & BC, 24*59599516SKenneth E. Jansen & ifath, velbar) 25*59599516SKenneth E. Jansenc 26*59599516SKenneth E. Jansenc.... enforce the new BC for SA variable 27*59599516SKenneth E. Jansenc 28*59599516SKenneth E. Jansen isclr = 1 29*59599516SKenneth E. Jansen if (iRANS.eq.-1) then ! S-A RANS 30*59599516SKenneth E. Jansen call itrBCSclr (y, ac, iBC, BC, iper, ilwork) 31*59599516SKenneth E. Jansen endif 32*59599516SKenneth E. Jansenc 33*59599516SKenneth E. Jansen return 34*59599516SKenneth E. Jansen end 35