Index: /issm/trunk-jpl/src/m/classes/gia.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/gia.m	(revision 14827)
+++ /issm/trunk-jpl/src/m/classes/gia.m	(revision 14828)
@@ -37,8 +37,8 @@
 					%recover the furthest time "in time": 
 					if(thickness(end,end)~=md.timestepping.start_time),
-						md = checkmessage(md,'if isprognostic is on, transient thickness forcing',...
-							' for the gia model should not be provided in the future.',...
-							' Synchronize your start_time to correspond to the most recent transient',...
-							' thickness forcing timestep');
+						md = checkmessage(md,['if isprognostic is on, transient thickness forcing'...
+							' for the gia model should not be provided in the future.'...
+							' Synchronize your start_time to correspond to the most recent transient'...
+							' thickness forcing timestep']);
 					end
 				end
@@ -55,5 +55,5 @@
 		end % }}}
 		function marshall(obj,fid) % {{{
-			WriteData(fid,'data',obj.lithosphere_thickness,'format','DoubleMat','mattype',1,'enum',GiaLithosphereThicknessEnum());
+			WriteData(fid,'object',obj,'fieldname','lithosphere_thickness','format','DoubleMat','mattype',1);
 			WriteData(fid,'object',obj,'fieldname','output_rates','format','Integer');
 			WriteData(fid,'object',obj,'fieldname','cross_section_shape','format','Integer');
Index: /issm/trunk-jpl/src/m/classes/gia.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/gia.py	(revision 14827)
+++ /issm/trunk-jpl/src/m/classes/gia.py	(revision 14828)
@@ -43,5 +43,5 @@
 			return md 
 		
-		md = checkfield(md,'gia.lithosphere_thickness','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0)
+		md = checkfield(md,'gia.lithosphere_thickness','NaN',1,'size',[md.mesh.numberofvertices,1],'>',0)
 		md = checkfield(md,'gia.output_rates','numel',[1],'values',[1,2])
 		md = checkfield(md,'gia.cross_section_shape','numel',[1],'values',[1,2])
@@ -53,5 +53,5 @@
 	# }}}
 	def marshall(self,fid):    # {{{
-		WriteData(fid,'data',self.lithosphere_thickness,'format','DoubleMat','mattype',1,'enum',GiaLithosphereThicknessEnum());
+		WriteData(fid,'object',self,'fieldname','lithosphere_thickness','format','DoubleMat','mattype',1);
 		WriteData(fid,'object',self,'fieldname','output_rates','format','Integer');
 		WriteData(fid,'object',self,'fieldname','cross_section_shape','format','Integer');
Index: /issm/trunk-jpl/src/m/classes/model/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model/model.py	(revision 14827)
+++ /issm/trunk-jpl/src/m/classes/model/model.py	(revision 14828)
@@ -118,4 +118,5 @@
 		        'steadystate',\
 		        'transient',\
+				  'gia',\
 		        'autodiff',\
 		        'flaim',\