INVALIDATES: primitive(everywhere)
  INVALIDATES: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)

      INVALIDATES: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)

  INVALIDATES: primitive(everywhere)
  INVALIDATES: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)

  INVALIDATES: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)

  INVALIDATES: primitive(everywhere)
  # Note: We need to keep the fluxes valid for the flux registers 

  # Refluxing is over, we can now mark the fluxes as invalid
  INVALIDATES: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)