Here is another bug:
for(int i = 0; i < mesh_dim_z; i++) {
printf("%12.5e ", mesh(0,0)->f(i));
}
after pC++ becomes
for (int i = 0 ; i < mesh_dim_z ; (i)++)
{
pcxx_printf56("%12.5e ", (mesh)(0, 0)->f);
} ^^^^
where f(i) is an overloaded operator for a class.
Shelby