Index: /issm/trunk/test/Exp/CrossLineEISMINT.exp
===================================================================
--- /issm/trunk/test/Exp/CrossLineEISMINT.exp	(revision 5106)
+++ /issm/trunk/test/Exp/CrossLineEISMINT.exp	(revision 5106)
@@ -0,0 +1,7 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+2 1.
+# X pos Y pos
+10000 199999 
+10000 1
Index: /issm/trunk/test/NightlyRun/test1106.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1106.m	(revision 5105)
+++ /issm/trunk/test/NightlyRun/test1106.m	(revision 5106)
@@ -3,5 +3,5 @@
 
 L_list={5000,10000,20000,40000,80000,160000};
-results{};
+results={};
 
 for i=1:length(L_list),
Index: /issm/trunk/test/NightlyRun/test1201.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1201.m	(revision 5106)
+++ /issm/trunk/test/NightlyRun/test1201.m	(revision 5106)
@@ -0,0 +1,58 @@
+%This test is a test from the EISMINT for Ice shelves 
+% Vincent Rommelaere 1996
+
+%The goal is to test the prognostic model
+%md=bamg(model,'domain','DomainOutline.exp','hmax',4550);
+md=bamg(model,'domain','../Exp/SquareEISMINT.exp','hmax',3000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareEISMINT.par');
+md=setelementstype(md,'macayeal','all');
+
+disp('      initial velocity');
+md.vx=zeros(md.numberofgrids,1);
+md.vy=-400*ones(md.numberofgrids,1);
+md.vz=zeros(md.numberofgrids,1);
+
+%analytical test
+md.spcvelocity(:,1:3)=1;
+md.spcvelocity(:,4)=0;
+md.spcvelocity(:,5)=-400;
+md.spcvelocity(1,1)=0;
+md.prognostic_DG=1;
+
+%Launch many prognostic as we don't really want a transient
+%FOR NOW:
+i=0;
+time=0;
+md.dummy=struct();
+pos=find(md.y>199999.9);
+connectivity=full(sparse(md.elements(:),1,1));
+while(time<500),
+	disp(['step ' num2str(i) '/' num2str(500/md.dt)]);
+	i=i+1;
+	time=time+md.dt;
+	md=solve(md,'analysis_type',PrognosticSolutionEnum);
+	thickness   =zeros(md.numberofgrids,1);
+	connectivity=sparse(md.results.PrognosticSolution.Thickness.index(:),1,1);
+	thickness   =sparse(md.results.PrognosticSolution.Thickness.index(:),1,md.results.PrognosticSolution.Thickness.value(:));
+	md.thickness=full(thickness./connectivity);
+	md.thickness(pos)=500+100*sin(2*pi*time/200)*ones(size(pos,1),1);
+	md.surface=md.bed+md.thickness;
+	md.dummy(i).thickness=md.thickness;
+	%plotmodel(md,'data',md.thickness,'sectionvalue','CrossLine.exp')
+end
+
+%plot results
+plotmodel(md,'data',md.thickness)
+%plotmodel(md,'data',[md.thickness 500+100*sin(2*pi/200*(500-md.y/400))],'sectionvalue','../Exp/CrossLineEISMINT.exp')
+
+%Fields and tolerances to track changes
+field_names     ={ ...
+	'Thickness', ...
+};
+field_tolerances={...
+	1e-13, ...
+};
+field_values={
+	md.thickness, ...
+};