Changeset 21458

Timestamp:

01/02/17 20:17:16 (8 years ago)

Author:

adhikari

Message:

BUG: rotational parameters are introduced for python

File:

• issm/trunk-jpl/src/m/classes/slr.py (modified) (4 diffs)

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

@@ -25 +25 @@
                 self.tide_love_h       = 0
                 self.tide_love_k       = 0
+                self.fluid_love        = 0; 
+                self.equatorial_moi    = 0; 
+                self.polar_moi         = 0; 
                 self.angular_velocity  = 0;
                 self.rigid             = 0
@@ -77 +80 @@
                 self.tide_love_h=0.6149; #degree 2
                 self.tide_love_k=0.3055; #degree 2
+                
+                #secular fluid love number: 
+                self.fluid_love=0.942; 
+                
+                #moment of inertia: 
+                self.equatorial_moi=8.0077*10^37; % [kg m^2] 
+                self.polar_moi     =8.0345*10^37; % [kg m^2] 
 
                 #mean rotational velocity of earth 
@@ -107 +117 @@
                 md = checkfield(md,'fieldname','slr.tide_love_h','NaN',1,'Inf',1)
                 md = checkfield(md,'fieldname','slr.tide_love_k','NaN',1,'Inf',1)
+                md = checkfield(md,'fieldname','slr.fluid_love','NaN',1,'Inf',1)
+                md = checkfield(md,'fieldname','slr.equatorial_moi','NaN',1,'Inf',1)
+                md = checkfield(md,'fieldname','slr.polar_moi','NaN',1,'Inf',1)
                 md = checkfield(md,'fieldname','slr.angular_velocity','NaN',1,'Inf',1)
                 md = checkfield(md,'fieldname','slr.reltol','size',[1,1])
@@ -133 +146 @@
                 WriteData(fid,prefix,'object',self,'fieldname','tide_love_h','format','Double');
                 WriteData(fid,prefix,'object',self,'fieldname','tide_love_k','format','Double');
+                WriteData(fid,prefix,'object',self,'fieldname','fluid_love','format','Double');
+                WriteData(fid,prefix,'object',self,'fieldname','equatorial_moi','format','Double');
+                WriteData(fid,prefix,'object',self,'fieldname','polar_moi','format','Double');
                 WriteData(fid,prefix,'object',self,'fieldname','angular_velocity','format','Double');
                 WriteData(fid,prefix,'object',self,'fieldname','rigid','format','Boolean')