Dear all,

We found a little bit serious bug in the present version.

The problem is on Line 359 of tephra_calc.c.

Present;

temp1 = GRAV_SQRD_x_4 * part_density * part_density / AIR_VISCOSITY_x_225 * rho; 

Correct;

temp1 = GRAV_SQRD_x_4 * part_density * part_density / AIR_VISCOSITY_x_225 / rho; 

The difference is * rho or / rho.

This problem causes slower fall of lapilli-ash size particles (-3 to +2 phi), especially in high altitude.

I recommend you to edit the line before you run your projects.

All the best,

Kaz