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