          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_vel, rhs_addto_buoy         real, dimension(ele_ngi(psi,ele))                      :: 