Index: /issm/trunk/src/m/model/extrude.m
===================================================================
--- /issm/trunk/src/m/model/extrude.m	(revision 6238)
+++ /issm/trunk/src/m/model/extrude.m	(revision 6239)
@@ -232,5 +232,7 @@
 md.melting_rate=project3d(md,md.melting_rate,'node',1); 
 md.observed_temperature=project3d(md,md.observed_temperature,'node'); 
-md.geothermalflux=project3d(md,md.geothermalflux,'node',1); %bedrock only gets geothermal flux
+if ~isnan(md.geothermalflux)
+	md.geothermalflux=project3d(md,md.geothermalflux,'node',1); %bedrock only gets geothermal flux
+end
 
 %increase connectivity if less than 25:
Index: /issm/trunk/src/m/model/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 6238)
+++ /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 6239)
@@ -276,5 +276,5 @@
 
 	%VELOCITIES AND PRESSURE
-	fields={'vx','vy','vz','pressure'};
+	fields={'vx','vy','vz','pressure','geothermalflux'};
 	checksize(md,fields,[md.numberofgrids 1]);
 	checknan(md,fields);