<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hey everyone,<br><div><br>  I'm pretty sure there's a bug in the unsplit_Rad Hydro_computeDT: <br>source/physics/Hydro/HydroMain/unsplit_rad/Hydro_computeDt.F90<br><br>Line 241 reads:<br>241:                 if (NDIM > 2) dt_ltemp = max(dt_ltemp,(abs(U(VELZ_VAR,i,j,k)-uzgrid(j))+sqrt(cfz2))*delzinv)<br><br>But should be:<br>241:                 if (NDIM > 2) dt_ltemp = max(dt_ltemp,(abs(U(VELZ_VAR,i,j,k)-uzgrid(k))+sqrt(cfz2))*delzinv)<br><br><br>However, I was getting the error:<br>***********************************************************<br>  Warning: The initial timestep is too large.<br>    initial timestep =   1.000000000000000E-010<br>    CFL timestep     =   5.215016855876254E-303<br>  Resetting dtinit to dr_tstepSlowStartFactor*dtcfl.<br> ***********************************************************<br>both with and without the fix above. With the fix, I dug deeper into why the dt was going to zero, and found that it related to guard cells. I then also initialized my guard cells in Simulation_initBlock and the small dt problem went away.<br><br>Oddly, a collaborator took my SimulationMain/Code/ and ran it, and got the small dt error again, even with the guardCells initialized. I forgot to give him the above fix to Hydro_computeDT. He made that fix, and now his small dt problem is gone as well. <br><br></div><div>I'm using the hypre unsplit solver, with radiation and mgd on. Any thoughts on why an intiial FillGuardCells call isn't fixing the problem?<br><br></div><div>Cheers,</div><div>  -Mark</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>-- <br></div><div><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span style="font-size:12.8px">Mark Richardson</span><br style="font-size:12.8px"><div style="font-size:12.8px">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><a href="tel:%2B44%201865%20283020" value="+441865283020" target="_blank">+44 1865 283020</a></div></div></div></div></div></div></div></div></div></div></div></div>
</div></font></span></div>
</blockquote></div><br>
</div></div>