Re: [vmtk-users] ReadTecplot surface data into VMTK

[Luca A. <luc...@gm...>]

Hello Yong,
 I've seen the file, there are a couple of issues.

First, the current Tecplot reader in vmtk is quite basic. It expects the header to be similar to the following:

VARIABLES = X,Y,Z,Something,SomethingElse
ZONE N=388,E=564,F=FEPOINT,ET=TRIANGLE

You have a more elaborate header:

TITLE     = "fluent12.0.16"
VARIABLES = "X"
"Y"
"Z"
"WSS dyne/cm2"
DATASETAUXDATA Common.PressureVar="4"
DATASETAUXDATA Common.UVar="5"
DATASETAUXDATA Common.VectorVarsAreVelocity="TRUE"
DATASETAUXDATA Common.VVar="6"
DATASETAUXDATA Common.WVar="7"
ZONE T="ew"
 STRANDID=1, SOLUTIONTIME=0.47499998
 Nodes=5850, Elements=5750, ZONETYPE=FEQuadrilateral
 DATAPACKING=POINT
 AUXDATA Common.BoundaryCondition="Wall"
 AUXDATA Common.IsBoundaryZone="TRUE"
 DT=(DOUBLE DOUBLE DOUBLE DOUBLE )

I reduced it to
VARIABLES = "X", "Y", "Z", "WSS dyne/cm2"
ZONE T="ew", N=5850, E=5750, ZONETYPE=FEQuadrilateral, DATAPACKING=POINT, AUXDATA Common.BoundaryCondition="Wall", AUXDATA Common.IsBoundaryZone="TRUE", DT=(DOUBLE DOUBLE DOUBLE DOUBLE )
 1.200498104E-02 8.789545441E-02 5.132280884E-01 1.406299010E+01
...

just to make the reader happy and I got something (see screenshot), but the current reader only supports triangles, not quads, so the surface isn't correct. Probably if you convert the quads to triangles within Tecplot (I don't know Tecplot, so I can't tell you how) the process will work.

Just so you know, Paraview (www.paraview.org) reads tec files. So I tried that too, but there was an error message saying that STRANDID is not currently supported. Can you save the file with an earlier Tecplot format version?

Luca




On May 1, 2013, at 5:09 PM, Yong He wrote:

> Hi VMTK users,
> 
> I want to use the VMTK mapping and patching scripts based on the ascii data file that is output from Tecplot 360 2011 (file format version 2009). I tried to use vmtksurfacereader to read the Tecplot .dat file and vmtksurfacewriter to save as a .vtp file, but I got the following error,
> 
> File "/usr/local/lib/vmtk/vmtk/vmtksurfacereader.py", line 98, in ReadTecplotSurfaceFile
>     self.PrintLog("Reading " + str(numberOfNodes)+" nodes.") UnboundLocalError: local variable 'numberOfNodes' referenced before assignment.
> 
> I know similar issues about reading Tecplot surface data into VMTK have been discussed before, but I still don't know how to proceed due to my limited knowledge of VTK, Python and VMTK. I attached a sample Tecplot data file, which contains the X, Y, Z and WSS data. I also tried to use Paraview and VisIt to convert to a vtp file, but neither worked. Hope you could guide me. Thank you very much.
> 
> Yong
> <tryvmtkpt2009.dat>------------------------------------------------------------------------------
> Introducing AppDynamics Lite, a free troubleshooting tool for Java/.NET
> Get 100% visibility into your production application - at no cost.
> Code-level diagnostics for performance bottlenecks with <2% overhead
> Download for free and get started troubleshooting in minutes.
> http://p.sf.net/sfu/appdyn_d2d_ap1_______________________________________________
> vmtk-users mailing list
> vmt...@li...
> https://lists.sourceforge.net/lists/listinfo/vmtk-users