Changeset 22970


Ignore:
Timestamp:
07/17/18 14:44:49 (7 years ago)
Author:
Mathieu Morlighem
Message:

CHG: fixing solid earth

Location:
issm/trunk-jpl
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/classes/slr.m

    r22955 r22970  
    66classdef slr
    77        properties (SetAccess=public)
    8                 deltathickness = NaN;
    9                 sealevel       = NaN;
    10                 spcthickness   = NaN;
    11                 maxiter        = 0;
    12                 reltol         = 0;
    13                 abstol         = 0;
    14                 love_h         = 0; %provided by PREM model
    15                 love_k         = 0; %ideam
    16                 love_l         = 0; %ideam
    17                 tide_love_k    = 0; %ideam
    18                 tide_love_h    = 0; %ideam
    19                 fluid_love     = 0;
    20                 equatorial_moi = 0;
    21                 polar_moi               = 0;
    22                 angular_velocity = 0;
    23                 rigid          = 0;
    24                 elastic        = 0;
    25                 rotation       = 0;
    26                 ocean_area_scaling = 0;
    27                 steric_rate    = 0; %rate of ocean expansion from steric effects.
    28                 geodetic_run_frequency  = 1; %how many time steps we skip before we run the geodetic part of the solver during transient
    29                 geodetic       = 0; %compute geodetic SLR? (in addition to steric?)
    30                 degacc         = 0;
    31                 loop_increment = 0;
    32                 horiz = 0;
    33                 Ngia  = 0;
    34                 Ugia  = 0;
     8                deltathickness         = NaN;
     9                sealevel               = NaN;
     10                spcthickness           = NaN;
     11                maxiter                = 0;
     12                reltol                 = 0;
     13                abstol                 = 0;
     14                love_h                 = 0; %provided by PREM model
     15                love_k                 = 0; %ideam
     16                love_l                 = 0; %ideam
     17                tide_love_k            = 0; %ideam
     18                tide_love_h            = 0; %ideam
     19                fluid_love             = 0;
     20                equatorial_moi         = 0;
     21                polar_moi              = 0;
     22                angular_velocity       = 0;
     23                rigid                  = 0;
     24                elastic                = 0;
     25                rotation               = 0;
     26                ocean_area_scaling     = 0;
     27                steric_rate            = 0; %rate of ocean expansion from steric effects.
     28                geodetic_run_frequency = 1; %how many time steps we skip before we run the geodetic part of the solver during transient
     29                geodetic               = 0; %compute geodetic SLR? (in addition to steric?)
     30                degacc                 = 0;
     31                loop_increment         = 0;
     32                horiz                  = 0;
     33                Ngia                   = NaN;
     34                Ugia                   = NaN;
    3535                requested_outputs      = {};
    36                 transitions    = {};
     36                transitions            = {};
    3737        end
    3838        methods

  • issm/trunk-jpl/test/NightlyRun/test2002.m

    r22109 r22970  
    5959md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
    6060
     61%New stuff
     62md.slr.spcthickness = NaN(md.mesh.numberofvertices,1);
     63md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
     64md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
     65
    6166%Miscellaneous
    6267md.miscellaneous.name='test2002';

  • issm/trunk-jpl/test/NightlyRun/test2003.m

    r21753 r22970  
    6060md.miscellaneous.name='test2003';
    6161
     62%New stuff
     63md.slr.spcthickness = NaN(md.mesh.numberofvertices,1);
     64md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
     65md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
     66
    6267%Solution parameters
    6368md.slr.reltol=NaN;

  • issm/trunk-jpl/test/NightlyRun/test2010.m

    r21947 r22970  
    6666% }}}
    6767
     68%New stuff
     69md.slr.spcthickness = NaN(md.mesh.numberofvertices,1);
     70md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
     71md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
     72
    6873%eustatic + rigid + elastic run:
    6974md.slr.rigid=1; md.slr.elastic=1; md.slr.rotation=1;

  • issm/trunk-jpl/test/NightlyRun/test2111.m

    r22350 r22970  
    4747% }}}
    4848
     49%New stuff
     50md.slr.spcthickness = NaN(md.mesh.numberofvertices,1);
     51md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
     52md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
     53
    4954%solve esa: {{{
    5055md.esa.requested_outputs = {'EsaUmotion','EsaNmotion','EsaEmotion','EsaXmotion','EsaYmotion'};

  • issm/trunk-jpl/test/NightlyRun/test2112.m

    r22350 r22970  
    4747% }}}
    4848
     49%New stuff
     50md.slr.spcthickness = NaN(md.mesh.numberofvertices,1);
     51md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
     52md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
     53
    4954%solve esa: {{{
    5055md.esa.requested_outputs = {'EsaUmotion','EsaNmotion','EsaEmotion','EsaXmotion','EsaYmotion'};

  • issm/trunk-jpl/test/NightlyRun/test2113.m

    r22350 r22970  
    4848% }}}
    4949
     50%New stuff
     51md.slr.spcthickness = NaN(md.mesh.numberofvertices,1);
     52md.slr.Ngia = zeros(md.mesh.numberofvertices,1);
     53md.slr.Ugia = zeros(md.mesh.numberofvertices,1);
     54
    5055%solve esa: {{{
    5156md.esa.requested_outputs = {'EsaUmotion','EsaNmotion','EsaEmotion'};
Note: See TracChangeset for help on using the changeset viewer.