[FLASH-BUGS] Bug in quadratic_cartesian interpolation scheme?

[Robi Banerjee]

Hi,

  I ran into a problem (non convergence in the Riemann solver)
using the new 'quadratic_cartesian' interpolation scheme in a
3D simulation including selfgravity. 
  Looking into the details after a refinement step I figured 
that the data in the newly refined blocks develop stripe-like
structures. This leads, for instance, to overshooting in 
velocities, resp. mach numbers. 
  Using the old interpolation scheme (second_order_old) 
everything is fine and the variables are smoothly prolongated 
onto the newly refined blocks.

  Is this a known problem, or is there a problem with the
new (quadratic_cartesian) interpolation scheme?

  If you need more information, please don't hesitate to
ask for (e.g. 2D slices of the simulation data with the
two different interpolations schemes)

  Thank you,
Robi

====================================
Robi Banerjee
Department of Physics and Astronomy
ABB-320, McMaster University
Hamilton, ON L8S 4M1
CANADA
e-mail: banerjee at physics.mcmaster.ca
phone : (905) 525-9140 x 23189
fax   : (905) 546-1252