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

Go to the documentation of this file.

subroutine preprocessor


use global

!use preprocessorParameters

implicit none

integer i,j,k,t

integer,pointer :: IM(:,:,:,:)

integer,pointer :: Node(:,:,:)

integer,pointer :: mytotalNode

real(8) dx,dy


integer startIndex,endIndex

! mesh will be called
! models initialization will be called - i.e. number of unk, grid, indexmatrix bla bla
! also bc's will be flagged here


integer,pointer :: allBC

integer dummy

integer,pointer :: minI,minJ,maxI,maxJ

modelDof=>myModel%totalDof

modelDof=0


modelNode=>myModel%totalNode

modelNode=0



do i=1,numberOfPhysics !myModel%numberOfPhysics

        currentPhysics=i

        numberOfDomains=>myModel%physics(currentPhysics)%numberOfDomains

        physicsDof=>myModel%physics(currentPhysics)%totalDof
        
        do j=1,numberOfDomains

                currentDomain=j

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

                allBC=>myModel%physics(currentPhysics)%domain(currentDomain)%allBC

                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
        

                !print*,size(im,dim=1)
                !stop


                call generateIndexMatrix()
        

                !do k=myModel%physics(currentPhysics)%domain(currentDomain)%maxJ,myModel%physics(currentPhysics)%domain(currentDomain)%minJ,-1

                        !print*,IM(:,k,3,absolute)

                !end do
                
                !allocate

                !stop



                call generateIndexVector() ! might be activated only set is on

                
                mytotalNode=>myModel%physics(currentPhysics)%domain(currentDomain)%totalNode


                allocate(myModel%physics(currentPhysics)%domain(currentDomain)%xCoor(mytotalNode))
                allocate(myModel%physics(currentPhysics)%domain(currentDomain)%yCoor(mytotalNode))

                

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


                dx=Lx/(myModel%physics(currentPhysics)%domain(currentDomain)%nx-1.d0)
                dy=Ly/(myModel%physics(currentPhysics)%domain(currentDomain)%ny-1.d0)

                !print*,dx,dy,Lx,Ly


                call generateCoordinates(dx,dy) !<------------ should be modified if not simple structured - probably user subroutine should be implemented
                                                                                !i.e. generateCoordinates will be renamed as generateRectangular Map, and user defined generateCoordinates will call
                                                                                !either generateRectangular Map as a demona function or userdefined subroutine



                call addBC()    !<------------ should be modified to cover all physics and ddm option


                !< print*,'at the initialization step we need full set'
                !IM(:,:,:,flag)=1


                !FORALL(I = 1:N, J = 1:N, A(I, J) .NE. 0.0)
                !FORALL(i = minI:maxI, j = minJ:maxJ, t=1:physicsDof, IM(i,j,t,absolute).NE.0) IM(i,j,t,flag)=1
                FORALL(i = minI:maxI, j = minJ:maxJ, t=1:myModel%physics(currentPhysics)%dof, IM(i,j,t,absolute).ne.0) IM(i,j,t,flag)=1

                !print*,minj,maxj
                !stop


                !do k=maxJ,minJ,-1

                !       print*,IM(:,k,1,flag)

                !end do

                !stop


!print*,globalCorrection

                call set

!               print*,globalCorrection
!pause



        end do

        startIndex=modelDof+1

        endIndex=startIndex+physicsDof-1

        !modelDof=physicsDof+modelDof



end do
!< at this point, totalDof of the domains as well as the physics are calculated.


!< correct absolute values of the 


!do i=1,18

!print*,mymodel%physics(2)%domain(1)%xcoor(i),mymodel%physics(2)%domain(1)%ycoor(i)

!end do

!stop






!< allocate x and tempx
allocate(myModel%x(modelDof),myModel%tempx(modelDof),myModel%xPTO(modelDof),myModel%xRTO(modelDof))

y=>myModel%x

y=0.d0

yPTO=>myModel%xPTO

yPTO=0.d0

yRTO=>myModel%xRTO

yRTO=0.d0



tempy=>myModel%tempx

tempy=0.d0

discreteY=>tempy


dummy=0

!< assign x of physics
do i=1,numberOfPhysics !myModel%numberOfPhysics

        currentPhysics=i

        physicsDof=>myModel%physics(currentPhysics)%totalDof

        startIndex=dummy+1

        endIndex=startIndex+physicsDof-1

        dummy=physicsDof+dummy

        myModel%physics(currentPhysics)%x=>discreteY(startIndex:endIndex)

        myModel%physics(currentPhysics)%xPTO=>yPTO(startIndex:endIndex)

        myModel%physics(currentPhysics)%xRTO=>yRTO(startIndex:endIndex)


end do


!discreteY=>myModel%tempx

!< assign x of domains
do i=1,numberOfPhysics !myModel%numberOfPhysics

        currentPhysics=i

        physicsDof=>myModel%physics(currentPhysics)%totalDof

        numberOfDomains=>myModel%physics(currentPhysics)%numberOfDomains


        dummy=0

        do j=1,numberOfDomains

                currentDomain=j

                startIndex=dummy+1

                domainDof=>myModel%physics(currentPhysics)%domain(currentDomain)%totalDof

                endIndex=startIndex+domainDof-1
                
                
                myModel%physics(currentPhysics)%domain(currentDomain)%x=>myModel%physics(currentPhysics)%x(startIndex:endIndex)

                myModel%physics(currentPhysics)%domain(currentDomain)%xPTO=>myModel%physics(currentPhysics)%xPTO(startIndex:endIndex)

                myModel%physics(currentPhysics)%domain(currentDomain)%xRTO=>myModel%physics(currentPhysics)%xRTO(startIndex:endIndex)

                dummy=domainDof+dummy


        end do

end do

discreteY=>y


deallocate(mymodel%pVec)
nullify(pVEc)







end subroutine preprocessor

---