Changeset 14862

Timestamp:

05/02/13 15:21:57 (12 years ago)

Author:

bdef

Message:

NEW: export md.mesh to gmsh format

File:

• issm/trunk-jpl/src/m/mesh/ExportGmsh.m (modified) (1 diff)

issm/trunk-jpl/src/m/mesh/ExportGmsh.m ¶

@@ -9 +9 @@
 fid=fopen(filename,'w');
 
-%initialiation
-fprintf(fid,'\n%s\n%i\n','MeshVersionFormatted',1);
+%initialiaion
+fprintf(fid,'$MeshFormat \n');
+fprintf(fid,'2.2 0 8 \n');
+fprintf(fid,'$EndMeshFormat \n');
+fprintf(fid,'$Nodes \n');
+fprintf(fid,'%i \n',md.mesh.numberofvertices);
+np=0;
+%printing point positions
+for ii=1:md.mesh.numberofvertices
+        np = np+1;
+        fprintf(fid,'%g %14.7e %14.7e 0.0 \n',np,md.mesh.x(np),md.mesh.y(np));
+end
 
-%dimension
-fprintf(fid,'\n%s\n%i\n','Dimension',2);
+fprintf(fid,'$EndNodes \n');
+fprintf(fid,'$Elements \n');
+fprintf(fid,'%i \n',md.mesh.numberofelements+size(md.mesh.segments,1));
+np=0;
 
-%Vertices
-fprintf(fid,'\n%s\n%i\n\n','Vertices',md.mesh.numberofvertices);
-fprintf(fid,'%8g %8g %i\n',[md.mesh.x md.mesh.y zeros(md.mesh.numberofvertices,1)]');
+%printing elements caracteristics for boundaries
 
-%Triangles
-fprintf(fid,'\n\n%s\n%i\n\n','Triangles',md.mesh.numberofelements);
-fprintf(fid,'%i %i %i %i\n',[md.mesh.elements zeros(md.mesh.numberofelements,1)]');
-
+for ii=1:size(md.mesh.segments,1)
+        np = np+1;
+        if(md.mesh.x(md.mesh.segments(np,1))==max(md.mesh.x(:))&&md.mesh.x(md.mesh.segments(np,2))==max(md.mesh.x(:))),
+                bc_id=1;
+        elseif(md.mesh.y(md.mesh.segments(np,1))==max(md.mesh.y(:))&&md.mesh.y(md.mesh.segments(np,2))==max(md.mesh.y(:))),
+                bc_id=2;
+        elseif(md.mesh.x(md.mesh.segments(np,1))==min(md.mesh.x(:))&&md.mesh.x(md.mesh.segments(np,2))==min(md.mesh.x(:))),
+                bc_id=3;
+        elseif(md.mesh.y(md.mesh.segments(np,1))==min(md.mesh.y(:))&&md.mesh.y(md.mesh.segments(np,2))==min(md.mesh.y(:))),
+                bc_id=4;
+  end
+                fprintf(fid,'%g 1 2 %g 1 %g %g \n',np,bc_id,md.mesh.segments(np,1),md.mesh.segments(np,2));
+end
+%and for the body
+body_id=1;
+for ii=1:md.mesh.numberofelements
+  np = np+1;
+  fprintf(fid,'%g 2 2 %g 3 %g %g %g \n',np,body_id,md.mesh.elements(ii,1),md.mesh.elements(ii,2),md.mesh.elements(ii,3));
+end
+fprintf(fid,'$EndElements \n');
 %close
 fclose(fid);
 t2=clock;fprintf('%s\n',[' done (' num2str(etime(t2,t1)) ' seconds)']);
+