Adding msg_stack to BoutMesh::communicate

Helps with debugging, as shows when something goes
wrong in the communication. Common causes are uninitialised
data being communicated, so getData call fails.
Tokamak geometry test case

o Advection, Delp2 and Laplace_par operators seem to work
  in ballooning coordinates or shifted metric.

o Test generates grid files from python script as needed