Changeset 21345

Timestamp:

11/04/16 17:43:29 (8 years ago)

Author:

adhikari

Message:

CHG: rotational feedback matlab parameterization

File:

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

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

@@ -16 +16 @@
                 tide_love_k    = 0; %ideam
                 tide_love_h    = 0; %ideam
+                fluid_love     = 0; 
+                equatorial_moi = 0; 
+                polar_moi               = 0; 
                 angular_velocity = 0;
                 rigid          = 0;
@@ -52 +55 @@
                 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 
                 self.angular_velocity=7.2921*10^-5; % [s^-1] 
@@ -76 +86 @@
                         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]);
@@ -113 +126 @@
                         fielddisplay(self,'tide_love_k','tidal load Love number (deg 2)');
                         fielddisplay(self,'tide_love_h','tidal load Love number (deg 2)');
+                        fielddisplay(self,'fluid_love','secular fluid Love number');
+                        fielddisplay(self,'equatorial_moi','mean equatorial moment of inertia [kg m^2]');
+                        fielddisplay(self,'polar_moi','polar moment of inertia [kg m^2]');
                         fielddisplay(self,'angular_velocity','mean rotational velocity of earth [per second]'); 
                         fielddisplay(self,'rigid','rigid earth graviational potential perturbation');
@@ -134 +150 @@
                         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');
@@ -164 +183 @@
                         writejsdouble(fid,[modelname '.slr.tide_love_k'],self.tide_love_k);
                         writejsdouble(fid,[modelname '.slr.tide_love_h'],self.tide_love_h);
+                        writejsdouble(fid,[modelname '.slr.fluid_love'],self.fluid_love);
+                        writejsdouble(fid,[modelname '.slr.equatorial_moi'],self.equatorial_moi);
+                        writejsdouble(fid,[modelname '.slr.polar_moi'],self.polar_moi);
                         writejsdouble(fid,[modelname '.slr.angular_velocity'],self.angular_velocity);
                         writejsdouble(fid,[modelname '.slr.rigid'],self.rigid);