, dim)          integer, intent(in)               :: ele, dim         type(scalar_field), intent(inout) :: psi          real, dimension(ele_loc(psi,ele),ele_ngi(psi,ele),dim) :: dshape_psi         real, dimension(ele_ngi(psi,ele))                      :: detwei         real, dimension(ele_loc(psi,ele))                      :: rhs_addto         type(element_type), pointer                            :: shape_psi         integer, pointer, dimension(:)                         :: nodes_psi              nodes_psi => ele_nodes(psi, ele)         shape_psi => ele_shape(psi, ele)         call transform_to_physical( positions, ele, shape_psi, dshape=dshape_psi, detwei=detwei )          