VTKPlot example 2

This example shows how to use plot function from [[VTKPlot_]] class. We plot a cartesian grid and function over it.

Use association

• [[VTKPlot_]]

Usage

MODULE func_module
  USE easifemBase
  CONTAINS
  PURE FUNCTION f( x ) RESULT( ans )
    REAL(DFP), INTENT( IN ) :: x( : )
    REAL(DFP) :: ans
    ans = x(1)*(x(1)-1) + x(2)*(x(2)-1)
  END FUNCTION
END MODULE func_module

PROGRAM MAIN
  USE easifemBase
  USE easifemClasses
  USE func_module
  IMPLICIT NONE
  TYPE( VTKPlot_ ) :: obj
  REAL( DFP ), ALLOCATABLE :: x( : ), y( : )
  PROCEDURE(f), POINTER :: func => NULL()

!!! note "" Plotting a two dimensional structured grid.

  x = linspace( 0.0_DFP, 1.0_DFP, 5)
  y = linspace( 0.0_DFP, 1.0_DFP, 5)
  CALL obj%initiate( )
  func => f
  CALL obj%plot( x, y, func, "./test_2.vts")
  CALL obj%deallocate()
  func => NULL()

end program