/* diffheav.h */
EXTERN struct t_DiffHeav {
	float DiffusHeavy[LIMELM - 1][LIMELM], 
	  xLyaHeavy[LIMELM - 1][LIMELM];
	}	DiffHeav;

/*
 *DiffusHeavy is total rec rate to grnd for nelem,ion element
 *xLyaHeavy is Lya for nelem,ion
 *do not do hydrogenic
 */