/*PhotoIonize fill array PhotoRate with photoionization rates for heavy elements */
void PhotoIonize(
	/* nlem is atomic number on C scale, 0 for H */
	long int nelem , 
	/* debugging flag to turn on print */
	int lgPrintIt );