Index: /issm/trunk/src/c/objects/Elements/Beam.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Beam.cpp	(revision 3706)
+++ /issm/trunk/src/c/objects/Elements/Beam.cpp	(revision 3707)
@@ -142,5 +142,4 @@
 	}	
 	if (iomodel->gridonbed) this->inputs->AddInput(new BoolInput(ElementOnBedEnum,(IssmBool)iomodel->gridonbed[index]));
-	if (iomodel->dt) this->inputs->AddInput(new DoubleInput(DtEnum,iomodel->dt));
 
 	//this->parameters: we still can't point to it, it may not even exist. Configure will handle this.
Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3706)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3707)
@@ -126,4 +126,8 @@
 		this->inputs->AddInput(new PentaVertexInput(TemperatureEnum,nodeinputs));
 	}
+	if (iomodel->dhdt) {
+		for(i=0;i<6;i++)nodeinputs[i]=iomodel->dhdt[penta_node_ids[i]-1];
+		this->inputs->AddInput(new PentaVertexInput(DhdtEnum,nodeinputs));
+	}
 
 	/*vx,vy and vz: */
@@ -161,5 +165,4 @@
 	if (iomodel->elementonwater) this->inputs->AddInput(new BoolInput(ElementOnWaterEnum,(IssmBool)iomodel->elementonwater[index]));
 	if (iomodel->elementonsurface) this->inputs->AddInput(new BoolInput(ElementOnSurfaceEnum,(IssmBool)iomodel->elementonsurface[index]));
-	if (iomodel->dt) this->inputs->AddInput(new DoubleInput(DtEnum,iomodel->dt));
 
 	//elements of type 3 are macayeal type penta. we collapse the formulation on their base.
Index: /issm/trunk/src/c/objects/Elements/Sing.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Sing.cpp	(revision 3706)
+++ /issm/trunk/src/c/objects/Elements/Sing.cpp	(revision 3707)
@@ -94,6 +94,4 @@
 	if (iomodel->vz_obs) this->inputs->AddInput(new SingVertexInput(VzEnum,iomodel->vz_obs[i]));
 	if (iomodel->weights) this->inputs->AddInput(new SingVertexInput(WeightsEnum,iomodel->weights[i]));
-
-	if (iomodel->dt) this->inputs->AddInput(new DoubleInput(DtEnum,iomodel->dt));
 
 	//this->parameters: we still can't point to it, it may not even exist. Configure will handle this.
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 3706)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 3707)
@@ -117,5 +117,5 @@
 		if (iomodel->drag_q) this->inputs->AddInput(new DoubleInput(DragQEnum,iomodel->drag_q[index]));
 		this->inputs->AddInput(new IntInput(DragTypeEnum,iomodel->drag_type));
-
+s
 	}
 	if (iomodel->melting_rate) {
@@ -130,4 +130,8 @@
 		for(i=0;i<3;i++)nodeinputs[i]=iomodel->geothermalflux[tria_node_ids[i]-1];
 		this->inputs->AddInput(new TriaVertexInput(GeothermalFluxEnum,nodeinputs));
+	}
+	if (iomodel->dhdt) {
+		for(i=0;i<3;i++)nodeinputs[i]=iomodel->dhdt[tria_node_ids[i]-1];
+		this->inputs->AddInput(new TriaVertexInput(DhDtEnum,nodeinputs));
 	}
 	if (iomodel->pressure) {
@@ -174,5 +178,4 @@
 	if (iomodel->elementonwater) this->inputs->AddInput(new BoolInput(ElementOnWaterEnum,(IssmBool)iomodel->elementonwater[index]));
 	if (iomodel->elementonsurface) this->inputs->AddInput(new BoolInput(ElementOnSurfaceEnum,(IssmBool)iomodel->elementonsurface[index]));
-	if (iomodel->dt) this->inputs->AddInput(new DoubleInput(DtEnum,iomodel->dt));
 
 	//this->parameters: we still can't point to it, it may not even exist. Configure will handle this.
