PETSc Global Vector to local process vector

[Ignacio Martínez]

Hello,

I'm a beginner user, and I have a problem:

I have a routine in Fortran and I just want to parallelize few lines
with PETSc, I'm using DA in order to gather the results, from local
vectors to a global vector.
It works fine. But I need to return those values to the Fortran
routine, and I only get access to local process values.

I try to create a sequential vectors with an scatter to get that
values. But, PETSc returns an error, something like: "Cannot create
VECSEQ on more than one process!" (I try to create the vector only in
one process).

Exactly in using DAGlobalToNaturalAllCreate and  VecScatterBegin to
get global values in to a process.

And my question is if someone has tried to do something like this, and
how has done it. Or if could you give me some orientation?

Thank you in advance,
Ignacio