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

real(8),pointer :: relXcoor,relYcoor

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

relXcoor=>myModel%physics(currentPhysics)%domain(currentDomain)%relXcoor
relYcoor=>myModel%physics(currentPhysics)%domain(currentDomain)%relYcoor

Lx=>myModel%physics(currentPhysics)%domain(currentDomain)%Lx
Ly=>myModel%physics(currentPhysics)%domain(currentDomain)%Ly


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

        case(1)

                
                select case(currentDomain)


                        case(1)

                                nx=33
                                ny=33
                                !nx=17
                                !ny=17
                                !nx=16
                                !ny=11
                                allBC=7

                                relXcoor=0.d0
                                relYcoor=0.d0

                                !Lx=1.d0
                                !Ly=1.d0

                                !or 

                                Lx=myModel%physics(currentPhysics)%Lx
                                Ly=myModel%physics(currentPhysics)%Ly


                                !numberOfGhostCells=2


                

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


                end select




        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

---