<div dir="ltr">Dear FLASH users,<div><br></div><div>TL;DR: Setting xl_boundary_type = "user", cells are filled with zero density although Grid_bcApplyToRegionSpecialized was not yet called. There is no such problem with all periodic BCs. If someone familiar with user BCs could mention which additional routines are executed that I might want to look into to debug this problem I'd very much appreciate it.</div><div><br></div><div>=========================================================</div><div><br></div><div>I am trying to modify the Torus problem included in FLASH4.5 to include a magnetized inflow from the lower x-boundary. </div><div><br></div><div>So, I modified flash.par to turn off gravity, comment out gr_prmp* variables, set geometry to cartesian, lrefine_min = 3, and set all BCs to periodic except xl_boundary_type = "user".</div><div><br></div><div>Additionally, I modified Simulation_initBlock so the domain is initially at a uniform density and temperature with zero velocity and zero magnetic field.</div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><br></div><div>When I run the simulation with all BCs set to periodic, I receive many warning messages about the density being the value I set (modulo roundoff error) with a lot of output (possibly printing the entire density array?); however, it runs to completion fine and plots show the domain remained uniform (as to be expected).</div><div><br></div><div>However, when I set xl_boundary_type = "user" then the warning messages say minimum density is -1.0 and it crashes with:</div><div>





<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures">ERROR After calling Eos, eosData(EOS_EINT) or eosData(EOS_DENS) are zero</span></p></div><div dir="ltr"><br></div><div>Compiling with -debug I receive a "floating invalid" error with traceback to eos_idealGamma.F90 with a line in which there is a division by density (which is zero, triggering the floating invalid error). However, I have write statements at the start of Grid_bcApplyToRegionSpecialized that are never executed suggesting that routine is not called yet.</div><div><br></div><div>I am seeking to eliminate this zero density issue; I'd be very grateful for any pointers those familiar with user BCs could point out.</div><div><br></div><div>Best,</div><div>--------</div><div dir="ltr"><div>Ryan Farber<br></div><div>PhD Candidate, University of Michigan</div></div></div></div></div></div></div>