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

Go to the documentation of this file.

subroutine tecFD


! instead of global or any other use methods
! add only use postModule and postModuel will use
! necessary modules
use global
use myFunctions


implicit none


integer d

integer i,j,k,t

integer myNodes,myDof,myNz,myNy,myNx,myNzCells,myNyCells,myNxCells,myLocalDof

integer myPhysics

integer numberOfElements,numberOfElementalNodes

integer gridInfo(2)

integer sum

integer axis

real(8), allocatable :: resultVec(:,:)   ! what if one variable is stored onlt at the centers, then results will vary.


integer, allocatable :: ConnectivityList(:,:)

character(len=36) :: file

logical,pointer :: realTime


!# dummy variables to be deleted

real(8) dum1,dum2




! strand id
! zones  (same physics domains with zone designation and shared x y data)
! if coordiates change then we need to add new x and y



! later a loop will be used as such:
! do myPhysics=1,totalPhysics  

! end do






myPhysics=currentPhysics




TecplotVariables=myModel%physics(myPhysics)%tecVariables


!print*,TecplotVariables

!stop




realTime=>myModel%physics(myPhysics)%realTime



do d=1,myModel%physics(currentPhysics)%numberOfDomains

        !print*,d



        ! what if we call tecFD for sets, and how to manage them
        ! normally we assume closed zones, ie. if the problem is 2d and we have a lined region
        ! in the set, then it cannot be represented. similarly a point in 1d or a surface in 3d cannot
        ! be represented. thats why we either call the routine for the same domain
        ! or when we are sure about closed zones, then we can call it. 
        ! however, this property will be implemented after the success of "set idea"


        !myNodes=myModel%physics(myPhysics)%domain(d)%dNodes
        myNodes=myModel%physics(myPhysics)%domain(d)%totalNode
        
        !myLocalDof=myModel%physics(myPhysics)%domain(d)%localDof
        myLocalDof=myModel%physics(myPhysics)%dof

        !myDof=myModel%physics(myPhysics)%domain(d)%dDof
        myDof=myModel%physics(myPhysics)%domain(d)%totalDof

        !myNz=myModel%physics(myPhysics)%domain(d)%Nz
        
        myNy=myModel%physics(myPhysics)%domain(d)%Ny

        myNx=myModel%physics(myPhysics)%domain(d)%Nx

        !myNzCells=myModel%physics(myPhysics)%domain(d)%NzCells
        
        !myNyCells=myModel%physics(myPhysics)%domain(d)%NyCells
        myNyCells=myModel%physics(myPhysics)%domain(d)%Ny-1

        !myNxCells=myModel%physics(myPhysics)%domain(d)%NxCells
        myNxCells=myModel%physics(myPhysics)%domain(d)%Nx-1


        
        
        IF (allocated(resultVec)) deallocate (resultVec)


        !print*,mydof
        !# correct this with interior zone dof definition for finite volume
        !if (dFX.eq.dFV) myDof=myNzCells*myNxCells*myNyCells*myLocalDof
        if (dFX.eq.dFV) myDof=1*myNxCells*myNyCells*myLocalDof
        !print*,mydof
        !stop


        allocate(resultVec(myDof,myLocalDof))

        resultVec=0.d0



        gridInfo=getElementNodes(myNx,myNy,1)

        numberOfElementalNodes=gridInfo(1)

        axis=gridInfo(2)


        NumberOfElements=0

        sum=1

        if (dFX.ne.dFV) then



                        do j=1,myNy

                                do i=1,myNx

                                        !if (myModel%physics(myPhysics)%domain(d)%indexMatrix(i,j,k,2,myDof+1).eq.0 ) cycle

                                        do t=1,myLocalDof

                                                !resultVec(myModel%physics(myPhysics)%domain(d)%dofMatrix(i,j,k,2,myLocalDof),t)=myModel%physics(myPhysics)%x(myModel%physics(myPhysics)%domain(d)%dofMatrix(i,j,k,1,t))
                                                !resultVec(sum,t)=myModel%physics(myPhysics)%x(myModel%physics(myPhysics)%domain(d)%dofMatrix(i,j,t,1))
                                                resultVec(sum,t)=y(myModel%physics(myPhysics)%domain(d)%dofMatrix(i,j,t,1))

                                                

                                        end do
                                        
                                        sum=sum+1

                                        ! here determine the number of element
                                        ! with another loop complete the connectivity list
                                        ! or with one loop we can prefer a reshape on the connectivity
                                        ! but I believe two loops are faster than one loop and continous pause and go s

                                        
                                end do
                
                        end do
                                

        else


                !do k=1,myNzCells

                        do j=1,myNyCells

                                do i=1,myNxCells

                                        !if (myModel%physics(myPhysics)%domain(d)%indexMatrix(i,j,k,2,myDof+1).eq.0 ) cycle


                                        !print*,i,j

                                        do t=1,myLocalDof

                                                !resultVec(myModel%physics(myPhysics)%domain(d)%dofMatrix(i,j,k,2,myLocalDof),t)=myModel%physics(myPhysics)%x(myModel%physics(myPhysics)%domain(d)%dofMatrix(i,j,k,1,t))
                                                resultVec(sum,t)=myModel%physics(myPhysics)%x(myModel%physics(myPhysics)%domain(d)%dofMatrix(i,j,t,1))

                                                

                                        end do
                                        
                                        sum=sum+1


                                        ! here determine the number of element
                                        ! with another loop complete the connectivity list
                                        ! or with one loop we can prefer a reshape on the connectivity
                                        ! but I believe two loops are faster than one loop and continous pause and go s

                                        
                                end do
                
                        end do
                                
                !end do 



        
        end if

        !print*,resultVec(:,1)

        !print*,sum-1
        
        !stop
        
        
        





        !allocate(ConnectivityList(NumberOfElements,numberOfElementalNodes))


        !ConnectivityList=0



        select case(axis)




        case(1)


                
                k=1
                j=1

                do i=1,myNx-1

                        if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i,j,2).eq.0 ) cycle

                        if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i+1,j,2).eq.0 ) cycle

                        NumberOfElements=NumberOfElements+1


                end do



        case(2)


                k=1

                do j=1,myNy-1

                        do i=1,myNx-1

                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i,j,2).eq.0 ) cycle

                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i,j+1,2).eq.0 ) cycle

                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i+1,j,2).eq.0 ) cycle
                                
                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i+1,j+1,2).eq.0 ) cycle
                                
                                NumberOfElements=NumberOfElements+1


                        end do

                end do


        !NumberOfElements=(myNx-1)*(myNy-1)


        !print*,NumberOfElements

        IF (allocated(ConnectivityList)) deallocate (ConnectivityList)
        allocate(ConnectivityList(NumberOfElements,numberOfElementalNodes))


        ConnectivityList=0

        sum=1


                do j=1,myNy-1

                        do i=1,myNx-1

                                !print*,i,j

                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i,j,2).eq.0 ) cycle   ! no node cycle
                                
                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i,j+1,2).eq.0 ) cycle

                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i+1,j,2).eq.0 ) cycle

                                if (myModel%physics(myPhysics)%domain(d)%nodeMatrix(i+1,j+1,2).eq.0 ) cycle



                                ConnectivityList(sum,1)=myModel%physics(myPhysics)%domain(d)%nodeMatrix(i,j,2)
                                ConnectivityList(sum,2)=myModel%physics(myPhysics)%domain(d)%nodeMatrix(i+1,j,2)
                                ConnectivityList(sum,3)=myModel%physics(myPhysics)%domain(d)%nodeMatrix(i+1,j+1,2)
                                ConnectivityList(sum,4)=myModel%physics(myPhysics)%domain(d)%nodeMatrix(i,j+1,2)


                                sum=sum+1


                        end do

                end do





        


        end select





        



        call fileNamer(myPhysics,d,file,timeIter)

        
        Open (4400,File=file,Status='unknown')  !<--- physics and domain and subdomain based numbering
        
        
        
        910 FORMAT((E18.10))
        911 FORMAT(4(I6,1X))

        WRITE(4400,*) 'TITLE="results"'


        !WRITE(4400,*) 'VARIABLES="x","y","z","SF","VOR"'
        !WRITE(4400,*) 'VARIABLES="x","y","z","SF","VOR","Temp"'
        !WRITE(4400,*) 'VARIABLES="x","y","z","p","u","v"'
        !WRITE(4400,*) 'VARIABLES="x","y","z","p","u","v","T"'
        !WRITE(4400,*) 'VARIABLES="x","y","z","p","u","v","F"'
        write(4400,*) TecplotVariables  !<--- define for the physics and point to it if valid for the domain else define other


        if (dFX.NE.dFV) then
        
                WRITE(4400,*) 'ZONE n=', myNodes, ', e=', NumberOfElements , ', DATAPACKING=BLOCK, ZONETYPE=FEQUADRILATERAL'
        
        else
        
                !WRITE(4400,*) 'ZONE n=', myNodes, ', e=', NumberOfElements , ', DATAPACKING=BLOCK, ZONETYPE=FEQUADRILATERAL,VARLOCATION=([4,5,6]=CELLCENTERED)'
                WRITE(4400,*) 'ZONE n=', myNodes, ', e=', NumberOfElements , ',DATAPACKING=BLOCK,ZONETYPE=FEQUADRILATERAL,VARLOCATION=([3,4,5,6]=CELLCENTERED)'

        end if
        

        if (realTime.EQV.YES) WRITE(4400,*) 'STRANDID=',myPhysics,' SOLUTIONTIME=',  dt*timeIter




        ! write the coordinates of the mesh nodes
        WRITE(4400,910) myModel%physics(myPhysics)%domain(d)%xCoor
        WRITE(4400,910) myModel%physics(myPhysics)%domain(d)%yCoor
        !WRITE(4400,910) myModel%physics(myPhysics)%domain(d)%zCoor


        do t=1,myLocalDof
                do i=1,myDof/myLocalDof
                        WRITE(4400,910) ResultVec(i,t)
                end do
        end do

        do i=1,NumberOfElements
                do j=1,numberOfElementalNodes
                WRITE(4400,911) ConnectivityList(i,j)
                end do
        end do

        close(4400)


        !dum1=myModel%physics(myPhysics)%x(myModel%physics(myPhysics)%domain(d)%dofmatrix(1,2,1,1,2))
        !dum2=myModel%physics(myPhysics)%x(myModel%physics(myPhysics)%domain(d)%dofmatrix(1,3,1,1,2))

        !print*,dum1,(dum1+dum2)/2.d0
        

        !print*,size(myModel%physics(myPhysics)%domain(d)%dofmatrix(0,:,1,1,1))
        !pause





end do




end subroutine tecFD

---