Index: /issm/trunk-jpl/examples/Inversion/Square.par
===================================================================
--- /issm/trunk-jpl/examples/Inversion/Square.par	(revision 27021)
+++ /issm/trunk-jpl/examples/Inversion/Square.par	(revision 27022)
@@ -11,5 +11,5 @@
 
 disp('      creating drag');
-md.friction.coefficient=200*ones(md.mesh.numberofvertices,1);
+md.friction.coefficient=50*ones(md.mesh.numberofvertices,1);
 md.friction.coefficient(find(md.mask.ocean_levelset<0.))=0.;
 md.friction.p = ones(md.mesh.numberofelements,1);
Index: /issm/trunk-jpl/examples/Inversion/runme.m
===================================================================
--- /issm/trunk-jpl/examples/Inversion/runme.m	(revision 27021)
+++ /issm/trunk-jpl/examples/Inversion/runme.m	(revision 27022)
@@ -1,3 +1,3 @@
-steps=[1];
+steps=[4];
 
 if any(steps==1) 
Index: /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 27021)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 27022)
@@ -382,4 +382,5 @@
 		dpressure_water[1] = rho_water*g*(dh[1] - dbed[1]);
 		PMPheat=CT*CW*conductivity*rho_water*(dh[0]*dpressure_water[0]+dh[1]*dpressure_water[1]);
+		PMPheat=0; /*** TEST no PMPheat***/
 
    	meltrate = 1/latentheat*(G+frictionheat+rho_water*g*conductivity*(dh[0]*dh[0]+dh[1]*dh[1])-PMPheat);
@@ -629,4 +630,5 @@
 		dpressure_water[1] = rho_water*g*(dh[1] - dbed[1]);
 		PMPheat=CT*CW*conductivity*rho_water*(dh[0]*dpressure_water[0]+dh[1]*dpressure_water[1]);
+PMPheat=0; /*** TEST no PMPheat***/
 		dissipation=rho_water*g*conductivity*(dh[0]*dh[0]+dh[1]*dh[1]);