subroutine FEM_Set_Node_Data_r(data)
REAL*8, intent(in),
        dimension(doublePerNode,nNode) :: data
subroutine FEM_Get_Node_Data_r(data)
REAL*8, intent(out),
        dimension(doublePerNode,nNode) :: data