[FLASH-USERS] Fix for Checkerboard of Magnetic Monopoles

[dongwook at flash.uchicago.edu]

Dear Aaron,

I am very glad that you identified and fixed the issue. We will make your
change in Cool.F90 right away. Thank you for your bug fixing and reporting
this to us.

Best,
Dongwook
=========================================
Dongwook Lee, Ph.D., Research Scientist
The Flash Center for Computational Science
The University of Chicago
5747 S. Ellis Ave., Room 319
Chicago, IL 60637
(773) 834-6830

> 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))