[FLASH-USERS] Fix for Checkerboard of Magnetic Monopoles

[Aaron Froese]

There is a bug in the distributed version of the Radloss cooling module in which the block is identified incorrectly when solving the equation of state.  The isothermal and SutherlandandDopita cooling modules do not suffer from this problem.  If implementing source terms based on Radloss, a lattice of magnetic monopoles will appear when using the unsplit MHD solver.  The monopoles surround all of the smallest blocks that have a blockID larger than the number of smallest blocks.

File: source/physics/sourceTerms/Cool/CoolMain/Radloss/Cool.F90

Existing Line 147:
call Eos_wrapped(MODE_DENS_EI,blkLimits,blockID)

Corrected Line 147:
call Eos_wrapped(MODE_DENS_EI,blkLimits,blockList(blockID))