     integer                          :: i, stat, ele     character(len=FIELD_NAME_LEN)    :: bc_type     type(scalar_field), pointer      :: scalar_surface     type(vector_field), pointer      :: positions     type(scalar_field), pointer      :: lumped_mass     type(scalar_field)               :: inverse_lumped_mass