Hey all,

Baiscally, I want to calculate the pressure for square obstacles and I did this sub-routine (using MATLAB) but the results were not correct

P_vector=reshape(rho,Nx,Ny)/3; % Pressure contour

for j=1:Ny

for i=1:Nx

P(i,j)=P_vector((i-1)*Ny+j);

if P(i,j)==0

P(i,j)=nan;

end

end

end

Could anyone helps to figure out what is the mistake or have another idea or I might do mistake in this routine.

Cheers,

William