Index: /issm/trunk/test/NightlyRun/test1201.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1201.m	(revision 9222)
+++ /issm/trunk/test/NightlyRun/test1201.m	(revision 9223)
@@ -4,5 +4,5 @@
 results={};
 
-for type=1:3;
+for stabilization=1:3;
 	%The goal is to test the prognostic model
 	md=bamg(model,'domain','../Exp/SquareEISMINT.exp','hmax',3000);
@@ -16,40 +16,29 @@
 	md.vy=-400*ones(md.numberofnodes,1);
 	md.vz=zeros(md.numberofnodes,1);
+	md.vel= 400*ones(md.numberofnodes,1);
 
-	%analytical test
-	md.spcvx=zeros(md.numberofnodes,1);
-	md.spcvy=zeros(md.numberofnodes,1);
-	md.spcvz=zeros(md.numberofnodes,1);
-	md.spcvx(:)=0;
-	md.spcvy(:)=-400;
-	md.spcvx(1)=0;
-	if type==1,
+	%Stabilization
+	if stabilization==1,
 		md.artificial_diffusivity=1;
 		md.prognostic_DG=0;
-	elseif type==2,
+	elseif stabilization==2,
 		md.artificial_diffusivity=0;
 		md.prognostic_DG=0;
-	elseif type==3,
+	elseif stabilization==3,
 		md.prognostic_DG=1;
 	end
 
-	%Launch many prognostic as we don't really want a transient
-	%FOR NOW:
-	i=0;
-	time=0;
-	md.dummy=struct();
+	%spc thickness
 	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,PrognosticSolutionEnum);
-		md.thickness=PatchToVec(md.results.PrognosticSolution.Thickness);
-		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;
-	end
-	results{type}=md.thickness;
+	times=0:1:500;
+	md.spcthickness=NaN*ones(md.numberofnodes+1,length(times));
+	md.spcthickness(end,:)=times;
+	md.spcthickness(pos,:)=repmat(500+100*sin(2*pi*times/200),length(pos),1);
+
+	%solve
+	md.isdiagnostic=0;
+	md.output_frequency=500; %keep only last step
+	md=solve(md,TransientSolutionEnum);
+	results{stabilization}=PatchToVec(md.results.TransientSolution(end).Thickness);
 end