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

Go to the documentation of this file.

subroutine domainInfo


use global


implicit none


integer,pointer :: minI,minJ,maxI,maxJ

integer,pointer :: nx,ny

integer,pointer :: allBC

minI=>myModel%physics(currentPhysics)%domain(currentDomain)%minI
maxI=>myModel%physics(currentPhysics)%domain(currentDomain)%maxI
minJ=>myModel%physics(currentPhysics)%domain(currentDomain)%minJ
maxJ=>myModel%physics(currentPhysics)%domain(currentDomain)%maxJ


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

numberOfGhostCells=>myModel%physics(currentPhysics)%domain(currentDomain)%numberOfGhostCells

allBC=>myModel%physics(currentPhysics)%domain(currentDomain)%allBC


!#####################################################################
!user supplied variables
!#####################################################################
select case(currentPhysics)

        case(1)

                
                select case(currentDomain)


                        case(1)

                                nx=5
                                ny=5
                                allBC=3


                                !numberOfGhostCells=2


                        case(2)


                                nx=3
                                ny=3
                                allBC=9


                                !numberOfGhostCells=2


                        case default
                        
                        print*,'!!!!!!!!!!!!!!!!!!!!!!!!!'
                        print*,'invalid domain definition'
                        print*,'demona stop:domainInfo'
                        stop


                end select



        case(2)


                select case(currentDomain)

                        case(1)

                                nx=3
                                ny=6
                                allBC=3


                        !case(2)
                        
                        
                        case default
                        
                        print*,'!!!!!!!!!!!!!!!!!!!!!!!!!'
                        print*,'invalid domain definition'
                        print*,'demona stop:domainInfo'
                        stop




                end select





        !< the format below could also been used
        !case(1:2)
        !       numberOfDomains=1
        !       dof=1


        case(3)


                select case(currentDomain)

                        case(1)

                                nx=4
                                ny=4
                                allBC=2


                        !case(2)
                        
                        
                        case default
                        
                        print*,'!!!!!!!!!!!!!!!!!!!!!!!!!'
                        print*,'invalid domain definition'
                        print*,'demona stop:domainInfo'
                        stop


                end select


        !< the format below could also been used
        !case(1:2)
        !       numberOfDomains=1
        !       dof=1


        case default

        print*,'!!!!!!!!!!!!!!!!!!!!!!!!!'
        print*,'invalid physics definition'
        print*,'demona stop:domainInfo'
        stop

end select
!#####################################################################

minI=1
minJ=1
maxI=nx
maxJ=ny


if (dFX.eq.dFV) then

        minI=minI-numberOfGhostCells
        minJ=minJ-numberOfGhostCells

        maxI=maxI+numberOfGhostCells
        maxJ=maxJ+numberOfGhostCells

end if














end subroutine domainInfo

---