Index: /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/configuration.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/configuration.m	(revision 2147)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/configuration.m	(revision 2147)
@@ -0,0 +1,31 @@
+%name
+testname='IceSheetNoIceFrontM2d_16';
+
+%meshing
+resolution=50000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''macayeal'',''all''';
+
+%extrusion
+numlayers=0;
+extrusionexponent=NaN;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'diagnostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'transient',   'none',      0 ,     0,      'nan'        0    };...
+				 {'diagnostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'transient',   'none',      0 ,     0,      'nan'        1    };...
+	};
Index: sm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/runme.m	(revision 2146)
+++ 	(revision )
@@ -1,103 +1,0 @@
-function varargout=runme(varargin)
-%RUNME - test deck for ISSM nightly runs
-%
-%   Usage:
-%      varargout=runme(varargin);
-%
-%   Examples:
-%      runme;
-%      runme({'ice'});
-%      runme({'ice','cielo_serial'},{'diagnostic'});
-%      md=runme({'cielo_parallel'});
-%
-%   See Also: UPDATEARCHIVE
-
-% This file can be run to check that the current version of issm is giving 
-% coherent results. This test deals with an icesheet without icefront for a 2d model. The geometry 
-% is square. Just run this file in Matlab, with a properly setup ISSM code. 
-% The results of this test will indicate if there is a difference between current computations 
-% and archived results.
-
-% Errors  between archived results and the current version will get flagged if they are not within 
-% a certain tolerance. The current tolerance is 10^-12. If you have good reasons to believe this 
-% tolerance should be lowered (for example, if you are running single precision compilers?), feel 
-% free to tweak the tolerance variable.
-
-%check arguments
-if (nargin>2 | nargout>1)
-	help runme
-	error('runme error message: bad usage');
-end
-
-%packages and analysis_types to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic','prognostic','transient'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'macayeal','ice','cielo_serial'};
-end
-
-%Initialize log message for nightly runs.
-testname='IceSheetNoIceFrontM2d_16';
-
-%go through the solutions requested
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check solution requested
-		if (strcmpi(package,'macayeal') & ~strcmpi(solution,'diagnostic')),
-			disp(sprintf(['\nsolution: ' solution  ', with package: ' package  ', in test: ' testname  ', not supported yet.\n']));
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',50000);
-		md=geography(md,'','');
-		md=parameterize(md,'Square.par');
-		md=setelementstype(md,'macayeal','all');
-		if md.numberofgrids==340
-			load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
-		end
-
-		%compute solution
-		[md analysis_type sub_analysis_type packagefinal]=testsparameterization(md,solution,package);
-		if ~testschecksolution(md.type,solution), continue; end
-		md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',packagefinal);
-
-		%compute fields to be checked
-		fields=testsgetfields(md,solution);
-
-		%load archive
-		eval(['load Archive' package solution ]);
-
-		for k=1:length(fields),
-
-			%Get field and tolerance
-			field=fields{k};
-			tolerance=testsgettolerance(md,package,solution,field);
-
-			%compare to archive
-			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=full(max(abs(Archive-md.results.' field '))/(max(abs(Archive))+eps));']);
-
-			%disp test result
-			if (error_diff>tolerance);
-				disp(sprintf(['\n\nERROR   difference: %-7.2g > %7.2g test: %-25s solution: %-16s package: %-14s field: ' field  '.\n\n'],error_diff,tolerance,testname,solution,package));
-			else
-				disp(sprintf(['\n\nSUCCESS difference: %-7.2g < %7.2g test: %-25s solution: %-16s package: %-14s field: ' field  '.\n\n'],error_diff,tolerance,testname,solution,package));
-			end
-
-		end
-	end
-end
-if nargout==1, 
-	varargout{1}=md; 
-end
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/testpresolve.m	(revision 2147)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/testpresolve.m	(revision 2147)
@@ -0,0 +1,3 @@
+if md.numberofgrids==340
+	load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
+end
Index: sm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/updatearchive.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/updatearchive.m	(revision 2146)
+++ 	(revision )
@@ -1,82 +1,0 @@
-function varargout=updatearchive(varargin)
-%updatearchive - test deck for ISSM nightly runs
-%
-%   Usage:
-%      varargout=updatearchive(varargin);
-%
-%   Examples:
-%      updatearchive;
-%      updatearchive({'ice'});
-%      updatearchive({'ice','cielo_serial'},{'diagnostic'});
-%      md=updatearchive({'cielo_parallel'});
-%
-%   See also: RUNME
-
-% This file can be run to update the velocity archives  of the test1.
-% This test deals with an icesheet without icefront for a 2d model. The geometry 
-% is square. Just run this file in Matlab, with a properly setup Ice code. 
-
-% The archive files will be saved in this directory but will not commited to ice1. 
-% Just commit the result if you want to.
-
-%check arguments
-if (nargin>2 | nargout>1)
-	help updatearchive
-	error('updatearchive error message: bad usage');
-end
-
-%packages and analysis_types to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic','prognostic','transient'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'macayeal','ice','cielo_serial','cielo_parallel'};
-end
-
-%go through the solutions requested
-testname='IceSheetNoIceFrontM2d_16';
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check solution requested
-		if (strcmpi(package,'macayeal') & ~strcmpi(solution,'diagnostic')),
-			disp(sprintf(['\nsolution: ' solution  ', with package: ' package  ', in test: ' testname  ', not supported yet.\n']));
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',50000);
-		md=geography(md,'','');
-		md=parameterize(md,'Square.par');
-		md=setelementstype(md,'macayeal','all');
-		if md.numberofgrids==340
-			load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
-		end
-
-		%compute solution
-		[md analysis_type sub_analysis_type packagefinal]=testsparameterization(md,solution,package);
-		if ~testschecksolution(md.type,solution), continue; end
-		md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',packagefinal);
-
-		%compute fields to be checked
-		fields=testsgetfields(md,solution);
-
-		%save new archive
-		for k=1:length(fields),
-			field=fields{k};
-			eval(['Archive' package solution '_field' num2str(k) '=md.results. ' field  ';']);
-		end
-		eval(['save Archive' package solution ' Archive' package solution '_field*']);
-	end
-end
-if nargout==1, 
-	varargout{1}=md; 
-end