<div dir="ltr">Hi Eddie :<div><br></div><div>   It's possible to interpolate the particles into 3D and then angle-average over azimuthal angle to get your desired r-z 2.5D mesh. However, in a paper I wrote with my students and colleagues, we developed an elegant way to interpolate directly into 2D axisymmetry. The method is described in an appendix. I am going to quote myself here :</div><div><br></div>"However, a more elegant and efficient method, which we adopt, is to break up the angle-dependent quantities, which involve only the smoothing kernel, from the particle data. This reduces the angle-averaging procedure to simple quadrature, with the weights written as integrals to be pre-tabulated once. The averaged quantities are then rapidly and efficiently calculated by weighting the particle data at each point on the r-z mesh over lookups of the pre-tabulated weights."<div><br></div><div>The paper is <br><div><br></div><div><a href="http://iopscience.iop.org/article/10.1088/0004-637X/773/2/136/meta#apj478831app1">http://iopscience.iop.org/article/10.1088/0004-637X/773/2/136/meta#apj478831app1</a><br></div><div><br></div><div>You can also write to me offline and we can share the relevant source code.</div><div><br></div><div>  Best wishes,</div><div>   </div><div>  Bob</div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 7, 2016 at 4:44 PM, Hix, William Raphael <span dir="ltr"><<a href="mailto:raph@ornl.gov" target="_blank">raph@ornl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Eddie-<br>
<br>
There is a tool SPLASH, intended for visualizing SPH data that may help, since mapping to pixels is similar to mapping to grids.<br>
<br>
<a href="http://users.monash.edu.au/~dprice/splash/" rel="noreferrer" target="_blank">http://users.monash.edu.au/~<wbr>dprice/splash/</a><br>
<br>
                                Raph<br>
<div class="HOEnZb"><div class="h5"><br>
> On Dec 7, 2016, at 4:40 PM, Baron, Edward A. <<a href="mailto:baron@ou.edu">baron@ou.edu</a>> wrote:<br>
><br>
><br>
><br>
> I am trying to take the output of a model calculated with SPH and use it as an initial model for a  FLASH run.<br>
><br>
> Can anyone point me to documentation on how to define initial physical values at each grid point?<br>
><br>
> This must be common (the results of 1-D hydrostatic equilibrium runs are often input), etc. But somehow I’m not finding the docs on how to set the initial values in FLASH.<br>
><br>
> I need to use the MHD solver for so I can run 2.5D (including rotation), these are pure hydro, just need the 2.5D capability (and also treatment of viscosity).<br>
><br>
> Any help with this would be greatly appreciated.<br>
><br>
> -eddie<br>
><br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Dr. Robert Fisher<br>Associate Professor / Graduate Program Director<br>University of Massachusetts/Dartmouth<br>Department of Physics<br>285 Old Westport Road<br>North Dartmouth, Massachusetts 02747<br><a href="mailto:robert.fisher@umassd.edu" target="_blank">robert.fisher@umassd.edu</a><br><a href="http://www.novastella.org" target="_blank">http://www.novastella.org</a></div></div>
</div>