Hello,

I am a FEniCS beginner. I am implementing Phase Field PDEs in 1-D. The initial condition is profile for each phase. Following would be my initial conditions:

eta_Cu = np.zeros([256]) # Cu

eta_Cu3Sn = np.zeros([256]) # Cu3Sn

eta_Cu6Sn5 = np.zeros([256]) # Cu6Sn5

eta_Sn = np.zeros([256]) # Sn

eta = [eta_Cu, eta_Cu3Sn, eta_Cu6Sn5, eta_Sn]

eta[0][0:108] = 1

eta[1][108:128] = 1

eta[2][128:148] = 1

eta[3][148:256] = 1

eta = np.reshape(eta, (4, 256))

So, I want these initial conditions to be imbeded in, mesh = IntervalMesh(256, 0, 256e-6).

I request some insight about how to do it. Thank you.