Index: ../trunk-jpl/src/m/classes/inversion.py
===================================================================
--- ../trunk-jpl/src/m/classes/inversion.py	(revision 18608)
+++ ../trunk-jpl/src/m/classes/inversion.py	(revision 18609)
@@ -32,6 +32,7 @@
 		self.vz_obs                      = float('NaN')
 		self.vel_obs                     = float('NaN')
 		self.thickness_obs               = float('NaN')
+		self.surface_obs                 = float('NaN')
 
 		#set defaults
 		self.setdefaultparameters()
@@ -55,6 +56,7 @@
 		string="%s\n%s"%(string,fielddisplay(self,'vy_obs','observed velocity y component [m/yr]'))
 		string="%s\n%s"%(string,fielddisplay(self,'vel_obs','observed velocity magnitude [m/yr]'))
 		string="%s\n%s"%(string,fielddisplay(self,'thickness_obs','observed thickness [m]'))
+		string="%s\n%s"%(string,fielddisplay(self,'surface_obs','observed surface elevation [m]'))
 		string="%s\n%s"%(string,'Available cost functions:')
 		string="%s\n%s"%(string,'   101: SurfaceAbsVelMisfit')
 		string="%s\n%s"%(string,'   102: SurfaceRelVelMisfit')
@@ -160,6 +162,7 @@
 		WriteData(fid,'object',self,'fieldname','vy_obs','format','DoubleMat','mattype',1,'scale',1./yts)
 		WriteData(fid,'object',self,'fieldname','vz_obs','format','DoubleMat','mattype',1,'scale',1./yts)
 		WriteData(fid,'object',self,'fieldname','thickness_obs','format','DoubleMat','mattype',1)
+		WriteData(fid,'object',self,'fieldname','surface_obs','format','DoubleMat','mattype',1)
 
 		#process control parameters
 		num_control_parameters=len(self.control_parameters)
