Index: /issm/trunk/test/NightlyRun/test1109.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1109.m	(revision 5178)
+++ /issm/trunk/test/NightlyRun/test1109.m	(revision 5178)
@@ -0,0 +1,30 @@
+%This test is a test from the ISMP-HOM Intercomparison project
+%TestE 
+Lx=10; %in m
+Ly=5000; %in m
+nx=3; %number of nodes in x direction
+ny=51;
+md=model;
+md=squaremesh(md,Lx,Ly,nx,ny);
+md=geography(md,'',''); %ice sheet test
+md=parameterize(md,'../Par/ISMIPE.par');
+md=extrude(md,10,1);
+
+md=setelementstype(md,'pattyn','all');
+%md=setelementstype(md,'pattyn','all','stokes','all'); 
+
+%md.gridonstokes(find(md.x==0 | md.x==Lx))=0;
+%%Create MPCs to have periodic boundary conditions
+%posx=find(md.x==0);
+%posx2=find(md.x==max(md.x));
+%
+%md.penalties=[posx,posx2];
+
+%Compute the diagnostic
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
+
+vx=PatchToVec(md.results.DiagnosticSolution.Vx);
+vy=PatchToVec(md.results.DiagnosticSolution.Vy);
+vz=PatchToVec(md.results.DiagnosticSolution.Vz);
+
+plotmodel(md,'data',vx,'data',vy,'data',vz,'layer#all',md.numlayers)
Index: /issm/trunk/test/Par/ISMIPE.par
===================================================================
--- /issm/trunk/test/Par/ISMIPE.par	(revision 5178)
+++ /issm/trunk/test/Par/ISMIPE.par	(revision 5178)
@@ -0,0 +1,39 @@
+%Ok, start defining model parameters here
+
+disp('      creating thickness');
+data=load('../Data/ISMIPE.data','-mat');
+data=data.data;
+md.surface=zeros(md.numberofgrids,1);
+md.bed=zeros(md.numberofgrids,1);
+for i=1:md.numberofgrids
+	y=md.y(i);
+	point1=floor(y/100)+1;
+	point2=min(point1+1,51);
+	coeff=(y-(point1-1)*100)/100;
+	md.bed(i)=(1-coeff)*data(point1,2)+coeff*data(point2,2);
+	md.surface(i)=(1-coeff)*data(point1,3)+coeff*data(point2,3);
+end
+md.thickness=md.surface-md.bed;
+md.firn_layer=0*ones(md.numberofgrids,1);
+md.thickness(find(~md.thickness))=0.01;
+md.bed=md.surface-md.thickness;
+
+disp('      creating drag');
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag_coefficient=zeros(md.numberofgrids,1);
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+
+disp('      creating flow law paramter');
+md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
+md.rheology_n=3*ones(md.numberofelements,1);
+
+disp('      boundary conditions for diagnostic model: ');
+%Create grid on boundary fist (because wi can not use mesh)
+md=SetIceSheetBC(md);
+
+md.np=3;
+md.cluster=oshostname();