Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15669)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15670)
@@ -6208,5 +6208,5 @@
 			switch(name){
 				case ThicknessEnum:
-					/*Update thickness + surface: assume bed is constant. On ice shelves, takes hydrostatic equilibrium {{{*/
+					/*Update thickness + surface: assume bed is constant. On ice shelves, takes hydrostatic equilibrium*/
 					IssmDouble  thickness[6];
 					IssmDouble  thickness_init[6];
@@ -6274,5 +6274,4 @@
 					this->inputs->AddInput(new PentaInput(SurfaceEnum,surface,P1Enum));
 
-					/*}}}*/
 					break;
 				default:
@@ -6300,10 +6299,8 @@
 void  Penta::InputUpdateFromMatrixDakota(IssmDouble* matrix, int nrows, int ncols, int name, int type){
 
-	int i,j,t;
-	TransientInput* transientinput=NULL;
-	IssmDouble values[6];
-	IssmDouble time;
-	int row;
-	IssmDouble yts;
+	int             i,t,row;
+	IssmDouble      time;
+	TransientInput *transientinput = NULL;
+	IssmDouble      values[6];
 
 	/*Check that name is an element input*/
@@ -6313,9 +6310,5 @@
 
 		case VertexEnum:
-
 			/*Create transient input: */
-
-			parameters->FindParam(&yts,ConstantsYtsEnum);
-
 			for(t=0;t<ncols;t++){ //ncols is the number of times
 
@@ -6323,9 +6316,9 @@
 				for(i=0;i<6;i++){
 					row=this->vertices[i]->Sid();
-					values[i]=(IssmDouble)matrix[ncols*row+t];
+					values[i]=matrix[ncols*row+t];
 				}
 
-				/*time? :*/
-				time=(IssmDouble)matrix[(nrows-1)*ncols+t]*yts;
+				/*time:*/
+				time=matrix[(nrows-1)*ncols+t];
 
 				if(t==0) transientinput=new TransientInput(name);
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15669)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15670)
@@ -6819,10 +6819,8 @@
 void  Tria::InputUpdateFromMatrixDakota(IssmDouble* matrix, int nrows, int ncols, int name, int type){
 
-	int i,j,t;
-	TransientInput* transientinput=NULL;
-	IssmDouble values[3];
-	IssmDouble time;
-	int row;
-	IssmDouble yts;
+	int             i,t,row;
+	IssmDouble      time;
+	TransientInput *transientinput = NULL;
+	IssmDouble      values[3];
 
 	/*Check that name is an element input*/
@@ -6832,8 +6830,5 @@
 
 		case VertexEnum:
-
 			/*Create transient input: */
-
-			parameters->FindParam(&yts,ConstantsYtsEnum);
 			for(t=0;t<ncols;t++){ //ncols is the number of times
 
@@ -6841,9 +6836,9 @@
 				for(i=0;i<3;i++){
 					row=this->vertices[i]->Sid();
-					values[i]=(IssmDouble)matrix[ncols*row+t];
+					values[i]=matrix[ncols*row+t];
 				}
 
-				/*time? :*/
-				time=(IssmDouble)matrix[(nrows-1)*ncols+t]*yts;
+				/*time:*/
+				time=matrix[(nrows-1)*ncols+t];
 
 				if(t==0) transientinput=new TransientInput(name);