Changeset 24793 for issm/trunk-jpl/src/m/classes/SMBsemic.m

Timestamp:

05/05/20 06:26:00 (5 years ago)

Author:

bdef

Message:

NEW: adding several averaging methods to pass from substep to main step

File:

• issm/trunk-jpl/src/m/classes/SMBsemic.m (modified) (4 diffs)

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

@@ -20 +20 @@
                 s0gcm                                   = NaN;
                 steps_per_step = 1;
+                averaging = 0;
                 requested_outputs = {};
         end
@@ -80 +81 @@
                         end
                         md = checkfield(md,'fieldname','smb.steps_per_step','>=',1,'numel',[1]);
+                        md = checkfield(md,'fieldname','smb.averaging',,'numel',[1],'values',[0 1 2]);
                         md = checkfield(md,'fieldname','smb.requested_outputs','stringrow',1);
 
@@ -104 +106 @@
                         fielddisplay(self,'s0gcm','GCM reference elevation; (default is 0) [m]');
                         fielddisplay(self, 'steps_per_step', 'number of smb steps per time step');
+                        fielddisplay(self, 'averaging', 'averaging methods from short to long steps');
+                        disp(sprintf('%55s  0: Arithmetic (default)');
+                        disp(sprintf('%55s  0: Geometric');
+                        disp(sprintf('%55s  0: Harmonic');
                         fielddisplay(self,'requested_outputs','additional outputs requested');
                 end % }}}
@@ -124 +130 @@
                         WriteData(fid,prefix,'object',self,'class','smb','fieldname','dailytemperature','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
                         WriteData(fid, prefix, 'object', self, 'fieldname', 'steps_per_step', 'format', 'Integer');
+                        WriteData(fid, prefix, 'object', self, 'fieldname', 'averaging', 'format', 'Integer');
                         %process requested outputs
                         outputs = self.requested_outputs;