Index: /issm/trunk/src/m/classes/initialization.m
===================================================================
--- /issm/trunk/src/m/classes/initialization.m	(revision 9685)
+++ /issm/trunk/src/m/classes/initialization.m	(revision 9685)
@@ -0,0 +1,44 @@
+%INITIALIZATION class definition
+%
+%   Usage:
+%      initialization=initialization();
+
+classdef initialization
+	properties (SetAccess=public) 
+		vx            = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1,'enum',VxEnum);
+		vy            = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1,'enum',VyEnum);
+		vz            = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1,'enum',VzEnum);
+		vel           = modelfield('default',NaN,'marshall',false);
+		pressure      = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1,'enum',PressureEnum);
+		temperature   = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1,'enum',TemperatureEnum);
+		watercolumn   = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1,'enum',WatercolumnEnum);
+		waterfraction = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1,'enum',WaterfractionEnum);
+	end
+	methods
+		function obj = initialization(varargin) % {{{
+			switch nargin
+				case 0
+					obj=setdefaultparameters(obj);
+				case 1
+					in=varargin{1};
+					if (isa(in,'numeric') & in==0),
+						% requesting templates do nothing
+					else
+						error('constructor not supported');
+					end 
+				otherwise
+					error('constructor not supported');
+			end
+		end % }}}
+		function obj = setdefaultparameters(obj) % {{{
+
+			%first, use the defaults provided by the properties definition above. 
+			fieldnames=fields(obj);
+			for i=1:length(fieldnames),
+				fieldname=fieldnames{i};
+				obj.(fieldname)=obj.(fieldname).default;
+			end
+
+		end % }}}
+	end
+end
Index: /issm/trunk/src/m/classes/model/model.m
===================================================================
--- /issm/trunk/src/m/classes/model/model.m	(revision 9684)
+++ /issm/trunk/src/m/classes/model/model.m	(revision 9685)
@@ -346,4 +346,5 @@
 			 if isfield(structmd,'vy'), md.initialization.vy=structmd.vy; end
 			 if isfield(structmd,'vz'), md.initialization.vz=structmd.vz; end
+			 if isfield(structmd,'vel'), md.initialization.vel=structmd.vel; end
 			 if isfield(structmd,'pressure'), md.initialization.pressure=structmd.pressure; end
 			 if isfield(structmd,'temperature'), md.initialization.temperature=structmd.temperature; end