Index: /issm/trunk/test/NightlyRun/Id2Name.m
===================================================================
--- /issm/trunk/test/NightlyRun/Id2Name.m	(revision 5061)
+++ /issm/trunk/test/NightlyRun/Id2Name.m	(revision 5062)
@@ -41,4 +41,5 @@
 elseif (id==135), name='SquareSheetConstrainedBalThic2dDGSerial';
 elseif (id==136), name='SquareSheetConstrainedBalThic2dDGParallel';
+elseif (id==137), name='SquareBamgMesh';
 elseif (id==201), name='SquareShelfDiagM2dSerial';
 elseif (id==202), name='SquareShelfDiagM2dParallel';
@@ -175,4 +176,5 @@
 elseif (id==525), name='PigCMDragSteaM3dSerial';
 elseif (id==526), name='PigCMDragSteaM3dparallel';
+elseif (id==527), name='PigBamgMesh';
 elseif (id==601), name='79NorthProg2dSerial';
 elseif (id==602), name='79NorthProg2dParallel';
Index: /issm/trunk/test/NightlyRun/test137.m
===================================================================
--- /issm/trunk/test/NightlyRun/test137.m	(revision 5062)
+++ /issm/trunk/test/NightlyRun/test137.m	(revision 5062)
@@ -0,0 +1,16 @@
+
+%Simple mesh
+md=bamg(model,'domain','../Exp/Square.exp','hmax',100000);
+x1=md.x;
+y1=md.y;
+
+%hVertices
+md=bamg(model,'domain','../Exp/Square.exp','hmax',300000,'hvertices',[10000 100000 400000 100000]');
+x2=md.x;
+y2=md.y;
+
+%big mesh
+t0=clock;
+md=bamg(model,'domain','../Exp/Square.exp','hmax',3000);
+nbelements=md.numberofelements;
+elapsedtime=etime(clock,t0);
Index: /issm/trunk/test/NightlyRun/test137_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test137_nightly.m	(revision 5062)
+++ /issm/trunk/test/NightlyRun/test137_nightly.m	(revision 5062)
@@ -0,0 +1,7 @@
+field_names     ={'x1' 'y1' 'x2' 'y2' 'nbelements' 'elapsed time'};
+field_tolerances={1e-13 1e-13 1e-13 1e-13 1e-13 0.5};
+field_values={...
+	x1, y1,...
+	x2, y2,...
+	nbelements,elapsedtime...
+	};
Index: /issm/trunk/test/NightlyRun/test527.m
===================================================================
--- /issm/trunk/test/NightlyRun/test527.m	(revision 5062)
+++ /issm/trunk/test/NightlyRun/test527.m	(revision 5062)
@@ -0,0 +1,43 @@
+%Simple mesh 1
+hVertices=10000*ones(27,1);
+hVertices(1:5)=1000;
+md=bamg(model,'domain','../Exp/Pig.exp','hmax',20000,'hVertices',hVertices,'gradation',3,'geometricalmetric',1);
+x1=md.x;
+y1=md.y;
+
+%Simple mesh 2
+md=bamg(model,'domain','../Exp/Pig.exp','hmax',10000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+x2=md.x;
+y2=md.y;
+
+%refine existing mesh 1
+hessian=ComputeHessian(md.elements,md.x,md.y,md.vy_obs,'node');
+metric=ComputeMetric(hessian,2/9,1,1000,25*10^3,[]);
+md.dummy=metric;
+md2=bamg(md,'metric',md.dummy,'hmin',1000,'hmax',20000,'gradation',3,'geometricalmetric',1);
+x3=md2.x;
+y3=md2.y;
+
+%refine existing mesh 2
+md2=bamg(md,'metric',md.dummy,'hmin',1000,'hmax',20000,'gradation',3,'geometricalmetric',1,'anisomax',1);
+x4=md2.x;
+y4=md2.y;
+
+%refine existing mesh 3
+hVertices=NaN*ones(md.numberofgrids,1);
+hVertices(find(md.gridonboundary))=500;
+md2=bamg(md,'metric',md.dummy,'hmin',1000,'hmax',20000,'gradation',3,'geometricalmetric',1,'anisomax',1,'hVertices',hVertices);
+x5=md2.x;
+y5=md2.y;
+
+%refine existing mesh 4
+md2=bamg(md,'field',md.vy_obs,'hmin',1000,'hmax',20000,'gradation',3,'geometricalmetric',1,'Hessiantype',0,'err',1);
+x6=md2.x;
+y6=md2.y;
+
+%refine existing mesh 5
+md2=bamg(md,'field',[md.vy_obs md.thickness],'hmin',1000,'hmax',20000,'gradation',3,'geometricalmetric',1,'Hessiantype',1,'err',[10 100]); 
+x7=md2.x;
+y7=md2.y;
Index: /issm/trunk/test/NightlyRun/test527_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test527_nightly.m	(revision 5062)
+++ /issm/trunk/test/NightlyRun/test527_nightly.m	(revision 5062)
@@ -0,0 +1,11 @@
+field_names     ={'x1' 'y1' 'x2' 'y2' 'x3' 'y3' 'x4' 'y4' 'x5' 'y5' 'x6' 'y6'  'x7' 'y7' };
+field_tolerances={1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13};
+field_values={...
+	x1, y1,...
+	y2, y2,...
+	y3, y3,...
+	y4, y4,...
+	y5, y5,...
+	y6, y6,...
+	y7, y7,...
+	};