Ignore:
Timestamp:
07/27/18 11:58:17 (7 years ago)
Author:
seroussi
Message:

NEW: missing features in slr.py

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r22955 r23019  
    1515       
    1616        def __init__(self): # {{{
    17                 self.deltathickness    = float('NaN')
    18                 self.sealevel          = float('NaN')
    19                 self.maxiter           = 0
    20                 self.reltol            = 0
    21                 self.abstol            = 0
    22                 self.love_h            = 0 #provided by PREM model()
    23                 self.love_k            = 0 #ideam
    24                 self.love_l            = 0 #ideam
    25                 self.tide_love_h       = 0
    26                 self.tide_love_k       = 0
    27                 self.fluid_love        = 0;
    28                 self.equatorial_moi    = 0;
    29                 self.polar_moi         = 0;
    30                 self.angular_velocity  = 0;
    31                 self.rigid             = 0
    32                 self.elastic           = 0
    33                 self.rotation          = 0
    34                 self.ocean_area_scaling = 0;
    35                 steric_rate                                     = 0; #rate of ocean expansion from steric effects.
    36                 self.degacc            = 0
    37                 self.requested_outputs = []
    38                 self.transitions       = []
     17                self.deltathickness         = float('NaN')
     18                self.sealevel               = float('NaN')
     19                self.maxiter                = 0
     20                self.reltol                 = 0
     21                self.abstol                 = 0
     22                self.love_h                 = 0 #provided by PREM model()
     23                self.love_k                 = 0 #ideam
     24                self.love_l                 = 0 #ideam
     25                self.tide_love_h            = 0
     26                self.tide_love_k            = 0
     27                self.fluid_love             = 0;
     28                self.equatorial_moi         = 0;
     29                self.polar_moi               = 0;
     30                self.angular_velocity       = 0;
     31                self.rigid                  = 0
     32                self.elastic                = 0
     33                self.rotation               = 0
     34                self.ocean_area_scaling     = 0;
     35                steric_rate                                         = 0; #rate of ocean expansion from steric effects.
     36                self.geodetic_run_frequency = 1; #how many time steps we skip before we run the geodetic part of the solver during transient
     37                self.geodetic               = 0; #compute geodetic SLR? (in addition to steric?)
     38                self.degacc                 = 0
     39                self.requested_outputs      = []
     40                self.transitions            = []
    3941               
    4042                #set defaults
    … …  
    7880
    7981                #computational flags:
     82                self.geodetic=0
    8083                self.rigid=1
    8184                self.elastic=1
    … …  
    8790                self.tide_love_k=0.3055; #degree 2
    8891               
    89                 #secular fluid love number:
     92      #secular fluid love number:
    9093                self.fluid_love=0.942;
    9194               
    9295                #moment of inertia:
    9396                self.equatorial_moi=8.0077*10**37; # [kg m^2]
    94                 self.polar_moi     =8.0345*10**37; # [kg m^2]
     97                self.polar_moi      =8.0345*10**37; # [kg m^2]
    9598               
    9699                #mean rotational velocity of earth
    … …  
    102105                #steric:
    103106                self.steric_rate=0;
     107
     108                #how many time steps we skip before we run SLR solver during transient
     109                self.geodetic_run_frequency=1;
    104110               
    105111                #output default:
    … …  
    133139                md = checkfield(md,'fieldname','slr.abstol','size',[1,1])
    134140                md = checkfield(md,'fieldname','slr.maxiter','size',[1,1],'>=',1)
     141                md = checkfield(md,'fieldname','slr.geodetic_run_frequency','size',[1,1],'>=',1)
    135142                md = checkfield(md,'fieldname','slr.steric_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices]);
    136143                md = checkfield(md,'fieldname','slr.degacc','size',[1,1],'>=',1e-10)
    … …  
    164171                WriteData(fid,prefix,'object',self,'fieldname','rotation','format','Boolean')
    165172                WriteData(fid,prefix,'object',self,'fieldname','ocean_area_scaling','format','Boolean')
     173                WriteData(fid,prefix,'object',self,'fieldname','geodetic_run_frequency','format','Integer')
    166174                WriteData(fid,prefix,'object',self,'fieldname','steric_rate','format','DoubleMat','mattype',1,'scale',1e-3/md.constants.yts);
    167175                WriteData(fid,prefix,'object',self,'fieldname','degacc','format','Double')
    168176                WriteData(fid,prefix,'object',self,'fieldname','transitions','format','MatArray')
     177                WriteData(fid,prefix,'object',self,'fieldname','geodetic','format','Integer')
    169178       
    170179                #process requested outputs
Note: See TracChangeset for help on using the changeset viewer.