Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/79North.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/79North.par	(revision 13232)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/79North.par	(revision 13233)
@@ -1,7 +1,13 @@
 %Start defining model parameters here
 
 %Geometry and observation
-load('../Data/79North.data','-mat');
+x=transpose(ncread('../Data/79North.nc','x'));
+y=transpose(ncread('../Data/79North.nc','y'));
+vx=transpose(ncread('../Data/79North.nc','vx'));
+vy=transpose(ncread('../Data/79North.nc','vy'));
+index=transpose(ncread('../Data/79North.nc','index'));
+surface=transpose(ncread('../Data/79North.nc','surface'));
+thickness=transpose(ncread('../Data/79North.nc','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);
 md.geometry.surface  =InterpFromMeshToMesh2d(index,x,y,surface,md.mesh.x,md.mesh.y);
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetShelf.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetShelf.par	(revision 13232)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetShelf.par	(revision 13233)
@@ -13,7 +13,11 @@
 md.geometry.surface=md.geometry.bed+md.geometry.thickness;
 
 %Initial velocity 
-load('../Data/SquareSheetShelf.data','-mat');
+x=transpose(ncread('../Data/SquareSheetShelf.nc','x'));
+y=transpose(ncread('../Data/SquareSheetShelf.nc','y'));
+vx=transpose(ncread('../Data/SquareSheetShelf.nc','vx'));
+vy=transpose(ncread('../Data/SquareSheetShelf.nc','vy'));
+index=transpose(ncread('../Data/SquareSheetShelf.nc','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);
 clear vx vy x y index;
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelfConstrained.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 13232)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelfConstrained.par	(revision 13233)
@@ -10,7 +10,11 @@
 md.geometry.surface=md.geometry.bed+md.geometry.thickness;
 
 %Initial velocity 
-load('../Data/SquareShelfConstrained.data','-mat');
+x=transpose(ncread('../Data/SquareShelfConstrained.nc','x'));
+y=transpose(ncread('../Data/SquareShelfConstrained.nc','y'));
+vx=transpose(ncread('../Data/SquareShelfConstrained.nc','vx'));
+vy=transpose(ncread('../Data/SquareShelfConstrained.nc','vy'));
+index=transpose(ncread('../Data/SquareShelfConstrained.nc','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);
 clear vx vy x y index;
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/Pig.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/Pig.par	(revision 13232)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/Pig.par	(revision 13233)
@@ -1,7 +1,13 @@
 %Start defining model parameters here
 
 %Geometry and observation
-load('../Data/Pig.data','-mat');
+x=transpose(ncread('../Data/Pig.nc','x'));
+y=transpose(ncread('../Data/Pig.nc','y'));
+vx_obs=transpose(ncread('../Data/Pig.nc','vx_obs'));
+vy_obs=transpose(ncread('../Data/Pig.nc','vy_obs'));
+index=transpose(ncread('../Data/Pig.nc','index'));
+surface=transpose(ncread('../Data/Pig.nc','surface'));
+thickness=transpose(ncread('../Data/Pig.nc','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);
 md.geometry.surface  =InterpFromMeshToMesh2d(index,x,y,surface,md.mesh.x,md.mesh.y);
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelf.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelf.par	(revision 13232)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareShelf.par	(revision 13233)
@@ -10,7 +10,11 @@
 md.geometry.surface=md.geometry.bed+md.geometry.thickness;
 
 %Initial velocity and pressure
-load('../Data/SquareShelf.data','-mat');
+x=transpose(ncread('../Data/SquareShelf.nc','x'));
+y=transpose(ncread('../Data/SquareShelf.nc','y'));
+vx=transpose(ncread('../Data/SquareShelf.nc','vx'));
+vy=transpose(ncread('../Data/SquareShelf.nc','vy'));
+index=transpose(ncread('../Data/SquareShelf.nc','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);
 clear vx vy x y index;
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetConstrained.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 13232)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/SquareSheetConstrained.par	(revision 13233)
@@ -10,7 +10,11 @@
 md.geometry.surface=md.geometry.bed+md.geometry.thickness;
 
 %Initial velocity 
-load('../Data/SquareSheetConstrained.data','-mat');
+x=transpose(ncread('../Data/SquareSheetConstrained.nc','x'));
+y=transpose(ncread('../Data/SquareSheetConstrained.nc','y'));
+vx=transpose(ncread('../Data/SquareSheetConstrained.nc','vx'));
+vy=transpose(ncread('../Data/SquareSheetConstrained.nc','vy'));
+index=transpose(ncread('../Data/SquareSheetConstrained.nc','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);
 clear vx vy x y index;
Index: /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPE.par
===================================================================
--- /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPE.par	(revision 13232)
+++ /u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/test/Par/ISMIPE.par	(revision 13233)
@@ -1,8 +1,7 @@
 %Ok, start defining model parameters here
 
 disp('      creating thickness');
-data=load('../Data/ISMIPE.data','-mat');
-data=data.data;
+data=transpose(ncread('../Data/ISMIPE.nc','data'));
 md.geometry.surface=zeros(md.mesh.numberofvertices,1);
 md.geometry.bed=zeros(md.mesh.numberofvertices,1);
 for i=1:md.mesh.numberofvertices