Index: /issm/trunk-jpl/src/m/classes/matdamageice.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/matdamageice.m	(revision 19987)
+++ /issm/trunk-jpl/src/m/classes/matdamageice.m	(revision 19988)
@@ -27,4 +27,7 @@
 		mantle_shear_modulus       = 0.;
 		mantle_density             = 0.;
+
+		%slr
+		earth_density              = 0;
 
 	end
@@ -100,4 +103,7 @@
 			self.mantle_density             = 3.34;       % (g/cm^-3)
 
+			%SLR
+			self.earth_density= 5512;  % average density of the Earth, (kg/m^3)
+
 		end % }}}
 		function md = checkconsistency(self,md,solution,analyses) % {{{
@@ -116,4 +122,8 @@
 				md = checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',1);
 			end
+			if ismember(SealevelriseAnalysisEnum(),analyses),
+				md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',1);
+			end
+
 
 		end % }}}
@@ -140,4 +150,5 @@
 			fielddisplay(self,'mantle_shear_modulus','Mantle shear modulus [Pa]');
 			fielddisplay(self,'mantle_density','Mantle density [g/cm^-3]');
+			fielddisplay(self,'earth_density','Mantle density [kg/m^-3]');
 		end % }}}
 		function marshall(self,md,fid) % {{{
@@ -163,4 +174,6 @@
 			WriteData(fid,'object',self,'class','materials','fieldname','mantle_shear_modulus','format','Double');
 			WriteData(fid,'object',self,'class','materials','fieldname','mantle_density','format','Double','scale',10^3);
+			WriteData(fid,'object',self,'class','materials','fieldname','earth_density','format','Double');
+
 		end % }}}
 	end
Index: /issm/trunk-jpl/src/m/classes/matdamageice.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/matdamageice.py	(revision 19987)
+++ /issm/trunk-jpl/src/m/classes/matdamageice.py	(revision 19988)
@@ -36,4 +36,8 @@
 		self.mantle_shear_modulus       = 0.
 		self.mantle_density             = 0.
+		
+		#SLR
+		self.earth_density= 5512;  # average density of the Earth, (kg/m^3)
+
 
 		self.setdefaultparameters()
@@ -61,4 +65,6 @@
 		string="%s\n%s"%(string,fielddisplay(self,"mantle_shear_modulus","Mantle shear modulus [Pa]"))
 		string="%s\n%s"%(string,fielddisplay(self,"mantle_density","Mantle density [g/cm^-3]"))
+		string="%s\n%s"%(string,fielddisplay(self,"earth_density","Mantle density [kg/m^-3]"))
+
 
 		return string
@@ -115,4 +121,8 @@
 		self.mantle_shear_modulus       = 1.45*10**11 # (Pa)
 		self.mantle_density             = 3.34        # (g/cm^-3)
+		
+		#SLR
+		self.earth_density= 5512;  #average density of the Earth, (kg/m^3)
+
 
 		return self
@@ -130,4 +140,6 @@
 		md = checkfield(md,'fieldname','materials.mantle_shear_modulus','>',0,'numel',[1]);
 		md = checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',[1]);
+		md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',[1]);
+
 		return md
 	# }}}
@@ -154,3 +166,5 @@
 		WriteData(fid,'object',self,'class','materials','fieldname','mantle_shear_modulus','format','Double');
 		WriteData(fid,'object',self,'class','materials','fieldname','mantle_density','format','Double','scale',10.**3.);
+		WriteData(fid,'object',self,'class','materials','fieldname','earth_density','format','Double');
+
 	# }}}
Index: /issm/trunk-jpl/src/m/classes/matice.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/matice.js	(revision 19987)
+++ /issm/trunk-jpl/src/m/classes/matice.js	(revision 19988)
@@ -53,4 +53,8 @@
 			this.mantle_shear_modulus       = 1.45*Math.pow(10,11); // (Pa)
 			this.mantle_density             = 3.34;       // (g/cm^-3)
+			
+			//SLR
+			this.earth_density= 5512;  // average density of the Earth, (kg/m^3)
+
 
 		} //}}}
@@ -77,4 +81,6 @@
 			fielddisplay(this,'mantle_shear_modulus','Mantle shear modulus [Pa]');
 			fielddisplay(this,'mantle_density','Mantle density [g/cm^-3]');
+			fielddisplay(this,'earth_density','Mantle density [kg/m^-3]');
+
 		} // }}}
 		this.classname = function() {// {{{
@@ -96,4 +102,8 @@
 				checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',1);
 			}
+			if ArrayAnyEqual(ArrayIsMember(SealevelriseAnalysisEnum(),analyses),1)){
+				checkfield(md,'fieldname','materials.earth_density','>',0,'numel',1);
+			}
+
 
 		} // }}}
@@ -119,4 +129,6 @@
 			WriteData(fid,'object',this,'class','materials','fieldname','mantle_shear_modulus','format','Double');
 			WriteData(fid,'object',this,'class','materials','fieldname','mantle_density','format','Double','scale',10^3);
+			WriteData(fid,'object',this,'class','materials','fieldname','earth_density','format','Double');
+
 		}//}}}
 		this.fix=function() { //{{{
@@ -145,4 +157,8 @@
 		this.mantle_shear_modulus       = 0.;
 		this.mantle_density             = 0.;
+
+		//SLR
+		this.earth_density= 5512;  // average density of the Earth, (kg/m^3)
+
 		this.setdefaultparameters();
 		//}}}
Index: /issm/trunk-jpl/src/m/classes/matice.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/matice.py	(revision 19987)
+++ /issm/trunk-jpl/src/m/classes/matice.py	(revision 19988)
@@ -35,5 +35,10 @@
 		self.lithosphere_density        = 0.
 		self.mantle_shear_modulus       = 0.
-		self.mantle_density             = 0.
+		self.mantle_density             = 0.  
+		
+		#SLR
+		self.earth_density= 5512;  
+
+
 
 		self.setdefaultparameters()
@@ -61,4 +66,6 @@
 		string="%s\n%s"%(string,fielddisplay(self,"mantle_shear_modulus","Mantle shear modulus [Pa]"))
 		string="%s\n%s"%(string,fielddisplay(self,"mantle_density","Mantle density [g/cm^-3]"))
+		string="%s\n%s"%(string,fielddisplay(self,"earth_density","Mantle density [kg/m^-3]"))
+
 
 		return string
@@ -115,4 +122,8 @@
 		self.mantle_shear_modulus       = 1.45*10**11 # (Pa)
 		self.mantle_density             = 3.34        # (g/cm^-3)
+		
+		#SLR
+		self.earth_density= 5512;  # average density of the Earth, (kg/m^3)
+
 
 		return self
@@ -130,4 +141,6 @@
 		md = checkfield(md,'fieldname','materials.mantle_shear_modulus','>',0,'numel',[1]);
 		md = checkfield(md,'fieldname','materials.mantle_density','>',0,'numel',[1]);
+		md = checkfield(md,'fieldname','materials.earth_density','>',0,'numel',[1]);
+
 		return md
 	# }}}
@@ -155,3 +168,5 @@
 		WriteData(fid,'object',self,'class','materials','fieldname','mantle_shear_modulus','format','Double');
 		WriteData(fid,'object',self,'class','materials','fieldname','mantle_density','format','Double','scale',10.**3.);
+		WriteData(fid,'object',self,'class','materials','fieldname','earth_density','format','Double');
+
 	# }}}