/*dcharg compute grains charge */
void dcharg(void);