D:/demona/2.0/manual/doc/2-31set/problems/test/2/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

integer bcValue


!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



IM(:,:,:,boundary)=0


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


select case(currentPhysics)


        case(1)

                select case(currentDomain)

                
                        case(1)

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

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

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

                        case(2)

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

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

                                !3-left
                                bcValue=3
                                IM(1,:,:,boundary)=bcValue

                                !4-right
                                bcValue=4
                                IM(nx,:,:,boundary)=bcValue





                        case default

                        print*,'no bc definition'
                        print*,'demona stop in addBC'
                        stop

                end select


        case(2)


                select case(currentDomain)

                
                        case(1)

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

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

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

                                !4-right
                                bcValue=4
                                IM(nx,2:ny-1,:,boundary)=bcValue


                        case default

                        print*,'no bc definition'
                        print*,'demona stop in addBC'
                        stop

                end select



        case(3)

                select case(currentDomain)

                
                        case(1)

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

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

                                !3-left
                                bcValue=3
                                IM(1,:,:,boundary)=bcValue

                                !4-right
                                bcValue=4
                                IM(nx,:,:,boundary)=bcValue

                        case(2)

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

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

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



                        case default

                        print*,'no bc definition'
                        print*,'demona stop in addBC'
                        stop

                end select


end select



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




end subroutine addBC

---