Index: /issm/trunk-jpl/src/m/classes/constants.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/constants.m	(revision 26224)
+++ /issm/trunk-jpl/src/m/classes/constants.m	(revision 26225)
@@ -6,8 +6,8 @@
 classdef constants
 	properties (SetAccess=public) 
-		g                    = 0.;
-		omega                = 0.;
-		yts                  = 0.;
-		referencetemperature = 0.;
+		g                      = 0.;
+		omega                  = 0.;
+		yts                    = 0.;
+		referencetemperature   = 0.;
 		gravitational_constant = 0.;
 	end
@@ -35,5 +35,5 @@
 			self.referencetemperature=223.15;
 		
-			%gravity: 
+			%gravitational constant: 
 			self.gravitational_constant = 6.67259e-11;
 
Index: /issm/trunk-jpl/src/m/classes/constants.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/constants.py	(revision 26224)
+++ /issm/trunk-jpl/src/m/classes/constants.py	(revision 26225)
@@ -13,12 +13,12 @@
 
     def __init__(self):  # {{{
-        self.g = 0.
-        self.omega = 0.
-        self.yts = 0.
-        self.referencetemperature = 0.
+        self.g                      = 0.
+        self.omega                  = 0.
+        self.yts                    = 0.
+        self.referencetemperature   = 0.
+        self.gravitational_constant = 0.;
 
-    #set defaults
+        #set defaults
         self.setdefaultparameters()
-
     #}}}
     def __repr__(self):  # {{{
@@ -29,8 +29,8 @@
         string = "%s\n%s" % (string, fielddisplay(self, "yts", "number of seconds in a year [s / yr]"))
         string = "%s\n%s" % (string, fielddisplay(self, "referencetemperature", "reference temperature used in the enthalpy model [K]"))
+        string = "%s\n%s" % (string, fielddisplay(self, "gravitational_constant", "Newtonian constant of gravitation [m^3/kg/s^2]"))
 
         return string
     #}}}
-
     def setdefaultparameters(self):  # {{{
         #acceleration due to gravity (m / s^2)
@@ -46,7 +46,9 @@
         self.referencetemperature = 223.15
 
+        #gravitational constant:
+        self.gravitational_constant = 6.67259e-11
+
         return self
     #}}}
-
     def checkconsistency(self, md, solution, analyses):  # {{{
         md = checkfield(md, 'fieldname', 'constants.g', '>=', 0, 'size', [1])
@@ -54,11 +56,12 @@
         md = checkfield(md, 'fieldname', 'constants.yts', '>', 0, 'size', [1])
         md = checkfield(md, 'fieldname', 'constants.referencetemperature', 'size', [1])
+        md = checkfield(md, 'fieldname', 'constants.gravitational_constant','size',[1]);
 
         return md
     # }}}
-
     def marshall(self, prefix, md, fid):  # {{{
         WriteData(fid, prefix, 'object', self, 'fieldname', 'g', 'format', 'Double')
         WriteData(fid, prefix, 'object', self, 'fieldname', 'yts', 'format', 'Double')
         WriteData(fid, prefix, 'object', self, 'fieldname', 'referencetemperature', 'format', 'Double')
+        WriteData(fid, prefix, 'object', self, 'fieldname', 'gravitational_constant', 'format', 'Double');
     # }}}