Index: sm/trunk/src/m/model/setelementstype.m
===================================================================
--- /issm/trunk/src/m/model/setelementstype.m	(revision 9663)
+++ 	(revision )
@@ -1,237 +1,0 @@
-function md=setelementstype(md,varargin)
-%SETELEMENTSTYPE - associate a solution type to each element
-%
-%   This routine works like plotmodel: it works with an even number of inputs
-%   'hutter','macayeal','pattyn','stokes' and 'fill' are the possible options
-%   that must be followed by the corresponding exp file or flags list
-%   It can either be a domain file (argus type, .exp extension), or an array of element flags. 
-%   If user wants every element outside the domain to be 
-%   setelementstyped, add '~' to the name of the domain file (ex: '~Pattyn.exp');
-%   an empty string '' will be considered as an empty domain
-%   a string 'all' will be considered as the entire domain
-%   You can specify the type of coupling, 'penalties' or 'tiling', to use with the input 'coupling'
-%
-%   Usage:
-%      md=setelementstype(md,varargin)
-%
-%   Example:
-%      md=setelementstype(md,'pattyn','Pattyn.exp','macayeal',md.mask.elementonfloatingice,'fill','hutter');
-%      md=setelementstype(md,'pattyn','Pattyn.exp',fill','hutter','coupling','tiling');
-
-%some checks on list of arguments
-if ((nargin<2) | (nargout~=1)),
-	error('setelementstype error message');
-end
-
-%Find_out what kind of coupling to use
-options=pairoptions(varargin{:});
-coupling_method=getfieldvalue(options,'coupling','tiling');
-if (~strcmpi(coupling_method,'tiling') & ~strcmpi(coupling_method,'penalties')),
-	error('coupling type can only be: tiling or penalties');
-end
-
-[hutterflag macayealflag pattynflag stokesflag filltype]=recover_areas(md,varargin{:});
-
-%Flag the elements that has not been flagged as filltype
-if strcmpi(filltype,'hutter'),
-	hutterflag(find(~macayealflag & ~pattynflag))=1;
-elseif strcmpi(filltype,'macayeal'),
-	macayealflag(find(~hutterflag & ~pattynflag & ~stokesflag))=1;
-elseif strcmpi(filltype,'pattyn'),
-	pattynflag(find(~hutterflag & ~macayealflag & ~stokesflag))=1;
-end
-
-%check that each element has at least one flag
-if any(hutterflag+ macayealflag+pattynflag+stokesflag==0),
-	error('setelementstype error message: elements type not assigned, must be specified')
-end
-
-%check that each element has only one flag
-if any(hutterflag+ macayealflag+pattynflag+stokesflag>1),
-	disp('setelementstype warning message: some elements have several types, higher order type is used for them')
-	hutterflag(find(hutterflag & macayealflag))=0;
-	hutterflag(find(hutterflag & pattynflag))=0;
-	macayealflag(find(macayealflag & pattynflag))=0;
-end
-
-%Check that no pattyn or stokes for 2d mesh
-if (md.dim==2),
-	if any(stokesflag | pattynflag)
-		error('setelementstype error message: stokes and pattyn elements no allowed in 2d mesh, extrude it first')
-	end
-end
-
-%Stokes can only be used alone for now:
-if any(stokesflag) &any(hutterflag),
-	error('setelementstype error message: stokes cannot be used with any other model for now, put stokes everywhere')
-end
-
-%add in model who is who
-md.flowequation.element_equation=zeros(md.numberofelements,1);
-
-%1: Hutter elements
-nodeonhutter=zeros(md.numberofnodes,1);
-nodeonhutter(md.elements(find(hutterflag),:))=1;
-md.flowequation.element_equation(find(hutterflag))=1;
-
-%2: MacAyeal elements
-nodeonmacayeal=zeros(md.numberofnodes,1);
-nodeonmacayeal(md.elements(find(macayealflag),:))=1;
-md.flowequation.bordermacayeal=nodeonmacayeal;
-md.flowequation.element_equation(find(macayealflag))=2;
-
-%3: Pattyn elements
-nodeonpattyn=zeros(md.numberofnodes,1);
-nodeonpattyn(md.elements(find(pattynflag),:))=1;
-md.flowequation.borderpattyn=nodeonpattyn;
-md.flowequation.element_equation(find(pattynflag))=3;
-
-%4: Stokes elements
-%First modify stokesflag to get rid of elements contrained everywhere (spc + border with pattyn or macayeal)
-if any(stokesflag),
-	nodeonstokes=zeros(md.numberofnodes,1);
-	nodeonstokes(md.elements(find(stokesflag),:))=1;
-	fullspcnodes=double((~isnan(md.spcvx)+~isnan(md.spcvy)+~isnan(md.spcvz))==3 | (nodeonpattyn & nodeonstokes));         %find all the nodes on the boundary of the domain without icefront
-	fullspcelems=double(sum(fullspcnodes(md.elements),2)==6);         %find all the nodes on the boundary of the domain without icefront
-	stokesflag(find(fullspcelems))=0;
-end
-nodeonstokes=zeros(md.numberofnodes,1);
-nodeonstokes(md.elements(find(stokesflag),:))=1;
-md.flowequation.borderstokes=nodeonstokes;
-md.flowequation.element_equation(find(stokesflag))=4;
-
-%Then complete with NoneApproximation or the other model used if there is no stokes
-if any(stokesflag), 
-	if any(pattynflag), %fill with pattyn
-		pattynflag(~stokesflag)=1;
-		nodeonpattyn(md.elements(find(pattynflag),:))=1;
-		md.flowequation.borderpattyn=nodeonpattyn;
-		md.flowequation.element_equation(find(~stokesflag))=3;
-	elseif any(macayealflag), %fill with macayeal
-		macayealflag(~stokesflag)=1;
-		nodeonmacayeal(md.elements(find(macayealflag),:))=1;
-		md.flowequation.bordermacayeal=nodeonmacayeal;
-		md.flowequation.element_equation(find(~stokesflag))=2;
-	else %fill with none 
-	%5: None elements (non Stokes)
-		md.flowequation.element_equation(find(~stokesflag))=0;
-	end
-end
-
-%Now take care of the coupling between MacAyeal and Pattyn
-md.penalties=[];
-nodeonmacayealpattyn=zeros(md.numberofnodes,1);
-nodeonpattynstokes=zeros(md.numberofnodes,1);
-nodeonmacayealstokes=zeros(md.numberofnodes,1);
-if strcmpi(coupling_method,'penalties'),
-	%Create the border nodes between Pattyn and MacAyeal and extrude them
-	numnodes2d=md.numberofnodes2d;
-	numlayers=md.numlayers;
-	bordernodes2d=find(nodeonpattyn(1:numnodes2d) & nodeonmacayeal(1:numnodes2d)); %Nodes connected to two different types of elements
-
-	%initialize and fill in penalties structure
-	if ~isnan(bordernodes2d),
-		penalties=[];
-		for	i=1:numlayers-1,
-			penalties=[penalties; [bordernodes2d bordernodes2d+md.numberofnodes2d*(i)]];
-		end
-		md.penalties=penalties;
-	end
-elseif strcmpi(coupling_method,'tiling'),
-	if any(macayealflag) & any(pattynflag), %coupling macayeal pattyn
-		%Find node at the border
-		nodeonmacayealpattyn(find(nodeonmacayeal & nodeonpattyn))=1;
-		%Macayeal elements in contact with this layer become MacAyealPattyn elements
-		matrixelements=ismember(md.elements,find(nodeonmacayealpattyn));
-		commonelements=sum(matrixelements,2)~=0;
-		commonelements(find(pattynflag))=0; %only one layer: the elements previously in macayeal
-		macayealflag(find(commonelements))=0; %these elements are now macayealpattynelements
-		macayealpattynflag=zeros(md.numberofelements,1);
-		macayealpattynflag(find(commonelements))=1;
-		nodeonmacayeal=zeros(md.numberofnodes,1);
-		nodeonmacayeal(md.elements(find(macayealflag),:))=1;
-		md.flowequation.bordermacayeal=nodeonmacayeal;
-
-		%Create MacaAyealPattynApproximation where needed
-		md.flowequation.element_equation(find(macayealpattynflag))=5;
-
-		%Now recreate nodeonmacayeal and nodeonmacayealpattyn
-		nodeonmacayealpattyn(md.elements(find(macayealpattynflag),:))=1;
-	elseif any(pattynflag) & any(stokesflag), %coupling pattyn stokes
-		%Find node at the border
-		nodeonpattynstokes(find(nodeonpattyn & nodeonstokes))=1;
-		%Stokes elements in contact with this layer become PattynStokes elements
-		matrixelements=ismember(md.elements,find(nodeonpattynstokes));
-		commonelements=sum(matrixelements,2)~=0;
-		commonelements(find(pattynflag))=0; %only one layer: the elements previously in macayeal
-		stokesflag(find(commonelements))=0; %these elements are now macayealpattynelements
-		pattynstokesflag=zeros(md.numberofelements,1);
-		pattynstokesflag(find(commonelements))=1;
-		nodeonstokes=zeros(md.numberofnodes,1);
-		nodeonstokes(md.elements(find(stokesflag),:))=1;
-		md.flowequation.borderstokes=nodeonstokes;
-
-		%Create MacaAyealPattynApproximation where needed
-		md.flowequation.element_equation(find(pattynstokesflag))=7;
-
-		%Now recreate nodeonpattynstokes
-		nodeonpattynstokes=zeros(md.numberofnodes,1);
-		nodeonpattynstokes(md.elements(find(pattynstokesflag),:))=1;
-	elseif any(stokesflag) & any(macayealflag),
-		%Find node at the border
-		nodeonmacayealstokes(find(nodeonmacayeal & nodeonstokes))=1;
-		%Stokes elements in contact with this layer become MacAyealStokes elements
-		matrixelements=ismember(md.elements,find(nodeonmacayealstokes));
-		commonelements=sum(matrixelements,2)~=0;
-		commonelements(find(macayealflag))=0; %only one layer: the elements previously in macayeal
-		stokesflag(find(commonelements))=0; %these elements are now macayealmacayealelements
-		macayealstokesflag=zeros(md.numberofelements,1);
-		macayealstokesflag(find(commonelements))=1;
-		nodeonstokes=zeros(md.numberofnodes,1);
-		nodeonstokes(md.elements(find(stokesflag),:))=1;
-		md.flowequation.borderstokes=nodeonstokes;
-
-		%Create MacaAyeal Approximation where needed
-		md.flowequation.element_equation(find(macayealstokesflag))=6;
-
-		%Now recreate nodeonmacayealstokes
-		nodeonmacayealstokes=zeros(md.numberofnodes,1);
-		nodeonmacayealstokes(md.elements(find(macayealstokesflag),:))=1;
-	elseif any(stokesflag) & any(hutterflag),
-		error('type of coupling not supported yet');
-	end
-end
-
-%Create vertices_type
-md.flowequation.vertex_equation=zeros(md.numberofnodes,1);
-pos=find(nodeonhutter);
-md.flowequation.vertex_equation(pos)=1;
-pos=find(nodeonmacayeal);
-md.flowequation.vertex_equation(pos)=2;
-pos=find(nodeonpattyn);
-md.flowequation.vertex_equation(pos)=3;
-pos=find(nodeonhutter);
-md.flowequation.vertex_equation(pos)=1;
-pos=find(nodeonpattyn & nodeonmacayeal);
-md.flowequation.vertex_equation(pos)=3;
-pos=find(nodeonmacayealpattyn);
-md.flowequation.vertex_equation(pos)=5;
-pos=find(nodeonstokes);
-md.flowequation.vertex_equation(pos)=4;
-if any(stokesflag),
-	pos=find(~nodeonstokes);
-	if(~any(pattynflag) & ~any(macayealflag)),
-		md.flowequation.vertex_equation(pos)=0;
-	end
-end
-pos=find(nodeonpattynstokes);
-md.flowequation.vertex_equation(pos)=7;
-pos=find(nodeonmacayealstokes);
-md.flowequation.vertex_equation(pos)=6;
-
-%figure out solution types
-md.flowequation.ishutter=double(any(md.flowequation.element_equation==1));
-md.flowequation.ismacayealpattyn=double(any(md.flowequation.element_equation==2 | md.flowequation.element_equation==3));
-md.flowequation.isstokes=double(any(md.flowequation.element_equation==4));
-
-end
Index: /issm/trunk/src/m/model/setflowequation.m
===================================================================
--- /issm/trunk/src/m/model/setflowequation.m	(revision 9664)
+++ /issm/trunk/src/m/model/setflowequation.m	(revision 9664)
@@ -0,0 +1,237 @@
+function md=setflowequation(md,varargin)
+%SETELEMENTSTYPE - associate a solution type to each element
+%
+%   This routine works like plotmodel: it works with an even number of inputs
+%   'hutter','macayeal','pattyn','stokes' and 'fill' are the possible options
+%   that must be followed by the corresponding exp file or flags list
+%   It can either be a domain file (argus type, .exp extension), or an array of element flags. 
+%   If user wants every element outside the domain to be 
+%   setflowequationd, add '~' to the name of the domain file (ex: '~Pattyn.exp');
+%   an empty string '' will be considered as an empty domain
+%   a string 'all' will be considered as the entire domain
+%   You can specify the type of coupling, 'penalties' or 'tiling', to use with the input 'coupling'
+%
+%   Usage:
+%      md=setflowequation(md,varargin)
+%
+%   Example:
+%      md=setflowequation(md,'pattyn','Pattyn.exp','macayeal',md.mask.elementonfloatingice,'fill','hutter');
+%      md=setflowequation(md,'pattyn','Pattyn.exp',fill','hutter','coupling','tiling');
+
+%some checks on list of arguments
+if ((nargin<2) | (nargout~=1)),
+	error('setflowequation error message');
+end
+
+%Find_out what kind of coupling to use
+options=pairoptions(varargin{:});
+coupling_method=getfieldvalue(options,'coupling','tiling');
+if (~strcmpi(coupling_method,'tiling') & ~strcmpi(coupling_method,'penalties')),
+	error('coupling type can only be: tiling or penalties');
+end
+
+[hutterflag macayealflag pattynflag stokesflag filltype]=recover_areas(md,varargin{:});
+
+%Flag the elements that has not been flagged as filltype
+if strcmpi(filltype,'hutter'),
+	hutterflag(find(~macayealflag & ~pattynflag))=1;
+elseif strcmpi(filltype,'macayeal'),
+	macayealflag(find(~hutterflag & ~pattynflag & ~stokesflag))=1;
+elseif strcmpi(filltype,'pattyn'),
+	pattynflag(find(~hutterflag & ~macayealflag & ~stokesflag))=1;
+end
+
+%check that each element has at least one flag
+if any(hutterflag+ macayealflag+pattynflag+stokesflag==0),
+	error('setflowequation error message: elements type not assigned, must be specified')
+end
+
+%check that each element has only one flag
+if any(hutterflag+ macayealflag+pattynflag+stokesflag>1),
+	disp('setflowequation warning message: some elements have several types, higher order type is used for them')
+	hutterflag(find(hutterflag & macayealflag))=0;
+	hutterflag(find(hutterflag & pattynflag))=0;
+	macayealflag(find(macayealflag & pattynflag))=0;
+end
+
+%Check that no pattyn or stokes for 2d mesh
+if (md.dim==2),
+	if any(stokesflag | pattynflag)
+		error('setflowequation error message: stokes and pattyn elements no allowed in 2d mesh, extrude it first')
+	end
+end
+
+%Stokes can only be used alone for now:
+if any(stokesflag) &any(hutterflag),
+	error('setflowequation error message: stokes cannot be used with any other model for now, put stokes everywhere')
+end
+
+%add in model who is who
+md.flowequation.element_equation=zeros(md.numberofelements,1);
+
+%1: Hutter elements
+nodeonhutter=zeros(md.numberofnodes,1);
+nodeonhutter(md.elements(find(hutterflag),:))=1;
+md.flowequation.element_equation(find(hutterflag))=1;
+
+%2: MacAyeal elements
+nodeonmacayeal=zeros(md.numberofnodes,1);
+nodeonmacayeal(md.elements(find(macayealflag),:))=1;
+md.flowequation.bordermacayeal=nodeonmacayeal;
+md.flowequation.element_equation(find(macayealflag))=2;
+
+%3: Pattyn elements
+nodeonpattyn=zeros(md.numberofnodes,1);
+nodeonpattyn(md.elements(find(pattynflag),:))=1;
+md.flowequation.borderpattyn=nodeonpattyn;
+md.flowequation.element_equation(find(pattynflag))=3;
+
+%4: Stokes elements
+%First modify stokesflag to get rid of elements contrained everywhere (spc + border with pattyn or macayeal)
+if any(stokesflag),
+	nodeonstokes=zeros(md.numberofnodes,1);
+	nodeonstokes(md.elements(find(stokesflag),:))=1;
+	fullspcnodes=double((~isnan(md.spcvx)+~isnan(md.spcvy)+~isnan(md.spcvz))==3 | (nodeonpattyn & nodeonstokes));         %find all the nodes on the boundary of the domain without icefront
+	fullspcelems=double(sum(fullspcnodes(md.elements),2)==6);         %find all the nodes on the boundary of the domain without icefront
+	stokesflag(find(fullspcelems))=0;
+end
+nodeonstokes=zeros(md.numberofnodes,1);
+nodeonstokes(md.elements(find(stokesflag),:))=1;
+md.flowequation.borderstokes=nodeonstokes;
+md.flowequation.element_equation(find(stokesflag))=4;
+
+%Then complete with NoneApproximation or the other model used if there is no stokes
+if any(stokesflag), 
+	if any(pattynflag), %fill with pattyn
+		pattynflag(~stokesflag)=1;
+		nodeonpattyn(md.elements(find(pattynflag),:))=1;
+		md.flowequation.borderpattyn=nodeonpattyn;
+		md.flowequation.element_equation(find(~stokesflag))=3;
+	elseif any(macayealflag), %fill with macayeal
+		macayealflag(~stokesflag)=1;
+		nodeonmacayeal(md.elements(find(macayealflag),:))=1;
+		md.flowequation.bordermacayeal=nodeonmacayeal;
+		md.flowequation.element_equation(find(~stokesflag))=2;
+	else %fill with none 
+	%5: None elements (non Stokes)
+		md.flowequation.element_equation(find(~stokesflag))=0;
+	end
+end
+
+%Now take care of the coupling between MacAyeal and Pattyn
+md.penalties=[];
+nodeonmacayealpattyn=zeros(md.numberofnodes,1);
+nodeonpattynstokes=zeros(md.numberofnodes,1);
+nodeonmacayealstokes=zeros(md.numberofnodes,1);
+if strcmpi(coupling_method,'penalties'),
+	%Create the border nodes between Pattyn and MacAyeal and extrude them
+	numnodes2d=md.numberofnodes2d;
+	numlayers=md.numlayers;
+	bordernodes2d=find(nodeonpattyn(1:numnodes2d) & nodeonmacayeal(1:numnodes2d)); %Nodes connected to two different types of elements
+
+	%initialize and fill in penalties structure
+	if ~isnan(bordernodes2d),
+		penalties=[];
+		for	i=1:numlayers-1,
+			penalties=[penalties; [bordernodes2d bordernodes2d+md.numberofnodes2d*(i)]];
+		end
+		md.penalties=penalties;
+	end
+elseif strcmpi(coupling_method,'tiling'),
+	if any(macayealflag) & any(pattynflag), %coupling macayeal pattyn
+		%Find node at the border
+		nodeonmacayealpattyn(find(nodeonmacayeal & nodeonpattyn))=1;
+		%Macayeal elements in contact with this layer become MacAyealPattyn elements
+		matrixelements=ismember(md.elements,find(nodeonmacayealpattyn));
+		commonelements=sum(matrixelements,2)~=0;
+		commonelements(find(pattynflag))=0; %only one layer: the elements previously in macayeal
+		macayealflag(find(commonelements))=0; %these elements are now macayealpattynelements
+		macayealpattynflag=zeros(md.numberofelements,1);
+		macayealpattynflag(find(commonelements))=1;
+		nodeonmacayeal=zeros(md.numberofnodes,1);
+		nodeonmacayeal(md.elements(find(macayealflag),:))=1;
+		md.flowequation.bordermacayeal=nodeonmacayeal;
+
+		%Create MacaAyealPattynApproximation where needed
+		md.flowequation.element_equation(find(macayealpattynflag))=5;
+
+		%Now recreate nodeonmacayeal and nodeonmacayealpattyn
+		nodeonmacayealpattyn(md.elements(find(macayealpattynflag),:))=1;
+	elseif any(pattynflag) & any(stokesflag), %coupling pattyn stokes
+		%Find node at the border
+		nodeonpattynstokes(find(nodeonpattyn & nodeonstokes))=1;
+		%Stokes elements in contact with this layer become PattynStokes elements
+		matrixelements=ismember(md.elements,find(nodeonpattynstokes));
+		commonelements=sum(matrixelements,2)~=0;
+		commonelements(find(pattynflag))=0; %only one layer: the elements previously in macayeal
+		stokesflag(find(commonelements))=0; %these elements are now macayealpattynelements
+		pattynstokesflag=zeros(md.numberofelements,1);
+		pattynstokesflag(find(commonelements))=1;
+		nodeonstokes=zeros(md.numberofnodes,1);
+		nodeonstokes(md.elements(find(stokesflag),:))=1;
+		md.flowequation.borderstokes=nodeonstokes;
+
+		%Create MacaAyealPattynApproximation where needed
+		md.flowequation.element_equation(find(pattynstokesflag))=7;
+
+		%Now recreate nodeonpattynstokes
+		nodeonpattynstokes=zeros(md.numberofnodes,1);
+		nodeonpattynstokes(md.elements(find(pattynstokesflag),:))=1;
+	elseif any(stokesflag) & any(macayealflag),
+		%Find node at the border
+		nodeonmacayealstokes(find(nodeonmacayeal & nodeonstokes))=1;
+		%Stokes elements in contact with this layer become MacAyealStokes elements
+		matrixelements=ismember(md.elements,find(nodeonmacayealstokes));
+		commonelements=sum(matrixelements,2)~=0;
+		commonelements(find(macayealflag))=0; %only one layer: the elements previously in macayeal
+		stokesflag(find(commonelements))=0; %these elements are now macayealmacayealelements
+		macayealstokesflag=zeros(md.numberofelements,1);
+		macayealstokesflag(find(commonelements))=1;
+		nodeonstokes=zeros(md.numberofnodes,1);
+		nodeonstokes(md.elements(find(stokesflag),:))=1;
+		md.flowequation.borderstokes=nodeonstokes;
+
+		%Create MacaAyeal Approximation where needed
+		md.flowequation.element_equation(find(macayealstokesflag))=6;
+
+		%Now recreate nodeonmacayealstokes
+		nodeonmacayealstokes=zeros(md.numberofnodes,1);
+		nodeonmacayealstokes(md.elements(find(macayealstokesflag),:))=1;
+	elseif any(stokesflag) & any(hutterflag),
+		error('type of coupling not supported yet');
+	end
+end
+
+%Create vertices_type
+md.flowequation.vertex_equation=zeros(md.numberofnodes,1);
+pos=find(nodeonhutter);
+md.flowequation.vertex_equation(pos)=1;
+pos=find(nodeonmacayeal);
+md.flowequation.vertex_equation(pos)=2;
+pos=find(nodeonpattyn);
+md.flowequation.vertex_equation(pos)=3;
+pos=find(nodeonhutter);
+md.flowequation.vertex_equation(pos)=1;
+pos=find(nodeonpattyn & nodeonmacayeal);
+md.flowequation.vertex_equation(pos)=3;
+pos=find(nodeonmacayealpattyn);
+md.flowequation.vertex_equation(pos)=5;
+pos=find(nodeonstokes);
+md.flowequation.vertex_equation(pos)=4;
+if any(stokesflag),
+	pos=find(~nodeonstokes);
+	if(~any(pattynflag) & ~any(macayealflag)),
+		md.flowequation.vertex_equation(pos)=0;
+	end
+end
+pos=find(nodeonpattynstokes);
+md.flowequation.vertex_equation(pos)=7;
+pos=find(nodeonmacayealstokes);
+md.flowequation.vertex_equation(pos)=6;
+
+%figure out solution types
+md.flowequation.ishutter=double(any(md.flowequation.element_equation==1));
+md.flowequation.ismacayealpattyn=double(any(md.flowequation.element_equation==2 | md.flowequation.element_equation==3));
+md.flowequation.isstokes=double(any(md.flowequation.element_equation==4));
+
+end
Index: /issm/trunk/test/Miscellaneous/Bump/runme.m
===================================================================
--- /issm/trunk/test/Miscellaneous/Bump/runme.m	(revision 9663)
+++ /issm/trunk/test/Miscellaneous/Bump/runme.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'Bump.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=astrid;
 md.cluster.np=13;
Index: /issm/trunk/test/Miscellaneous/GJM_test1/test1.m
===================================================================
--- /issm/trunk/test/Miscellaneous/GJM_test1/test1.m	(revision 9663)
+++ /issm/trunk/test/Miscellaneous/GJM_test1/test1.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/Miscellaneous/connectivity/runme.m
===================================================================
--- /issm/trunk/test/Miscellaneous/connectivity/runme.m	(revision 9663)
+++ /issm/trunk/test/Miscellaneous/connectivity/runme.m	(revision 9664)
@@ -40,7 +40,7 @@
 		if is3d,
 			md=extrude(md,4,5);
-			md=setelementstype(md,'pattyn','all');
+			md=setflowequation(md,'pattyn','all');
 		else
-			md=setelementstype(md,'macayeal','all');
+			md=setflowequation(md,'macayeal','all');
 		end
 
Index: /issm/trunk/test/NightlyRun/test101.m
===================================================================
--- /issm/trunk/test/NightlyRun/test101.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test101.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test102.m
===================================================================
--- /issm/trunk/test/NightlyRun/test102.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test102.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test103.m
===================================================================
--- /issm/trunk/test/NightlyRun/test103.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test103.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test104.m
===================================================================
--- /issm/trunk/test/NightlyRun/test104.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test104.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test105.m
===================================================================
--- /issm/trunk/test/NightlyRun/test105.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test105.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test106.m
===================================================================
--- /issm/trunk/test/NightlyRun/test106.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test106.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test107.m
===================================================================
--- /issm/trunk/test/NightlyRun/test107.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test107.m	(revision 9664)
@@ -4,5 +4,5 @@
 md=extrude(md,3,2);
 md.cluster=none;
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md=solve(md,DiagnosticSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test108.m
===================================================================
--- /issm/trunk/test/NightlyRun/test108.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test108.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test109.m
===================================================================
--- /issm/trunk/test/NightlyRun/test109.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test109.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test110.m
===================================================================
--- /issm/trunk/test/NightlyRun/test110.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test110.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test1101.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1101.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1101.m	(revision 9664)
@@ -18,5 +18,5 @@
 	md=extrude(md,9,1);
 
-	md=setelementstype(md,'pattyn','all');
+	md=setflowequation(md,'pattyn','all');
 
 	%Create dirichlet on the bed only
Index: /issm/trunk/test/NightlyRun/test1102.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1102.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1102.m	(revision 9664)
@@ -25,5 +25,5 @@
 	md=parameterize(md,'../Par/ISMIPA.par');
 	md=extrude(md,10,1);
-	md=setelementstype(md,'stokes','all');
+	md=setflowequation(md,'stokes','all');
 
 	%Create dirichlet on the bed only
Index: /issm/trunk/test/NightlyRun/test1103.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1103.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1103.m	(revision 9664)
@@ -16,5 +16,5 @@
 	md=extrude(md,10,1);
 
-	md=setelementstype(md,'pattyn','all');
+	md=setflowequation(md,'pattyn','all');
 
 	%Create dirichlet on the bed only
Index: /issm/trunk/test/NightlyRun/test1104.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1104.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1104.m	(revision 9664)
@@ -15,5 +15,5 @@
 	md=extrude(md,10,1);
 
-	md=setelementstype(md,'pattyn','all','stokes','all');
+	md=setflowequation(md,'pattyn','all','stokes','all');
 
 	%Create dirichlet on the bed only
Index: /issm/trunk/test/NightlyRun/test1105.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1105.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1105.m	(revision 9664)
@@ -18,5 +18,5 @@
 	md=extrude(md,10,1);
 
-	md=setelementstype(md,'pattyn','all'); 
+	md=setflowequation(md,'pattyn','all'); 
 
 	%Create MPCs to have periodic boundary conditions
Index: /issm/trunk/test/NightlyRun/test1106.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1106.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1106.m	(revision 9664)
@@ -15,5 +15,5 @@
 	md=extrude(md,6,1);
 
-	md=setelementstype(md,'pattyn','all','stokes','all'); 
+	md=setflowequation(md,'pattyn','all','stokes','all'); 
 
 	%Create MPCs to have periodic boundary conditions
Index: /issm/trunk/test/NightlyRun/test1107.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1107.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1107.m	(revision 9664)
@@ -18,5 +18,5 @@
 	md=extrude(md,10,1);
 
-	md=setelementstype(md,'pattyn','all');
+	md=setflowequation(md,'pattyn','all');
 
 	%We need one grd on dirichlet: the 4 corners are set to zero
Index: /issm/trunk/test/NightlyRun/test1108.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1108.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1108.m	(revision 9664)
@@ -15,5 +15,5 @@
 	md=extrude(md,6,1);
 
-	md=setelementstype(md,'pattyn','all','stokes','all');
+	md=setflowequation(md,'pattyn','all','stokes','all');
 
 	%We need one grd on dirichlet: the 4 corners are set to zero
Index: /issm/trunk/test/NightlyRun/test1109.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1109.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1109.m	(revision 9664)
@@ -20,7 +20,7 @@
 
 	if i==1 | i==3,
-		md=setelementstype(md,'pattyn','all');
+		md=setflowequation(md,'pattyn','all');
 	elseif i==2 | i==4,
-		md=setelementstype(md,'stokes','all');
+		md=setflowequation(md,'stokes','all');
 	end
 
Index: /issm/trunk/test/NightlyRun/test111.m
===================================================================
--- /issm/trunk/test/NightlyRun/test111.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test111.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.prognostic.stabilization=3;
Index: /issm/trunk/test/NightlyRun/test1110.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1110.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1110.m	(revision 9664)
@@ -15,7 +15,7 @@
 
 	if (i==1 | i==2),
-		md=setelementstype(md,'pattyn','all');
+		md=setflowequation(md,'pattyn','all');
 	else
-		md=setelementstype(md,'stokes','all');
+		md=setflowequation(md,'stokes','all');
 	end
 
Index: /issm/trunk/test/NightlyRun/test112.m
===================================================================
--- /issm/trunk/test/NightlyRun/test112.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test112.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.prognostic.stabilization=3;
Index: /issm/trunk/test/NightlyRun/test113.m
===================================================================
--- /issm/trunk/test/NightlyRun/test113.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test113.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md=extrude(md,5,3);
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test114.m
===================================================================
--- /issm/trunk/test/NightlyRun/test114.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test114.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md=extrude(md,5,3);
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test115.m
===================================================================
--- /issm/trunk/test/NightlyRun/test115.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test115.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.timestepping.time_step=0;
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test116.m
===================================================================
--- /issm/trunk/test/NightlyRun/test116.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test116.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.timestepping.time_step=0;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test117.m
===================================================================
--- /issm/trunk/test/NightlyRun/test117.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test117.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test118.m
===================================================================
--- /issm/trunk/test/NightlyRun/test118.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test118.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test119.m
===================================================================
--- /issm/trunk/test/NightlyRun/test119.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test119.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test120.m
===================================================================
--- /issm/trunk/test/NightlyRun/test120.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test120.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test1201.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1201.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1201.m	(revision 9664)
@@ -10,5 +10,5 @@
 	md=parameterize(md,'../Par/SquareEISMINT.par');
 	md.surfaceforcings.mass_balance(:)=0;
-	md=setelementstype(md,'macayeal','all');
+	md=setflowequation(md,'macayeal','all');
 	md.cluster=generic('name',oshostname(),'np',8);
 
Index: /issm/trunk/test/NightlyRun/test1202.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1202.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1202.m	(revision 9664)
@@ -7,5 +7,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareEISMINT.par');
-md=setelementstype(md,'macayeal','all'); %MacAyeal's model and 2d
+md=setflowequation(md,'macayeal','all'); %MacAyeal's model and 2d
 
 %Compute solution for MacAyeal's model 
Index: /issm/trunk/test/NightlyRun/test1203.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1203.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1203.m	(revision 9664)
@@ -7,5 +7,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareEISMINT.par');
-md=setelementstype(md,'macayeal','all'); %MacAyeal's model and 2d
+md=setflowequation(md,'macayeal','all'); %MacAyeal's model and 2d
 
 %Impose a non zero velocity on the upper boundary condition (y=max(y))
Index: /issm/trunk/test/NightlyRun/test1204.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1204.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1204.m	(revision 9664)
@@ -7,5 +7,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareEISMINT.par');
-md=setelementstype(md,'macayeal','all'); %MacAyeal's model and 2d
+md=setflowequation(md,'macayeal','all'); %MacAyeal's model and 2d
 
 %Impose a non zero velocity on the upper boundary condition (y=max(y))
Index: /issm/trunk/test/NightlyRun/test1205.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1205.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1205.m	(revision 9664)
@@ -19,5 +19,5 @@
 %We extrude the model to have a 3d model
 md=extrude(md,numlayers,1);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 
 %Spc the nodes on the bed
Index: /issm/trunk/test/NightlyRun/test1206.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1206.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1206.m	(revision 9664)
@@ -19,5 +19,5 @@
 %We extrude the model to have a 3d model
 md=extrude(md,numlayers,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %Spc the nodes on the bed
Index: /issm/trunk/test/NightlyRun/test1207.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1207.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1207.m	(revision 9664)
@@ -19,5 +19,5 @@
 %We extrude the model to have a 3d model
 md=extrude(md,numlayers,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 
 %Spc the nodes on the bed
Index: /issm/trunk/test/NightlyRun/test1208.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1208.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1208.m	(revision 9664)
@@ -10,5 +10,5 @@
 %We extrude the model to have a 3d model
 md=extrude(md,numlayers,1);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 
 %Spc the nodes on the bed
Index: /issm/trunk/test/NightlyRun/test121.m
===================================================================
--- /issm/trunk/test/NightlyRun/test121.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test121.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test122.m
===================================================================
--- /issm/trunk/test/NightlyRun/test122.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test122.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test123.m
===================================================================
--- /issm/trunk/test/NightlyRun/test123.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test123.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test124.m
===================================================================
--- /issm/trunk/test/NightlyRun/test124.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test124.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test125.m
===================================================================
--- /issm/trunk/test/NightlyRun/test125.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test125.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test126.m
===================================================================
--- /issm/trunk/test/NightlyRun/test126.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test126.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test127.m
===================================================================
--- /issm/trunk/test/NightlyRun/test127.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test127.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test128.m
===================================================================
--- /issm/trunk/test/NightlyRun/test128.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test128.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test129.m
===================================================================
--- /issm/trunk/test/NightlyRun/test129.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test129.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test130.m
===================================================================
--- /issm/trunk/test/NightlyRun/test130.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test130.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test1301.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1301.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1301.m	(revision 9664)
@@ -10,5 +10,5 @@
 md=parameterize(md,'../Par/SquareThermal.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'Pattyn','all');
+md=setflowequation(md,'Pattyn','all');
 
 %Some conditions specific to melting test
Index: /issm/trunk/test/NightlyRun/test1302.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1302.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1302.m	(revision 9664)
@@ -9,5 +9,5 @@
 md=parameterize(md,'../Par/SquareThermal.par');
 md=extrude(md,30,1);   %NB: the more one extrudes, the better (10-> relative~0.35%, 20->0.1%, 30->0.05%)
-md=setelementstype(md,'Pattyn','all'); 
+md=setflowequation(md,'Pattyn','all'); 
 
 %Thermal boundary conditions
Index: /issm/trunk/test/NightlyRun/test1303.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1303.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1303.m	(revision 9664)
@@ -10,5 +10,5 @@
 md=parameterize(md,'../Par/SquareThermal.par');
 md=extrude(md,11,2);
-md=setelementstype(md,'Pattyn','all');
+md=setflowequation(md,'Pattyn','all');
 pos1=find(md.elementonbed);     md.thermal.spctemperature(md.elements(pos1,1:3))=10;
 pos2=find(md.elementonsurface); md.thermal.spctemperature(md.elements(pos2,4:6))=0;
Index: /issm/trunk/test/NightlyRun/test1304.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1304.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1304.m	(revision 9664)
@@ -10,5 +10,5 @@
 md=parameterize(md,'../Par/SquareThermal.par');
 md=extrude(md,11,1);
-md=setelementstype(md,'Pattyn','all');
+md=setflowequation(md,'Pattyn','all');
 
 pos2=find(md.elementonsurface); md.thermal.spctemperature(md.elements(pos2,4:6))=0;
Index: /issm/trunk/test/NightlyRun/test131.m
===================================================================
--- /issm/trunk/test/NightlyRun/test131.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test131.m	(revision 9664)
@@ -5,5 +5,5 @@
 pos=find(md.nodeonboundary);
 md.balancethickness.spcthickness(pos)=md.thickness(pos);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test132.m
===================================================================
--- /issm/trunk/test/NightlyRun/test132.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test132.m	(revision 9664)
@@ -5,5 +5,5 @@
 pos=find(md.nodeonboundary);
 md.balancethickness.spcthickness(pos)=md.thickness(pos);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test133.m
===================================================================
--- /issm/trunk/test/NightlyRun/test133.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test133.m	(revision 9664)
@@ -6,5 +6,5 @@
 pos=find(md.nodeonboundary);
 md.balancethickness.spcthickness(pos)=md.thickness(pos);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test134.m
===================================================================
--- /issm/trunk/test/NightlyRun/test134.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test134.m	(revision 9664)
@@ -6,5 +6,5 @@
 pos=find(md.nodeonboundary);
 md.balancethickness.spcthickness(pos)=md.thickness(pos);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test135.m
===================================================================
--- /issm/trunk/test/NightlyRun/test135.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test135.m	(revision 9664)
@@ -5,5 +5,5 @@
 md.prognostic.stabilization=3;
 md.vy=md.vy+400;
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test136.m
===================================================================
--- /issm/trunk/test/NightlyRun/test136.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test136.m	(revision 9664)
@@ -5,5 +5,5 @@
 md.prognostic.stabilization=3;
 md.vy=md.vy+400;
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test1501.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1501.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1501.m	(revision 9664)
@@ -4,5 +4,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.isthermal=0;
Index: /issm/trunk/test/NightlyRun/test1502.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1502.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test1502.m	(revision 9664)
@@ -4,5 +4,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md=extrude(md,3,1);
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test201.m
===================================================================
--- /issm/trunk/test/NightlyRun/test201.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test201.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test202.m
===================================================================
--- /issm/trunk/test/NightlyRun/test202.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test202.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test203.m
===================================================================
--- /issm/trunk/test/NightlyRun/test203.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test203.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test204.m
===================================================================
--- /issm/trunk/test/NightlyRun/test204.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test204.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test205.m
===================================================================
--- /issm/trunk/test/NightlyRun/test205.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test205.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test206.m
===================================================================
--- /issm/trunk/test/NightlyRun/test206.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test206.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test207.m
===================================================================
--- /issm/trunk/test/NightlyRun/test207.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test207.m	(revision 9664)
@@ -4,5 +4,5 @@
 md=extrude(md,3,2);
 md.cluster=none;
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md=solve(md,DiagnosticSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test208.m
===================================================================
--- /issm/trunk/test/NightlyRun/test208.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test208.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test209.m
===================================================================
--- /issm/trunk/test/NightlyRun/test209.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test209.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal','coupling','penalties');
+md=setflowequation(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal','coupling','penalties');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test210.m
===================================================================
--- /issm/trunk/test/NightlyRun/test210.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test210.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal','coupling','penalties');
+md=setflowequation(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal','coupling','penalties');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test211.m
===================================================================
--- /issm/trunk/test/NightlyRun/test211.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test211.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.timestepping.time_step=0;
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test212.m
===================================================================
--- /issm/trunk/test/NightlyRun/test212.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test212.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.timestepping.time_step=0;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test213.m
===================================================================
--- /issm/trunk/test/NightlyRun/test213.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test213.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test214.m
===================================================================
--- /issm/trunk/test/NightlyRun/test214.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test214.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test215.m
===================================================================
--- /issm/trunk/test/NightlyRun/test215.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test215.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test216.m
===================================================================
--- /issm/trunk/test/NightlyRun/test216.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test216.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test217.m
===================================================================
--- /issm/trunk/test/NightlyRun/test217.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test217.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test218.m
===================================================================
--- /issm/trunk/test/NightlyRun/test218.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test218.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test219.m
===================================================================
--- /issm/trunk/test/NightlyRun/test219.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test219.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test220.m
===================================================================
--- /issm/trunk/test/NightlyRun/test220.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test220.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test221.m
===================================================================
--- /issm/trunk/test/NightlyRun/test221.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test221.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.eps_rel=NaN;
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test222.m
===================================================================
--- /issm/trunk/test/NightlyRun/test222.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test222.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.eps_rel=NaN;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test223.m
===================================================================
--- /issm/trunk/test/NightlyRun/test223.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test223.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test224.m
===================================================================
--- /issm/trunk/test/NightlyRun/test224.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test224.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test225.m
===================================================================
--- /issm/trunk/test/NightlyRun/test225.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test225.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test226.m
===================================================================
--- /issm/trunk/test/NightlyRun/test226.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test226.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test227.m
===================================================================
--- /issm/trunk/test/NightlyRun/test227.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test227.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test228.m
===================================================================
--- /issm/trunk/test/NightlyRun/test228.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test228.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test229.m
===================================================================
--- /issm/trunk/test/NightlyRun/test229.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test229.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test230.m
===================================================================
--- /issm/trunk/test/NightlyRun/test230.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test230.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test232.m
===================================================================
--- /issm/trunk/test/NightlyRun/test232.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test232.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test233.m
===================================================================
--- /issm/trunk/test/NightlyRun/test233.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test233.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test234.m
===================================================================
--- /issm/trunk/test/NightlyRun/test234.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test234.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test235.m
===================================================================
--- /issm/trunk/test/NightlyRun/test235.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test235.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test236.m
===================================================================
--- /issm/trunk/test/NightlyRun/test236.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test236.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test237.m
===================================================================
--- /issm/trunk/test/NightlyRun/test237.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test237.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test238.m
===================================================================
--- /issm/trunk/test/NightlyRun/test238.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test238.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'pattyn','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test239.m
===================================================================
--- /issm/trunk/test/NightlyRun/test239.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test239.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn');
 md=solve(md,DiagnosticSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test240.m
===================================================================
--- /issm/trunk/test/NightlyRun/test240.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test240.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test241.m
===================================================================
--- /issm/trunk/test/NightlyRun/test241.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test241.m	(revision 9664)
@@ -4,5 +4,5 @@
 md=extrude(md,2,1);
 md.viscosity_overshoot=0.0;
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal');
 md=solve(md,DiagnosticSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test242.m
===================================================================
--- /issm/trunk/test/NightlyRun/test242.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test242.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=generic('name',oshostname(),'np',3);
 md.viscosity_overshoot=0;
Index: /issm/trunk/test/NightlyRun/test243.m
===================================================================
--- /issm/trunk/test/NightlyRun/test243.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test243.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.prognostic.hydrostatic_adjustment='Incremental';
Index: /issm/trunk/test/NightlyRun/test244.m
===================================================================
--- /issm/trunk/test/NightlyRun/test244.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test244.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.prognostic.hydrostatic_adjustment='Incremental';
Index: /issm/trunk/test/NightlyRun/test245.m
===================================================================
--- /issm/trunk/test/NightlyRun/test245.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test245.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.bed=md.bed+50; md.surface=md.surface+50;
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test246.m
===================================================================
--- /issm/trunk/test/NightlyRun/test246.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test246.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.bed=md.bed+50; md.surface=md.surface+50;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test247.m
===================================================================
--- /issm/trunk/test/NightlyRun/test247.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test247.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md.prognostic.hydrostatic_adjustment='Incremental';
Index: /issm/trunk/test/NightlyRun/test248.m
===================================================================
--- /issm/trunk/test/NightlyRun/test248.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test248.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.prognostic.hydrostatic_adjustment='Incremental';
Index: /issm/trunk/test/NightlyRun/test249.m
===================================================================
--- /issm/trunk/test/NightlyRun/test249.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test249.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.bed=md.bed+50; md.surface=md.surface+50;
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test250.m
===================================================================
--- /issm/trunk/test/NightlyRun/test250.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test250.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.bed=md.bed+50; md.surface=md.surface+50;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test251.m
===================================================================
--- /issm/trunk/test/NightlyRun/test251.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test251.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test252.m
===================================================================
--- /issm/trunk/test/NightlyRun/test252.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test252.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test253.m
===================================================================
--- /issm/trunk/test/NightlyRun/test253.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test253.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test254.m
===================================================================
--- /issm/trunk/test/NightlyRun/test254.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test254.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test255.m
===================================================================
--- /issm/trunk/test/NightlyRun/test255.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test255.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test256.m
===================================================================
--- /issm/trunk/test/NightlyRun/test256.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test256.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test257.m
===================================================================
--- /issm/trunk/test/NightlyRun/test257.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test257.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test258.m
===================================================================
--- /issm/trunk/test/NightlyRun/test258.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test258.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test259.m
===================================================================
--- /issm/trunk/test/NightlyRun/test259.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test259.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test260.m
===================================================================
--- /issm/trunk/test/NightlyRun/test260.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test260.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test261.m
===================================================================
--- /issm/trunk/test/NightlyRun/test261.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test261.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test262.m
===================================================================
--- /issm/trunk/test/NightlyRun/test262.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test262.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test263.m
===================================================================
--- /issm/trunk/test/NightlyRun/test263.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test263.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5, md.thermal.spctemperature+10, md.thermal.spctemperature+15; 1.5 2.5 3.5 4];
Index: /issm/trunk/test/NightlyRun/test264.m
===================================================================
--- /issm/trunk/test/NightlyRun/test264.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test264.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5, md.thermal.spctemperature+10, md.thermal.spctemperature+15; 1.5 2.5 3.5 4];
Index: /issm/trunk/test/NightlyRun/test265.m
===================================================================
--- /issm/trunk/test/NightlyRun/test265.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test265.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5; 1 2];
Index: /issm/trunk/test/NightlyRun/test266.m
===================================================================
--- /issm/trunk/test/NightlyRun/test266.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test266.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.thermal.spctemperature=[md.thermal.spctemperature, md.thermal.spctemperature+5; 1 2];
Index: /issm/trunk/test/NightlyRun/test301.m
===================================================================
--- /issm/trunk/test/NightlyRun/test301.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test301.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test302.m
===================================================================
--- /issm/trunk/test/NightlyRun/test302.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test302.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test303.m
===================================================================
--- /issm/trunk/test/NightlyRun/test303.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test303.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test304.m
===================================================================
--- /issm/trunk/test/NightlyRun/test304.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test304.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test305.m
===================================================================
--- /issm/trunk/test/NightlyRun/test305.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test305.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test306.m
===================================================================
--- /issm/trunk/test/NightlyRun/test306.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test306.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test307.m
===================================================================
--- /issm/trunk/test/NightlyRun/test307.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test307.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test308.m
===================================================================
--- /issm/trunk/test/NightlyRun/test308.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test308.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test309.m
===================================================================
--- /issm/trunk/test/NightlyRun/test309.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test309.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,5,2);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test310.m
===================================================================
--- /issm/trunk/test/NightlyRun/test310.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test310.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,5,2);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test311.m
===================================================================
--- /issm/trunk/test/NightlyRun/test311.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test311.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,4,2);
-md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test312.m
===================================================================
--- /issm/trunk/test/NightlyRun/test312.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test312.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,4,2);
-md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test313.m
===================================================================
--- /issm/trunk/test/NightlyRun/test313.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test313.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test314.m
===================================================================
--- /issm/trunk/test/NightlyRun/test314.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test314.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test315.m
===================================================================
--- /issm/trunk/test/NightlyRun/test315.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test315.m	(revision 9664)
@@ -4,5 +4,5 @@
 md=extrude(md,3,1);
 md.cluster=none;
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md=solve(md,DiagnosticSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test316.m
===================================================================
--- /issm/trunk/test/NightlyRun/test316.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test316.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test317.m
===================================================================
--- /issm/trunk/test/NightlyRun/test317.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test317.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test318.m
===================================================================
--- /issm/trunk/test/NightlyRun/test318.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test318.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test319.m
===================================================================
--- /issm/trunk/test/NightlyRun/test319.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test319.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.prognostic.stabilization=3;
 md.prognostic.spcthickness=md.thickness;
Index: /issm/trunk/test/NightlyRun/test320.m
===================================================================
--- /issm/trunk/test/NightlyRun/test320.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test320.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.prognostic.stabilization=3;
 md.prognostic.spcthickness=md.thickness;
Index: /issm/trunk/test/NightlyRun/test321.m
===================================================================
--- /issm/trunk/test/NightlyRun/test321.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test321.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md=extrude(md,5,0.5);
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test322.m
===================================================================
--- /issm/trunk/test/NightlyRun/test322.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test322.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md=extrude(md,5,0.5);
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test323.m
===================================================================
--- /issm/trunk/test/NightlyRun/test323.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test323.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.timestepping.time_step=0;
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test324.m
===================================================================
--- /issm/trunk/test/NightlyRun/test324.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test324.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.timestepping.time_step=0;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test325.m
===================================================================
--- /issm/trunk/test/NightlyRun/test325.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test325.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test326.m
===================================================================
--- /issm/trunk/test/NightlyRun/test326.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test326.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test327.m
===================================================================
--- /issm/trunk/test/NightlyRun/test327.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test327.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test328.m
===================================================================
--- /issm/trunk/test/NightlyRun/test328.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test328.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test329.m
===================================================================
--- /issm/trunk/test/NightlyRun/test329.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test329.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,5,1.2);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test330.m
===================================================================
--- /issm/trunk/test/NightlyRun/test330.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test330.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,5,1.2);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test331.m
===================================================================
--- /issm/trunk/test/NightlyRun/test331.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test331.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test332.m
===================================================================
--- /issm/trunk/test/NightlyRun/test332.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test332.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test333.m
===================================================================
--- /issm/trunk/test/NightlyRun/test333.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test333.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test334.m
===================================================================
--- /issm/trunk/test/NightlyRun/test334.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test334.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test335.m
===================================================================
--- /issm/trunk/test/NightlyRun/test335.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test335.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=none;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test336.m
===================================================================
--- /issm/trunk/test/NightlyRun/test336.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test336.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test337.m
===================================================================
--- /issm/trunk/test/NightlyRun/test337.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test337.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test338.m
===================================================================
--- /issm/trunk/test/NightlyRun/test338.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test338.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test339.m
===================================================================
--- /issm/trunk/test/NightlyRun/test339.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test339.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test340.m
===================================================================
--- /issm/trunk/test/NightlyRun/test340.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test340.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test341.m
===================================================================
--- /issm/trunk/test/NightlyRun/test341.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test341.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test342.m
===================================================================
--- /issm/trunk/test/NightlyRun/test342.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test342.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test343.m
===================================================================
--- /issm/trunk/test/NightlyRun/test343.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test343.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test344.m
===================================================================
--- /issm/trunk/test/NightlyRun/test344.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test344.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test345.m
===================================================================
--- /issm/trunk/test/NightlyRun/test345.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test345.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test346.m
===================================================================
--- /issm/trunk/test/NightlyRun/test346.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test346.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test347.m
===================================================================
--- /issm/trunk/test/NightlyRun/test347.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test347.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,4,1.2);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=none;
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test348.m
===================================================================
--- /issm/trunk/test/NightlyRun/test348.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test348.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=extrude(md,5,1.2);
-md=setelementstype(md,'hutter','all');
+md=setflowequation(md,'hutter','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_adapt=1;
Index: /issm/trunk/test/NightlyRun/test401.m
===================================================================
--- /issm/trunk/test/NightlyRun/test401.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test401.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test402.m
===================================================================
--- /issm/trunk/test/NightlyRun/test402.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test402.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test403.m
===================================================================
--- /issm/trunk/test/NightlyRun/test403.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test403.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test404.m
===================================================================
--- /issm/trunk/test/NightlyRun/test404.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test404.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test405.m
===================================================================
--- /issm/trunk/test/NightlyRun/test405.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test405.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test406.m
===================================================================
--- /issm/trunk/test/NightlyRun/test406.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test406.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test407.m
===================================================================
--- /issm/trunk/test/NightlyRun/test407.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test407.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md=solve(md,DiagnosticSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test408.m
===================================================================
--- /issm/trunk/test/NightlyRun/test408.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test408.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test409.m
===================================================================
--- /issm/trunk/test/NightlyRun/test409.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test409.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
+md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test410.m
===================================================================
--- /issm/trunk/test/NightlyRun/test410.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test410.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
+md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test411.m
===================================================================
--- /issm/trunk/test/NightlyRun/test411.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test411.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test412.m
===================================================================
--- /issm/trunk/test/NightlyRun/test412.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test412.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test413.m
===================================================================
--- /issm/trunk/test/NightlyRun/test413.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test413.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test414.m
===================================================================
--- /issm/trunk/test/NightlyRun/test414.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test414.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,4,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test415.m
===================================================================
--- /issm/trunk/test/NightlyRun/test415.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test415.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test416.m
===================================================================
--- /issm/trunk/test/NightlyRun/test416.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test416.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test417.m
===================================================================
--- /issm/trunk/test/NightlyRun/test417.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test417.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
+md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test418.m
===================================================================
--- /issm/trunk/test/NightlyRun/test418.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test418.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
+md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test419.m
===================================================================
--- /issm/trunk/test/NightlyRun/test419.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test419.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test420.m
===================================================================
--- /issm/trunk/test/NightlyRun/test420.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test420.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test421.m
===================================================================
--- /issm/trunk/test/NightlyRun/test421.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test421.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test422.m
===================================================================
--- /issm/trunk/test/NightlyRun/test422.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test422.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,2);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test423.m
===================================================================
--- /issm/trunk/test/NightlyRun/test423.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test423.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test424.m
===================================================================
--- /issm/trunk/test/NightlyRun/test424.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test424.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test425.m
===================================================================
--- /issm/trunk/test/NightlyRun/test425.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test425.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test426.m
===================================================================
--- /issm/trunk/test/NightlyRun/test426.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test426.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test427.m
===================================================================
--- /issm/trunk/test/NightlyRun/test427.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test427.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.materials.rho_ice=10^7; %involved in the mass flux, make it easy
Index: /issm/trunk/test/NightlyRun/test428.m
===================================================================
--- /issm/trunk/test/NightlyRun/test428.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test428.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.materials.rho_ice=10^7; %involved in the mass flux, make it easy
Index: /issm/trunk/test/NightlyRun/test429.m
===================================================================
--- /issm/trunk/test/NightlyRun/test429.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test429.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test430.m
===================================================================
--- /issm/trunk/test/NightlyRun/test430.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test430.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test431.m
===================================================================
--- /issm/trunk/test/NightlyRun/test431.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test431.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test432.m
===================================================================
--- /issm/trunk/test/NightlyRun/test432.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test432.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test434.m
===================================================================
--- /issm/trunk/test/NightlyRun/test434.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test434.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.materials.rho_ice=10^7; %involved in the mass flux, make it easy
Index: /issm/trunk/test/NightlyRun/test435.m
===================================================================
--- /issm/trunk/test/NightlyRun/test435.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test435.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test437.m
===================================================================
--- /issm/trunk/test/NightlyRun/test437.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test437.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','tiling');
+md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','tiling');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test438.m
===================================================================
--- /issm/trunk/test/NightlyRun/test438.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test438.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn');
+md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test439.m
===================================================================
--- /issm/trunk/test/NightlyRun/test439.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test439.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 
Index: /issm/trunk/test/NightlyRun/test440.m
===================================================================
--- /issm/trunk/test/NightlyRun/test440.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test440.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/SquareShelf.exp','');
 md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk/test/NightlyRun/test441.m
===================================================================
--- /issm/trunk/test/NightlyRun/test441.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test441.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','tiling');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','tiling');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test442.m
===================================================================
--- /issm/trunk/test/NightlyRun/test442.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test442.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn');
 md.cluster=generic('name',oshostname(),'np',7);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test443.m
===================================================================
--- /issm/trunk/test/NightlyRun/test443.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test443.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal','coupling','tiling');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal','coupling','tiling');
 md.cluster=none;
 md.eps_rel=0.4;
Index: /issm/trunk/test/NightlyRun/test444.m
===================================================================
--- /issm/trunk/test/NightlyRun/test444.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test444.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/SquareSheetShelf.par');
 md=extrude(md,5,1);
-md=setelementstype(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal');
+md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','macayeal');
 md.cluster=generic('name',oshostname(),'np',3);
 md.eps_rel=0.4;
Index: /issm/trunk/test/NightlyRun/test446.m
===================================================================
--- /issm/trunk/test/NightlyRun/test446.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test446.m	(revision 9664)
@@ -20,5 +20,5 @@
 %}}}
 md=parameterize(md,'../Par/RoundSheetShelf.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 %plug holes into the ice sheet, to test for grounding line migration. {{{1
Index: /issm/trunk/test/NightlyRun/test501.m
===================================================================
--- /issm/trunk/test/NightlyRun/test501.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test501.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
 md=parameterize(md,'../Par/Pig.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test502.m
===================================================================
--- /issm/trunk/test/NightlyRun/test502.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test502.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
 md=parameterize(md,'../Par/Pig.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test503.m
===================================================================
--- /issm/trunk/test/NightlyRun/test503.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test503.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,0.9);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test504.m
===================================================================
--- /issm/trunk/test/NightlyRun/test504.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test504.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,0.9);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test505.m
===================================================================
--- /issm/trunk/test/NightlyRun/test505.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test505.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,0.9);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md=solve(md,DiagnosticSolutionEnum);
 
Index: /issm/trunk/test/NightlyRun/test506.m
===================================================================
--- /issm/trunk/test/NightlyRun/test506.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test506.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,0.9);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test507.m
===================================================================
--- /issm/trunk/test/NightlyRun/test507.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test507.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
 md=parameterize(md,'../Par/Pig.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test508.m
===================================================================
--- /issm/trunk/test/NightlyRun/test508.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test508.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
 md=parameterize(md,'../Par/Pig.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test509.m
===================================================================
--- /issm/trunk/test/NightlyRun/test509.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test509.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test510.m
===================================================================
--- /issm/trunk/test/NightlyRun/test510.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test510.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test511.m
===================================================================
--- /issm/trunk/test/NightlyRun/test511.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test511.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test512.m
===================================================================
--- /issm/trunk/test/NightlyRun/test512.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test512.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test513.m
===================================================================
--- /issm/trunk/test/NightlyRun/test513.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test513.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=none;
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test514.m
===================================================================
--- /issm/trunk/test/NightlyRun/test514.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test514.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test515.m
===================================================================
--- /issm/trunk/test/NightlyRun/test515.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test515.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,4,1.1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test516.m
===================================================================
--- /issm/trunk/test/NightlyRun/test516.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test516.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,4,1.1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test517.m
===================================================================
--- /issm/trunk/test/NightlyRun/test517.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test517.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test518.m
===================================================================
--- /issm/trunk/test/NightlyRun/test518.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test518.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test519.m
===================================================================
--- /issm/trunk/test/NightlyRun/test519.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test519.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=none;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test520.m
===================================================================
--- /issm/trunk/test/NightlyRun/test520.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test520.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test521.m
===================================================================
--- /issm/trunk/test/NightlyRun/test521.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test521.m	(revision 9664)
@@ -7,5 +7,5 @@
 md.surface=md.bed+md.thickness;
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md=modelextract(md,md.mask.elementonfloatingice);
 
Index: /issm/trunk/test/NightlyRun/test522.m
===================================================================
--- /issm/trunk/test/NightlyRun/test522.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test522.m	(revision 9664)
@@ -7,5 +7,5 @@
 md.surface=md.bed+md.thickness;
 md=extrude(md,3,1);
-md=setelementstype(md,'stokes','all');
+md=setflowequation(md,'stokes','all');
 md=modelextract(md,md.mask.elementonfloatingice);
 
Index: /issm/trunk/test/NightlyRun/test523.m
===================================================================
--- /issm/trunk/test/NightlyRun/test523.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test523.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test524.m
===================================================================
--- /issm/trunk/test/NightlyRun/test524.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test524.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test525.m
===================================================================
--- /issm/trunk/test/NightlyRun/test525.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test525.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test526.m
===================================================================
--- /issm/trunk/test/NightlyRun/test526.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test526.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test529.m
===================================================================
--- /issm/trunk/test/NightlyRun/test529.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test529.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.thermal.stabilization=2;
 md=solve(md,ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test530.m
===================================================================
--- /issm/trunk/test/NightlyRun/test530.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test530.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.thermal.stabilization=2;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test531.m
===================================================================
--- /issm/trunk/test/NightlyRun/test531.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test531.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.thermal.stabilization=2;
 md.timestepping.time_step=0;
Index: /issm/trunk/test/NightlyRun/test532.m
===================================================================
--- /issm/trunk/test/NightlyRun/test532.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test532.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/Pig.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.thermal.stabilization=2;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test601.m
===================================================================
--- /issm/trunk/test/NightlyRun/test601.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test601.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test602.m
===================================================================
--- /issm/trunk/test/NightlyRun/test602.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test602.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test603.m
===================================================================
--- /issm/trunk/test/NightlyRun/test603.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test603.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.prognostic.stabilization=3;
 md.prognostic.spcthickness=md.thickness;
Index: /issm/trunk/test/NightlyRun/test604.m
===================================================================
--- /issm/trunk/test/NightlyRun/test604.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test604.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.prognostic.stabilization=3;
 md.prognostic.spcthickness=md.thickness;
Index: /issm/trunk/test/NightlyRun/test605.m
===================================================================
--- /issm/trunk/test/NightlyRun/test605.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test605.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,6,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test606.m
===================================================================
--- /issm/trunk/test/NightlyRun/test606.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test606.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,6,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test607.m
===================================================================
--- /issm/trunk/test/NightlyRun/test607.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test607.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test608.m
===================================================================
--- /issm/trunk/test/NightlyRun/test608.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test608.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test609.m
===================================================================
--- /issm/trunk/test/NightlyRun/test609.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test609.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,5,1.5);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=none;
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test610.m
===================================================================
--- /issm/trunk/test/NightlyRun/test610.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test610.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,5,1.5);
-md=setelementstype(md,'pattyn','all');
+md=setflowequation(md,'pattyn','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test611.m
===================================================================
--- /issm/trunk/test/NightlyRun/test611.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test611.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test612.m
===================================================================
--- /issm/trunk/test/NightlyRun/test612.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test612.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test613.m
===================================================================
--- /issm/trunk/test/NightlyRun/test613.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test613.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test614.m
===================================================================
--- /issm/trunk/test/NightlyRun/test614.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test614.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,2,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test615.m
===================================================================
--- /issm/trunk/test/NightlyRun/test615.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test615.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test616.m
===================================================================
--- /issm/trunk/test/NightlyRun/test616.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test616.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test617.m
===================================================================
--- /issm/trunk/test/NightlyRun/test617.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test617.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.prognostic.stabilization=3;
 md.cluster=none;
Index: /issm/trunk/test/NightlyRun/test618.m
===================================================================
--- /issm/trunk/test/NightlyRun/test618.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test618.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.balancethickness.stabilization=3;
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk/test/NightlyRun/test619.m
===================================================================
--- /issm/trunk/test/NightlyRun/test619.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test619.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=none;
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test620.m
===================================================================
--- /issm/trunk/test/NightlyRun/test620.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test620.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=parameterize(md,'../Par/79North.par');
 md=extrude(md,3,1);
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md=solve(md,BalancethicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test621.m
===================================================================
--- /issm/trunk/test/NightlyRun/test621.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test621.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test622.m
===================================================================
--- /issm/trunk/test/NightlyRun/test622.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test622.m	(revision 9664)
@@ -2,5 +2,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test623.m
===================================================================
--- /issm/trunk/test/NightlyRun/test623.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test623.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test624.m
===================================================================
--- /issm/trunk/test/NightlyRun/test624.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test624.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %control parameters
Index: /issm/trunk/test/NightlyRun/test625.m
===================================================================
--- /issm/trunk/test/NightlyRun/test625.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test625.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %Ice sheet only
Index: /issm/trunk/test/NightlyRun/test626.m
===================================================================
--- /issm/trunk/test/NightlyRun/test626.m	(revision 9663)
+++ /issm/trunk/test/NightlyRun/test626.m	(revision 9664)
@@ -3,5 +3,5 @@
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=setelementstype(md,'macayeal','all');
+md=setflowequation(md,'macayeal','all');
 
 %Ice sheet only
