As you may know, I am obsessed with Wolf-Rayet stars. I'm a bit frustrated whenever I play around with WRs and see that they're so much cooler than they should be. For example, HIP 5100 (WR 2) in Celestia has a temperature of 50,000 K, while in reality that temperature is over 140,000 K.

I would like the Wolf-Rayet temperatures to be raised, to reflect their real-life temperatures. Note: This is for the hydrogen-free Wolf-Rayet stars. For the hydrogen-rich ones (WNh stars) I would like a separate class to be created that has temperatures that reflect the WNh stars.

Sometimes, I use the hottest white dwarf (D0) which has a temperature of 100,000 K, to depict hot Wolf-Rayet stars. This works to an extent, but not very well for WRs between 60kK and 100kK, and ones between 110kK and 130kK.

Possible temperatures for WRs:
(Note: The temperatures for early and late WCs are mostly for CSPNe with [WC] spectra)
static float tempWN[10] =
{
185000, 158000, 140000, 115000, 89000, 75000, 63000, 54000, 47000, 39000
};

static float tempWC[12] =
{
225000, 190000, 165000, 140000, 125000, 105000, 83000, 70000, 58000, 46000, 40000, 34000
};

static float tempWNh[12] =
{
133000, 106000, 92000, 72000, 61000, 51000, 42000, 36000, 31000, 27000, 22000, 15000
};

and because there are the 10 and 11 subclasses,

const char* SubclassNames[11] = {
"0", "1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", ""
};

const char* SpectralClassNames[StellarClass::NormalClassCount] = {
"O", "B", "A", "F", "G", "K", "M", "R",
"S", "N", "WC", "WN", "WNh", "WO", "?", "L", "T", "Y", "C",
};