Changeset 19988


Ignore:
Timestamp:
01/22/16 16:52:38 (9 years ago)
Author:
Eric.Larour
Message:

CHG: missing earth density

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

Legend:

Unmodified
Added
Removed

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

    r19958 r19988  
    2727                mantle_shear_modulus       = 0.;
    2828                mantle_density             = 0.;
     29
     30                %slr
     31                earth_density              = 0;
    2932
    3033        end
    … …  
    100103                        self.mantle_density             = 3.34;       % (g/cm^-3)
    101104
     105                        %SLR
     106                        self.earth_density= 5512;  % average density of the Earth, (kg/m^3)
     107
    102108                end % }}}
    103109                function md = checkconsistency(self,md,solution,analyses) % {{{
    … …  
    116122                                md = checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',1);
    117123                        end
     124                        if ismember(SealevelriseAnalysisEnum(),analyses),
     125                                md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',1);
     126                        end
     127
    118128
    119129                end % }}}
    … …  
    140150                        fielddisplay(self,'mantle_shear_modulus','Mantle shear modulus [Pa]');
    141151                        fielddisplay(self,'mantle_density','Mantle density [g/cm^-3]');
     152                        fielddisplay(self,'earth_density','Mantle density [kg/m^-3]');
    142153                end % }}}
    143154                function marshall(self,md,fid) % {{{
    … …  
    163174                        WriteData(fid,'object',self,'class','materials','fieldname','mantle_shear_modulus','format','Double');
    164175                        WriteData(fid,'object',self,'class','materials','fieldname','mantle_density','format','Double','scale',10^3);
     176                        WriteData(fid,'object',self,'class','materials','fieldname','earth_density','format','Double');
     177
    165178                end % }}}
    166179        end

  • issm/trunk-jpl/src/m/classes/matdamageice.py

    r19048 r19988  
    3636                self.mantle_shear_modulus       = 0.
    3737                self.mantle_density             = 0.
     38               
     39                #SLR
     40                self.earth_density= 5512;  # average density of the Earth, (kg/m^3)
     41
    3842
    3943                self.setdefaultparameters()
    … …  
    6165                string="%s\n%s"%(string,fielddisplay(self,"mantle_shear_modulus","Mantle shear modulus [Pa]"))
    6266                string="%s\n%s"%(string,fielddisplay(self,"mantle_density","Mantle density [g/cm^-3]"))
     67                string="%s\n%s"%(string,fielddisplay(self,"earth_density","Mantle density [kg/m^-3]"))
     68
    6369
    6470                return string
    … …  
    115121                self.mantle_shear_modulus       = 1.45*10**11 # (Pa)
    116122                self.mantle_density             = 3.34        # (g/cm^-3)
     123               
     124                #SLR
     125                self.earth_density= 5512;  #average density of the Earth, (kg/m^3)
     126
    117127
    118128                return self
    … …  
    130140                md = checkfield(md,'fieldname','materials.mantle_shear_modulus','>',0,'numel',[1]);
    131141                md = checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',[1]);
     142                md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',[1]);
     143
    132144                return md
    133145        # }}}
    … …  
    154166                WriteData(fid,'object',self,'class','materials','fieldname','mantle_shear_modulus','format','Double');
    155167                WriteData(fid,'object',self,'class','materials','fieldname','mantle_density','format','Double','scale',10.**3.);
     168                WriteData(fid,'object',self,'class','materials','fieldname','earth_density','format','Double');
     169
    156170        # }}}

  • issm/trunk-jpl/src/m/classes/matice.js

    r19901 r19988  
    5353                        this.mantle_shear_modulus       = 1.45*Math.pow(10,11); // (Pa)
    5454                        this.mantle_density             = 3.34;       // (g/cm^-3)
     55                       
     56                        //SLR
     57                        this.earth_density= 5512;  // average density of the Earth, (kg/m^3)
     58
    5559
    5660                } //}}}
    … …  
    7781                        fielddisplay(this,'mantle_shear_modulus','Mantle shear modulus [Pa]');
    7882                        fielddisplay(this,'mantle_density','Mantle density [g/cm^-3]');
     83                        fielddisplay(this,'earth_density','Mantle density [kg/m^-3]');
     84
    7985                } // }}}
    8086                this.classname = function() {// {{{
    … …  
    96102                                checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',1);
    97103                        }
     104                        if ArrayAnyEqual(ArrayIsMember(SealevelriseAnalysisEnum(),analyses),1)){
     105                                checkfield(md,'fieldname','materials.earth_density','>',0,'numel',1);
     106                        }
     107
    98108
    99109                } // }}}
    … …  
    119129                        WriteData(fid,'object',this,'class','materials','fieldname','mantle_shear_modulus','format','Double');
    120130                        WriteData(fid,'object',this,'class','materials','fieldname','mantle_density','format','Double','scale',10^3);
     131                        WriteData(fid,'object',this,'class','materials','fieldname','earth_density','format','Double');
     132
    121133                }//}}}
    122134                this.fix=function() { //{{{
    … …  
    145157                this.mantle_shear_modulus       = 0.;
    146158                this.mantle_density             = 0.;
     159
     160                //SLR
     161                this.earth_density= 5512;  // average density of the Earth, (kg/m^3)
     162
    147163                this.setdefaultparameters();
    148164                //}}}

  • issm/trunk-jpl/src/m/classes/matice.py

    r19897 r19988  
    3535                self.lithosphere_density        = 0.
    3636                self.mantle_shear_modulus       = 0.
    37                 self.mantle_density             = 0.
     37                self.mantle_density             = 0. 
     38               
     39                #SLR
     40                self.earth_density= 5512; 
     41
     42
    3843
    3944                self.setdefaultparameters()
    … …  
    6166                string="%s\n%s"%(string,fielddisplay(self,"mantle_shear_modulus","Mantle shear modulus [Pa]"))
    6267                string="%s\n%s"%(string,fielddisplay(self,"mantle_density","Mantle density [g/cm^-3]"))
     68                string="%s\n%s"%(string,fielddisplay(self,"earth_density","Mantle density [kg/m^-3]"))
     69
    6370
    6471                return string
    … …  
    115122                self.mantle_shear_modulus       = 1.45*10**11 # (Pa)
    116123                self.mantle_density             = 3.34        # (g/cm^-3)
     124               
     125                #SLR
     126                self.earth_density= 5512;  # average density of the Earth, (kg/m^3)
     127
    117128
    118129                return self
    … …  
    130141                md = checkfield(md,'fieldname','materials.mantle_shear_modulus','>',0,'numel',[1]);
    131142                md = checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',[1]);
     143                md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',[1]);
     144
    132145                return md
    133146        # }}}
    … …  
    155168                WriteData(fid,'object',self,'class','materials','fieldname','mantle_shear_modulus','format','Double');
    156169                WriteData(fid,'object',self,'class','materials','fieldname','mantle_density','format','Double','scale',10.**3.);
     170                WriteData(fid,'object',self,'class','materials','fieldname','earth_density','format','Double');
     171
    157172        # }}}
Note: See TracChangeset for help on using the changeset viewer.