Index: /issm/trunk-jpl/src/m/classes/linearbasalforcings.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/linearbasalforcings.m	(revision 21417)
+++ /issm/trunk-jpl/src/m/classes/linearbasalforcings.m	(revision 21418)
@@ -7,7 +7,7 @@
 	properties (SetAccess=public) 
 		groundedice_melting_rate  = NaN;
-		deepwater_melting_rate    = NaN;
-		deepwater_elevation       = NaN;
-		upperwater_elevation      = NaN;
+		deepwater_melting_rate    = 0.;
+		deepwater_elevation       = 0.;
+		upperwater_elevation      = 0.;
 		geothermalflux            = NaN;
 	end
@@ -22,4 +22,8 @@
 					error('constructor not supported');
 			end
+		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 = initialize(self,md) % {{{
Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 21417)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 21418)
@@ -248,4 +248,5 @@
 			if numel(md.levelset.spclevelset)>1, md.levelset.spclevelset=project2d(md,md.levelset.spclevelset,md.mesh.numberoflayers); end
 			md.thermal.spctemperature=project2d(md,md.thermal.spctemperature,md.mesh.numberoflayers);
+			md.damage.spcdamage=project2d(md,md.damage.spcdamage,md.mesh.numberoflayers);
 
 			% Hydrologydc variables