Output solution across surface of mesh with element coordinates to a .txt file

Ok, I have been able to successfully extract the approximation of the normals at the facets, but I am now trying to extract the individual components for this list.

from dolfin import *
import numpy as np
from mpi4py import MPI


mesh = UnitCubeMesh(1,1,1)
Space = FunctionSpace(mesh,'P', 1)

n = FacetNormal(mesh)
V = VectorFunctionSpace(mesh, "CG", 1)
u = TrialFunction(V)
v = TestFunction(V)
a = inner(u,v)*ds
l = inner(n, v)*ds
A = assemble(a, keep_diagonal=True)
L = assemble(l)

A.ident_zeros()
nh = Function(V)

solve(A, nh.vector(), L)

[nx,ny,nz] = nh.split(True)

print(nx, ny, nz)
File("nh.pvd") << nh

See: How to get the solution obatained by fenics as a vector? - #7 by dokken