found bug in NR absorption term solver-- updates to intensity(1:na-1) didnt scale the B coefficient by the new temp
Small error 10^-2 in the fluid momentum updates. It can be eliminated by dropping the InvCrat term in the dGasMomentum, and reverting the velocity estimate at t+1/2 to the velocity at t, but this is too hand wavey. How do I use the velocity estimate?