Changeset 22216 for issm/trunk-jpl/src/m/mesh

Timestamp:

11/02/17 16:07:10 (8 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing python

File:

• issm/trunk-jpl/src/m/mesh/bamg.py (modified) (2 diffs)

issm/trunk-jpl/src/m/mesh/bamg.py

@@ -72 +72 @@
         bamg_mesh=bamgmesh()
 
+        subdomain_ref = 1
+        hole_ref = 1
+
         # Bamg Geometry parameters {{{
         if options.exist('domain'):
@@ -99 +102 @@
                         bamg_geometry.Vertices=np.vstack((bamg_geometry.Vertices,np.vstack((domaini['x'][0:nods],domaini['y'][0:nods],np.ones((nods)))).T))
                         bamg_geometry.Edges   =np.vstack((bamg_geometry.Edges,np.vstack((np.arange(count+1,count+nods+1),np.hstack((np.arange(count+2,count+nods+1),count+1)),1.*np.ones((nods)))).T))
-                        if i:
-                                bamg_geometry.SubDomains=np.vstack((bamg_geometry.SubDomains,[2,count+1,1,1]))
+                        if i:
+                            bamg_geometry.SubDomains=np.vstack((bamg_geometry.SubDomains,[2,count+1,1,-subdomain_ref]))
+                            subdomain_ref = subdomain_ref+1;
+                        else:
+                            bamg_geometry.SubDomains=np.vstack((bamg_geometry.SubDomains,[2,count+1,1,0]))
 
                         # bamg_geometry.Vertices=np.hstack((bamg_geometry.Vertices,np.vstack((domaini['x'][0:nods].reshape(-1),domaini['y'][0:nods].reshape(-1),np.ones((nods))))))