Aw man... I hate to be giving you work, but I think I might've found another bug, and this one will probably not be a trivial fix...
Using system seed
020201 as an example, the innermost planet has Orbit#0.05, and is listed as being in the HZ. This shouldn't, however, be the case.
Page 43 of the WBH defines HZCO Deviation as
Orbit# - HZCO, naturally enough. But it then notes the following:
If either the HZCO or the Orbit# of the world lies below Orbit#1, the calculation becomes more complicated. The same basic subtraction occurs but the result is modified by dividing by the smaller of the HZCO or the worlds’ Orbit number and will result in a greater effective deviation.
And provides the following formula:
Effective HZCO Deviation = (Orbit# - HZCO)/Smaller of Orbit# or HZCO.
In the system generated by the above seed, the HZCO is 0.13, so using the above formula we find that the
Effective HZCO should be:
(0.05 - 0.13)/0.05 = -1.6, which is smaller than -1 and therefore inside the inner edge of the HZ.
Using another system as example, that generated via continuation method by seed
000, 0 Gas Giants, 0 Belts, 12 Terrestrials; the HZCO for this one is 0.04, the outermost planet has Orbit#0.14. Using the formula, we find:
(0.14 - 0.04)/0.04 = 2.5, most definitely outside the habitable zone, yet listed as HZ.
I think the 'fix' for this one would be adding an
IF logic block to the HZCO Deviation calculation that is triggered whenever either the object's orbit or HZCO is
< 1, and then using the smaller of the two values as the dividend for the alternate formula.