Dear all, I want to create a C0-continuous field from interpolated input data. The input field is obtained from a FVM code and was filtered with a gaussian in space and interpolated to the quadrature points. I store the field in a ContField2D m_bf = Array<OneD, MultiRegions::ExpListSharedPtr>(m_spacedim + 2); m_bf[i] = MemoryManager<MultiRegions::ContField2D>::AllocateSharedPtr(m_session, m_graph); copy the interpolated data into its m_phys array Vmath::Vcopy(nq, filteredField[i], 1, m_bf[i]->UpdatePhys(), 1); and then try to smooth the field. Until now, i just did: m_bf[i]->IProductWRTBase(m_bf[i]->GetPhys(), m_bf[i]->UpdateCoeffs()); m_bf[i]->MultiplyByElmtInvMass(m_bf[i]->GetCoeffs(), m_bf[i]->UpdateCoeffs()); m_bf[i]->BwdTrans(m_bf[i]->GetCoeffs(), m_bf[i]->UpdatePhys()); which yields the field depicted in DG.png. However, this does not ensure C0- continuity, so i tried to call ExpList::SmoothField afterward: m_bf[i]->SmoothField(m_bf[i]->UpdatePhys()); However, this produces a distorted field as shown in smoothed.png. I am using a 6th order modified expansion basis: <E COMPOSITE="C[10]" NUMMODES="6" FIELDS="p,u,v,w" TYPE="MODIFIED"/> and the mesh is depicted in smoothed_mesh.png. What is the cause of this problem and what can i do do fix it? Cheers, Kilian -- Kilian Lackhove, M.Sc. Fachgebiet für Energie- und Kraftwerkstechnik L1|08 114 Technische Universität Darmstadt Jovanka-Bontschits-Straße 2 D-64287 Darmstadt Germany Tel.: +49 6151 16 - 28915 Fax: +49 6151 16 - 6555 e-mail: lackhove@ekt.tu-darmstadt.de