Quick question about "surf" or "surfc"

Hi

The problem is, I need to plot a 3D surface in each loop of my computations, lets say water depth. the computational domain is rectangular but it is not square, so the number of elements in X direction is not the same as the number of elements in Y direction.

for square domains I used to use the command :

    xx1=linspace(0,Xlength,NX);
    yy1=linspace(0,Ylength,NY);
    [x1,y1]=meshgrid(xx1,yy1);
    surf(x1,y1,z);
    view(125,20);

but when managing to use a non-square domain it gives the error:

??? Error using
==> surf at 78
Data dimensions
must agree.

Error in ==>
test at 53
surf(x1,y1,z);

removing the x1,y1 from the surf(x1,y1,z) command gives a 3D with the wrong scale and it is not useful at all

, so please advice

Answer 1

[Sean de Wolski]

size(x1)
size(y1)
size(z)

Will show that not all three matrix sizes are the same. Remember meshgrid swaps the inputs unlike ndgrid

[seyedashraf]

sure the x1 and y1 sizes are not same, please check the link in my first post to see the 2D domain (x-y) I am managing to plot z on it

[Walter Roberson]

Yup, but that 2D domain was not produced by meshgrid.

Use ndgrid() instead of meshgrid() to do what you want.

[seyedashraf]

Thanks Sean de Wolski. Thanks Walter Roberson, that replacement made it work.