Index: /issm/trunk/test/Par/ISMIPF.par
===================================================================
--- /issm/trunk/test/Par/ISMIPF.par	(revision 9195)
+++ /issm/trunk/test/Par/ISMIPF.par	(revision 9196)
@@ -31,6 +31,6 @@
 md.temperature=255*ones(md.numberofnodes,1);
 pos=find(md.x==min(md.x) | md.x==max(md.x) | md.y==min(md.y) | md.y==max(md.y));
-md.spcthickness(pos,1)=1;
-md.spcthickness(pos,2)=md.thickness(pos);
+md.spcthickness=NaN*ones(md.numberofnodes,1);
+md.spcthickness(pos)=md.thickness(pos);
 md.spctemperature=255*ones(md.numberofnodes,1);
 md.geothermalflux=0.4*ones(md.numberofnodes,1);
Index: /issm/trunk/test/Par/SquareEISMINT.par
===================================================================
--- /issm/trunk/test/Par/SquareEISMINT.par	(revision 9195)
+++ /issm/trunk/test/Par/SquareEISMINT.par	(revision 9196)
@@ -35,6 +35,6 @@
 %Evolution of the ice shelf
 pos=find(md.y==200000); %nodes on the upper boundary condition
-md.spcthickness(pos,1)=1;
-md.spcthickness(pos,2)=500;
+md.spcthickness=NaN*ones(md.numberofnodes,1);
+md.spcthickness(pos)=500;
 md.ndt=500;
 md.dt=1;