Changeset 25144


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

CHG: resolved serious conflict in naming, all in the name of love

Location:
issm/trunk-jpl/src/m/classes
Files:
1 edited
1 moved

Legend:

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

    r25143 r25144  
    1 %LOVE numbers class definition
     1%LOVENUMBERS numbers class definition
    22%
    33%   Usage:
    4 %      love=love() %will setup love numbers deg 1001 by default
    5 %      love=love('maxdeg',10001);  %supply numbers of degrees required (here, 10001)
     4%      lovenumbers=lovenumbers() %will setup love numbers deg 1001 by default
     5%      lovenumbers=lovenumbers('maxdeg',10001);  %supply numbers of degrees required (here, 10001)
    66
    7 classdef love
     7classdef lovenumbers
    88        properties (SetAccess=public)
    99               
     
    2121        end
    2222        methods
    23                 function self = love(varargin) % {{{
     23                function self = lovenumbers(varargin) % {{{
    2424                        options=pairoptions(varargin{:});
    2525                        maxdeg=getfieldvalue(options,'maxdeg',1000);
     
    2929               
    3030                        %initialize love numbers:
    31                         self.h=love_numbers('type','loadingverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
    32                         self.k=love_numbers('type','loadinggravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
    33                         self.l=love_numbers('type','loadinghorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
    34                         self.th=love_numbers('type','tidalverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
    35                         self.tk=love_numbers('type','tidalgravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
    36                         self.tl=love_numbers('type','tidalhorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
     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);
    3737
    3838                        %secular fluid love number:
     
    4747                        end
    4848
    49                         md = checkfield(md,'fieldname','solidearth.love.h','NaN',1,'Inf',1);
    50                         md = checkfield(md,'fieldname','solidearth.love.k','NaN',1,'Inf',1);
    51                         md = checkfield(md,'fieldname','solidearth.love.l','NaN',1,'Inf',1);
     49                        md = checkfield(md,'fieldname','solidearth.lovenumbers.h','NaN',1,'Inf',1);
     50                        md = checkfield(md,'fieldname','solidearth.lovenumbers.k','NaN',1,'Inf',1);
     51                        md = checkfield(md,'fieldname','solidearth.lovenumbers.l','NaN',1,'Inf',1);
    5252                       
    5353                       
    54                         md = checkfield(md,'fieldname','solidearth.love.th','NaN',1,'Inf',1);
    55                         md = checkfield(md,'fieldname','solidearth.love.tk','NaN',1,'Inf',1);
    56                         md = checkfield(md,'fieldname','solidearth.love.tk2secular','NaN',1,'Inf',1);
     54                        md = checkfield(md,'fieldname','solidearth.lovenumbers.th','NaN',1,'Inf',1);
     55                        md = checkfield(md,'fieldname','solidearth.lovenumbers.tk','NaN',1,'Inf',1);
     56                        md = checkfield(md,'fieldname','solidearth.lovenumbers.tk2secular','NaN',1,'Inf',1);
    5757
    5858                        %check that love numbers are provided at the same level of accuracy:
    5959                        if (size(self.h,1)~=size(self.k,1) | size(self.h,1)~=size(self.l,1)),
    60                                 error('love error message: love numbers should be provided at the same level of accuracy');
     60                                error('lovenumbers error message: love numbers should be provided at the same level of accuracy');
    6161                        end
    6262
     
    6767                end % }}}
    6868                function disp(self) % {{{
    69                         disp(sprintf('   love parameters:'));
     69                        disp(sprintf('   lovenumbers parameters:'));
    7070
    7171                        fielddisplay(self,'h','load Love number for radial displacement');
     
    8080                function marshall(self,prefix,md,fid) % {{{
    8181
    82                         WriteData(fid,prefix,'object',self,'fieldname','h','name','md.solidearth.love.h','format','DoubleMat','mattype',1);
    83                         WriteData(fid,prefix,'object',self,'fieldname','k','name','md.solidearth.love.k','format','DoubleMat','mattype',1);
    84                         WriteData(fid,prefix,'object',self,'fieldname','l','name','md.solidearth.love.l','format','DoubleMat','mattype',1);
     82                        WriteData(fid,prefix,'object',self,'fieldname','h','name','md.solidearth.lovenumbers.h','format','DoubleMat','mattype',1);
     83                        WriteData(fid,prefix,'object',self,'fieldname','k','name','md.solidearth.lovenumbers.k','format','DoubleMat','mattype',1);
     84                        WriteData(fid,prefix,'object',self,'fieldname','l','name','md.solidearth.lovenumbers.l','format','DoubleMat','mattype',1);
    8585
    86                         WriteData(fid,prefix,'object',self,'fieldname','th','name','md.solidearth.love.th','format','DoubleMat','mattype',1);
    87                         WriteData(fid,prefix,'object',self,'fieldname','tk','name','md.solidearth.love.tk','format','DoubleMat','mattype',1);
    88                         WriteData(fid,prefix,'object',self,'fieldname','tl','name','md.solidearth.love.tl','format','DoubleMat','mattype',1);
    89                         WriteData(fid,prefix,'object',self,'data',self.tk2secular,'fieldname','love.tk2secular','format','Double');
     86                        WriteData(fid,prefix,'object',self,'fieldname','th','name','md.solidearth.lovenumbers.th','format','DoubleMat','mattype',1);
     87                        WriteData(fid,prefix,'object',self,'fieldname','tk','name','md.solidearth.lovenumbers.tk','format','DoubleMat','mattype',1);
     88                        WriteData(fid,prefix,'object',self,'fieldname','tl','name','md.solidearth.lovenumbers.tl','format','DoubleMat','mattype',1);
     89                        WriteData(fid,prefix,'object',self,'data',self.tk2secular,'fieldname','lovenumbers.tk2secular','format','Double');
    9090
    9191                end % }}}
    9292                function savemodeljs(self,fid,modelname) % {{{
    93                         writejs1Darray(fid,[modelname '.love.h'],self.h);
    94                         writejs1Darray(fid,[modelname '.love.k'],self.k);
    95                         writejs1Darray(fid,[modelname '.love.l'],self.l);
     93                        writejs1Darray(fid,[modelname '.lovenumbers.h'],self.h);
     94                        writejs1Darray(fid,[modelname '.lovenumbers.k'],self.k);
     95                        writejs1Darray(fid,[modelname '.lovenumbers.l'],self.l);
    9696                end % }}}
    9797                function self = extrude(self,md) % {{{
  • TabularUnified issm/trunk-jpl/src/m/classes/solidearth.m

    r25138 r25144  
    99                settings               = solidearthsettings();
    1010                surfaceload            = surfaceload();
    11                 love                   = love();
     11                lovenumbers            = lovenumbers();
    1212                rotational             = rotational();
    1313                planetradius           = planetradius('earth');
     
    4747                        self.settings.checkconsistency(md,solution,analyses);
    4848                        self.surfaceload.checkconsistency(md,solution,analyses);
    49                         self.love.checkconsistency(md,solution,analyses);
     49                        self.lovenumbers.checkconsistency(md,solution,analyses);
    5050                        self.rotational.checkconsistency(md,solution,analyses);
    5151
     
    6464                        self.settings.disp();
    6565                        self.surfaceload.disp();
    66                         self.love.disp();
     66                        self.lovenumbers.disp();
    6767                        self.rotational.disp();
    6868
     
    7676                        self.settings.marshall(prefix,md,fid);
    7777                        self.surfaceload.marshall(prefix,md,fid);
    78                         self.love.marshall(prefix,md,fid);
     78                        self.lovenumbers.marshall(prefix,md,fid);
    7979                        self.rotational.marshall(prefix,md,fid);
    8080                       
     
    9494                        self.settings.savemodeljs(fid,modelname);
    9595                        self.surfaceload.savemodeljs(fid,modelname);
    96                         self.love.savemodeljs(fid,modelname);
     96                        self.lovenumbers.savemodeljs(fid,modelname);
    9797                        self.rotational.savemodeljs(fid,modelname);
    9898                        writejscellstring(fid,[modelname '.solidearth.requested_outputs'],self.requested_outputs);
Note: See TracChangeset for help on using the changeset viewer.