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

Timestamp:

01/09/23 08:52:51 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added age to matlab classes

File:

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

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

@@ -26 +26 @@
                 sample              = NaN;
                 debris              = NaN;
+                age                 = NaN;
         end
         methods
@@ -127 +128 @@
                                 if ~isnan(md.initialization.debris)
                                         md = checkfield(md,'fieldname','initialization.debris','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+                                end
+                        end
+                        if ismember('AgeAnalysis',analyses),
+                                if ~isnan(md.initialization.age)
+                                        md = checkfield(md,'fieldname','initialization.age','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
                                 end
                         end
@@ -152 +158 @@
                         fielddisplay(self,'str','Steric sea level.');
                         fielddisplay(self,'debris','Surface debris layer [m]');
+                        fielddisplay(self,'age','Initial age [yr]');
                 end % }}}
                 function marshall(self,prefix,md,fid) % {{{
@@ -174 +181 @@
                         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);
+                        WriteData(fid,prefix,'object',self,'fieldname','debris','format','DoubleMat','mattype',1,'scale',yts);
 
                         if md.thermal.isenthalpy,
@@ -206 +213 @@
                         self.str=project3d(md,'vector',self.str,'type','node','layer',1);
                         self.str=project3d(md,'vector',self.debris,'type','node','layer',1);
+                        self.str=project3d(md,'vector',self.age,'type','node','layer',1);
 
                         %Lithostatic pressure by default
@@ -228 +236 @@
                         writejs1Darray(fid,[modelname '.initialization.sample'],self.sample);
                         writejs1Darray(fid,[modelname '.initialization.debris'],self.debris);
+                        writejs1Darray(fid,[modelname '.initialization.age'],self.age);
 
                 end % }}}