/*grngam compute grains photoionization rates */
double grngam(long int nd);