Index: /issm/trunk/test/NightlyRun/test445.m
===================================================================
--- /issm/trunk/test/NightlyRun/test445.m	(revision 10303)
+++ /issm/trunk/test/NightlyRun/test445.m	(revision 10304)
@@ -17,27 +17,11 @@
 md=setflowequation(md,'macayeal','all');
 md.cluster=none;
-%plug holes into the ice sheet, to test for grounding line migration. 
-di=md.materials.rho_ice/md.materials.rho_water;
-rad=sqrt(md.mesh.x.^2+md.mesh.y.^2);
-pos=find(rad<200000);
-md.geometry.thickness(pos)=100;
-md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
-md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
 
-pos=find(md.mesh.x<.2*1e6 & md.mesh.x>-.2*1e6 & md.mesh.y>0);
-md.geometry.thickness(pos)=100;
-md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
-md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
-
-pos=find(md.mesh.x<.1*1e6 & md.mesh.x>-.1*1e6 & md.mesh.y<-.5*1e6 & md.mesh.y>-.6*1e6);
-md.geometry.thickness(pos)=100;
-md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
-md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
-
-%test different grounding line dynamics.
 md.transient.isthermal=0;
 md.transient.isprognostic=0;
 md.transient.isdiagnostic=0;
 md.transient.isgroundingline=1;
+
+%test different grounding line dynamics
 md.groundingline.migration='AgressiveMigration';
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test446.m
===================================================================
--- /issm/trunk/test/NightlyRun/test446.m	(revision 10303)
+++ /issm/trunk/test/NightlyRun/test446.m	(revision 10304)
@@ -17,27 +17,11 @@
 md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
-%plug holes into the ice sheet, to test for grounding line migration. 
-di=md.materials.rho_ice/md.materials.rho_water;
-rad=sqrt(md.mesh.x.^2+md.mesh.y.^2);
-pos=find(rad<200000);
-md.geometry.thickness(pos)=100;
-md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
-md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
 
-pos=find(md.mesh.x<.2*1e6 & md.mesh.x>-.2*1e6 & md.mesh.y>0);
-md.geometry.thickness(pos)=100;
-md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
-md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
-
-pos=find(md.mesh.x<.1*1e6 & md.mesh.x>-.1*1e6 & md.mesh.y<-.5*1e6 & md.mesh.y>-.6*1e6);
-md.geometry.thickness(pos)=100;
-md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
-md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
-
-%test different grounding line dynamics.
 md.transient.isthermal=0;
 md.transient.isprognostic=0;
 md.transient.isdiagnostic=0;
 md.transient.isgroundingline=1;
+
+%test different grounding line dynamics.
 md.groundingline.migration='AgressiveMigration';
 md=solve(md,TransientSolutionEnum);
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 10303)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 10304)
@@ -33,4 +33,21 @@
 md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
 
+%plug holes into the ice sheet, to test for grounding line migration. 
+di=md.materials.rho_ice/md.materials.rho_water;
+rad=sqrt(md.mesh.x.^2+md.mesh.y.^2);
+pos=find(rad<200000);
+md.geometry.thickness(pos)=100;
+md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
+md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
+
+pos=find(md.mesh.x<.2*1e6 & md.mesh.x>-.2*1e6 & md.mesh.y>0);
+md.geometry.thickness(pos)=100;
+md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
+md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
+
+pos=find(md.mesh.x<.1*1e6 & md.mesh.x>-.1*1e6 & md.mesh.y<-.5*1e6 & md.mesh.y>-.6*1e6);
+md.geometry.thickness(pos)=100;
+md.geometry.bed(pos)=-di*md.geometry.thickness(pos)-20;
+md.geometry.surface(pos)=md.geometry.bed(pos)+md.geometry.thickness(pos);
 %Initial velocity 
 md.initialization.vx=zeros(md.mesh.numberofvertices,1);