Program test

  implicit none

#include "petsc/finclude/petsc.h"

  DM :: da3
  AO :: ao3
  PetscInt :: ierr,nx,ny,nz
  PetscInt, parameter :: i1=1
  PetscInt, parameter :: i2=2
  PetscInt, parameter :: i3=3

  call PetscInitialize(PETSC_NULL_CHARACTER,ierr)

  ! Create DMDA
  nx=10
  ny=10
  nz=10

  call DMDACreate3d(PETSC_COMM_WORLD,DM_BOUNDARY_NONE,DM_BOUNDARY_NONE,DM_BOUNDARY_NONE, &
    DMDA_STENCIL_BOX,nx,ny,nz,PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,i3,i2,PETSC_NULL_INTEGER, &
    PETSC_NULL_INTEGER,PETSC_NULL_INTEGER, da3,ierr)

  ! Set AO
  call DMDASetAOType(da3,AOBASIC,ierr)

  call PetscFinalize(ierr)


end program test