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

Use association​

Usage​

Use association

Usage