Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 24053)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 24054)
@@ -2330,6 +2330,6 @@
 							ISSM_MPI_Bcast(&numchannels,1,ISSM_MPI_INT,0,IssmComm::GetComm());
 
-							IssmDouble* values    = xNewZeroInit<IssmDouble>(numchannels);
-							IssmDouble* allvalues = xNew<IssmDouble>(numchannels);
+							IssmPDouble* values    = xNewZeroInit<IssmPDouble>(numchannels);
+							IssmPDouble* allvalues = xNew<IssmPDouble>(numchannels);
 
 							/*Fill-in vector*/
@@ -2343,8 +2343,8 @@
 							/*Gather from all cpus*/
 							ISSM_MPI_Allreduce((void*)values,(void*)allvalues,numchannels,ISSM_MPI_PDOUBLE,ISSM_MPI_SUM,IssmComm::GetComm());
-							xDelete<IssmDouble>(values);
-
-							if(save_results)results->AddResult(new GenericExternalResult<IssmDouble*>(results->Size()+1,output_enum,allvalues,numchannels,1,step,time));
-							xDelete<IssmDouble>(allvalues);
+							xDelete<IssmPDouble>(values);
+
+							if(save_results)results->AddResult(new GenericExternalResult<IssmPDouble*>(results->Size()+1,output_enum,allvalues,numchannels,1,step,time));
+							xDelete<IssmPDouble>(allvalues);
 
 							isvec = true;
Index: /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp	(revision 24053)
+++ /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp	(revision 24054)
@@ -683,8 +683,8 @@
 }
 /*}}}*/
-void           Channel::WriteChannelCrossSection(IssmDouble* values){/*{{{*/
+void           Channel::WriteChannelCrossSection(IssmPDouble* values){/*{{{*/
 
 	_assert_(values);
-	values[this->sid] = this->S;
-}
-/*}}}*/
+	values[this->sid] = reCast<IssmPDouble>(this->S);
+}
+/*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Channel.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Channel.h	(revision 24053)
+++ /issm/trunk-jpl/src/c/classes/Loads/Channel.h	(revision 24054)
@@ -78,5 +78,5 @@
 		ElementVector* CreatePVectorHydrologyGlaDS(void);
 		ElementMatrix* CreateKMatrixHydrologyGlaDS(void);
-		void           WriteChannelCrossSection(IssmDouble* values);
+		void           WriteChannelCrossSection(IssmPDouble* values);
 		/*}}}*/