On Wed, Jul 31, 2013 at 4:59 PM, Jed Brown <jedbrown@mcs.anl.gov> wrote:

Chris Eldred <chris.eldred@gmail.com> writes:
> if (rank==0) then
> call DMPlexCreateFromCellList(PETSC_COMM_WORLD, 2, 2, 4, 3, 1 , cell_list
> , 2 , vertex_coords , dm, ierr)

The problem here is that after this line, dm is garbage on rank 1. You
want to start by creating the DM on COMM_WORLD, but have all the cells
on rank 0. Try the attached.