E:/turkuaz/krylov/codes/newton/1.2b/physics/bratu/anaJ/postProcessor.f90

Go to the documentation of this file.

module postProcessor
use localParameters
implicit none


contains




!#####################################################


subroutine visualizer(X,n) 

implicit none

integer n

real(8) X(n)

real(8) dx,dy


integer i,j,k
integer index
real(8) CoorX(nx,ny),CoorY(nx,ny),ResultVectors(nx*ny,unk)  ! size of the each result may vary if staggered. rather use type!
real(8) ResultMatrices(nx,ny,unk)






dx=Lx/(nx-1.d0)
dy=Ly/(ny-1.d0)



do i=1,nx*ny
    index=(i-1)*(unk-1)+i
        do k=1,unk
                ResultVectors(i,k)=x(index+k-1)
    end do
end do


dx=Lx/(nx-1.d0)
dy=Ly/(ny-1.d0)


!fill xcoor matrix
!$OMP PARALLEL DO
do i=1,nx
        CoorX(i,:)=(i-1)*dx
end do
!$OMP PARALLEL DO

!fill ycoor matrix
!$OMP PARALLEL DO
do j=1,ny
        CoorY(:,j)=(j-1)*dy
end do
!$OMP END PARALLEL DO


do j=1,ny
        do i=1,nx
                index=(j-1)*nx+i
                do k=1,unk
                        ResultMatrices(i,j,k)=ResultVectors(index,k)
                end do
                

        end do
end do

Open (4400,File='results.dat',Status='unknown')
910 FORMAT(5(E18.10,1X))

WRITE(4400,*) 'TITLE="results"'
!WRITE(4400,*) 'VARIABLES="x","y","u","v","p"'
if (unk.eq.1) WRITE(4400,*) 'VARIABLES="x","y","T"'
if (unk.eq.2) WRITE(4400,*) 'VARIABLES="x","y","str","vor"'
if (unk.eq.3) WRITE(4400,*) 'VARIABLES="x","y","u","v","p"'
if (unk.eq.4) WRITE(4400,*) 'VARIABLES="x","y","u","v","p","f"'
WRITE(4400,*) 'ZONE I=', nx, ',J=', ny,', ZONETYPE = Ordered, DATAPACKING = BLOCK'


DO j= 1, ny
        DO i= 1, nx
                WRITE(4400,910) CoorX(i,j)
        END DO
END DO


DO j= 1, ny
        DO i= 1, nx
                WRITE(4400,910) CoorY(i,j)
        END DO
END DO


do k=1,unk

        DO j= 1, ny
                DO i= 1, nx
                        WRITE(4400,910) ResultMatrices(i,j,k)
                END DO
        END DO

end do





close(4400)


end subroutine visualizer



end module postProcessor

---