source: issm/oecreview/Archive/21337-21723/ISSM-21417-21418.diff@ 21726

../trunk-jpl/src/m/classes/model.m

@@ -247 +247 @@
                         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');

../trunk-jpl/src/m/classes/linearbasalforcings.m

@@ -6 +6 @@
 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 +22 @@
                                         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),