<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Thanks for your response Klaus.  That clears up a lot.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Jon</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 14, 2016 at 2:07 PM, Klaus Weide <span dir="ltr"><<a href="mailto:klaus@flash.uchicago.edu" target="_blank">klaus@flash.uchicago.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On Wed, 14 Dec 2016, Slavin, Jonathan wrote:<br>
<br>
> Hi all,<br>
><br>
> I'm wondering what the effective difference is for initiating a simulation<br>
> between 1) using Grid_putPointData and 2) using assignment to the elements<br>
> of solnData (where solnData is returned by Grid_getBlkPtr).  In the default<br>
> supplied Simulation_initBlock for the Sedov simulation, the latter is used<br>
> if FL_NON_PERMANENT_GUARDCELLS is defined and the former is used<br>
> otherwise.  In the Simulation_initBlock for the BlastBS simulation, there<br>
> is no test whether FL_NON_PERMANENT_GUARDCELLS is defined and method 2 is<br>
> used, despite the fact that (as far as I can tell) that symbol is not<br>
> defined.  Any help on this would be appreciated.  Also it'd be helpful if<br>
> someone could tell me under what circumstances FL_NON_PERMANENT_GUARDCELLS<br>
> gets defined.<br>
<br>
</span>Jon,<br>
<br>
The symbol FL_NON_PERMANENT_GUARDCELLS should only be defined if FLASH is<br>
set up to use the no_permanent_guardcells mode of PARAMESH 4.0 or<br>
PARAMESH4dev. You should be able to find more information on this mode in<br>
the PARAMESH documentation.<br>
<br>
We have experimented with this mode years ago, as a way to lower memory<br>
requirements for large simulations. There is some supporting code still<br>
in the FLASH distribution - including support for optionally running some<br>
of the tests in no_permanent_guardcells mode, as you have found - but to<br>
my knowledge this has not been exercised in several years.<br>
<br>
About the question which routine should be called for UNK variable<br>
initialization - both Grid_putPointData and<br>
Grid_getBlkPtr/Grid_<wbr>releaseBlockPtr work and can be used (as well<br>
as the further alternatives Grid_put{Blk,Plane,Row}Data).<br>
I consider this  matter of personal preference and convenience.<br>
<br>
There are probably differences in performance (I would expect the<br>
variant that works with pointers to be the fastest), but this should not<br>
matter much for Simulation_initBlock which is only called at<br>
initialization.<br>
<span class="HOEnZb"><font color="#888888"><br>
Klaus<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">________________________________________________________<br>Jonathan D. Slavin                 Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>       60 Garden Street, MS 83<br>phone: (617) 496-7981       Cambridge, MA 02138-1516<br>cell: (781) 363-0035             USA<br>________________________________________________________<br><br></div></div></div></div>
</div>