Index: /issm/trunk/test/NightlyRun/test446.m
===================================================================
--- /issm/trunk/test/NightlyRun/test446.m	(revision 7358)
+++ /issm/trunk/test/NightlyRun/test446.m	(revision 7359)
@@ -2,5 +2,4 @@
 shelfextent=2e5;
 resolution=.5e5;
-
 %mesh {{{1
 md=roundmesh(model,radius,resolution);
@@ -23,13 +22,35 @@
 md=setelementstype(md,'macayeal','all');
 md=SetParallel(md,3);
-md=solve(md,'analysis_type',Transient2DSolutionEnum);
+%plug holes into the ice sheet, to test for grounding line migration. {{{1
+di=md.rho_ice/md.rho_water;
+rad=sqrt((md.x).*md.x+(md.y).*md.y);
+pos=find(rad<200000);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+%}}}
+
+%test different grounding line dynamics.
+%grounding line: 
+md.gl_migration=AgressiveMigrationEnum;
+md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum);
+element_on_iceshelf_agressive=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf);
+
+md.gl_migration=SoftMigrationEnum;
+md=solve(md,'analysis_type',GroundingLineMigration2DSolutionEnum);
+element_on_iceshelf_soft=PatchToVec(md.results.GroundingLineMigration2DSolution.ElementOnIceShelf);
 
 %Fields and tolerances to track changes
-field_names     ={'ElementOnIceShelf','Vx50','Vy50','Vel50'};
-field_tolerances={1e-13,1e-13,1e-13,1e-13};
-field_values={...
-	PatchToVec(md.results.Transient2DSolution(50).ElementOnIceShelf),...
-	PatchToVec(md.results.Transient2DSolution(50).Vx),...
-	PatchToVec(md.results.Transient2DSolution(50).Vy),...
-	PatchToVec(md.results.Transient2DSolution(50).Vel),...
-	};
+field_names     ={'ElementOnIceShelfAgressive','ElementOnIceShelfSoft'};
+field_tolerances={1e-13,1e-13};
+field_values={element_on_iceshelf_agressive,element_on_iceshelf_soft};
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 7358)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 7359)
@@ -1,4 +1,5 @@
 %Start defining model parameters here
 
+di=md.rho_ice/md.rho_water;
 rad=1e6;
 shelfextent=2e5;
@@ -15,4 +16,20 @@
 md.bed(pos)=md.bed(pos)-300*(radius(pos)-(rad-shelfextent))/(rad-shelfextent);
 md.surface=md.bed+md.thickness;
+
+pos=find(radius<200000);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+
+pos=find(md.x<.2*1e6 & md.x>-.2*1e6 & md.y>0);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
+
+pos=find(md.x<.1*1e6 & md.x>-.1*1e6 & md.y<-.5*1e6 & md.y>-.6*1e6);
+md.thickness(pos)=100;
+md.bed(pos)=-di*md.thickness(pos)-20;
+md.surface(pos)=md.bed(pos)+md.thickness(pos);
 
 %Initial velocity