In the interpolators the last depth point was not saved, resulting in only 100 depth points beingused, instead of all 101 as in the actual model atmosphere. Not sure if that caused any issues, but in NLTE there were 101 departure coefficient computed. Not sure if that could have an effect (since # of departure coefficients =/= depth points). By default TS only can use up to 100 depth points, so had to change spectrum.inc as well to 101.
In the interpolators the last depth point was not saved, resulting in only 100 depth points beingused, instead of all 101 as in the actual model atmosphere. Not sure if that caused any issues, but in NLTE there were 101 departure coefficient computed. Not sure if that could have an effect (since # of departure coefficients =/= depth points). By default TS only can use up to 100 depth points, so had to change spectrum.inc as well to 101.
Cheers to SP for identifying the bug.