Index: /issm/trunk/examples/IceStream/DomainOutline.exp
===================================================================
--- /issm/trunk/examples/IceStream/DomainOutline.exp	(revision 2461)
+++ /issm/trunk/examples/IceStream/DomainOutline.exp	(revision 2461)
@@ -0,0 +1,10 @@
+## Name:DomainOutline
+## Icon:0
+# Points Count  Value
+5 1.000000
+# X pos Y pos
+20000 -5000
+20000 5000
+-10000 5000
+-10000 -5000
+20000 -5000
Index: /issm/trunk/examples/IceStream/Front.exp
===================================================================
--- /issm/trunk/examples/IceStream/Front.exp	(revision 2461)
+++ /issm/trunk/examples/IceStream/Front.exp	(revision 2461)
@@ -0,0 +1,10 @@
+## Name:DomainOutline
+## Icon:0
+# Points Count  Value
+5 1.000000
+# X pos Y pos
+-9000 -5000
+-9000 5000
+-11000 5000
+-11000 -5000
+-9000 -5000
Index: /issm/trunk/examples/IceStream/IceShelf.exp
===================================================================
--- /issm/trunk/examples/IceStream/IceShelf.exp	(revision 2461)
+++ /issm/trunk/examples/IceStream/IceShelf.exp	(revision 2461)
@@ -0,0 +1,18 @@
+## Name:IceShelf
+## Icon:0
+# Points Count  Value
+12 1.000000
+# X pos Y pos
+510.5667474181 -1726.4530739959
+470.2122912151 347.5327616485
+550.9212036211 1928.2462904369
+-7519.9700369756 1804.9282137229
+-11010.6304985337 2219.7253808518
+-10808.8582175188 257.8468876747
+-10849.2126737218 -1311.6559068670
+-11172.0483233457 -2410.3078630462
+-9335.9205661099 -1939.4570246837
+-6470.7541756981 -2107.6180383846
+-2475.6630116027 -1939.4570246837
+510.5667474181 -1726.4530739959
+
Index: /issm/trunk/examples/IceStream/Stream.par
===================================================================
--- /issm/trunk/examples/IceStream/Stream.par	(revision 2461)
+++ /issm/trunk/examples/IceStream/Stream.par	(revision 2461)
@@ -0,0 +1,64 @@
+%Start defining model parameters here
+
+disp('      creating surface and bed');
+
+Hm=1600;
+xt=15000;
+xr=7000;
+betat=2*10^-4;
+betar=-8*10^-5;
+mu=7*10^-4;
+
+md.bed=zeros(md.numberofgrids,1);
+pos=find(md.x<xr);
+md.bed(pos)=1+betar*(xr-md.x(pos));
+pos=find(md.x<xt & md.x>=xr);
+md.bed(pos)=1;
+pos=find(md.x>=xt);
+md.bed(pos)=exp(-(betat*(xt-md.x(pos))).^2);
+pos=find(md.x<0);
+md.bed(pos)=1+betar*xr;
+
+md.bed=-md.bed.*exp(-(mu*md.y).^2)*Hm;
+
+md.surface=zeros(md.numberofgrids,1);
+pos=find(md.x>=0);
+md.surface(pos)=50+4.5*sqrt(md.x(pos));
+pos=find(md.x<0);
+md.surface(pos)=50;
+
+md.thickness=md.surface-md.bed;
+
+disp('      creating drag');
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=150*ones(md.numberofgrids,1); %q=1.
+in=ContourToMesh(md.elements,md.x,md.y,expread('drag.exp',1),'node',1);
+md.drag(find(in))=10;
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+disp('      creating temperature');
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.temperature=(273-20)*ones(md.numberofgrids,1);
+
+disp('      creating flow law paramter');
+md.B=5.2738*10^7*ones(md.numberofgrids,1);
+in=ContourToMesh(md.elements,md.x,md.y,expread('WeakB.exp',1),'node',1);
+md.B(find(in))=md.B(find(in))/3;
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetIceShelfBC(md,'Front.exp');
+
+%Parallel options
+md.np=6;
+
+%dynamics
+md.dt=0.1; %1 year
+md.ndt=1; 
+md.artificial_diffusivity=1;
+
Index: /issm/trunk/examples/IceStream/WeakB.exp
===================================================================
--- /issm/trunk/examples/IceStream/WeakB.exp	(revision 2461)
+++ /issm/trunk/examples/IceStream/WeakB.exp	(revision 2461)
@@ -0,0 +1,29 @@
+## Name:WeakB
+## Icon:0
+# Points Count  Value
+8 1.000000
+# X pos Y pos
+-9799.9968124442 2522.4152055134
+-599.1807981640 2881.1587014087
+2124.7449955374 2163.6717096182
+1680.8459773046 986.5446137119
+-4089.8412597220 1154.7056274128
+-9255.2116537039 1188.3378301530
+-11958.9602193038 1939.4570246837
+-9799.9968124442 2522.4152055134
+
+## Name:WeakB
+## Icon:0
+# Points Count  Value
+9 1.000000
+# X pos Y pos
+-8327.0591610353 -1737.6638082426
+6.1360448808 -1479.8169205679
+2407.2261889583 -1647.9779342688
+1822.0865740150 -2668.1547507209
+-3545.0561009818 -2993.2660438759
+-7943.6918271070 -2869.9479671619
+-10788.6809894173 -2612.1010794872
+-11999.3146755068 -1816.1389479697
+-8327.0591610353 -1737.6638082426
+
Index: /issm/trunk/examples/IceStream/drag.exp
===================================================================
--- /issm/trunk/examples/IceStream/drag.exp	(revision 2461)
+++ /issm/trunk/examples/IceStream/drag.exp	(revision 2461)
@@ -0,0 +1,20 @@
+## Name:drag
+## Icon:0
+# Points Count  Value
+14 1.000000
+# X pos Y pos
+19921.0601810532 1345.2881096072
+14876.7531556802 1322.8666411137
+7007.6341960984 1390.1310465941
+-2152.8273619788 1558.2920602950
+-9416.6294785159 1334.0773753604
+-10586.9087084024 1210.7592986464
+-10385.1364273875 -1468.6061863211
+-5441.7155425220 -1647.9779342688
+450.0350631136 -1939.4570246837
+5938.2411067194 -1569.5027945417
+11870.3461685579 -1692.8208712557
+17358.5522121637 -1502.2383890613
+20829.0354456203 -1580.7135287884
+19921.0601810532 1345.2881096072
+
Index: /issm/trunk/examples/IceStream/runme.m
===================================================================
--- /issm/trunk/examples/IceStream/runme.m	(revision 2461)
+++ /issm/trunk/examples/IceStream/runme.m	(revision 2461)
@@ -0,0 +1,11 @@
+cluster='astrid';
+
+md=model;
+md=mesh(md,'DomainOutline.exp',400);
+md=geography(md,'IceShelf.exp','');
+md=parameterize(md,'Stream.par');
+%md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster=cluster;
+md=solve(md,'analysis_type','transient');
+%md=solve(md,'analysis_type','diagnostic');