In the diffusion creep database, the A parameter has units of Pa^(-1) * m^3 * s^(-1). It seems like the A parameter it is already multiplied with the grainsize (with exponent) here. However later, in the compute_εII and compute_τII functions, the grainsize dependence is accounted for again.