Unverified Commit e2e94af5 authored by Stefan K's avatar Stefan K Committed by GitHub
Browse files

add boundary for t-mesh (#640)

parent 8958b38e
......@@ -90,5 +90,5 @@ c
$ ,ifmseg,ifmscr,ifnskp
$ ,ifbcor
integer boundaryID(6,lelt)
common /cbbid/ boundaryID
integer boundaryID(6,lelv), boundaryIDt(6,lelt)
common /cbbid/ boundaryID, boundaryIDt
......@@ -2164,13 +2164,22 @@ c
$ 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.bid)
$ cbc(ifc,iel,ifld) = cbc_bmap(bid,ifld)
enddo
enddo
if (iftmsh(ifld)) then
do iel = 1,nelt
do ifc = 1,2*ndim
if (boundaryIDt(ifc,iel).eq.bid)
$ cbc(ifc,iel,ifld) = cbc_bmap(bid,ifld)
enddo
enddo
else
do iel = 1,nelv
do ifc = 1,2*ndim
if (boundaryID(ifc,iel).eq.bid)
$ cbc(ifc,iel,ifld) = cbc_bmap(bid,ifld)
enddo
enddo
endif
return
end
......@@ -107,11 +107,27 @@ C End of input data, close read file.
$ dnekclock()-etime0,' sec'
endif
99 do iel = 1,nelt
99 call izero(boundaryID, size(boundaryID))
call izero(boundaryIDt, size(boundaryIDt))
do iel = 1,nelv
do ifc = 1,2*ndim
boundaryID(ifc,iel) = bc(5,ifc,iel,1)
enddo
enddo
enddo
ntmsh = 0
do i=1,ldimt
if(iftmsh(1+i)) ntmsh = ntmsh + 1
enddo
if (ntmsh.gt.0) then
do iel = 1,nelt
do ifc = 1,2*ndim
boundaryIDt(ifc,iel) = bc(5,ifc,iel,2)
enddo
enddo
endif
return
END
......
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