[FLASH-USERS] Error triggering gr_sanitizeData

[Joshua Wall]

Dear Flash users,

     I'm currently testing a cooling function I've written using the
SinkRotatingCloudCore collapse problem from the provided simulations. I've
set the gamma=5/3 and swtiched to the ideal EOS. My cooling function seems
to be working correctly in testing outside of Flash, and cools properly for
a given cell at a certain temperature and density.

    However when I try and let this problem form a sink particle, the cloud
collapses and refines down to a high density, then begins to give the
gr_sanitizeDataAfterInterpolation warnings like:

WARNING after gc filling: min.
unk(EINT_VAR)=-2.4992047803499256E+82          PE=3
block=269                              type=2

This seems to happen close to sink formation, and usually (maybe always,
but I can't remember all of the runs!) after a refinement. I've tried
setting gr_sanitizeDataMode = 3, but this only delays the issue. Also, if I
increase the minimum temperature for the cooling, it delays the issue, but
it always still occurs right up until temperatures that prevent the cloud
from collapsing. I've also tested by switching off sink particles, and it
still occurs, so its not sink formation.

I saw in the flash users list history some mention of whether to apply
sourceTerms to all the blocks or just leaf blocks, but it now appears to me
in Driver_evolveFlash that all blocks are evolved, so I'm not sure if
there's a fix there.

I've also tested this with the split solver and the unsplit Roe and PPM
solvers, and see the same behavior.

Any suggestions would be most welcome and appreciated at this stage, if
anyone has any ideas of how I might tackle this problem.

King regards,

Joshua Wall
Doctorial Candidate
Drexel University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20150630/22c9beb4/attachment.htm>