Zero values using sampleSurface
When I execute sampleSurface (constantPlane, interpolatedPlane) on my parallel case and everything runs fine. After I try to look at the results in paraview 2.4.2 I already get an error message in my console:
Error or warning: There was a VTK Warning in file: ../../../../../paraview-2.4.2/VTK/IO/vtkDataReader.cxx (844) Error reading ascii data! The field to plot is zero. The "Data Range" in paraview also say that everything is zero. But when I apply foamToVTK to the same case I can do a cut with a plane at the same position and everything looks fine. Also paraFoam gives the correct picture. I use OpenFOAM-1.3 on FC5_X86_64 and here is my sampleSurfaceDict: =============================================== surfaceFormat vtk; interpolationScheme cellPointFace; surfaces ( constantPlane { name constantPlane; basePoint (0 0 0); normalVector (0 0 1); // Optional: whether to leave as faces or triangulate (=default) triangulate false; } interpolatedPlane { name interpolatedPlane; basePoint (0 0 0); normalVector (0 0 1); triangulate false; } ); ( p U F cI cD nuSgs ); ================================================ I experienced a similar problem using the probes utility, i.e. zero values eventhough they are not zero. They probably use the same interpolation method. I also checked running in serial and different interpolation methods. Any suggestions are welcome. |
Change your plane so it doesn'
Change your plane so it doesn't coincide with any faces/edges/points. Just perturb by small amount.
If problem persists please create small testcase and send to me. |
Thanks!
Using a slightly diff
Thanks!
Using a slightly different basepoint for each of 2 coordinates solved the problem, but nevertheless it seems to be very sensitive especially if you have a structured mesh. I had to play around several times to get sampleSurface the right output. |
All times are GMT -4. The time now is 20:48. |