[cig-commits] r18492 - seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab

[liuqy at geodynamics.org]

Author: liuqy
Date: 2011-05-30 14:03:22 -0700 (Mon, 30 May 2011)
New Revision: 18492

Modified:
   seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/read_cubit_basic.m
   seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/run_cubit2specfem2d.m
Log:

Fix iunit in run_cubit2specfem2d.py so that it indeed corresponds the number of material ids present.



Modified: seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/read_cubit_basic.m
===================================================================
--- seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/read_cubit_basic.m	2011-05-30 20:31:49 UTC (rev 18491)
+++ seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/read_cubit_basic.m	2011-05-30 21:03:22 UTC (rev 18492)
@@ -32,14 +32,15 @@
     ncol1 = 3; ncol2 = 8; icol1 = [2 3 4]; i2D = 0;
 end
 icol2 = 2:(ncol2+1);
+disp(sprintf('read_cubit_basic.m: %s', ifile));
 
 lines = textread(ifile,'%s','delimiter','\n');
 
 % line numbers for start and finish
 in1 = 10;
 in2 = in1 + nnod - 1;
-ie1 = in2 + 4;
-ie2 = ie1 + nele - 1 + (nunit-1);
+ie1 = in2 + 3;
+ie2 = ie1 + nele + nunit - 1;
 
 % read in nodes
 dnode = zeros(nnod,ncol1);
@@ -62,7 +63,7 @@
 
 % read in elements
 iele = zeros(nele,ncol2);
-mele = zeros(nele,1); im = 1;
+mele = zeros(nele,1); im = 0;
 ieleind = zeros(nele,1);  % index listed by cubit
 jj = 0;
 imax = ie2-ie1+1;

Modified: seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/run_cubit2specfem2d.m
===================================================================
--- seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/run_cubit2specfem2d.m	2011-05-30 20:31:49 UTC (rev 18491)
+++ seismo/2D/SPECFEM2D/trunk/UTILS/cubit2specfem2d/matlab/run_cubit2specfem2d.m	2011-05-30 21:03:22 UTC (rev 18492)
@@ -19,7 +19,7 @@
 sfac = 1;
 icol = [2 4];
 nbound = 2;
-nnod = 2673; nele = 2560; stag = 'Tromp2005'; iunit = [1 2];
+nnod = 2673; nele = 2560; stag = 'Tromp2005'; iunit = [1];
 dir0 = '/data/svn/seismo/2D/SPECFEM2D/UTILS/cubit2specfem2d/matlab/';
 
 ifile1 = [dir0 stag '_elements_nodes'];