/*phfit derive photoionization cross sectoins for first 30 elements */
void phfit(long int nz, 
	  long int ne, 
	  long int is, 
	  double e, 
	  float *s, 
	  int lghs);