function [da,err] = PetscDMDACreate2d(wrap,stentype,M,N,mm,nn,dof,s)
  da = PetscDM();
  err = da.SetType('da');
  da.SetPeriodicity(wrap);
  da.SetSizes([M N 1]);
  da.SetDof(dof);
  da.SetStencilType(stentype);
  da.SetStencilWidth(s);
  da.SetDim(2);
  da.SetFromOptions();
  err = da.SetUp();