Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 6348)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 6349)
@@ -575,4 +575,7 @@
 	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]));
+
+	/*time: */
+	this->inputs->AddInput(new DoubleInput(DtEnum,iomodel->dt));
 
 	/*Control Inputs*/
@@ -2881,5 +2884,5 @@
 	heatcapacity=matpar->GetHeatCapacity();
 	thermalconductivity=matpar->GetThermalConductivity();
-	this->parameters->FindParam(&dt,DtEnum);
+	this->inputs->GetParameterValue(&dt,DtEnum);
 	this->parameters->FindParam(&artdiff,ArtDiffEnum);
 	this->parameters->FindParam(&epsvel,EpsVelEnum);
@@ -3783,5 +3786,5 @@
 	heatcapacity=matpar->GetHeatCapacity();
 	thermalconductivity=matpar->GetThermalConductivity();
-	this->parameters->FindParam(&dt,DtEnum);
+	this->inputs->GetParameterValue(&dt,DtEnum);
 	this->parameters->FindParam(&artdiff,ArtDiffEnum);
 	Input* vx_input=inputs->GetInput(VxEnum); ISSMASSERT(vx_input);
@@ -5332,4 +5335,5 @@
 				name==SurfaceEnum ||
 				name==BedEnum ||
+				name==DtEnum ||
 				name==SurfaceSlopeXEnum ||
 				name==SurfaceSlopeYEnum ||
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 6348)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 6349)
@@ -642,4 +642,7 @@
 	if (iomodel->elementonsurface) this->inputs->AddInput(new BoolInput(ElementOnSurfaceEnum,(IssmBool)iomodel->elementonsurface[index]));
 
+	/*time: */
+	this->inputs->AddInput(new DoubleInput(DtEnum,iomodel->dt));
+
 	/*Control Inputs*/
 	if (iomodel->control_analysis && iomodel->control_type){
@@ -3257,5 +3260,5 @@
 	/*Retrieve all inputs and parameters*/
 	GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
-	this->parameters->FindParam(&dt,DtEnum);
+	this->inputs->GetParameterValue(&dt,DtEnum);
 	this->parameters->FindParam(&dim,DimEnum);
 	this->parameters->FindParam(&artdiff,ArtDiffEnum);
@@ -3376,5 +3379,5 @@
 	/*Retrieve all inputs and parameters*/
 	GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
-	this->parameters->FindParam(&dt,DtEnum);
+	this->inputs->GetParameterValue(&dt,DtEnum);
 	this->parameters->FindParam(&dim,DimEnum);
 	Input* vxaverage_input=NULL;
@@ -3498,5 +3501,5 @@
 	/*Retrieve all inputs and parameters*/
 	GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
-	this->parameters->FindParam(&dt,DtEnum);
+	this->inputs->GetParameterValue(&dt,DtEnum);
 	mixed_layer_capacity=matpar->GetMixedLayerCapacity();
 	thermal_exchange_velocity=matpar->GetThermalExchangeVelocity();
@@ -4503,5 +4506,5 @@
 	rho_ice=matpar->GetRhoIce();
 	heatcapacity=matpar->GetHeatCapacity();
-	this->parameters->FindParam(&dt,DtEnum);
+	this->inputs->GetParameterValue(&dt,DtEnum);
 	Input* vx_input=inputs->GetInput(VxEnum);                         ISSMASSERT(vx_input);
 	Input* vy_input=inputs->GetInput(VyEnum);                         ISSMASSERT(vy_input);
@@ -5092,4 +5095,5 @@
 				name==SurfaceSlopeYEnum ||
 				name==MeltingRateEnum ||
+				name==DtEnum ||
 				name==AccumulationRateEnum ||
 				name==SurfaceAreaEnum||