[FLASH-USERS] Bug in 3Thy_uhd_DataReconstructNormalDir_MH.F90

[Mark Richardson]

Hello,

  there appears to be a bug in hy_uhd_DataReconstructNormalDir_MH.F90, at
least for the 3T case.

  For the MUSCL-Hancock case, korder is set to 1 in
hy_uhd_dataReconstOneStep.F90, but in
hy_uhd_DataReconstructNormalDir_MH.F90 at line 227 Data1D (which points to
U(:,ix-korder:ix+korder,iy,iz) for the x sweep say, thus sizeof(Data1D) =
NVAR * (2korder+1) = NVAR*3 for MH) is indexed at 4+iBeg (iBeg=0), and is
thus indexing out of range.

 Where should this be pointing? Should Vm,Vc,Vp point to 1:3 in Data1D?

Thank you,
  -Mark

-- 

Mark Richardson
Beecroft Postdoctoral Fellow, Department of Physics
Denys Wilkinson Building - 555F
University of Oxford
Mark.Richardson at physics.ox.ac.uk
+44 1865 280793
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20160509/28abd6d1/attachment.htm>