Commit 7c4dc92c authored by Stefan's avatar Stefan
Browse files

add cbc_mp

parent ef9a5f8a
......@@ -104,9 +104,13 @@ c
integer igroup(lelt),object(maxobj,maxmbr,2)
common /input6/ igroup,object
integer lbid
parameter(lbid = 100)
character*1 ccurve(12,lelt),cdof(6,lelt)
character*3 cbc(6,lelt,0:ldimt1)
common /input8/ cbc,ccurve,cdof
character*3 cbc_bmap(lbid,ldimt1)
common /input8/ cbc,ccurve,cdof,cbc_bmap
integer ieact(lelt),neact
common /input9/ ieact,neact
......
......@@ -2149,7 +2149,7 @@ c write(6,1) iobj,mem,f,ieg,e,nid,' OBJ'
return
end
c-----------------------------------------------------------------------
subroutine setbc(sid,ifld,cbci)
subroutine setbc(bid,ifld,cbci)
c
c sets boundary condition for a given surface id and field
c
......@@ -2158,11 +2158,17 @@ c
include 'GEOM'
character*3 cbci
integer sid
integer bid
if (bid.lt.1 .or. bid.gt.lbid)
$ call exitti('invalid boundary id!$',bid)
cbc_bmap(bid,ifld) = cbci
do iel = 1,nelt
do ifc = 1,2*ndim
if (boundaryID(ifc,iel).eq.sid) cbc(ifc,iel,ifld) = cbci
if (boundaryID(ifc,iel).eq.bid)
$ cbc(ifc,iel,ifld) = cbc_bmap(bid,ifld)
enddo
enddo
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment