Originally Posted by lbmagis (Post 525592)

Hello,
Could any one tell us how to use this way to visualize isotherms and streamlines in 2Dimensions.
Thank you

Originally Posted by sarah440 (Post 410899)

Hi,

I have a DNS database which includes the coordinate matrices x, y, z as well as the velocity fields u, v and w. All matrices have the same size of nl*nc*ns. Then I use Matlab to compute some new parameters like swirling rate. At the end of my own script I use mat2tecplot.m to convert mat files to PLT format. Here is what I put at the end of my script:

PLTout = strcat(outputAddress,PLTout_name); %The output name and address
tdata=[];
tdata.Nvar=7; %number of variables
tdata.vformat = 2*ones(1,tdata.Nvar); %Double precision for each variable (optional)
tdata.varnames={'x','y','z','u','v','w','Lambda'};
tdata.cubes(1).zonename=eval('PLTout_name'); %any name
tdata.cubes(1).x=x(nl*nc*ns);
tdata.cubes(1).y=y(nl*nc*ns);
tdata.cubes(1).z=z(nl*nc*ns);
tdata.cubes(1).v(1,:,:,: )=u(nl*nc*ns);
tdata.cubes(1).v(2,:,:,: )=v(nl*nc*ns);
tdata.cubes(1).v(3,:,:,: )=w(nl*nc*ns);
tdata.cubes(1).v(4,:,:,: )=lambda(nl*nc*ns);
mat2tecplot(tdata,eval('PLTout'))

I can easily visulize this PLT file in Tecplot
Hope this helps

Originally Posted by shhe (Post 533732)

Hi,
I have tried to do as you but unfortunately, it failed to generate the tecplot file. My code as following:

tdata=[];
tdata.Nvar=5; %number of variables
tdata.vformat = 2*ones(1,tdata.Nvar); %Double precision for each variable (optional)
tdata.varnames={'X','Y','Z','TEMP','TEMP_GRADIENT_ NON_UNITY'};
tdata.cubes(1).zonename='TEMP_GRADIENT_NON_UNITY'; %any name
%tdata.cubes(1).zonename=eval('TEMP_GRADIENT_NON_U NITY'); %any name
tdata.cubes(1).x=X;
tdata.cubes(1).y=Y;
tdata.cubes(1).z=Z;
tdata.cubes(1).v(1,:,:,: )=TEMP;
tdata.cubes(1).v(2,:,:,: )=TEMP_GRADIENT_NON_UNITY;
mat2tecplot(tdata,'TEMP_GRADIENT_NON_UNITY.plt')


I have gotten the following error:
Error : cube 1 v 2nd, 3rd, 4th dimension
must be of the same size as x, y, z

To your knwledge, my domain is 201*201*201

Originally Posted by sarah440 (Post 533798)

Hi,

It seems that the sizes of your arrays are not the same. Here is the error display lines in the script:

if(Imax_v~=Imax || Jmax_v~=Jmax || Kmax_v ~= Kmax)
s=-1;
display(['Error : cube ' int2str(icube) ' v 2nd, 3rd, 4th dimension']);
display([' must be of the same size as x, y, z']);

Are you sure that your arrays (including x, y, z) are all 3D with the same size?

Originally Posted by sarah440 (Post 533855)

What I see here is that X, Y, Z are 201*1, you need them to be 201*201*201 like TEMP. First, use "meshgrid" function to convert X, Y, Z to rectangular grid. It should work then.

Originally Posted by Hoggs17 (Post 405985)

After reading through the mat2tecplot.m documentation, this is what I did in order to export my matlab data into a readable Tecplot file. Hope this will help any other users out there that aren't quite using it right. This is just my interpretation of Wen Long's code; all credit must be given to him for writing it.

1. You must keep the mat2tecplot.m in the same directory as your other files as your main program will call mat2tecplot.m (should be self-explanatory but you know never know).

2. In your program, you must call mat2tecplot.m and specify what data you are transferring.

3. Following the examples given in the documentation (I followed example 3 I believe), you can change a few lines of code to reflect which data you want to export. My code is as follows.

function name(inputs)

tdata=[];
tdata.Nvar=6;
tdata.varnames={'x','y','z','u','v','p_new'};
tdata.surfaces(1).zonename='mysurface zone';
tdata.surfaces(1).x=x;
tdata.surfaces(1).y=y;
tdata.surfaces(1).v(1,:,: )=u;
tdata.surfaces(1).v(2,:,: )=v;
tdata.surfaces(1).v(3,:,: )=p_new;
mat2tecplot(tdata,'SIMPLE Results.plt')

end

You have to make sure that the value specified in tdata.Nvar is equal to the number of variables you place in tdata.varnames. I specified 6 variables and then had 6 names given. I have 2-D flow (x,y,z=0) and want my data for u and v momentum and pressure exported. "tdata.surfaces(1).x= " is used to specify your domain. In the example, coordinates are given for x,y,z; in my code; i simply set the .x value equal to my x-values on my defined mesh. The same was done for the y-values. if you specify nothing (as with z), it is assumed 0.

The "tdata.surfaces(1).v(1,:,: )= " code specifies which data you want displayed as a result. i set it equal to my u momentum. To add multiple data to the same file, you change the .v(1,:,: ) to .v(2,:,: ) etc... for as many different data sets as you like. I had 3 data sets (u,v,pressure) so I have 3 lines.

Finally, the last line in the code calls mat2tecplot using the values you have specified and outputs an easily readable tecplot input file (called 'SIMPLE Results'). My contour plots for momentum and pressure look much better than they did on matlab plus you now have the flexibility of using tecplot to manipulate your visualization axis, contour # etc...

You just have to add a few lines of code to your matlab program and the mat2tecplot.m does all the work for you.

Now if I could just figure out how to transfer my streamline data I'd be all set - anyone have any tips?