Index: ../trunk-jpl/src/c/cores/controladm1qn3_core.cpp
===================================================================
--- ../trunk-jpl/src/c/cores/controladm1qn3_core.cpp	(revision 22793)
+++ ../trunk-jpl/src/c/cores/controladm1qn3_core.cpp	(revision 22794)
@@ -105,6 +105,7 @@
 			if(X[index]>XU[index]) X[index]=XU[index];
 			if(X[index]<XL[index]) X[index]=XL[index];
 		}
+		N_add+=N[c];
 	}
 
 	/*Start Tracing*/
@@ -319,6 +320,7 @@
 			X[index] = X[index]/reCast<double>(scaling_factors[c]);
 			Gnorm += G[index]*G[index];
 		}
+		N_add+=N[c];
 	}
 	Gnorm = sqrt(Gnorm);
 
Index: ../trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 22793)
+++ ../trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 22794)
@@ -1661,7 +1661,7 @@
             IssmDouble* times = xNew<IssmDouble>(N);
             for(int t=0;t<N;t++) times[t] = vector[(M-1)*N+t];
 				/*Create the three transient inputs for the control input*/
-            TransientInput* values_input=new TransientInput(ControlInputValuesEnum,times,N);
+            TransientInput* values_input=new TransientInput(input_enum,times,N);
 				TransientInput* mins_input = new TransientInput(ControlInputMinsEnum,times,N);
 				TransientInput* maxs_input = new TransientInput(ControlInputMaxsEnum,times,N);
 			   TransientInput* grad_input = new TransientInput(ControlInputGradEnum);
@@ -1673,7 +1673,7 @@
 					 } 
 					switch(this->ObjectEnum()){
                     case TriaEnum:
-									values_input->AddTimeInput(new TriaInput(ControlInputValuesEnum,values,P1Enum)); 
+									values_input->AddTimeInput(new TriaInput(input_enum,values,P1Enum)); 
 									mins_input->AddTimeInput(new TriaInput(ControlInputMinsEnum,values_min,P1Enum)); 
 									maxs_input->AddTimeInput(new TriaInput(ControlInputMaxsEnum,values_max,P1Enum));
 								break;
Index: ../trunk-jpl/src/c/classes/Inputs/ControlInput.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Inputs/ControlInput.cpp	(revision 22793)
+++ ../trunk-jpl/src/c/classes/Inputs/ControlInput.cpp	(revision 22794)
@@ -329,7 +329,7 @@
 	if(input->ObjectEnum()!=TransientInputEnum)_error_("cannot have timeoffset argument if not TransientInput Control");
 	TransientInput* transient_input = xDynamicCast<TransientInput*>(input);
 	IssmDouble time = transient_input->GetTimeByOffset(timeoffset);
-	TransientInput* new_trans_input = new TransientInput(ControlInputValuesEnum);
+	TransientInput* new_trans_input = new TransientInput(this->enum_type);
 	for(int i=0;i<transient_input->numtimesteps;i++){
 		if(transient_input->timesteps[i]==time) new_trans_input->AddTimeInput(xDynamicCast<TriaInput*>(in_input),time);
 		else {
Index: ../trunk-jpl/src/c/shared/io/Marshalling/IoCodeConversions.cpp
===================================================================
--- ../trunk-jpl/src/c/shared/io/Marshalling/IoCodeConversions.cpp	(revision 22793)
+++ ../trunk-jpl/src/c/shared/io/Marshalling/IoCodeConversions.cpp	(revision 22794)
@@ -142,6 +142,11 @@
 		input_enum        = CalvingStressThresholdGroundediceEnum;
 		fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1));
 	}
+	else if(strcmp(string_in,"DamageDbar")==0){
+		const char* field = "md.damage.D";
+		input_enum        = DamageDbarEnum;
+		fieldname=xNew<char>((strlen(field)+1)); xMemCpy<char>(fieldname,field,(strlen(field)+1));
+	}
 	else{
 		_error_("Field \""<<string_in<<"\" not supported yet");
 	}