Index: /issm/trunk-jpl/test/SandBox/test2004.m
===================================================================
--- /issm/trunk-jpl/test/SandBox/test2004.m	(revision 24763)
+++ /issm/trunk-jpl/test/SandBox/test2004.m	(revision 24764)
@@ -8,4 +8,5 @@
 
 skipmesh=1;
+skipparameterization=0;
 
 %create sealevel model to hold our information: 
@@ -91,4 +92,5 @@
 %}}}
 end
+if ~skipparameterization,
 %Parameterize ice sheets : {{{
 
@@ -180,4 +182,6 @@
 %}}}
 % ParameterizeContinents {{{
+
+sl.basinindx('continent',{'hemisphereeast','hemispherewest'})
 
 for ind=sl.basinindx('continent',{'hemisphereeast','hemispherewest'}),
@@ -240,4 +244,6 @@
 	%grounded ice: 
 	md.mask.groundedice_levelset=-ones(md.mesh.numberofvertices,1); 
+	
+	md.mask.glacier_levelset=zeros(md.mesh.numberofvertices,1);
 
 	% }}}
@@ -269,5 +275,5 @@
 end
 % }}}
-
+end
 %Assemble Earth in 3D {{{
 
@@ -279,6 +285,4 @@
 %create earth model by concatenating all the icecaps in 3d: 
 sl.caticecaps('tolerance',tolerance,'loneedgesdetect',loneedgesdetect);
-
-error;
 
 %figure out how each icecap's mesh connects to the larger earth mesh: 
@@ -299,10 +303,5 @@
 sl.transfer('mask.land_levelset');
 sl.transfer('mask.glacier_levelset');
-sl.transfer('geometry.thickness');
-sl.transfer('geometry.surface');
-sl.transfer('geometry.base');
 sl.transfer('geometry.bed');
-sl.transfer('initialization.temperature');
-sl.transfer('materials.rheology_B');
 sl.transfer('mesh.lat');
 sl.transfer('mesh.long');
@@ -310,4 +309,5 @@
 sl.transfer('slr.Ngia');
 sl.transfer('slr.Ugia');
+error;
 
 %radius: