Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14809)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14810)
@@ -3083,5 +3083,4 @@
 
 	/*thickness averages: */
-	IssmDouble current_he=NULL;
 	IssmDouble* hes=NULL;
 	IssmDouble* times=NULL;
@@ -3106,4 +3105,7 @@
 	IssmDouble rho_ice;
 
+	/*constants: */
+	IssmDouble yts;
+
 	/*output: */
 	IssmDouble  wi;
@@ -3115,4 +3117,5 @@
 	/*how many dofs are we working with here? */
 	this->parameters->FindParam(&gsize,MeshNumberofverticesEnum);
+	this->parameters->FindParam(&yts,ConstantsYtsEnum);
 
 	/*recover gia solution parameters: */
@@ -3134,7 +3137,6 @@
 	thickness_input=inputs->GetInput(ThicknessEnum); 
 	if (!thickness_input)_error_("thickness input needed to compute gia deflection!");
-	thickness_input->GetInputAllTimeAverages(&hes,&times,&numtimes);
-	thickness_input->GetInputAverage(&current_he);
-
+	thickness_input->GetInputUpToCurrentTimeAverages(&hes,&times,&numtimes,currenttime);
+	
 	/*recover lithosphere thickness: */
 	lithosphere_thickness_input=inputs->GetInput(GiaLithosphereThicknessEnum);
@@ -3155,5 +3157,4 @@
 	/*start loading GiaDeflectionCore arguments: */
 	arguments.re=re;
-	arguments.current_he=current_he;
 	arguments.hes=hes;
 	arguments.times=times;
@@ -3170,4 +3171,5 @@
 	arguments.irate=output_rates;
 	arguments.iedge=cross_section_shape;
+	arguments.yts=yts;
 
 	for(i=0;i<gsize;i++){