Skip to main content

BlockMatrixField example 6

This example shows the use of BlockMatrixField

Usage​

PROGRAM main
  USE easifemBase
  USE easifemClasses
  IMPLICIT NONE
  TYPE( BlockMatrixField_ ) :: obj
  TYPE( ParameterList_ ) :: param
  TYPE( Domain_ ), TARGET :: pressureDomain
  TYPE( Domain_ ), TARGET :: velocityDomain
  TYPE( DomainPointer_ ) :: domains( 2 )
  TYPE( HDF5File_ ) :: pressureMeshFile
  TYPE( HDF5File_ ) :: velocityMeshFile
  TYPE( HDF5File_ ) :: blockMatrixOutfile
  CHARACTER( LEN = * ), PARAMETER :: pressureMeshFilename="./mesh_tri3.h5"
  CHARACTER( LEN = * ), PARAMETER :: velocityMeshFilename="./mesh_tri6.h5"
  CHARACTER( LEN = * ), PARAMETER :: blockMatrixOutfileName= &
    & "./blockMatrixOutfile.h5"
  CALL pressureMeshfile%Initiate( FileName=pressureMeshFilename, MODE="READ" )
  CALL velocityMeshfile%Initiate( FileName=velocityMeshFilename, MODE="READ" )
  CALL pressureMeshfile%Open()
  CALL velocityMeshfile%Open()
  CALL pressureDomain%Initiate( pressureMeshfile, "")
  CALL velocityDomain%Initiate( velocityMeshfile, "")
  CALL pressureMeshfile%Deallocate()
  CALL velocityMeshfile%Deallocate()
  domains(1)%ptr => velocityDomain
  domains(2)%ptr => pressureDomain
  CALL blockMatrixOutfile%Initiate( blockMatrixOutfileName, "READWRITE" )
  CALL blockMatrixOutfile%open()
  CALL obj%Import(blockMatrixOutfile, "/BlockMatrixField", domains=domains)
  CALL obj%Export(blockMatrixOutfile, "/BlockMatrixFieldState" )
  CALL blockMatrixOutfile%Deallocate()
  CALL obj%Deallocate()
  !> cleanup
  NULLIFY(domains(1)%ptr, domains(2)%ptr)
  CALL pressureDomain%Deallocate()
  CALL velocityDomain%Deallocate()
  CALL obj%Deallocate()
END PROGRAM main