<div dir="ltr">Hello, <div><br></div><div>  there appears to be a bug in hy_uhd_DataReconstructNormalDir_MH.F90, at least for the 3T case. </div><div><br></div><div>  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. </div><div><br></div><div> Where should this be pointing? Should Vm,Vc,Vp point to 1:3 in Data1D? </div><div><br></div><div>Thank you,</div><div>  -Mark<br><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><br>Mark Richardson<br><div>Beecroft Postdoctoral Fellow, Department of Physics<div>Denys Wilkinson Building - 555F</div><div>University of Oxford<div><a href="mailto:Mark.Richardson@physics.ox.ac.uk" target="_blank">Mark.Richardson@physics.ox.ac.uk</a><br><div>+44 1865 280793</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div>