Changeset 25149

Timestamp:

06/24/20 19:53:41 (5 years ago)

Author:

Eric.Larour

Message:

CHG: transfer reference frame better.

File:

• issm/trunk-jpl/src/m/classes/lovenumbers.m (modified) (1 diff)

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

@@ -24 +24 @@
                         options=pairoptions(varargin{:});
                         maxdeg=getfieldvalue(options,'maxdeg',1000);
-                        self=setdefaultparameters(self,maxdeg);
+                        referenceframe=getfieldvalue(options,'referenceframe','CM');
+                        self=setdefaultparameters(self,maxdeg,referenceframe);
                 end % }}}
-                function self = setdefaultparameters(self,maxdeg) % {{{
+                function self = setdefaultparameters(self,maxdeg,referenceframe) % {{{
                 
                         %initialize love numbers:
-                        self.h=getlovenumbers('type','loadingverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
-                        self.k=getlovenumbers('type','loadinggravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
-                        self.l=getlovenumbers('type','loadinghorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
-                        self.th=getlovenumbers('type','tidalverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
-                        self.tk=getlovenumbers('type','tidalgravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
-                        self.tl=getlovenumbers('type','tidalhorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
+                        self.h=getlovenumbers('type','loadingverticaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
+                        self.k=getlovenumbers('type','loadinggravitationalpotential','referenceframe',referenceframe,'maxdeg',maxdeg);
+                        self.l=getlovenumbers('type','loadinghorizontaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
+                        self.th=getlovenumbers('type','tidalverticaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
+                        self.tk=getlovenumbers('type','tidalgravitationalpotential','referenceframe',referenceframe,'maxdeg',maxdeg);
+                        self.tl=getlovenumbers('type','tidalhorizontaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
 
                         %secular fluid love number: