Create_field Example
! 3D force is contained as fXYZ variable
! in a user-defined type node_type
TYPE(node_type), ALLOCATABLE, DIMENTION(:) :: nodes
INTEGER :: fid
...allocate nodes array as n_nodes...
fid = FEM_Create_Field(FEM_DOUBLE,3,
          foffsetof(nodes(1), nodes(1)%fXYZ),
          foffsetof(nodes(1), nodes(2)) )