Index: /issm/trunk-jpl/src/m/classes/basalforcings.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/basalforcings.m	(revision 26951)
+++ /issm/trunk-jpl/src/m/classes/basalforcings.m	(revision 26952)
@@ -15,4 +15,6 @@
 				case 0
 					self=setdefaultparameters(self);
+				case 1
+					self =structtoobj(basalforcings(),varargin{1});
 				otherwise
 					error('constructor not supported');
Index: /issm/trunk-jpl/src/m/classes/initialization.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.m	(revision 26951)
+++ /issm/trunk-jpl/src/m/classes/initialization.m	(revision 26952)
@@ -39,5 +39,5 @@
 		function md = checkconsistency(self,md,solution,analyses) % {{{
 			if ismember('StressbalanceAnalysis',analyses) & ~(strcmp(solution,'TransientSolution') & md.transient.isstressbalance == 0),
-				if ~(isnan(md.initialization.vx) | isnan(md.initialization.vy)),
+				if numel(md.initialization.vx)>1 | numel(md.initialization.vy)>1
 					md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
 					md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);