D:/demona/2.0/manual/doc/2-31set/problems/sf_vor/1/solver/physicsRule.f90

Go to the documentation of this file.

subroutine physicsRule

use global

use model

!use mf

implicit none


integer i

integer myP,myD

real(8) myEps

real(8),allocatable :: testVec(:)






print*,pVec

print*,'now decide on the phyics rule'




stop





do myP=1,size(pVec)

        currentPhysics=pVec(myP)

        numberOfDomains=>myModel%physics(currentPhysics)%numberOfDomains

        select case(currentPhysics)


                case default

                        allocate(myModel%pg(currentPhysics)%dg(1))

                        allocate(myModel%pg(currentPhysics)%dg(1)%dVec(numberOfDomains))

                        do i=1,numberOfDomains

                                myModel%pg(currentPhysics)%dg(1)%dVec(i)=i


                        end do


        end select


end do







maxDomainGroup=>myModel%physics(currentPhysics)%numberOfDomainGroups


maxDomainGroup=size(myModel%pg(currentPhysics)%dg)





end subroutine physicsRule

---