Index: ../trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 19031)
+++ ../trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 19032)
@@ -230,9 +230,10 @@
 						basis,1,numnodes,0,
 						&Ke->values[0],1);
 			
+			
 			/*Transfer EPL part*/
 			transfer=GetHydrologyKMatrixTransfer(basalelement,gauss,sed_head_input,epl_head_input,thick_input,base_input);
-			D_scalar=transfer*gauss->weight*Jdet;
+			D_scalar=dt*transfer*gauss->weight*Jdet;
 			TripleMultiply(basis,numnodes,1,0,
 										 &D_scalar,1,1,0,
 										 basis,1,numnodes,0,
@@ -327,7 +328,7 @@
 			
 			/*Dealing with the epl part of the transfer term*/
 			transfer=GetHydrologyPVectorTransfer(basalelement,gauss,sed_head_input,epl_head_input,thick_input,base_input);
-			scalar = Jdet*gauss->weight*((water_head*epl_specificstoring*epl_thickness)+(transfer));
+			scalar = Jdet*gauss->weight*((water_head*epl_specificstoring*epl_thickness)+(dt*transfer));
 			for(int i=0;i<numnodes;i++)pe->values[i]+=scalar*basis[i];
 		}
 	}
@@ -537,7 +538,7 @@
 		
 		switch(domaintype){
 		case Domain2DhorizontalEnum:
-			if(!element->IsOnBase()) return;			
+			if(!element->IsOnBase()) return;	
 			B = element->GetMaterialParameter(MaterialsRheologyBbarEnum);
 			break;
 		case Domain3DEnum:
@@ -596,7 +597,8 @@
 				EPLgrad2 = (epl_slopeX[i]*epl_slopeX[i])+(epl_slopeY[i]*epl_slopeY[i]);
 				
 				/*And proceed to the real thing*/
-				thickness[i] = old_thickness[i]*(1+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2-2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n)));
+				thickness[i] = old_thickness[i]*(1.0+((rho_water*gravity*dt)/(rho_ice*latentheat))*epl_conductivity*EPLgrad2-
+																				 (2.0*(A*dt/(pow(n,n)))*(pow(EPL_N,n))));
 				
 				/*Take care of otherthikening*/
 				if(thickness[i]>max_thick){
Index: ../trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 19031)
+++ ../trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 19032)
@@ -250,7 +250,7 @@
 				if(active_element){
 					transfer=GetHydrologyKMatrixTransfer(basalelement,gauss,sed_head_input,epl_head_input,thick_input,base_input);
 					basalelement->NodalFunctions(&basis[0],gauss);
-					D_scalar=transfer*gauss->weight*Jdet;
+					D_scalar=dt*transfer*gauss->weight*Jdet;
 					TripleMultiply(basis,numnodes,1,0,
 												 &D_scalar,1,1,0,
 												 basis,1,numnodes,0,
@@ -352,7 +352,7 @@
 				else{
 					transfer=0.0;
 				}
-				scalar = Jdet*gauss->weight*((water_head*sediment_storing)+(transfer));
+				scalar = Jdet*gauss->weight*((water_head*sediment_storing)+(dt*transfer));
 				for(int i=0;i<numnodes;i++)pe->values[i]+=scalar*basis[i];
 			}
 			else{
