Index: /issm/trunk/test/NightlyRun/test1208.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1208.m	(revision 5176)
+++ /issm/trunk/test/NightlyRun/test1208.m	(revision 5176)
@@ -0,0 +1,23 @@
+%EISMINT benchmark experiment A
+numlayers=10;
+resolution=30000;
+
+%To begin with the numerical model
+md=model;
+md=roundmesh(md,750000,resolution);
+md=geography(md,'',''); %We can not test iceshelves nor ice rises with this analytical solution
+md=parameterize(md,'../Par/RoundSheetEISMINT.par');
+
+%We extrude the model to have a 3d model
+md=extrude(md,numlayers,1);
+md=setelementstype(md,'hutter','all');
+
+%Spc the grids on the bed
+pos=find(md.gridonbed);
+md.spcvelocity(pos,1:3)=1;
+
+%Adapt the time steps to the resolution
+md.dt=25;
+
+%Now we can solve the problem 
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
Index: /issm/trunk/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 5176)
+++ /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 5176)
@@ -0,0 +1,62 @@
+%Ok, start defining model parameters here
+disp('      creating thickness');
+md.thickness=10*ones(md.numberofgrids,1);
+md.bed=zeros(md.numberofgrids,1);
+md.surface=md.bed+md.thickness;
+
+disp('      creating drag');
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag_coefficient=20*ones(md.numberofgrids,1); %q=1. %no drag is specified in the analytical solution
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+
+disp('      creating temperatures');
+tmin=238.15; %K
+st=1.67*10^-2/1000; %k/m;
+radius=sqrt(md.x.^2+md.y.^2);
+md.observed_temperature=(tmin+st*radius);
+md.temperature=md.observed_temperature;
+md.geothermalflux=4.2*10^-2*ones(md.numberofgrids,1);
+
+disp('      creating flow law paramter');
+md.rheology_B=6.81*10^(7)*ones(md.numberofgrids,1); %to have the same B as the analytical solution 
+md.rheology_n=3*ones(md.numberofelements,1);
+
+disp('      creating accumulation rates');
+acc_max=0.5; %m/yr
+sb=10^-2/1000; %m/yr/m
+rel=450*1000; %m
+md.accumulation_rate=min(acc_max,sb*(rel-radius));
+
+disp('      creating velocities');
+constant=0.3;
+md.vx_obs=constant/2*md.x.*(md.thickness).^-1;
+md.vy_obs=constant/2*md.y.*(md.thickness).^-1;
+md.vel_obs=(sqrt((md.vx_obs).^2+(md.vy_obs).^2));
+md.vx=zeros(md.numberofgrids,1);
+md.vy=zeros(md.numberofgrids,1);
+md.vz=zeros(md.numberofgrids,1);
+md.pressure=zeros(md.numberofgrids,1);
+
+%Deal with boundary conditions:
+disp('      boundary conditions for diagnostic model: ');
+md=SetMarineIceSheetBC(md,'../Exp/RoundFrontEISMINT.exp');
+
+radius=sqrt((md.x).*md.x+(md.y).*md.y);
+pos=find(radius==min(radius));
+md.x(pos)=0; md.y(pos)=0; %the closest node to the center is changed to be exactly at the center
+
+md.spcvelocity(pos,1:3)=1;
+md.spcvelocity(pos,4:6)=0;
+
+%parallel options
+md.cluster=oshostname();
+md.np=8;
+md.ndt=50000;
+
+%Constants
+md.rho_ice=910;
+md.thermalconductivity=2.1;
+md.latentheat=3.35*10^5;
+md.beta=8.66*10^4/(md.rho_ice*md.g); %conversion from K/m to K/Pa
+md.yts=31556926;