I want to export the boundary value of velocity, exact nodes on the boundary labrle 1 and 2. I wrote the following:
int Nbtriangle=th.nbe;
for (int k=0; k<Nbtriangle; k++)
{
int L=th.be(k).label;
if (L==1 || L==2)
{
real xx=th.be(k)[1].x;
real yy=th.be(k)[1].y;
file1 << xx << " " << yy << " " << Cr(xx,yy)<< endl;
}
}
when I compared the profile with the paper it was different. then I exported the vtk plot of my domain to tecplot and wanted to plot velocity profile at inlet (Labels 1 and 2) the plot was exactly the same as paper. how I can write it to give me the right value?
Thank you for your response. Sorry, I didn’t understand if(abs(bcLab1vec[i]-1)<1e-10) because its name presents the boundary label but when I use it is not that and doesn’t show me anything.