Index: /issm/trunk-jpl/test/Par/79North.par
===================================================================
--- /issm/trunk-jpl/test/Par/79North.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/79North.par	(revision 21167)
@@ -2,11 +2,11 @@
 
 %Geometry and observation
-x=archread('../Data/79North.nc','x');
-y=archread('../Data/79North.nc','y');
-vx=archread('../Data/79North.nc','vx');
-vy=archread('../Data/79North.nc','vy');
-index=archread('../Data/79North.nc','index');
-surface=archread('../Data/79North.nc','surface');
-thickness=archread('../Data/79North.nc','thickness');
+x = archread('../Data/79North.arch','x');
+y = archread('../Data/79North.arch','y');
+vx = archread('../Data/79North.arch','vx');
+vy = archread('../Data/79North.arch','vy');
+index = archread('../Data/79North.arch','index');
+surface = archread('../Data/79North.arch','surface');
+thickness = archread('../Data/79North.arch','thickness');
 md.initialization.vx =InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
 md.initialization.vy =InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
Index: /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par
===================================================================
--- /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/GiaBenchmarksAB.par	(revision 21167)
@@ -27,9 +27,9 @@
 
 %Initial velocity 
-x     = archread('../Data/SquareSheetConstrained.nc','x');
-y     = archread('../Data/SquareSheetConstrained.nc','y');
-vx    = archread('../Data/SquareSheetConstrained.nc','vx');
-vy    = archread('../Data/SquareSheetConstrained.nc','vy');
-index = archread('../Data/SquareSheetConstrained.nc','index');
+x     = archread('../Data/SquareSheetConstrained.arch','x');
+y     = archread('../Data/SquareSheetConstrained.arch','y');
+vx    = archread('../Data/SquareSheetConstrained.arch','vx');
+vy    = archread('../Data/SquareSheetConstrained.arch','vy');
+index = archread('../Data/SquareSheetConstrained.arch','index');
 
 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
Index: /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par
===================================================================
--- /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/GiaBenchmarksCD.par	(revision 21167)
@@ -26,9 +26,9 @@
 
 %Initial velocity 
-x     = archread('../Data/SquareSheetConstrained.nc','x');
-y     = archread('../Data/SquareSheetConstrained.nc','y');
-vx    = archread('../Data/SquareSheetConstrained.nc','vx');
-vy    = archread('../Data/SquareSheetConstrained.nc','vy');
-index = archread('../Data/SquareSheetConstrained.nc','index');
+x     = archread('../Data/SquareSheetConstrained.arch','x');
+y     = archread('../Data/SquareSheetConstrained.arch','y');
+vx    = archread('../Data/SquareSheetConstrained.arch','vx');
+vy    = archread('../Data/SquareSheetConstrained.arch','vy');
+index = archread('../Data/SquareSheetConstrained.arch','index');
 
 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
Index: /issm/trunk-jpl/test/Par/ISMIPE.par
===================================================================
--- /issm/trunk-jpl/test/Par/ISMIPE.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/ISMIPE.par	(revision 21167)
@@ -2,5 +2,5 @@
 
 disp('      creating thickness');
-data=archread('../Data/ISMIPE.nc','data');
+data=archread('../Data/ISMIPE.arch','data');
 md.geometry.surface=zeros(md.mesh.numberofvertices,1);
 md.geometry.base=zeros(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/Par/Pig.par
===================================================================
--- /issm/trunk-jpl/test/Par/Pig.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/Pig.par	(revision 21167)
@@ -2,11 +2,11 @@
 
 %Geometry and observation
-x         = archread('../Data/Pig.nc','x');
-y         = archread('../Data/Pig.nc','y');
-vx_obs    = archread('../Data/Pig.nc','vx_obs');
-vy_obs    = archread('../Data/Pig.nc','vy_obs');
-index     = archread('../Data/Pig.nc','index');
-surface   = archread('../Data/Pig.nc','surface');
-thickness = archread('../Data/Pig.nc','thickness');
+x         = archread('../Data/Pig.arch','x');
+y         = archread('../Data/Pig.arch','y');
+vx_obs    = archread('../Data/Pig.arch','vx_obs');
+vy_obs    = archread('../Data/Pig.arch','vy_obs');
+index     = archread('../Data/Pig.arch','index');
+surface   = archread('../Data/Pig.arch','surface');
+thickness = archread('../Data/Pig.arch','thickness');
 md.inversion.vx_obs   =InterpFromMeshToMesh2d(index,x,y,vx_obs,md.mesh.x,md.mesh.y);
 md.inversion.vy_obs   =InterpFromMeshToMesh2d(index,x,y,vy_obs,md.mesh.x,md.mesh.y);
Index: /issm/trunk-jpl/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 21167)
@@ -13,9 +13,9 @@
 
 %Initial velocity 
-x     = archread('../Data/SquareSheetConstrained.nc','x');
-y     = archread('../Data/SquareSheetConstrained.nc','y');
-vx    = archread('../Data/SquareSheetConstrained.nc','vx');
-vy    = archread('../Data/SquareSheetConstrained.nc','vy');
-index = archread('../Data/SquareSheetConstrained.nc','index');
+x     = archread('../Data/SquareSheetConstrained.arch','x');
+y     = archread('../Data/SquareSheetConstrained.arch','y');
+vx    = archread('../Data/SquareSheetConstrained.arch','vx');
+vy    = archread('../Data/SquareSheetConstrained.arch','vy');
+index = archread('../Data/SquareSheetConstrained.arch','index');
 
 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
@@ -41,13 +41,13 @@
 
 %Numerical parameters
-md.stressbalance.viscosity_overshoot=0.0;
+md.stressbal.arche.viscosity_overshoot=0.0;
 md.masstransport.stabilization=1.;
 md.thermal.stabilization=1.;
 md.verbose=verbose(0);
 md.settings.waitonlock=30;
-md.stressbalance.restol=0.05;
+md.stressbal.arche.restol=0.05;
 md.steadystate.reltol=0.05;
-md.stressbalance.reltol=0.05;
-md.stressbalance.abstol=NaN;
+md.stressbal.arche.reltol=0.05;
+md.stressbal.arche.abstol=NaN;
 md.timestepping.time_step=1.;
 md.timestepping.final_time=3.;
Index: /issm/trunk-jpl/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareSheetShelf.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/SquareSheetShelf.par	(revision 21167)
@@ -16,9 +16,9 @@
 
 %Initial velocity 
-x     = archread('../Data/SquareSheetShelf.nc','x');
-y     = archread('../Data/SquareSheetShelf.nc','y');
-vx    = archread('../Data/SquareSheetShelf.nc','vx');
-vy    = archread('../Data/SquareSheetShelf.nc','vy');
-index = archread('../Data/SquareSheetShelf.nc','index');
+x     = archread('../Data/SquareSheetShelf.arch','x');
+y     = archread('../Data/SquareSheetShelf.arch','y');
+vx    = archread('../Data/SquareSheetShelf.arch','vx');
+vy    = archread('../Data/SquareSheetShelf.arch','vy');
+index = archread('../Data/SquareSheetShelf.arch','index');
 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
Index: /issm/trunk-jpl/test/Par/SquareShelf.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/SquareShelf.par	(revision 21167)
@@ -13,9 +13,9 @@
 
 %Initial velocity and pressure
-x     = archread('../Data/SquareShelf.nc','x');
-y     = archread('../Data/SquareShelf.nc','y');
-vx    = archread('../Data/SquareShelf.nc','vx');
-vy    = archread('../Data/SquareShelf.nc','vy');
-index = archread('../Data/SquareShelf.nc','index');
+x     = archread('../Data/SquareShelf.arch','x');
+y     = archread('../Data/SquareShelf.arch','y');
+vx    = archread('../Data/SquareShelf.arch','vx');
+vy    = archread('../Data/SquareShelf.arch','vy');
+index = archread('../Data/SquareShelf.arch','index');
 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
Index: /issm/trunk-jpl/test/Par/SquareShelf2.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf2.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/SquareShelf2.par	(revision 21167)
@@ -13,9 +13,9 @@
 
 %Initial velocity and pressure
-x     = archread('../Data/SquareShelf.nc','x');
-y     = archread('../Data/SquareShelf.nc','y');
-vx    = archread('../Data/SquareShelf.nc','vx');
-vy    = archread('../Data/SquareShelf.nc','vy');
-index = archread('../Data/SquareShelf.nc','index');
+x     = archread('../Data/SquareShelf.arch','x');
+y     = archread('../Data/SquareShelf.arch','y');
+vx    = archread('../Data/SquareShelf.arch','vx');
+vy    = archread('../Data/SquareShelf.arch','vy');
+index = archread('../Data/SquareShelf.arch','index');
 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);
Index: /issm/trunk-jpl/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 21166)
+++ /issm/trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 21167)
@@ -14,9 +14,9 @@
 
 %Initial velocity 
-x     = archread([issmdir() '/test/Data/SquareShelfConstrained.nc'],'x');
-y     = archread([issmdir() '/test/Data/SquareShelfConstrained.nc'],'y');
-vx    = archread([issmdir() '/test/Data/SquareShelfConstrained.nc'],'vx');
-vy    = archread([issmdir() '/test/Data/SquareShelfConstrained.nc'],'vy');
-index = archread([issmdir() '/test/Data/SquareShelfConstrained.nc'],'index');
+x     = archread([issmdir() '/test/Data/SquareShelfConstrained.arch'],'x');
+y     = archread([issmdir() '/test/Data/SquareShelfConstrained.arch'],'y');
+vx    = archread([issmdir() '/test/Data/SquareShelfConstrained.arch'],'vx');
+vy    = archread([issmdir() '/test/Data/SquareShelfConstrained.arch'],'vy');
+index = archread([issmdir() '/test/Data/SquareShelfConstrained.arch'],'index');
 md.initialization.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.mesh.x,md.mesh.y);
 md.initialization.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.mesh.x,md.mesh.y);