/*dcolid compute grains collisional heating cooling */
void dcolid(float *dcheat, 
  float *dccool, 
  long int nd);