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

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),