<div dir="ltr"><div><div><div>Dear Slah,<br><br></div>  there is a default parameter that sets the minimum density (and e and p, etc) that can happen. The default value, I believe, is 1e-10. It's called sml_rho (there's also small_e, small_p I believe). You could set this to a lower value in your flash.par file. <br><br></div>Cheers,<br></div>  -Mark<br><div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On 2 July 2018 at 13:31, Slah Abdellaoui <span dir="ltr"><<a href="mailto:slah.abdella@gmail.com" target="_blank">slah.abdella@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear All</div><div>I am running a problem, where the initial density is proportional to 1/r^2, whenever I change the parameter the density profile reaches the of 1e-10 (figure dens.png) and it doesn't go further, as it suppose to vary as a function of 1/ r^2. I don't know what causing that. I am using 1d spherical coordinate problem with PPM hydro solver and including the gravity point mass.</div><div>I have set "convertToConsvdInMeshInterp =.True." and I saw the output in attached log file.</div><div>Some of the configuration I get non convergence in Riemann subroutine.</div><div><br></div><div>I look forward to receiving your response.</div><div>Best regards,</div><div>Slah​<div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:rgb(245,245,245);padding:5px;font-family:arial;font-weight:bold;border:1px solid rgb(221,221,221);line-height:1"><a href="https://drive.google.com/file/d/1mMbxAvAXxlFNEzMD8CwmuXPe37FmHah-/view?usp=drive_web" style="display:inline-block;max-width:366px;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0;border:none" target="_blank"><img style="vertical-align:bottom;border:none" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png"> <span dir="ltr" style="color:rgb(17,85,204);text-decoration:none;vertical-align:bottom">log</span></a><img src="" style="opacity:0.55;float:right;display:none"></div>​<br></div></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br><span style="font-size:12.8px">Mark Richardson</span><br style="font-size:12.8px"><div style="font-size:small"><span style="font-size:12.8px">MAT Postdoctoral Fellow</span><div style="font-size:12.8px">Department of Astrophysics</div><div><span style="font-size:12.8px">American Museum of Natural History</span><div style="font-size:12.8px"><a href="mailto:Mark.Richardson.Work@gmail.com" style="color:rgb(17,85,204)" target="_blank">Mark.Richardson.Work@gmail.com</a></div><div><a href="https://sites.google.com/site/marklarichardson/" style="color:rgb(17,85,204)" target="_blank">My Website</a><br><div style="font-size:12.8px">212 496 3432</div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div></div></div></div>