/*gbar0 compute g-bar gaunt factor for neutrals */
void gbar0(double ex, 
	  float *g);