Index: ../trunk-jpl/src/m/classes/model.m
===================================================================
--- ../trunk-jpl/src/m/classes/model.m	(revision 21417)
+++ ../trunk-jpl/src/m/classes/model.m	(revision 21418)
@@ -247,6 +247,7 @@
 			if ~isnan(md.damage.spcdamage), md.damage.spcdamage=project2d(md,md.damage.spcdamage,md.mesh.numberoflayers); end
 			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
 			if isa(md.hydrology,'hydrologydc');
Index: ../trunk-jpl/src/m/classes/linearbasalforcings.m
===================================================================
--- ../trunk-jpl/src/m/classes/linearbasalforcings.m	(revision 21417)
+++ ../trunk-jpl/src/m/classes/linearbasalforcings.m	(revision 21418)
@@ -6,9 +6,9 @@
 classdef linearbasalforcings
 	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
 	methods
@@ -22,6 +22,10 @@
 					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) % {{{
 
 			if isnan(self.groundedice_melting_rate),