Changeset 27297 for issm/trunk-jpl/src/m/classes/initialization.m

Timestamp:

10/06/22 23:34:44 (2 years ago)

Author:

rueckamp

Message:

NEW: SMB debris model based on Mayer & Licuilli 2021

File:

• issm/trunk-jpl/src/m/classes/initialization.m (modified) (6 diffs)

issm/trunk-jpl/src/m/classes/initialization.m

@@ -25 +25 @@
                 str                 = NaN;
                 sample              = NaN;
+                debris              = NaN;
         end
         methods
@@ -121 +122 @@
                                 if ~isnan(md.initialization.sample)
                                         md = checkfield(md,'fieldname','initialization.sample','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+                                end
+                        end
+                        if ismember('DebrisAnalysis',analyses),
+                                if ~isnan(md.initialization.debris)
+                                        md = checkfield(md,'fieldname','initialization.debris','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
                                 end
                         end
@@ -145 +151 @@
                         fielddisplay(self,'dsl','Dynamic sea level.');
                         fielddisplay(self,'str','Steric sea level.');
+                        fielddisplay(self,'debris','Surface debris layer [m]');
                 end % }}}
                 function marshall(self,prefix,md,fid) % {{{
@@ -167 +174 @@
                         WriteData(fid,prefix,'object',self,'fieldname','hydraulic_potential','format','DoubleMat','mattype',1);
                         WriteData(fid,prefix,'object',self,'fieldname','sample','format','DoubleMat','mattype',1);
+                        WriteData(fid,prefix,'object',self,'fieldname','debris','format','DoubleMat','mattype',1);
 
                         if md.thermal.isenthalpy,
@@ -197 +205 @@
                         self.dsl=project3d(md,'vector',self.dsl,'type','node','layer',1);
                         self.str=project3d(md,'vector',self.str,'type','node','layer',1);
+                        self.str=project3d(md,'vector',self.debris,'type','node','layer',1);
 
                         %Lithostatic pressure by default
@@ -218 +227 @@
                         writejs1Darray(fid,[modelname '.initialization.channel'],self.channelarea);
                         writejs1Darray(fid,[modelname '.initialization.sample'],self.sample);
+                        writejs1Darray(fid,[modelname '.initialization.debris'],self.debris);
 
                 end % }}}