Changeset 25149


Ignore:
Timestamp:
06/24/20 19:53:41 (5 years ago)
Author:
Eric.Larour
Message:

CHG: transfer reference frame better.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/classes/lovenumbers.m

    r25144 r25149  
    2424                        options=pairoptions(varargin{:});
    2525                        maxdeg=getfieldvalue(options,'maxdeg',1000);
    26                         self=setdefaultparameters(self,maxdeg);
     26                        referenceframe=getfieldvalue(options,'referenceframe','CM');
     27                        self=setdefaultparameters(self,maxdeg,referenceframe);
    2728                end % }}}
    28                 function self = setdefaultparameters(self,maxdeg) % {{{
     29                function self = setdefaultparameters(self,maxdeg,referenceframe) % {{{
    2930               
    3031                        %initialize love numbers:
    31                         self.h=getlovenumbers('type','loadingverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
    32                         self.k=getlovenumbers('type','loadinggravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
    33                         self.l=getlovenumbers('type','loadinghorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
    34                         self.th=getlovenumbers('type','tidalverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
    35                         self.tk=getlovenumbers('type','tidalgravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
    36                         self.tl=getlovenumbers('type','tidalhorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
     32                        self.h=getlovenumbers('type','loadingverticaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
     33                        self.k=getlovenumbers('type','loadinggravitationalpotential','referenceframe',referenceframe,'maxdeg',maxdeg);
     34                        self.l=getlovenumbers('type','loadinghorizontaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
     35                        self.th=getlovenumbers('type','tidalverticaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
     36                        self.tk=getlovenumbers('type','tidalgravitationalpotential','referenceframe',referenceframe,'maxdeg',maxdeg);
     37                        self.tl=getlovenumbers('type','tidalhorizontaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
    3738
    3839                        %secular fluid love number:
Note: See TracChangeset for help on using the changeset viewer.