Changeset 21676
- Timestamp:
- 04/18/17 10:08:39 (8 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/amr.py
r21675 r21676 4 4 5 5 class amr(object): 6 """7 AMR Class definition6 """ 7 AMR Class definition 8 8 9 10 11 """9 Usage: 10 amr=amr(); 11 """ 12 12 13 def __init__(self): # {{{14 self.level_max= 0.15 self.region_level_1= 0.16 self.region_level_max= 0.13 def __init__(self): # {{{ 14 self.level_max = 0. 15 self.region_level_1 = 0. 16 self.region_level_max = 0. 17 17 18 #set defaults 19 self.setdefaultparameters() 18 #set defaults 19 self.setdefaultparameters() 20 #}}} 21 def __repr__(self): # {{{ 22 string=" amr parameters:" 23 string="%s\n%s"%(string,fielddisplay(self,"level_max","maximum refinement level (1, 2, 3 or 4)")) 24 string="%s\n%s"%(string,fielddisplay(self,"region_level_1","region which will be refined once (level 1) [ m ]")) 25 string="%s\n%s"%(string,fielddisplay(self,"region_level_max","region which will be refined with level_max [ m ]")) 26 return string 27 #}}} 28 def setdefaultparameters(self): # {{{ 20 29 21 #}}} 22 def __repr__(self): # {{{ 23 string=" amr parameters:" 24 string="%s\n%s"%(string,fielddisplay(self,"level_max","maximum refinement level (1, 2, 3 or 4)")) 25 string="%s\n%s"%(string,fielddisplay(self,"region_level_1","region which will be refined once (level 1) [ m ]")) 26 string="%s\n%s"%(string,fielddisplay(self,"region_level_max","region which will be refined with level_max [ m ]")) 27 return string 28 #}}} 29 def setdefaultparameters(self): # {{{ 30 #level_max: 2 to 4 31 self.level_max=2 30 32 31 #level_max: 2 to 432 self.level_max=233 #region_level_1: region around (m) the discontinuity (grounding line or ice front) where the mesh will be refined once (h=1). 34 self.region_level_1=20000. 33 35 34 #region_level_1: region around (m) the discontinuity (grounding line or ice front) where the mesh will be refined once (h=1).35 self.region_level_1=20000.36 #region_level_max: region around (m) the discontinuity (grounding line or ice front) where the mesh will be refined with max level of refinement (h=level_max). 37 self.region_level_max=15000. 36 38 37 #region_level_max: region around (m) the discontinuity (grounding line or ice front) where the mesh will be refined with max level of refinement (h=level_max). 38 self.region_level_max=15000. 39 40 return self 41 #}}} 42 def checkconsistency(self,md,solution,analyses): # {{{ 39 return self 40 #}}} 41 def checkconsistency(self,md,solution,analyses): # {{{ 42 md = checkfield(md,'fieldname','amr.level_max','numel',[1],'>=',0,'<=',4) 43 md = checkfield(md,'fieldname','amr.region_level_1','numel',[1],'>',0,'NaN',1,'Inf',1) 44 md = checkfield(md,'fieldname','amr.region_level_max','numel',[1],'>',0,'NaN',1,'Inf',1) 45 #it was adopted 20% of the region_level_1 46 if self.region_level_1-self.region_level_max<0.2*self.region_level_1: 47 md.checkmessage("region_level_max should be lower than 80% of region_level_1") 43 48 44 md = checkfield(md,'fieldname','amr.level_max','numel',[1],'>=',0,'<=',4) 45 md = checkfield(md,'fieldname','amr.region_level_1','numel',[1],'>',0,'NaN',1,'Inf',1) 46 md = checkfield(md,'fieldname','amr.region_level_max','numel',[1],'>',0,'NaN',1,'Inf',1) 47 #it was adopted 20% of the region_level_1 48 if self.region_level_1-self.region_level_max<0.2*self.region_level_1: 49 md.checkmessage("region_level_max should be lower than 80% of region_level_1") 50 51 return md 52 # }}} 53 def marshall(self,prefix,md,fid): # {{{ 54 55 WriteData(fid,prefix,'object',self,'fieldname','level_max','format','Integer') 56 WriteData(fid,prefix,'object',self,'fieldname','region_level_1','format','Double') 57 WriteData(fid,prefix,'object',self,'fieldname','region_level_max','format','Double') 58 # }}} 49 return md 50 # }}} 51 def marshall(self,prefix,md,fid): # {{{ 52 WriteData(fid,prefix,'object',self,'fieldname','level_max','format','Integer') 53 WriteData(fid,prefix,'object',self,'fieldname','region_level_1','format','Double') 54 WriteData(fid,prefix,'object',self,'fieldname','region_level_max','format','Double') 55 # }}} -
issm/trunk-jpl/src/m/classes/bamgmesh.m
r18455 r21676 10 10 Edges=[]; 11 11 Triangles=[]; 12 Quadrilaterals=[];13 12 IssmEdges=[]; 14 13 IssmSegments=[]; -
issm/trunk-jpl/src/m/classes/bamgmesh.py
r21303 r21676 13 13 self.Edges=np.empty((0,3)) 14 14 self.Triangles=np.empty((0,0)) 15 self.Quadrilaterals=np.empty((0,0))16 15 self.IssmEdges=np.empty((0,0)) 17 16 self.IssmSegments=np.empty((0,0)) … … 45 44 s+=" Edges: %s\n" % str(self.Edges) 46 45 s+=" Triangles: %s\n" % str(self.Triangles) 47 s+=" Quadrilaterals: %s\n" % str(self.Quadrilaterals)48 46 s+=" IssmEdges: %s\n" % str(self.IssmEdges) 49 47 s+=" IssmSegments: %s\n" % str(self.IssmSegments)
Note:
See TracChangeset
for help on using the changeset viewer.