D:/demona/2.0/manual/doc/2-31set/preprocessor/generateIndexVector.f90

Go to the documentation of this file.

subroutine generateIndexVector

use global

implicit none

integer sumDof

integer i,j,t

integer,pointer :: IM(:,:,:,:)
integer,pointer :: myDofVectorArray(:,:)
integer,pointer :: minI,minJ,maxI,maxJ

IM=>mymodel%physics(currentPhysics)%domain(currentDomain)%dofMatrix




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


allocate(mymodel%physics(currentPhysics)%domain(currentDomain)%dofVector(size(IM)/4,4))

myDofVectorArray=>mymodel%physics(currentPhysics)%domain(currentDomain)%dofVector



myDofVectorArray=0

sumDof=1


dof=>myModel%physics(currentPhysics)%dof


!print*,dof,currentPhysics


do j=minJ,maxJ

        do i=minI,maxI

                ! skip function here
                ! if skip_is_yes cycle
                                
                do t=1,dof

                        if (IM(i,j,t,relative).ne.0) then
                                        
                                myDofVectorArray(sumDof,1)=IM(i,j,t,absolute)
                                myDofVectorArray(sumDof,2)=i
                                myDofVectorArray(sumDof,3)=j
                                myDofVectorArray(sumDof,4)=t

                                sumDof=sumDof+1

                        end if

                end do

        end do

end do


end subroutine generateIndexVector

---