Index: /issm/trunk/examples/IceStream/Stream.par
===================================================================
--- /issm/trunk/examples/IceStream/Stream.par	(revision 2896)
+++ /issm/trunk/examples/IceStream/Stream.par	(revision 2897)
@@ -24,5 +24,6 @@
 md.surface=zeros(md.numberofgrids,1);
 pos=find(md.x>=0);
-md.surface(pos)=50+4.5*sqrt(md.x(pos));
+md.surface(pos)=50+0.5*sqrt(md.x(pos));
+%md.surface(pos)=50+4.5*sqrt(md.x(pos));
 pos=find(md.x<0);
 md.surface(pos)=50;
@@ -30,7 +31,12 @@
 md.thickness=md.surface-md.bed;
 
+%Add extra slope 
+pos=find(md.x>0);
+md.bed(pos)=md.bed(pos)+500*betat*md.x(pos);
+md.surface(pos)=md.surface(pos)+500*betat*md.x(pos);
+
 disp('      creating drag');
 md.drag_type=2; %0 none 1 plastic 2 viscous
-md.drag=150*ones(md.numberofgrids,1); %q=1.
+md.drag=5000*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;
@@ -45,9 +51,10 @@
 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
 md.temperature=(273-20)*ones(md.numberofgrids,1);
+md.pressure=md.g*md.rho_ice*(md.surface-md.z); %lithostatique
 
 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.B(find(in))=md.B(find(in))/3;
 md.n=3*ones(md.numberofelements,1);
 
@@ -56,8 +63,11 @@
 
 %Parallel options
-md.np=6;
+md.cluster='cosmos';
+md.time=60;
+md.queue='shortq';
+md.np=10;
 
 %dynamics
-md.dt=0.1; %1 year
+md.dt=0; %1 year
 md.ndt=1; 
 md.artificial_diffusivity=1;
Index: /issm/trunk/examples/IceStream/runme.m
===================================================================
--- /issm/trunk/examples/IceStream/runme.m	(revision 2896)
+++ /issm/trunk/examples/IceStream/runme.m	(revision 2897)
@@ -1,11 +1,11 @@
-cluster='astrid';
 
 md=model;
-md=mesh(md,'DomainOutline.exp',400);
+md=mesh(md,'DomainOutline.exp',500);
 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');
+md=extrude(md,10,1);
+md=setelementstype(md,'pattyn','all');
+%md=setelementstype(md,'macayeal','all');
+%md=solve(md,'analysis_type','transient');
+md=solve(md,'analysis_type','diagnostic');
+%md=solve(md,'analysis_type','steadystate');