Looks like it was another scale issue. My density function:
1013.25 * pow(e,-.029*9.8*(pow(x*x+y*y+z*z,0.5)-earth_radius)/8.134/298)
with earth_radius a little less than 1000, may have pushed some
floating-point limits. Shrinking the entire scene helped, for now.