Index: /issm/trunk/test/Validation/EISMINT/MassConservation/CrossLine2.exp
===================================================================
--- /issm/trunk/test/Validation/EISMINT/MassConservation/CrossLine2.exp	(revision 3580)
+++ /issm/trunk/test/Validation/EISMINT/MassConservation/CrossLine2.exp	(revision 3581)
@@ -4,4 +4,4 @@
 2 1.
 # X pos Y pos
-10000  199999
-110000 0
+0 0
+200000 100000
Index: /issm/trunk/test/Validation/EISMINT/MassConservation/runme2.m
===================================================================
--- /issm/trunk/test/Validation/EISMINT/MassConservation/runme2.m	(revision 3580)
+++ /issm/trunk/test/Validation/EISMINT/MassConservation/runme2.m	(revision 3581)
@@ -4,5 +4,5 @@
 %The goal is to test the prognostic model
 %md=bamg(model,'domain','DomainOutline.exp','hmax',4550);
-md=bamg(model,'domain','DomainOutline.exp','hmax',3000);
+md=bamg(model,'domain','DomainOutline.exp','hmax',2000);
 md=geography(md,'all','');
 md=parameterize(md,'Square.par');
@@ -15,4 +15,6 @@
 md.vx=160*sqrt(5)*ones(md.numberofgrids,1);
 md.vy= 80*sqrt(5)*ones(md.numberofgrids,1);
+sinus=1/5*sqrt(5);
+cosinus=2/5*sqrt(5);
 
 %0launch transient solution
@@ -20,5 +22,5 @@
 %FOR NOW:
 md.cluster=oshostname();
-md.np=10;
+md.np=14;
 i=0;
 time=0;
@@ -27,6 +29,6 @@
 posLEFT=find(md.x<0.01);
 connectivity=full(sparse(md.elements(:),1,1));
-while(time<100),
-	disp(['step ' num2str(i) '/' num2str(500/md.dt)]);
+while(time<600),
+	disp(['step ' num2str(i) '/' num2str(600/md.dt)]);
 	i=i+1;
 	time=time+md.dt;
@@ -35,6 +37,6 @@
 	%md.thickness=thickness./connectivity;
 	md.thickness=md.results.prognostic2.thickness;
-	md.thickness(posLEFT)=500+100*sin((2*pi/200)*(time-md.y(posLEFT)/md.vy(1)));
-	md.thickness(posDOWN)=500+100*sin((2*pi/200)*(time-md.x(posDOWN)/md.vx(1)));
+	md.thickness(posLEFT)=500+100*sin((2*pi/200)*(time-sinus*md.y(posLEFT)/sqrt(md.vx(1)^2+md.vy(1)^2)));
+	md.thickness(posDOWN)=500+100*sin((2*pi/200)*(time-cosinus*md.x(posDOWN)/sqrt(md.vx(1)^2+md.vy(1)^2)));
 	md.surface=md.bed+md.thickness;
 	md.dummy(i).thickness=md.thickness;
@@ -43,5 +45,5 @@
 
 %plot results
-plotmodel(md,'data',[md.thickness 500+100*sin(2*pi/200*(500-md.y/400))],'sectionvalue','CrossLine2.exp')
+plotmodel(md,'data',[md.thickness 500+100*sin(2*pi/200*(600-(md.x*cosinus+md.y*sinus)/sqrt(md.vx(1)^2+md.vy(1)^2)))],'sectionvalue','CrossLine2.exp')
 
 %Don't forget to add these lines in icetransient2d.m, just before the computation of the thickness to change the thickness on the upper boundary condition