D:/demona/2.0/manual/doc/2-31set/problems/test/1/addBC.f90

Go to the documentation of this file.

subroutine addBC
!< bc is a misnomer. it is used to add flexibility in defining the boundary conditions
!< however it is also used to represent interior. Interior always starts with 1
!< addBC is a user supplied routine
!< future work is to create a gallery considering the examples below and mant more

! most of the time, 9 definitions are enough sample cases are

!2 - example: bratu problem

!       2       2       2       2       2 

!       2       1       1       1       2

!       2       1       1       1       2

!       2       1       1       1       2

!       2       2       2       2       2


!5 - example: laplace with 4 different bc's

!       3       3       3       3       3 

!       4       1       1       1       5

!       4       1       1       1       5

!       4       1       1       1       5

!       2       2       2       2       2


!9 - example: 5 + 4 corner discretization

!       8       3       3       3       9 

!       4       1       1       1       5

!       4       1       1       1       5

!       4       1       1       1       5

!       6       2       2       2       7


!9 - example: finite volume incompressible flow - quick with two layers ghost cells.

!       0       0       7       7       7       0       0 

!       0       0       3       3       3       0       0 

!       8       4       1       1       1       5       9

!       8       4       1       1       1       5       9

!       8       4       1       1       1       5       9

!       0       0       2       2       2       0       0

!       0       0       6       6       6       0       0




use global


implicit none

integer,pointer :: IM(:,:,:,:)
integer,pointer :: nx,ny




!should also be physics and domain dependent
! also consider artificial bc's.



IM=>myModel%physics(currentPhysics)%domain(currentDomain)%DofMatrix

nx=>myModel%physics(currentPhysics)%domain(currentDomain)%nx
ny=>myModel%physics(currentPhysics)%domain(currentDomain)%ny




!#####################################################################
!user supplied variables
!#####################################################################



select case(currentPhysics)


        case(1)

                !1-interior
                IM(:,:,:,boundary)=1

                !2-all boundary except right
                IM(1,1:ny,:,boundary)=2
                IM(2:nx-1,1,:,boundary)=2
                IM(2:nx-1,ny,:,boundary)=2

                !3-right
                IM(nx,1:ny,:,boundary)=3

        case(2)

                !1-interior
                IM(:,:,:,boundary)=1

                !2-all boundary except right
                IM(nx,1:ny,:,boundary)=2
                IM(1:nx-1,1,:,boundary)=2
                IM(1:nx-1,ny,:,boundary)=2

                !3-left
                IM(1,2:ny-1,absolute,boundary)=3


end select



!#####################################################################




end subroutine addBC

---