Index: /issm/trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 25511)
+++ /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 25512)
@@ -357,5 +357,5 @@
 
 	for(int i=0;i<Xsize;i++) Xb[i] += X[i].gradient();
-	InitAdjointFunctor* marshallhandle3 = new InitAdjointFunctor(Yb);
+	SetAdjointFunctor* marshallhandle3 = new SetAdjointFunctor(Yb);
 	femmodel->Marshall(marshallhandle3);
 	delete marshallhandle3;
Index: /issm/trunk-jpl/src/c/shared/io/Marshalling/Marshalling.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/io/Marshalling/Marshalling.h	(revision 25511)
+++ /issm/trunk-jpl/src/c/shared/io/Marshalling/Marshalling.h	(revision 25512)
@@ -20,5 +20,4 @@
 	AD_REGISTERINPUT,
 	AD_REGISTEROUTPUT,
-	AD_INITADJOINT,
 	AD_SETADJOINT,
 #endif
@@ -212,40 +211,4 @@
 		}
 }; /*}}}*/
-class InitAdjointFunctor:     public MarshallHandle{ /*{{{*/
-
-	private:
-		int                   double_count;
-		IssmDouble::TapeType* tape_codi;
-		double*               adjoint;
-
-	public:
-		InitAdjointFunctor(double* adjoint_in) : MarshallHandle(AD_INITADJOINT){
-			this->double_count = 0;
-			this->tape_codi    = &(IssmDouble::getGlobalTape());
-			this->adjoint      = adjoint_in;
-		}
-		void Echo(void){
-			printf("InitAdjointFunctor Echo:\n");
-			printf("   double_count: %i\n",double_count);
-		}
-		template<typename T> void call(T & value){
-			/*General case: do nothing*/
-		}
-		template<typename T> void call(T* & value,int size){
-			/*General case: do nothing*/
-		}
-		void call(IssmDouble value){
-			value.gradient() = adjoint[this->double_count];
-			this->double_count++;
-		}
-		void call(IssmDouble* value,int size){
-			if(value){
-				for(int i=0;i<size;i++){
-					value[i].gradient() = adjoint[this->double_count];
-					this->double_count++;
-				}
-			}
-		}
-}; /*}}}*/
 class SetAdjointFunctor:      public MarshallHandle{ /*{{{*/
 
@@ -295,5 +258,4 @@
 		case AD_REGISTERINPUT: {RegisterInputFunctor*  temp = xDynamicCast<RegisterInputFunctor*>(this);   temp->call(value); break;}
 		case AD_REGISTEROUTPUT:{RegisterOutputFunctor* temp = xDynamicCast<RegisterOutputFunctor*>(this);  temp->call(value); break;}
-		case AD_INITADJOINT:   {InitAdjointFunctor*    temp = xDynamicCast<InitAdjointFunctor*>(this);     temp->call(value); break;}
 		case AD_SETADJOINT:    {SetAdjointFunctor*     temp = xDynamicCast<SetAdjointFunctor*>(this);      temp->call(value); break;}
 #endif
@@ -310,5 +272,4 @@
 		case AD_REGISTERINPUT: {RegisterInputFunctor*  temp = xDynamicCast<RegisterInputFunctor*>(this);   temp->call(value,size); break;}
 		case AD_REGISTEROUTPUT:{RegisterOutputFunctor* temp = xDynamicCast<RegisterOutputFunctor*>(this);  temp->call(value,size); break;}
-		case AD_INITADJOINT:   {InitAdjointFunctor*    temp = xDynamicCast<InitAdjointFunctor*>(this);     temp->call(value,size); break;}
 		case AD_SETADJOINT:    {SetAdjointFunctor*     temp = xDynamicCast<SetAdjointFunctor*>(this);      temp->call(value,size); break;}
 #endif