Index: /issm/trunk/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk/test/NightlyRun/IdToName.m	(revision 7074)
+++ /issm/trunk/test/NightlyRun/IdToName.m	(revision 7075)
@@ -188,4 +188,5 @@
 	case 443, name='SquareSheetShelfDiagMS3dTilingSerial';
 	case 444, name='SquareSheetShelfDiagMS3dTilingParallel';
+	case 446, name='RoundSheetShelfDiagM2dParallel';
 	case 501, name='PigDiagM2dSerial';
 	case 502, name='PigDiagM2dParallel';
Index: /issm/trunk/test/NightlyRun/test446.m
===================================================================
--- /issm/trunk/test/NightlyRun/test446.m	(revision 7075)
+++ /issm/trunk/test/NightlyRun/test446.m	(revision 7075)
@@ -0,0 +1,35 @@
+radius=1e6;
+shelfextent=2e5;
+resolution=.5e5;
+
+%mesh {{{1
+md=roundmesh(model,radius,resolution);
+%fix center grid to 0,0
+rad=sqrt((md.x).*md.x+(md.y).*md.y);
+pos=find(rad==min(rad));
+md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
+%}}}
+%put our grounding line 'shelfextent' meters from the icefront {{{1
+xelem=md.x(md.elements)*[1;1;1]/3;
+yelem=md.y(md.elements)*[1;1;1]/3;
+rad=sqrt((xelem).*xelem+(yelem).*yelem);
+flags=zeros(md.numberofelements,1);
+pos=find(rad>=(radius-shelfextent));
+flags(pos)=1;
+flags2=zeros(md.numberofelements,1);
+md=geography(md,flags,flags2); 
+%}}}
+md=parameterize(md,'../Par/RoundSheetShelf.par');
+md=setelementstype(md,'macayeal','all');
+md=SetParallel(md,3);
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
+
+%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),...
+	};
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 7075)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 7075)
@@ -0,0 +1,72 @@
+%Start defining model parameters here
+
+rad=1e6;
+shelfextent=2e5;
+%Geometry
+hmin=300;
+hmax=1000;
+radius=sqrt((md.x).*md.x+(md.y).*md.y);
+ymin=min(radius);
+ymax=max(radius);
+md.thickness=hmax+(hmin-hmax)*(radius-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+
+pos=find(md.gridonicesheet);
+md.bed(pos)=md.bed(pos)-300*(radius(pos)-(rad-shelfextent))/(rad-shelfextent);
+md.surface=md.bed+md.thickness;
+
+%Initial velocity 
+md.vx=zeros(md.numberofgrids,1);
+md.vy=zeros(md.numberofgrids,1);
+md.vz=zeros(md.numberofgrids,1);
+md.pressure=zeros(md.numberofgrids,1);
+
+%Materials
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+md.temperature=md.observed_temperature;
+
+%Accumulation and melting
+md.accumulation_rate=-10*ones(md.numberofgrids,1);
+md.melting_rate=zeros(md.numberofgrids,1);
+pos=find(md.gridoniceshelf);md.melting_rate(pos)=10;
+
+%Friction
+pos=find(md.elementoniceshelf);
+md.drag_type=2;
+md.drag_coefficient=20*ones(md.numberofgrids,1);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+
+%Numerical parameters
+md.viscosity_overshoot=0.0;
+md.artificial_diffusivity=1;
+md.verbose=verbose(0);
+md.waitonlock=30;
+md.eps_res=0.05;
+md.eps_rel=0.05;
+md.eps_abs=NaN;
+md.dt=.1;
+md.ndt=5;
+
+%bathymetry and grounding line migration:
+md.gl_migration=1;
+md.bathymetry=md.bed;
+pos=find(md.gridoniceshelf); md.bathymetry(pos)=md.bed(pos)-900;
+
+%Deal with boundary conditions:
+md.spcvelocity=zeros(md.numberofgrids,6);
+pos=find(md.x==0 & md.y==0);
+md.spcvelocity(pos,1:2)=1;
+pressureload=md.segments;
+pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))+AirEnum*md.elementonicesheet(pressureload(:,end))];
+md.pressureload=pressureload;
+md.spcthickness=zeros(md.numberofgrids,2);
+
+
+
+%Change name so that no test have the same name
+A=dbstack;
+if (length(A)>2), md.name=A(3).file(1:end-2); end