Skip to main content

VTKFile example 5

This example shows how to construct structuredGrid (ASCII) vtkfile. We show how to create two piece.

Use association​

  • [[VTKFile_]]

Usage​

PROGRAM MAIN
  USE easifemBase
  USE easifemClasses
  Implicit None
  TYPE(VTKFile_) :: aVTKfile
  REAL(DFP), ALLOCATABLE :: xv(:), yv(:), zv(:), x(:, :, :), &
    & y(:, :, :), z(:, :, :)
  INTEGER(I4B) :: nx1, nx2, ny1, ny2, nz1, nz2
  CHARACTER( LEN=* ), PARAMETER :: filename="./example5.vts"

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

  nx1 = 0
  nx2 = 4
  ny1 = 0
  ny2 = 1
  nz1 = 0
  nz2 = 1
  CALL aVTKfile%InitiateVTKFile( &
    & filename=filename, &
    & 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.

  nx1 = 0
  nx2 = 2
  ny1 = 0
  ny2 = 1
  nz1 = 0
  nz2 = 1
  CALL aVTKfile%WritePiece(extent=[nx1, nx2, ny1, ny2, nz1, nz2])
  xv = linspace(0.0_DFP, 10.0_DFP, (nx2 - nx1) + 1)
  yv = linspace(0.0_DFP, 1.0_DFP, (ny2 - ny1) + 1)
  zv = linspace(0.0_DFP, 1.0_DFP, (nz2 - nz1) + 1)

!!! note "" Now we create meshgrid

  CALL MeshGrid(x=x, y=y, z=z, xgv=xv, ygv=yv, zgv=zv)
  CALL aVTKfile%WritePoints(x=x, y=y, z=z)

!!! note "" Closing the piece.

  CALL aVTKfile%WritePiece()

!!! note "Piece 2"

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

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

  CALL aVTKfile%Deallocate()
END PROGRAM main