SetGaussPoint.f90 Source File


Contents

Source Code


Source Code

if(allocated(obj%gzi) ) then
    deallocate(obj%gzi)
endif
allocate(obj%gzi(obj%NumOfDim) )

if(obj%NumOfDim /= size(obj%GaussPoint,1) )then
    print *, "ERROR::SetGaussPoint",obj%NumOfDim, size(obj%GaussPoint,1)
    obj%ErrorMsg="ERROR::SetGaussPoint"
    obj%ierr=1
else
    obj%gzi(:) = obj%GaussPoint(:,obj%GpID )
    obj%ErrorMsg="Succeed::SetGaussPoint"
    obj%ierr=0
endif