/*dintg compute total radiative cooling due to large grains */
double dintg(double tdust, 
  long int nd, 
  long int ip);