Ok, coma disappearance isn't implemented.
Something odd seems to be going on here with PSR B1257+12...
Celestia reports the following parameters for the pulsar:
r=10 km
T=5,000,000 K
L=0.00000000136 times solar
Now as I understand it the displayed value corresponds to visual luminosity, but we probably want the bolometric luminosity.
By calculation, we have:
L=4?€r????T^4
which works out as 116 times solar -- a discrepancy of 11 orders of magnitude!!!
Now we take the description in the thread:
t00fri @ thread 8453 wrote:For our solar system, the threshold irradiance where the tail starts fading, is (empirically) known to be
Irr_crit =L_sol/x0_sol^2, with x0_sol ~ 5 AU and L_sol being the known solar luminosity.
It is now very sensible to assume that comet tails start vanishing at a universal threshold irradiance
Irr_crit that we know already from our solar system.
Hence for any other system (with 1 sun), we immediately obtain the scaling law for the fading start x0
+++++++++++++++
x0 = x0_sol sqrt(L/L_sol).
+++++++++++++++
This gives x0_pulsar=x0_sol*sqrt(116)=54 AU!
However the tail is not being drawn!
Does the comet tails code use the bolometric luminosity (as used in the temperature calculation) or the displayed (visual) luminosity?
If it is using visual luminosity, we get x0_pulsar=x0_sol*sqrt(0.00000000136)=0.00018 AU
Lo and behold, when the comet goes inside this radius, the tail gets drawn!
So it looks like the visual luminosity is being used for the comet tail calculation! Shouldn't the code really be using the bolometric luminosity though?