Ignore:
Timestamp:
02/12/19 06:10:51 (6 years ago)
Author:
bdef
Message:

CHG: shifting to py3 version of python interface (py2 compatible)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/classes/slr.py

    r23088 r23716  
    99        """
    1010        SLR class definition
    11        
     11
    1212                Usage:
    1313                  slr=slr()
    1414        """
    15        
    1615        def __init__(self): # {{{
    1716                self.deltathickness         = float('NaN')
    1817                self.sealevel               = float('NaN')
    19                 self.spcthickness           = float('NaN')
     18                self.spcthickness                                               = float('NaN')
    2019                self.maxiter                = 0
    2120                self.reltol                 = 0
    … …  
    2625                self.tide_love_k            = 0 #ideam
    2726                self.tide_love_h            = 0 #ideam
    28                 self.fluid_love             = 0 
    29                 self.equatorial_moi         = 0 
    30                 self.polar_moi              = 0 
     27                self.fluid_love             = 0
     28                self.equatorial_moi         = 0
     29                self.polar_moi              = 0
    3130                self.angular_velocity       = 0
    3231                self.rigid                  = 0
    … …  
    4443                self.requested_outputs      = []
    4544                self.transitions            = []
    46                
     45
    4746                #set defaults
    4847                self.setdefaultparameters()
    4948                #}}}
     49
    5050        def __repr__(self): # {{{
    5151                        string='   slr parameters:'
    52                         string="%s\n%s"%(string,fielddisplay(self,'deltathickness','thickness change: ice height equivalent [m]'))
     52                        string="%s\n%s"%(string,fielddisplay(self,'deltathickness','thickness change: ice height equivalent [m]'))
    5353                        string="%s\n%s"%(string,fielddisplay(self,'sealevel','current sea level (prior to computation) [m]'))
    5454                        string="%s\n%s"%(string,fielddisplay(self,'spcthickness','thickness constraints (NaN means no constraint) [m]'))
    … …  
    6464                        string="%s\n%s"%(string,fielddisplay(self,'equatorial_moi','mean equatorial moment of inertia [kg m^2]'))
    6565                        string="%s\n%s"%(string,fielddisplay(self,'polar_moi','polar moment of inertia [kg m^2]'))
    66                         string="%s\n%s"%(string,fielddisplay(self,'angular_velocity','mean rotational velocity of earth [per second]')) 
     66                        string="%s\n%s"%(string,fielddisplay(self,'angular_velocity','mean rotational velocity of earth [per second]'))
    6767                        string="%s\n%s"%(string,fielddisplay(self,'ocean_area_scaling','correction for model representation of ocean area [default: No correction]'))
    6868                        string="%s\n%s"%(string,fielddisplay(self,'steric_rate','rate of steric ocean expansion [mm/yr]'))
    69                         string="%s\n%s"%(string,fielddisplay(self,'Ngia','rate of viscous (GIA) geoid expansion (in mm/yr)')) 
     69                        string="%s\n%s"%(string,fielddisplay(self,'Ngia','rate of viscous (GIA) geoid expansion (in mm/yr)'))
    7070                        string="%s\n%s"%(string,fielddisplay(self,'Ugia','rate of viscous (GIA) bedrock uplift (in mm/yr)'))
    7171                        string="%s\n%s"%(string,fielddisplay(self,'loop_increment','vector assembly (in the convolution) framentation'))
    … …  
    8181                        return string
    8282                # }}}
     83
    8384        def setdefaultparameters(self): # {{{
    84                
    8585                #Convergence criterion: absolute, relative and residual
    86                 self.reltol=0.01 #default
    87                 self.abstol=float('NaN') #1 mm of sea level rise
     86                self.reltol     =       0.01 #default
     87                self.abstol     =       float('NaN') #1 mm of sea level rise
    8888
    8989                #maximum of non-linear iterations.
    90                 self.maxiter=5
    91                 self.loop_increment=200
    92 
    93                 #computational flags: 
    94                 self.geodetic=0
    95                 self.rigid=1
    96                 self.elastic=1
    97                 self.ocean_area_scaling=0
    98                 self.rotation=1
    99 
    100                 #tidal love numbers: 
    101                 self.tide_love_h=0.6149 #degree 2
    102                 self.tide_love_k=0.3055 #degree 2
    103                
    104       #secular fluid love number: 
    105                 self.fluid_love=0.942
    106                
    107                 #moment of inertia: 
    108                 self.equatorial_moi=8.0077*10**37 # [kg m^2]
    109                 self.polar_moi      =8.0345*10**37 # [kg m^2]
    110                
    111                 #mean rotational velocity of earth 
    112                 self.angular_velocity=7.2921*10**-5 # [s^-1]
     90                self.maxiter                            =       5
     91                self.loop_increment     =       200
     92
     93                #computational flags:
     94                self.geodetic                                           =       0
     95                self.rigid                                                      =       1
     96                self.elastic                                            =       1
     97                self.ocean_area_scaling =       0
     98                self.rotation                                           =       1
     99
     100                #tidal love numbers:
     101                self.tide_love_h = 0.6149 #degree 2
     102                self.tide_love_k = 0.3055 #degree 2
     103
     104      #secular fluid love number:
     105                self.fluid_love =       0.942
     106
     107                #moment of inertia:
     108                self.equatorial_moi     =       8.0077*10**37 # [kg m^2]
     109                self.polar_moi      =   8.0345*10**37 # [kg m^2]
     110
     111                #mean rotational velocity of earth
     112                self.angular_velocity   =       7.2921*10**-5 # [s^-1]
    113113
    114114                #numerical discretization accuracy
    115                 self.degacc=.01
     115                self.degacc     =       .01
    116116
    117117                #steric:
    118                 self.steric_rate=0
     118                self.steric_rate = 0
    119119
    120120                #how many time steps we skip before we run SLR solver during transient
    121                 self.geodetic_run_frequency=1
    122                
     121                self.geodetic_run_frequency     =       1
     122
    123123                #output default:
    124                 self.requested_outputs=['default']
    125 
    126                 #transitions should be a cell array of vectors: 
    127                 self.transitions=[]
     124                self.requested_outputs = ['default']
     125
     126                #transitions should be a cell array of vectors:
     127                self.transitions = []
    128128
    129129                #horizontal displacement?  (not by default)
    130                 self.horiz=0
     130                self.horiz = 0
    131131
    132132                return self
    133133                #}}}
     134
    134135        def checkconsistency(self,md,solution,analyses):    # {{{
    135 
    136136                #Early return
    137137                if (solution!='SealevelriseAnalysis'):
    … …  
    162162                md = checkfield(md,'fieldname','slr.Ugia','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    163163
    164                 #check that love numbers are provided at the same level of accuracy: 
     164                #check that love numbers are provided at the same level of accuracy:
    165165                if (size(self.love_h,0) != size(self.love_k,0) | size(self.love_h,0) != size(self.love_l,0)):
    166166                        error('slr error message: love numbers should be provided at the same level of accuracy')
    167167
    168                 #cross check that whereever we have an ice load, the mask is <0 on each vertex: 
     168                #cross check that whereever we have an ice load, the mask is <0 on each vertex:
    169169                pos=np.where(self.deltathickness)
    170                 maskpos=md.mask.ice_levelset[md.mesh.elements[pos,:]] 
     170                maskpos=md.mask.ice_levelset[md.mesh.elements[pos,:]]
    171171                els=np.where(maskpos>0)
    172172                if len(els[0])>0:
    173173                        warnings.warn('slr checkconsistency fail: there are elements with ice loads where some vertices are not on the ice!')
    174                
    175                 #check that  if geodetic is requested, we are a mesh3dsurface model (planet), or if we are not, 
    176                 #a coupler to a planet model is provided. 
     174
     175                #check that  if geodetic is requested, we are a mesh3dsurface model (planet), or if we are not,
     176                #a coupler to a planet model is provided.
    177177                if self.geodetic and not md.transient.iscoupler and domaintype(md.mesh)!='mesh3dsurface':
    178178                        error('model is requesting geodetic computations without being a mesh3dsurface, or being coupled to one!')
    179179                return md
    180180        # }}}
     181
    181182        def defaultoutputs(self,md): # {{{
    182183                return ['Sealevel']
    183184        # }}}
     185
    184186        def marshall(self,prefix,md,fid): # {{{
    185187                WriteData(fid,prefix,'object',self,'fieldname','deltathickness','format','DoubleMat','mattype',2)
    … …  
    211213                WriteData(fid,prefix,'object',self,'fieldname','horiz','format','Integer')
    212214                WriteData(fid,prefix,'object',self,'fieldname','geodetic','format','Integer')
    213        
     215
    214216                #process requested outputs
    215217                outputs = self.requested_outputs
Note: See TracChangeset for help on using the changeset viewer.