/*csphot returns photoionization cross section from opacity stage using std pointers */
double csphot(long int inu, 
  long int ithr, 
  long int iofset);