/*plankf evaluate Planck function for any cell at current electron temperature */
double plankf(long int ip);