Changeset 24805

Timestamp:

05/05/20 11:21:02 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing 245

File:

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

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

@@ -6 +6 @@
 classdef SMBpddSicopolis
         properties (SetAccess=public)
-                precipitation                                   = NaN;
-                monthlytemperatures                     = NaN;
-                temperature_anomaly                     = NaN;
-                precipitation_anomaly           = NaN;
-                smb_corr                                                        = NaN;
-                desfac                                                  = 0;
-                s0p                                                             = NaN;
-                s0t                                                             = NaN;
-                rlaps                                                           = 0;
-                isfirnwarming                                   = 0;
-                steps_per_step=1
-                requested_outputs                               = {};
+                precipitation         = NaN;
+                monthlytemperatures   = NaN;
+                temperature_anomaly   = NaN;
+                precipitation_anomaly = NaN;
+                smb_corr              = NaN;
+                desfac                = 0;
+                s0p                   = NaN;
+                s0t                   = NaN;
+                rlaps                 = 0;
+                isfirnwarming         = 0;
+                steps_per_step        = 1
+                averaging             = 0
+                requested_outputs     = {};
         end
         methods
@@ -86 +87 @@
                         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,'isfirnwarming','is firnwarming (Reeh 1991) activated (0 or 1, default is 1)');
                         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'));
@@ -127 +130 @@
                         WriteData(fid,prefix,'object',self,'class','smb','fieldname','smb_corr','format','DoubleMat','mattype',1,'scale',1./yts,'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