/*ddrift compute grains drift velocity */
void ddrift(void);