Index: /issm/trunk-jpl/test/NightlyRun/runme.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/runme.m	(revision 20181)
+++ /issm/trunk-jpl/test/NightlyRun/runme.m	(revision 20182)
@@ -118,5 +118,5 @@
 	test_ids=intersect(test_ids,[1601:1602]);
 elseif strcmpi(benchmark,'slr'),
-	test_ids=intersect(test_ids,[2001:2100]);
+	test_ids=intersect(test_ids,[2001:2500]);
 elseif strcmpi(benchmark,'adolc'),
 	test_ids=intersect(test_ids,[3001:3200]);
Index: /issm/trunk-jpl/test/NightlyRun/runme.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/runme.py	(revision 20181)
+++ /issm/trunk-jpl/test/NightlyRun/runme.py	(revision 20182)
@@ -117,5 +117,5 @@
 		test_ids=test_ids.intersection(set(range(1601,1603)))
 	elif strcmpi(benchmark,'slr'):
-		test_ids=test_ids.intersection(set(range(2001,2100)))
+		test_ids=test_ids.intersection(set(range(2001,2500)))
 	elif strcmpi(benchmark,'adolc'):
 		test_ids=test_ids.intersection(set(range(3001,3200)))
Index: /issm/trunk-jpl/test/NightlyRun/test2424.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2424.m	(revision 20182)
+++ /issm/trunk-jpl/test/NightlyRun/test2424.m	(revision 20182)
@@ -0,0 +1,38 @@
+%Test Name: SquareSheetShelfGroundingLine2dAggressive. From test424, with sea level increasing.
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=setflowequation(md,'SSA','all');
+md.initialization.vx(:)=0.;
+md.initialization.vy(:)=0.;
+md.smb.mass_balance(:)=0.;
+
+md.geometry.base=-700.-abs(md.mesh.y-500000.)/1000.;
+md.geometry.bed =-700.-abs(md.mesh.y-500000.)/1000.;
+md.geometry.thickness(:)=1000.;
+md.geometry.surface=md.geometry.base+md.geometry.thickness;
+
+md.transient.isstressbalance=0;
+md.transient.isgroundingline=1;
+md.transient.isthermal=0;
+md.groundingline.migration='AggressiveMigration';
+md.transient.requested_outputs={'IceVolume','IceVolumeAboveFloatation','Sealevel'};
+
+md.timestepping.time_step=.1;
+md.slr.sealevel=newforcing(md.timestepping.start_time, md.timestepping.final_time, md.timestepping.time_step,-200,200,md.mesh.numberofvertices);
+
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,TransientSolutionEnum());
+
+%we are checking that the grounding line position is near the theorical one, which is the 0 contour level 
+%of surface - sealevel - (1-di)* thickness; 
+
+nsteps=length(md.results.TransientSolution);
+field_names={};
+field_tolerances={};
+field_values={};
+for i=1:nsteps,
+	field_names{end+1}=sprintf('Time-%g-yr-ice_levelset-S-sl-(1-di)*H',md.results.TransientSolution(i).time);
+	field_tolerances{end+1}=1e-13;
+	field_values{end+1}= md.results.TransientSolution(i).MaskGroundediceLevelset- (md.geometry.surface-md.results.TransientSolution(i).Sealevel-(1-md.materials.rho_ice/md.materials.rho_water)*md.geometry.thickness);
+end