Index: /issm/trunk-jpl/test/NightlyRun/test332.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test332.m	(revision 14747)
+++ /issm/trunk-jpl/test/NightlyRun/test332.m	(revision 14748)
@@ -5,4 +5,5 @@
 md.cluster=generic('name',oshostname(),'np',2);
 md.hydrology=(hydrologydc);
+md.hydrology.isefficientlayer=0;
 md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/NightlyRun/tmp332.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/tmp332.m	(revision 14747)
+++ /issm/trunk-jpl/test/NightlyRun/tmp332.m	(revision 14748)
@@ -3,9 +3,16 @@
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 md=setflowequation(md,'macayeal','all');
+md.cluster=generic('name',oshostname(),'np',2);
 md.hydrology=(hydrologydc);
-md.cluster=generic('name',oshostname(),'np',3);
 md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
-pos=find(md.mesh.y==50000.);
+pos=find(md.mesh.x==0);
 md.hydrology.spcsediment_head(pos)=0.0;
+md.basalforcings.melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
+md.timestepping.time_step=0;
+md.timestepping.final_time=1.0;
 md=solve(md,HydrologySolutionEnum());
+
+analitic=(md.mesh.x.^2-2*md.mesh.x*10^6)*(-2.0/(365*3600*24*2.0*md.hydrology.sediment_transmitivity));
+plotmodel(md,'data',md.results.HydrologySolution.SedimentHead-analitic)
+%plotmodel(md,'data',analitic)