PRO read_sphercube, filename, r, theta, phi, rho

; Read the data into a structure

  s = H5_PARSE(filename, /READ_DATA)

; Extract the relevant quantities

  r = s.x_min._data + (s.x_max._data - s.x_min._data)*(FINDGEN(s.n_x._data) + 0.5)/s.n_x._data

  theta = s.y_min._data + (s.y_max._data - s.y_min._data)*(FINDGEN(s.n_y._data) + 0.5)/s.n_y._data

  phi = s.z_min._data + (s.z_max._data - s.z_min._data)*(FINDGEN(s.n_z._data) + 0.5)/s.n_z._data

  rho = s.data._data

; Finish

  return

END