Changeset 21676


Ignore:
Timestamp:
04/18/17 10:08:39 (8 years ago)
Author:
Mathieu Morlighem
Message:

CHG: fixing amr indentation

Location:
issm/trunk-jpl/src/m/classes
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/classes/amr.py

    r21675 r21676  
    44
    55class amr(object):
    6    """
    7    AMR Class definition
     6    """
     7    AMR Class definition
    88
    9       Usage:
    10          amr=amr();
    11    """
     9    Usage:
     10        amr=amr();
     11    """
    1212
    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.
    1717
    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): # {{{
    2029
    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
    3032
    31                 #level_max: 2 to 4
    32       self.level_max=2
     33        #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.
    3335
    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.
    3638
    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")
    4348
    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  
    1010                Edges=[];
    1111                Triangles=[];
    12                 Quadrilaterals=[];
    1312                IssmEdges=[];
    1413                IssmSegments=[];
  • issm/trunk-jpl/src/m/classes/bamgmesh.py

    r21303 r21676  
    1313                self.Edges=np.empty((0,3))
    1414                self.Triangles=np.empty((0,0))
    15                 self.Quadrilaterals=np.empty((0,0))
    1615                self.IssmEdges=np.empty((0,0))
    1716                self.IssmSegments=np.empty((0,0))
     
    4544                s+="    Edges: %s\n" % str(self.Edges)
    4645                s+="    Triangles: %s\n" % str(self.Triangles)
    47                 s+="    Quadrilaterals: %s\n" % str(self.Quadrilaterals)
    4846                s+="    IssmEdges: %s\n" % str(self.IssmEdges)
    4947                s+="    IssmSegments: %s\n" % str(self.IssmSegments)
Note: See TracChangeset for help on using the changeset viewer.