Assistance with Peclet Number Calculation in FreeFem++

Hi everyone,

I would like to ask if there is any previous code where the calculation of the Peclet number was addressed. I tried to calculate it myself, but I am not getting reasonable results. I think I might have made a mistake.

Thank you,

according to what I know

	//Peclet Number
	VNe hK = hTriangle ;
	VNe PeK = We*hK/2./D ;
Where We is the drift velocity and D is the diffusion coefficient 

I would wait for feedback