Skip to main content

VTKFile example 7

This example shows how to construct structuredGrid (ASCII) vtkfile. We make 2D planar grid. We also plot a function on this structuredGrid.

Use association​

  • [[VTKFile_]]

Usage​

PROGRAM MAIN
  USE easifemBase
  USE easifemClasses
  Implicit None
  TYPE(VTKFile_) :: obj
  REAL(DFP), ALLOCATABLE :: xv(:), yv(:), zv(:), x(:, :, :), &
    & y(:, :, :), z(:, :, :), func(:,:,:)
  INTEGER(I4B) :: nx1 = 0, nx2 = 9, ny1 = 0, ny2 = 5, nz1 = 0, nz2 = 0

!!! note "" Let us initiate a vtkfile for structured grid. We need to specify the WholeExtent.

  CALL obj%InitiateVTKFile( &
    & filename="./example7.vts", &
    & mode="NEW", &
    & DataFormat=VTK_ASCII, &
    & DataStructureType=VTK_StructuredGrid, &
    & WholeExtent=[nx1, nx2, ny1, ny2, nz1, nz2])

!!! note "Piece 1" In the following we add a piece (portion of geometry defined by the extent). Note that we can specify as many piece as we want.

  CALL obj%WritePiece(extent=[nx1, nx2, ny1, ny2, nz1, nz2])
  !!
  xv = linspace(0.0_DFP, 1.0_DFP, (nx2 - nx1) + 1)
  yv = linspace(0.0_DFP, 1.0_DFP, (ny2 - ny1) + 1)
  zv = linspace(1.0_DFP, 1.0_DFP, (nz2 - nz1) + 1)
  CALL MeshGrid(x=x, y=y, z=z, xgv=xv, ygv=yv, zgv=zv)
  !!
  CALL obj%WritePoints(x=x, y=y, z=z)
  !!

!!! note "PointData"

  CALL obj%WriteDataArray( &
    & location=String("node"), &
    & action=String("open") )
  func = x**2 + y**2
  CALL obj%WriteDataArray( &
    & name=String("func"), &
    & x=func, &
    & numberOfComponents=1 )
  CALL obj%WriteDataArray( &
    & location=String("node"), &
    & action=String("close") )

!!! note "Close Piece"

    CALL obj%WritePiece()

!!! note "cleanup" Closing the file, a simple Deallocate is enough.

  CALL obj%Deallocate()
END PROGRAM main