Another change to calculate the lumped mass matrix via a par_loop.
Only support spherical meshes

Cartesian meshes are not supported any longer.
Reformatted comments and code for better readability.