Index: sm/trunk/test/Verification/PigControlMethodBM2d_21/runme.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodBM2d_21/runme.m	(revision 2167)
+++ 	(revision )
@@ -1,101 +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'},{'absolute'});
-%      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 Pine Island Glacierfor a 2d model.
-% 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'macayeal','cielo_parallel','cielo_serial'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodBM2d_21';
-
-%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_absolute')) | strcmpi(package,'ice')),
-			disp(['package: ' package  ' and solution: ' solution  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=setelementstype(md,'macayeal','all');
-		md=modelextract(md,md.elementoniceshelf);
-
-		%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: sm/trunk/test/Verification/PigControlMethodBM2d_21/updatearchive.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodBM2d_21/updatearchive.m	(revision 2167)
+++ 	(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 with 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'macayeal','cielo_parallel','cielo_serial'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodBM2d_21';
-
-%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_absolute')) | strcmpi(package,'ice')),
-			disp(['package: ' package  ' and solution: ' solution  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=setelementstype(md,'macayeal','all');
-		md=modelextract(md,md.elementoniceshelf);
-
-		%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
Index: sm/trunk/test/Verification/PigControlMethodDragM2d_20/runme.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodDragM2d_20/runme.m	(revision 2167)
+++ 	(revision )
@@ -1,100 +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'},{'absolute'});
-%      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 Pine Island Glacierfor a 2d model.
-% 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'cielo_serial','cielo_parallel'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodDragM2d_20';
-
-%go through the solutions requested
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check package
-		if ~(strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel'));
-			disp(['package: ' package  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=setelementstype(md,'macayeal','all');
-
-		%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: sm/trunk/test/Verification/PigControlMethodDragM2d_20/updatearchive.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodDragM2d_20/updatearchive.m	(revision 2167)
+++ 	(revision )
@@ -1,81 +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'},{'absolute'});
-%      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 with 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'cielo_parallel','cielo_serial'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodDragM2d_20';
-
-%go through the solutions requested
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check package
-		if ~(strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel'));
-			disp(['package: ' package  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=setelementstype(md,'macayeal','all');
-
-		%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
Index: sm/trunk/test/Verification/PigControlMethodDragP3d_22/runme.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodDragP3d_22/runme.m	(revision 2167)
+++ 	(revision )
@@ -1,101 +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'},{'absolute'});
-%      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 Pine Island Glacierfor a 2d model.
-% 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'cielo_serial','cielo_parallel'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodDragP3d_22';
-
-%go through the solutions requested
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check package
-		if ~(strcmpi(package,'cielo_parallel'));
-			disp(['package: ' package  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=extrude(md,4,2);
-		md=setelementstype(md,'pattyn','all');
-
-		%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: sm/trunk/test/Verification/PigControlMethodDragP3d_22/updatearchive.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodDragP3d_22/updatearchive.m	(revision 2167)
+++ 	(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'},{'absolute'});
-%      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 with 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'cielo_parallel','cielo_serial'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodDragM2d_20';
-
-%go through the solutions requested
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check package
-		if ~(strcmpi(package,'cielo_parallel'));
-			disp(['package: ' package  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=extrude(md,4,2);
-		md=setelementstype(md,'pattyn','all');
-
-		%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
Index: sm/trunk/test/Verification/PigControlMethodDragS3d_23/runme.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodDragS3d_23/runme.m	(revision 2167)
+++ 	(revision )
@@ -1,101 +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'},{'absolute'});
-%      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 Pine Island Glacierfor a 2d model.
-% 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'cielo_serial','cielo_parallel'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodDragS3d_23';
-
-%go through the solutions requested
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check package
-		if ~(strcmpi(package,'cielo_parallel'));
-			disp(['package: ' package  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=extrude(md,4,3);
-		md=setelementstype(md,'pattyn','all','stokes','all');
-
-		%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: sm/trunk/test/Verification/PigControlMethodDragS3d_23/updatearchive.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodDragS3d_23/updatearchive.m	(revision 2167)
+++ 	(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'},{'absolute'});
-%      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 with 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 solutions to be tested
-if nargin==2,
-	solutions=varargin{2};
-else
-	solutions={'diagnostic_absolute','diagnostic_relative','diagnostic_logarithmic'};
-end
-if nargin,
-	packages=varargin{1};
-else
-	packages={'cielo_parallel','cielo_serial'};
-end
-
-%Initialize log message for nightly runs.
-testname='PigControlMethodDragS3d_23';
-
-%go through the solutions requested
-for i=1:length(packages),
-	package=packages{i};
-
-	for j=1:length(solutions),
-		solution=solutions{j};
-
-		%check package
-		if ~(strcmpi(package,'cielo_parallel'));
-			disp(['package: ' package  ' in test: ' testname  ' not supported yet']);
-			continue
-		end
-
-		%initialize model
-		md=model;
-		md=mesh(md,'DomainOutline.exp',9000);
-		md=geography(md,'Iceshelves.exp','Islands.exp');
-		md=parameterize(md,'Pig.par');
-		md=extrude(md,4,3);
-		md=setelementstype(md,'pattyn','all','stokes','all');
-
-		%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
Index: /issm/trunk/test/Verification/list
===================================================================
--- /issm/trunk/test/Verification/list	(revision 2168)
+++ /issm/trunk/test/Verification/list	(revision 2168)
@@ -0,0 +1,174 @@
+IceSheetNoIceFrontH2d_12/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontH2d_12/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontH2d_12/Archive_transient_none_0_0_nan_0.mat
+IceSheetNoIceFrontH2d_12/Archive_transient_none_0_0_nan_1.mat
+IceSheetNoIceFrontH2d_12/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontH2d_12/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetIceFrontM2dDakota_25/Archive_diagnostic_none_1_0_nan_1.mat
+IceSheetIceFrontM2dDakota_25/Archive_diagnostic_none_1_0_nan_0.mat
+IceShelfIceFrontMP3d_3/Archive_thermal_steady_0_0_nan_0.mat
+IceShelfIceFrontMP3d_3/Archive_thermal_steady_0_0_nan_1.mat
+IceShelfIceFrontMP3d_3/Archive_prognostic_none_0_0_nan_0.mat
+IceShelfIceFrontMP3d_3/Archive_prognostic_none_0_0_nan_1.mat
+IceShelfIceFrontMP3d_3/Archive_thermal_transient_0_0_nan_0.mat
+IceShelfIceFrontMP3d_3/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetNoIceFrontH3d_13/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontH3d_13/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontH3d_13/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetNoIceFrontH3d_13/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetNoIceFrontH3d_13/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontH3d_13/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontH3d_13/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetNoIceFrontH3d_13/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetNoIceFrontH3d_13/Archive_transient_none_0_0_nan_0.mat
+IceSheetNoIceFrontH3d_13/Archive_transient_none_0_0_nan_1.mat
+IceSheetNoIceFrontM2d_16/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontM2d_16/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontM2d_16/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontM2d_16/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontM2d_16/Archive_transient_none_0_0_nan_0.mat
+IceSheetNoIceFrontM2d_16/Archive_transient_none_0_0_nan_1.mat
+IceSheetIceFrontP3d_10/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetIceFrontP3d_10/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetIceFrontP3d_10/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetIceFrontP3d_10/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetIceFrontP3d_10/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetIceFrontP3d_10/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetIceFrontP3d_10/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetIceFrontP3d_10/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetIceFrontP3d_10/Archive_transient_none_0_0_nan_0.mat
+IceSheetIceFrontP3d_10/Archive_transient_none_0_0_nan_1.mat
+IceSheetNoIceFrontM3d_17/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontM3d_17/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontM3d_17/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetNoIceFrontM3d_17/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetNoIceFrontM3d_17/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontM3d_17/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontM3d_17/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetNoIceFrontM3d_17/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetNoIceFrontM3d_17/Archive_transient_none_0_0_nan_0.mat
+IceSheetNoIceFrontM3d_17/Archive_transient_none_0_0_nan_1.mat
+IceShelfIceFrontM2dRifts_24/Archive_diagnostic_none_0_0_nan_0.mat
+IceShelfIceFrontM2dRifts_24/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetIceFrontS3d_11/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetIceFrontS3d_11/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetIceFrontS3d_11/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetIceFrontS3d_11/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetIceFrontS3d_11/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetIceFrontS3d_11/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetIceFrontS3d_11/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetIceFrontS3d_11/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetNoIceFrontHM2d_14/Archive_transient_none_0_0_nan_0.mat
+IceSheetNoIceFrontHM2d_14/Archive_transient_none_0_0_nan_1.mat
+IceSheetNoIceFrontHM2d_14/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontHM2d_14/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontHM2d_14/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontHM2d_14/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontP3d_18/Archive_transient_none_0_0_nan_0.mat
+IceSheetNoIceFrontP3d_18/Archive_transient_none_0_0_nan_1.mat
+IceSheetNoIceFrontP3d_18/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontP3d_18/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontP3d_18/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetNoIceFrontP3d_18/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetNoIceFrontP3d_18/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontP3d_18/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontP3d_18/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetNoIceFrontP3d_18/Archive_thermal_transient_0_0_nan_1.mat
+IceShelfIceFrontM2d_1/Archive_transient_none_0_0_nan_0.mat
+IceShelfIceFrontM2d_1/Archive_transient_none_0_0_nan_1.mat
+IceShelfIceFrontM2d_1/Archive_diagnostic_none_0_0_nan_0.mat
+IceShelfIceFrontM2d_1/Archive_diagnostic_none_0_0_nan_1.mat
+IceShelfIceFrontM2d_1/Archive_prognostic_none_0_0_nan_0.mat
+IceShelfIceFrontM2d_1/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetIceFrontMP3d_9/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetIceFrontMP3d_9/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetIceFrontMP3d_9/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetIceFrontMP3d_9/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetIceFrontMP3d_9/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetIceFrontMP3d_9/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetNoIceFrontHM3d_15/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetNoIceFrontHM3d_15/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontHM3d_15/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontHM3d_15/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetNoIceFrontHM3d_15/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetNoIceFrontHM3d_15/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontHM3d_15/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontHM3d_15/Archive_thermal_transient_0_0_nan_0.mat
+IceShelfIceFrontM3d_2/Archive_thermal_transient_0_0_nan_1.mat
+IceShelfIceFrontM3d_2/Archive_transient_none_0_0_nan_0.mat
+IceShelfIceFrontM3d_2/Archive_transient_none_0_0_nan_1.mat
+IceShelfIceFrontM3d_2/Archive_diagnostic_none_0_0_nan_0.mat
+IceShelfIceFrontM3d_2/Archive_diagnostic_none_0_0_nan_1.mat
+IceShelfIceFrontM3d_2/Archive_thermal_steady_0_0_nan_0.mat
+IceShelfIceFrontM3d_2/Archive_thermal_steady_0_0_nan_1.mat
+IceShelfIceFrontM3d_2/Archive_prognostic_none_0_0_nan_0.mat
+IceShelfIceFrontM3d_2/Archive_prognostic_none_0_0_nan_1.mat
+IceShelfIceFrontM3d_2/Archive_steadystate_none_0_0_nan_0.mat
+IceShelfIceFrontM3d_2/Archive_steadystate_none_0_0_nan_1.mat
+IceShelfIceFrontM3d_2/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetNoIceFrontS3d_19/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetNoIceFrontS3d_19/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetNoIceFrontS3d_19/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontS3d_19/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetNoIceFrontS3d_19/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetNoIceFrontS3d_19/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetNoIceFrontS3d_19/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetNoIceFrontS3d_19/Archive_prognostic_none_0_0_nan_1.mat
+IceShelfIceFrontMPS3d_4/Archive_thermal_transient_0_0_nan_0.mat
+IceShelfIceFrontMPS3d_4/Archive_thermal_transient_0_0_nan_1.mat
+IceShelfIceFrontMPS3d_4/Archive_thermal_steady_0_0_nan_0.mat
+IceShelfIceFrontMPS3d_4/Archive_thermal_steady_0_0_nan_1.mat
+IceShelfIceFrontMPS3d_4/Archive_prognostic_none_0_0_nan_0.mat
+IceShelfIceFrontMPS3d_4/Archive_prognostic_none_0_0_nan_1.mat
+IceShelfIceFrontP3d_5/Archive_thermal_transient_0_0_nan_0.mat
+IceShelfIceFrontP3d_5/Archive_thermal_transient_0_0_nan_1.mat
+IceShelfIceFrontP3d_5/Archive_transient_none_0_0_nan_0.mat
+IceShelfIceFrontP3d_5/Archive_transient_none_0_0_nan_1.mat
+IceShelfIceFrontP3d_5/Archive_diagnostic_none_0_0_nan_0.mat
+IceShelfIceFrontP3d_5/Archive_diagnostic_none_0_0_nan_1.mat
+IceShelfIceFrontP3d_5/Archive_thermal_steady_0_0_nan_0.mat
+IceShelfIceFrontP3d_5/Archive_thermal_steady_0_0_nan_1.mat
+IceShelfIceFrontP3d_5/Archive_prognostic_none_0_0_nan_0.mat
+IceShelfIceFrontP3d_5/Archive_prognostic_none_0_0_nan_1.mat
+IceShelfIceFrontS3d_6/Archive_steadystate_none_0_0_nan_1.mat
+IceShelfIceFrontS3d_6/Archive_thermal_transient_0_0_nan_0.mat
+IceShelfIceFrontS3d_6/Archive_thermal_transient_0_0_nan_1.mat
+IceShelfIceFrontS3d_6/Archive_diagnostic_none_0_0_nan_0.mat
+IceShelfIceFrontS3d_6/Archive_diagnostic_none_0_0_nan_1.mat
+IceShelfIceFrontS3d_6/Archive_thermal_steady_0_0_nan_0.mat
+IceShelfIceFrontS3d_6/Archive_thermal_steady_0_0_nan_1.mat
+IceShelfIceFrontS3d_6/Archive_prognostic_none_0_0_nan_0.mat
+IceShelfIceFrontS3d_6/Archive_prognostic_none_0_0_nan_1.mat
+IceShelfIceFrontS3d_6/Archive_steadystate_none_0_0_nan_0.mat
+IceSheetIceFrontM2d_7/Archive_transient_none_0_0_nan_0.mat
+IceSheetIceFrontM2d_7/Archive_transient_none_0_0_nan_1.mat
+IceSheetIceFrontM2d_7/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetIceFrontM2d_7/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetIceFrontM2d_7/Archive_prognostic_none_0_0_nan_0.mat
+IceSheetIceFrontM2d_7/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetIceFrontM3d_8/Archive_prognostic_none_0_0_nan_1.mat
+IceSheetIceFrontM3d_8/Archive_thermal_transient_0_0_nan_0.mat
+IceSheetIceFrontM3d_8/Archive_thermal_transient_0_0_nan_1.mat
+IceSheetIceFrontM3d_8/Archive_transient_none_0_0_nan_0.mat
+IceSheetIceFrontM3d_8/Archive_transient_none_0_0_nan_1.mat
+IceSheetIceFrontM3d_8/Archive_diagnostic_none_0_0_nan_0.mat
+IceSheetIceFrontM3d_8/Archive_diagnostic_none_0_0_nan_1.mat
+IceSheetIceFrontM3d_8/Archive_thermal_steady_0_0_nan_0.mat
+IceSheetIceFrontM3d_8/Archive_thermal_steady_0_0_nan_1.mat
+IceSheetIceFrontM3d_8/Archive_prognostic_none_0_0_nan_0.mat
+PigControlMethodDragM2d_20/Archive_diagnostic_none_0_1_logarithmic_0.mat
+PigControlMethodDragM2d_20/Archive_diagnostic_none_0_1_logarithmic_1.mat
+PigControlMethodDragM2d_20/Archive_diagnostic_none_0_1_relative_0.mat
+PigControlMethodDragM2d_20/Archive_diagnostic_none_0_1_relative_1.mat
+PigControlMethodDragM2d_20/Archive_diagnostic_none_0_1_absolute_0.mat
+PigControlMethodDragM2d_20/Archive_diagnostic_none_0_1_absolute_1.mat
+PigControlMethodBM2d_21/Archive_diagnostic_none_0_1_logarithmic_0.mat
+PigControlMethodBM2d_21/Archive_diagnostic_none_0_1_logarithmic_1.mat
+PigControlMethodBM2d_21/Archive_diagnostic_none_0_1_relative_0.mat
+PigControlMethodBM2d_21/Archive_diagnostic_none_0_1_relative_1.mat
+PigControlMethodBM2d_21/Archive_diagnostic_none_0_1_absolute_0.mat
+PigControlMethodBM2d_21/Archive_diagnostic_none_0_1_absolute_1.mat
+PigControlMethodDragP3d_22/Archive_diagnostic_none_0_1_absolute_1.mat
+PigControlMethodDragP3d_22/Archive_diagnostic_none_0_1_logarithmic_1.mat
+PigControlMethodDragP3d_22/Archive_diagnostic_none_0_1_relative_1.mat
+PigControlMethodDragS3d_23/Archive_diagnostic_none_0_1_absolute_1.mat
