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

CHG: materials hydro class was missing fresh water density for hydrological applications.
Also added hydrological fingerprints capabilities for sealevelchange core.

File:
1 edited

Legend:

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

    r24756 r25150  
    5858                                        self.addprop('rho_ice');
    5959                                        self.addprop('rho_water');
     60                                        self.addprop('rho_freshwater');
    6061                                        self.addprop('earth_density');
    6162                                otherwise
     
    138139                                        % average density of the Earth (kg/m^3)
    139140                                        self.earth_density=5512;
     141
     142                                        %fresh water density (kg/m^3)
     143                                        self.rho_freshwater=1000.;
    140144
    141145                                otherwise
     
    184188                                        fielddisplay(self,'rho_water','ocean water density [kg/m^3]');
    185189                                        fielddisplay(self,'earth_density','mantle density [kg/m^3]');
     190                                        fielddisplay(self,'rho_freshwater','fresh water density [kg/m^3]');
    186191
    187192                                otherwise
     
    233238                                        md = checkfield(md,'fieldname','materials.rho_water','>',0);
    234239                                        md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',1);
     240                                        md = checkfield(md,'fieldname','materials.rho_freshwater','>',0);
    235241
    236242                                otherwise
     
    280286                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','rho_water','format','Double');
    281287                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','earth_density','format','Double');
     288                                        WriteData(fid,prefix,'object',self,'class','materials','fieldname','rho_freshwater','format','Double');
    282289
    283290                                otherwise
     
    334341                                        writejsdouble(fid,[modelname '.materials.rho_water'],self.rho_water);
    335342                                        writejsdouble(fid,[modelname '.materials.earth_density'],self.earth_density);
     343                                        writejsdouble(fid,[modelname '.materials.rho_freshwater'],self.rho_freshwater);
    336344
    337345                                otherwise
Note: See TracChangeset for help on using the changeset viewer.