Changeset 25166 for issm/trunk-jpl/src/m/boundaryconditions

Timestamp:

06/26/20 20:45:42 (5 years ago)

Author:

jdquinn

Message:

CHG: Saving MATLAB -> Python changes; still hunting for implementation bug

Location:

issm/trunk-jpl/src/m/boundaryconditions

Files:

• getlovenumbers.m (modified) (1 diff)
• getlovenumbers.py (modified) (1 diff)

issm/trunk-jpl/src/m/boundaryconditions/getlovenumbers.m

@@ -10062 +10062 @@
         if strcmpi(frame,'CM'), 
                 return; 
-        elseif strcmpi(frame,'CF'), % from Blewitt, 2003, JGR 
+        elseif strcmpi(frame,'CF'), % from Blewitt, 2003, JGR
                 if strcmpi(type,'loadingverticaldisplacement'),
                         series(2,1) = -0.269; 

issm/trunk-jpl/src/m/boundaryconditions/getlovenumbers.py

@@ -10069 +10069 @@
     elif frame == 'CF': #from Blewitt, 2003, JGR
         if type == 'loadingverticaldisplacement':
-            series[1, 0] = -0.269
+            series[0] = -0.269
         elif type == 'loadinggravitationalpotential':
-            series[1, 0] = 0.021
+            series[0] = 0.021
         elif type == 'loadinghorizontaldisplacement':
-            series[1, 0] = 0.134
+            series[0] = 0.134
     else:
         raise Exception("love_numbers error message: unknown reference frame: {}".format(frame))