Changeset 22353

Timestamp:

01/14/18 21:56:59 (7 years ago)

Author:

adhikari

Message:

CHG: center of frame reference frame implemented

File:

• issm/trunk-jpl/src/m/boundaryconditions/love_numbers.py (modified) (3 diffs)

issm/trunk-jpl/src/m/boundaryconditions/love_numbers.py ¶

@@ -7 +7 @@
 #                          retrieved from: http://www.srosat.com/iag-jsg/loveNb.php
 #    Usage:   series=love_numbers(value) 
+#             series=love_numbers(value,reference_frame) 
+# 
 #             where value is one of 'h','k','l','gamma' and 'lambda'. 
+#             reference_frame = one of 'CM' (default) and 'CF'.
 #
 #    Example:  
 #          love_k=love_numbers('k');
+#          love_k=love_numbers('k','CF');
 # 
 
-        #some checks:
-        if len(varargin)>0:
-                raise RuntimeError('love_numbers error message: wrong usage')
-
-        if value not in ['h','k','l','gamma','lambda']:
-                raise RuntimeError('value should be one of ''h'',''k'',''l'',''gamma'' and ''lambda''')
-
-        love_numbers=np.array([[    0         , 0          ,0          ,0          ,0          ,0          ,0          ],
+        # some checks:
+        if varargin==1:
+            frame='CM'; 
+        elif varargin==2: 
+            frame=reference_frame;
+        else:
+            raise RuntimeError('love_numbers error message: bad usage') 
+        
+        if value not in ['h','k','l','gamma','lambda']:
+            raise RuntimeError('value should be one of ''h'',''k'',''l'',''gamma'' and ''lambda''') 
+        
+        love_numbers=np.array([[    0         , 0          ,0          ,0          ,0          ,0          ,0          ],
                                                                                                  [      -1.28740059,-1.00000000,-0.89858519,1.28740059, 0.42519882  ,0.89858519 ,0.00000000 ],
                                                                                                  [      -1.00025365, -0.30922675, 0.02060926, 1.69102690, 0.46358648, 0.67016399, 0.61829668],
@@ -10021 +10029 @@
                                                                                                  [      -6.27342771, -0.00030948, 0.00018957, 7.27311822, 0.99905470, 0.99950094, 0.49327205],
                                                                                                  [      -6.27342778, -0.00030945, 0.00018956, 7.27311833, 0.99905480, 0.99950099, 0.49327194]]);
-
-
-        if value=='h':
+                                                                                                 
+        if value=='h':
                 series=love_numbers[:,0];
         elif value=='k':
@@ -10035 +10042 @@
         else:
                 raise RuntimeError(['love_numbers error message: unknow value:', value])
-                
-        return series
+        
+        % choose degree 1 term for CF reference system 
+        if frame=='CF': # from Blewitt, 2003, JGR 
+            if value=='h':
+                    series[2,1] = -0.269; 
+            elif value=='k':
+                    series[2,1] = 0.021;  
+            elif value=='l':
+                    series[2,1] = 0.134; 
+        
+        return series
+