Changeset 21460

Timestamp:

01/02/17 21:41:30 (8 years ago)

Author:

adhikari

Message:

CHG: rotational parameters introduced

File:

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

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

@@ -25 +25 @@
                 this.tide_love_k=0.3055; //degree 2
         
+                //secular fluid love number: 
+                this.fluid_love=0.942; 
+                
+                //moment of inertia: 
+                this.equatorial_moi=8.0077*10^37; // [kg m^2] 
+                this.polar_moi           =8.0345*10^37; // [kg m^2] 
+
                 // mean rotational velocity of earth:
                 this.angular_velocity=7.2921*10^-5; // [s^-1] 
@@ -50 +57 @@
                         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]);
@@ -83 +93 @@
                 fielddisplay(this,'tide_love_h','tidal love number (degree 2)');
                 fielddisplay(this,'tide_love_k','tidal love number (degree 2)');
+                fielddisplay(this,'fluid_love','secular fluid Love number');
+                fielddisplay(this,'equatorial_moi','mean equatorial moment of inertia [kg m^2]');
+                fielddisplay(this,'polar_moi','polar moment of inertia [kg m^2]');
                 fielddisplay(this,'angular_velocity','mean rotational velocity of earth [per second]'); 
                 fielddisplay(this,'rigid','rigid earth graviational potential perturbation');
@@ -104 +117 @@
                         WriteData(fid,prefix,'object',this,'fieldname','tide_love_h','format','Double');
                         WriteData(fid,prefix,'object',this,'fieldname','tide_love_k','format','Double');
+                        WriteData(fid,prefix,'object',this,'fieldname','fluid_love','format','Double');
+                        WriteData(fid,prefix,'object',this,'fieldname','equatorial_moi','format','Double');
+                        WriteData(fid,prefix,'object',this,'fieldname','polar_moi','format','Double');
                         WriteData(fid,prefix,'object',this,'fieldname','angular_velocity','format','Double');
                         WriteData(fid,prefix,'object',this,'fieldname','rigid','format','Boolean');
@@ -134 +150 @@
                         this.tide_love_h=NullFix(this.tide_love_h,NaN);
                         this.tide_love_k=NullFix(this.tide_love_k,NaN);
+                        this.fluid_love=NullFix(this.fluid_love,NaN);
+                        this.equatorial_moi=NullFix(this.equatorial_moi,NaN);
+                        this.polar_moi=NullFix(this.polar_moi,NaN);
                         this.angular_velocity=NullFix(this.angular_velocity,NaN);
                         this.rigid=NullFix(this.rigid,NaN);
@@ -153 +172 @@
         this.tide_love_h    = 0; 
         this.tide_love_k    = 0; 
+        this.fluid_love   = 0; 
+        this.equatorial_moi     = 0; 
+        this.polar_moi                  = 0; 
         this.angular_velocity = 0; 
         this.rigid          = 0;