We guess what the right edge is in parse_parameter_file, based on the usual setup that Athena right edge = - left edge (maybe not so usual!). This gets corrected during the hierachy build, but the damage to the units is already done in that unitary is 1./(dre-dle), and when dle = 0.0, dre = -0.0, and we get -inf.
Hi Sam, can you open an issue? And, I think we should accept once we have Athena answer testing in place -- which should be pretty easy to do. Ping @jzuhone .
Issue #513 is has been opened. I think, however, that this can be accepted before answer testing...
Alright, I suppose I am inclined to agree. This looks good to me.