D:/demona/2.0/manual/doc/2-31set/processor/activatePhysics.f90

Go to the documentation of this file.

subroutine activatePhysics(vec,maxVec)

use global


implicit none

integer i

integer maxVec

integer vec(maxVec)

integer myP,p

integer startIndex,endIndex

integer currentXLength

currentXLength=0


do myP=1,maxVec
        
        p=Vec(myP)
        
        startIndex=currentXLength+1

        endIndex=startIndex+myModel%physics(p)%totalDof-1

        currentXLength=endIndex

end do


print*,currentXlength


allocate(setY(currentXlength))


currentXLength=0


!print*,myModel%physics(1)%totalDof
!print*,myModel%physics(2)%totalDof


!stop




do myP=1,maxVec
        
        p=Vec(myP)
        
        startIndex=currentXLength+1

        endIndex=startIndex+myModel%physics(p)%totalDof-1

        do i=startIndex,endIndex


                setY(i)%y=>myModel%physics(p)%x(i-currentXLength)

        end do

        currentXLength=endIndex

end do





print*,setY(1)%y,myModel%physics(p)%x(i-currentXLength)

stop



call printY(1,currentXLength)






end subroutine activatePhysics

---