Index: /issm/trunk-jpl/src/m/classes/mismipbasalforcings.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/mismipbasalforcings.m	(revision 19910)
+++ /issm/trunk-jpl/src/m/classes/mismipbasalforcings.m	(revision 19911)
@@ -20,4 +20,8 @@
              
         end % }}}
+		function self = extrude(self,md) % {{{
+			self.groundedice_melting_rate=project3d(md,'vector',self.groundedice_melting_rate,'type','node','layer',1); 
+			self.geothermalflux=project3d(md,'vector',self.geothermalflux,'type','node','layer',1); %bedrock only gets geothermal flux
+		end % }}}
 		function self = mismipbasalforcings(varargin) % {{{
 			switch nargin
Index: /issm/trunk-jpl/src/m/classes/mismipbasalforcings.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/mismipbasalforcings.py	(revision 19910)
+++ /issm/trunk-jpl/src/m/classes/mismipbasalforcings.py	(revision 19911)
@@ -38,8 +38,7 @@
 	return string
     #}}}
-    def extrude(self,md): # {{{
-	self.coefficient=project3d(md,'vector',self.coefficient,'type','node','layer',1)
-	self.p=project3d(md,'vector',self.p,'type','element')
-	self.q=project3d(md,'vector',self.q,'type','element')
+	 def extrude(self,md): # {{{
+		 self.groundedice_melting_rate=project3d(md,'vector',self.groundedice_melting_rate,'type','node','layer',1)
+		 self.geothermalflux=project3d(md,'vector',self.geothermalflux,'type','node','layer',1)    #bedrock only gets geothermal flux
 	return self
     #}}}