Index: /issm/trunk/src/c/objects/ElementResults/NodalValuesUnitConversion.cpp
===================================================================
--- /issm/trunk/src/c/objects/ElementResults/NodalValuesUnitConversion.cpp	(revision 5291)
+++ /issm/trunk/src/c/objects/ElementResults/NodalValuesUnitConversion.cpp	(revision 5292)
@@ -27,8 +27,9 @@
 	/*Get scaling factor: */
 	switch(enum_type){
-		case VxEnum:  scale=yts;break; //m/yr
-		case VyEnum:  scale=yts;break; //m/yr
-		case VzEnum:  scale=yts;break; //m/yr
-		case VelEnum: scale=yts;break; //m/yr
+		case VxEnum:   scale=yts;break; //m/yr
+		case VyEnum:   scale=yts;break; //m/yr
+		case VzEnum:   scale=yts;break; //m/yr
+		case VelEnum:  scale=yts;break; //m/yr
+		case DhDtEnum: scale=yts;break; //m/yr
 		case MeltingRateEnum: scale=yts;break; //m/yr
 		case MisfitEnum:      scale=pow(yts,2);break; //(m/yr)^2
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 5291)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 5292)
@@ -1278,5 +1278,5 @@
 		for (i=0;i<numgrids;i++){
 			//standard gradient dJ/dki
-			grade_g_gaussian[i]=lambda*Jdet*gauss_weight*l1l2l3[i]; 
+			grade_g_gaussian[i]=-lambda*Jdet*gauss_weight*l1l2l3[i]; 
 		}
 
@@ -4401,8 +4401,10 @@
 	double  accumulation_g;
 	double  melting_g;
+	double  dhdt_g;
 
 	/*inputs: */
 	Input* accumulation_input=NULL;
 	Input* melting_input=NULL;
+	Input* dhdt_input=NULL;
 
 	/* Get node coordinates and dof list: */
@@ -4414,6 +4416,7 @@
 
 	/*retrieve inputs :*/
-	accumulation_input=inputs->GetInput(AccumulationRateEnum);
-	melting_input=inputs->GetInput(MeltingRateEnum);
+	accumulation_input=inputs->GetInput(AccumulationRateEnum); ISSMASSERT(accumulation_input);
+	melting_input=inputs->GetInput(MeltingRateEnum);           ISSMASSERT(melting_input);
+	dhdt_input=inputs->GetInput(DhDtEnum);                     ISSMASSERT(dhdt_input);
 	
 	/* Start  looping on the number of gaussian points: */
@@ -4434,7 +4437,8 @@
 		accumulation_input->GetParameterValue(&accumulation_g, &gauss_l1l2l3[0]);
 		melting_input->GetParameterValue(&melting_g, &gauss_l1l2l3[0]);
+		dhdt_input->GetParameterValue(&dhdt_g, &gauss_l1l2l3[0]);
 
 		/* Add value into pe_g: */
-		for( i=0; i<numdof; i++) pe_g[i]+=Jdettria*gauss_weight*(accumulation_g-melting_g)*L[i];
+		for( i=0; i<numdof; i++) pe_g[i]+=Jdettria*gauss_weight*(accumulation_g-melting_g-dhdt_g)*L[i];
 
 	} // for (ig=0; ig<num_gauss; ig++)
@@ -4520,5 +4524,5 @@
 
 		/* Add value into pe_g: */
-		for( i=0; i<numdof; i++) pe_g[i]+=Jdettria*gauss_weight*(accumulation_g-melting_g+dhdt_g)*L[i];
+		for( i=0; i<numdof; i++) pe_g[i]+=Jdettria*gauss_weight*(accumulation_g-melting_g-dhdt_g)*L[i];
 
 	} // for (ig=0; ig<num_gauss; ig++)
Index: /issm/trunk/src/c/objects/Loads/Numericalflux.h
===================================================================
--- /issm/trunk/src/c/objects/Loads/Numericalflux.h	(revision 5291)
+++ /issm/trunk/src/c/objects/Loads/Numericalflux.h	(revision 5292)
@@ -45,9 +45,9 @@
 		/*}}}*/
 		/*Update virtual functions resolution: {{{1*/
-		void    InputUpdateFromVector(double* vector, int name, int type){ISSMERROR("Not implemented yet!");}
+		void    InputUpdateFromVector(double* vector, int name, int type){/*Do nothing*/}
 		void    InputUpdateFromVector(int* vector, int name, int type){ISSMERROR("Not implemented yet!");}
 		void    InputUpdateFromVector(bool* vector, int name, int type){ISSMERROR("Not implemented yet!");}
 		void    InputUpdateFromConstant(double constant, int name){ISSMERROR("Not implemented yet!");}
-		void    InputUpdateFromConstant(int constant, int name){ISSMERROR("Not implemented yet!");}
+		void    InputUpdateFromConstant(int constant, int name){/*Do nothing*/};
 		void    InputUpdateFromConstant(bool constant, int name){ISSMERROR("Not implemented yet!");}
 		void    InputUpdateFromSolution(double* solution){ISSMERROR("Not implemented yet!");}