Hi
I have a fortran code that I'm trying to get streamlines from. The code
generates a spherical grid of points and a velocity on each of those
points.
I have been using a structured grid but paraview doesn't seem to read the
vectors, instead it thinks they are more grid points and plots them as
such. From what I've read it looks like I need to use an unstructured grid
but I don't know how to set up the cell commands for that.
Here are the lines of my code that write the .vts files:
write(113,*) '# vtk DataFile Version 3.0'
write(113,*) 'This is the cartesian velocity field'
write(113,*) 'ASCII'
write(113,*) 'DATASET STRUCTURED_GRID'
write(113,*) 'DIMENSIONS',nx,ny,nz
do k = 1, nphi + 1
do j = 1, nth + 2
do i = 1, nr
write(113,*) xc(i,j,k), yc(i,j,k), zc(i,j,k)
enddo
enddo
enddo
write(113,*) 'POINT_DATA',nx
write(113,*) 'VECTORS velocityfield float'
write(113,*) 'LOOKUP_TABLE default'
do k = 1, nphi + 1
do j = 1, nth + 2
do i = 1, nr
write(113,*) ux(i,j,k), uy(i,j,k), uz(i,j,k)
enddo
enddo
enddo
Could anyone let me know how I should modify this for an unstructured grid?
Thanks
George
_______________________________________________
Powered by www.kitware.com
Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html
Please keep messages on-topic and check the ParaView Wiki at:
http://paraview.org/Wiki/ParaView
Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview
