Hi: i want write my results into VTK format, then i can use the paraview. but i have some problem. my compute model is 700KM*700KM, there are 71*71 nodes. it's a Quadrilateral mesh.
my matlab code as followed % Model size, m xsize=700000; ysize=700000; % Defining resolution xnum=71; ynum=71; NP=xnum*ynum; NC=xnum*ynum-2*xnum+1; % Defining gridline positions for regular basic grid % In horizontal direction gridx=zeros(xnum,1); for i=2:1:xnum gridx(i)=gridx(i-1)+xstp; end % Defining gridline positions for 51x51 irregular basic grid % In vertical direction gridy=zeros(ynum,1); for i=2:1:6 gridy(i)=gridy(i-1)+20000; end for i=7:1:46 gridy(i)=gridy(i-1)+1250; end for i=46:1:51 gridy(i)=gridy(i-1)+10000; end for i=51:1:71 gridy(i)=gridy(i-1)+25000; end gridy(ynum)=ysize; %write the points a=0; for i=1:1:ynum for j=1:1:xnum a=a+1; Points(1,a)=gridx(j); Points(2,a)=gridy(i); end end % write the cell nn=0; for j=1:1:ynum for i=1:1:xnum nn=nn+1; Cells(1,nn)=(j-1)*(xnum+1)+i; Cells(2,nn)=(j-1)*(xnum+1)+i+1; Cells(3,nn)=j*(xnum+1)+i+1; Cells(4,nn)=j*(xnum+1)+i; Cells(5,nn)=nn; end end % write VTK file fp=fopen('Resault.vtk','w'); fprintf(fp,'# vtk DataFile Version 3.0\n'); fprintf(fp,'vtk output\n'); fprintf(fp,'ASCII\n'); fprintf(fp,'DATASET STRUCTURED_GRID\n'); fprintf(fp,'POINTS %d float\n',NP); for i=1:NP fprintf(fp,'%e %e %e \n',Points(1,i),Points(2,i),0); end fprintf(fp,'\n'); fprintf(fp,'CELLS %d %d\n',NC,NC*4); aa=0; for i=1:NC aa=aa+1; fprintf(fp,'%d %d %d %d %d\n',aa,Cells(1,i)-1,Cells(2,i)-1,Cells(3,i)-1,Cells(4,i)-1); end fprintf(fp,'\n'); fprintf(fp,'CELL_TYPES %d\n',NC); for i=1:NC fprintf(fp,'%d\n',4); end fprintf(fp,'\n'); fprintf(fp,'POINT_DATA %d\n',NP); when i tried to open it with paraview. it sayed "ERROR: In C:\berk\ReleaseBuilds\ParaView3\VTK\IO\vtkDataReader.cxx, line 2078 ERROR: In C:\berk\ReleaseBuilds\ParaView3\VTK\IO\vtkDataReader.cxx, line 2078 ERROR: In C:\berk\ReleaseBuilds\ParaView3\VTK\Parallel\vtkPDataSetReader.cxx, line 705 vtkPDataSetReader (0AEA4D18): Expecting 'DIMENSIONS' insted of: POINTS ERROR: In C:\berk\ReleaseBuilds\ParaView3\VTK\Parallel\vtkPDataSetReader.cxx, line 705 vtkPDataSetReader (0AEA4D18): Expecting 'DIMENSIONS' insted of: POINTS" What's the problem? could you help me? thank you very much. Best regards! -- Yanchao Qiao Ph.D student of Geodynamics +86-136 010 926 41 College of Earth Science Graduate University of Chinese Academy of Sciences 19A, YuQuan Road, Beijing, China, 100049 Email:oliver.q...@163.com <email%3aoliver.q...@163.com>; oliver.q...@gmail.com;qiaoyancha...@mails.gucas.ac.cn
_______________________________________________ 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