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