subroutine DMGlobalToLocalSolve(a,b,c,z)
       import tDM,tVec
       DM a ! DM
       Vec b ! Vec
       Vec c ! Vec
       PetscErrorCode z
       end subroutine DMGlobalToLocalSolve