PetscSection feature added for dof=1 but has not worked yet. Uncomment line 68 to use PetscSection.

Currently works only with dmda without using petscsection. In this case, dof = no. of dimension of the model.