/*dtempr compute grains temperature, heating */
void dtempr(void);