Index: /issm/trunk-jpl/src/m/classes/SMBpddSicopolis.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBpddSicopolis.m	(revision 24804)
+++ /issm/trunk-jpl/src/m/classes/SMBpddSicopolis.m	(revision 24805)
@@ -6,16 +6,17 @@
 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,4 +87,5 @@
 			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,4 +106,5 @@
 			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,4 +130,5 @@
 			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