Hi,

I am trying to save the results of a 1D problem with mixed elements. After using dofmap().dofs() the order is still wrong when plotting. The code is shown below:

```
import matplotlib.pyplot as plt
from fenics import *
import numpy as np
rho_p = 3000.0
d_p = 100.0e-6;
m_p = rho_p*pi*d_p**3/6.0
l = 30.0e-3
numElems = 100
abc = 2.0
dfg = abc/60000.0
n = dfg/m_p
r_0 = 1.0e-3
h_d = 17.0e-3
U = 0.017
Dz_F = 10.0e-3;
PI = 3.1415926535
theta = 10.95*PI/180.0
# Define function space for system
deg = 2
mesh = IntervalMesh(numElems, 0.0,l/2.0)
P1 = FiniteElement('Lagrange', interval, deg)
element = MixedElement([P1, P1])
V = FunctionSpace(mesh, element)
R = Expression("r_0 + (h_d-Dz_F/2.0)*tan(theta)", degree=deg, r_0=r_0, h_d=h_d, Dz_F=Dz_F, theta=theta)
u_ic='tan(theta)*U*x[0]/R'
c_ic='2.0*n/U/pi/pow(R,2)*exp( -2.0*pow(x[0],2)/pow(R,2) )'
w_0 = Expression((u_ic, c_ic), degree=deg, theta=theta, U=U, R=R, n=n, pi=PI)
#Write values in text file
w_p = project(w_0, V)
dof_coordinates = V.tabulate_dof_coordinates()
n = V.dim()
d = mesh.geometry().dim()
dof_coordinates.resize((n, d))
u_dofs = V.sub(0).dofmap().dofs()
c_dofs = V.sub(1).dofmap().dofs()
plt.plot(dof_coordinates[c_dofs], w_p.vector()[c_dofs])
plt.show()
```

Thank you for your time,

Alex