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

Go to the documentation of this file.

subroutine physicsInfo


use global


implicit none



numberOfDomains=>myModel%physics(currentPhysics)%numberOfDomains
dof=>myModel%physics(currentPhysics)%dof
dFX=>myModel%physics(currentPhysics)%dFX


refXcoor=>myModel%physics(currentPhysics)%refXcoor
refYcoor=>myModel%physics(currentPhysics)%refYcoor

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

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

        case(1)

                numberOfDomains=2
                dof=2
                !dof=1
                dFX=1
                !TecplotVariables='VARIABLES="x","y","T"'
                TecplotVariables='VARIABLES="x","y","T","U"'
                !TecplotVariables='VARIABLES="x","y","str","vor"'
                refXcoor=0.d0
                refYcoor=0.d0
                Lx=1.d0
                Ly=2.5d0

                allocate(myModel%physics(currentPhysics)%damping(dof))
                damping=>myModel%physics(currentPhysics)%damping
                damping(1)=1.0
                damping(2)=1.0

                

        case(2)
                numberOfDomains=1
                dof=2
                dFX=1
                !TecplotVariables='VARIABLES="x","y","T"'
                TecplotVariables='VARIABLES="x","y","T","U"'
                refXcoor=2.5d0
                refYcoor=0.d0
                Lx=0.5d0
                Ly=1.25d0


                allocate(myModel%physics(currentPhysics)%damping(dof))
                damping=>myModel%physics(currentPhysics)%damping
                damping(1)=1.0
                damping(2)=1.0

        case(3)
                numberOfDomains=2
                dof=2
                dFX=1
                !TecplotVariables='VARIABLES="x","y","str","vor"'
                TecplotVariables='VARIABLES="x","y","T","U"'
                refXcoor=3.25d0
                refYcoor=0.25d0
                Lx=0.5d0
                Ly=1.25d0

                allocate(myModel%physics(currentPhysics)%damping(dof))
                damping=>myModel%physics(currentPhysics)%damping
                damping(1)=1.0
                damping(2)=1.0
        
        !< the format below could also been used
        !case(1:2)
        !       numberOfDomains=1
        !       dof=1
        !       dFX=1

        case default

        print*,'invalid physics definition'
        print*,'demona stop'

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

!#####################################################################
!user supplied variables
!#####################################################################
! if different pVec are to be used, something similar to the code lines
! below should be implemented.
!deallocate(myModel%pVec)
!allocate(myModel%pVec(*))
!pVec(1)=2
!pVec(2)=5
!#####################################################################


myModel%physics(currentPhysics)%tecVariables=TecplotVariables

allocate(myModel%physics(currentPhysics)%endIndex(numberOfDomains))
myModel%physics(currentPhysics)%endIndex=1




end subroutine physicsInfo

---