Get1 Method in STVectorField

[vickysharma0812]

Feature 1

We want to get x, y, z coordinate of a STVector field at a given time step by calling following method.

CALL Obj%Get(value(:), globalNode, timeCompo) 

Feature 2

Get time component of a given space component at a given node.

CALL Obj%Get(value(:), globalNode, spaceCompo) 

Feature 3

Get space-time component a t given node.

CALL obj%Get(value(:), globalNode, spaceCompo, timeCompo)

[vickysharma0812]

PROGRAM main
USE easifemBase
USE easifemClasses
TYPE(Domain_) :: dom
TYPE(STVectorField_) :: obj
TYPE(HDF5File_) :: meshfile
TYPE(ParameterList_) :: param
INTEGER(I4B) :: ierr
INTEGER(I4B), PARAMETER :: spaceCompo = 2_I4B
INTEGER(I4B), PARAMETER :: timeCompo = 2_I4B
REAL(DFP), ALLOCATABLE :: svalue(:, :, :)
REAL(DFP), ALLOCATABLE :: VALUE(:)
CHARACTER(*), PARAMETER :: engine = "NATIVE_SERIAL"
INTEGER(I4B) :: tnodes
LOGICAL(LGT) :: isok

! import domain
CALL meshfile%initiate(filename="./mesh.h5", mode="READ")
CALL meshfile%OPEN()
CALL dom%initiate(hdf5=meshfile, group="")
CALL meshfile%DEALLOCATE()
! set parameters for VectorField_
CALL FPL_INIT()
CALL param%initiate()
CALL SetSTVectorFieldParam( &
  & param=param, &
  & fieldType=FIELD_TYPE_NORMAL, &
  & name="U", &
  & engine=engine, &
  & spaceCompo=spaceCompo, &
  & timeCompo=timeCompo)
CALL obj%initiate(param, dom)
tnodes = dom%GetTotalNodes()
CALL Reallocate(svalue, spaceCompo, timeCompo, tnodes)
CALL RANDOM_NUMBER(svalue)
CALL obj%set(VALUE=svalue)

CALL obj%Get(VALUE=VALUE, spaceCompo=1, globalNode=3)
isok = ALL(svalue(1, :, 3) .approxeq.VALUE)
CALL ok(isok, "Get(1) ")

CALL obj%Get(VALUE=VALUE, timeCompo=2, globalNode=3)
isok = ALL(svalue(:, 2, 3) .approxeq.VALUE)
CALL ok(isok, "Get(2) ")

CALL obj%Get(VALUE=VALUE, timeCompo=2, spaceCompo=1, globalNode=3)
isok = ALL(svalue(1, 2, 3) .approxeq.VALUE)
CALL ok(isok, "Get(3) ")

CALL obj%DEALLOCATE(); CALL dom%DEALLOCATE()
CALL param%DEALLOCATE(); CALL FPL_FINALIZE()
END PROGRAM main