Don't!! initialise the domain to be a tiny RectangleMesh by default. When you run on
many processors, it errors because that mesh cannot be distributed among so many
cores (most of them wind up empty, and crash). This explains the mysterious

Process 0: Number of global vertices: 84
Process 0: Number of global cells: 120

we've been seeing on every run (no matter how big the domain).