Index: /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 24131)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyGlaDSAnalysis.cpp	(revision 24132)
@@ -29,6 +29,9 @@
 
 	/*Add channels?*/
+	int K,L;
 	bool ischannels;
+	IssmDouble* channelarea;
 	iomodel->FindConstant(&ischannels,"md.hydrology.ischannels");
+	iomodel->FetchData(&channelarea,&K,&L,"md.initialization.channelarea");
 	if(ischannels){
 		/*Get faces (edges in 2d)*/
@@ -42,6 +45,15 @@
 			if(!iomodel->my_elements[element]) continue;
 
-			/* Add load */
-			loads->AddObject(new Channel(i+1,i,i,iomodel));
+			/*Add channelarea from initialization if exists*/
+			if(K!=0 && K!=iomodel->numberoffaces){
+				_error_("Unknown dimension for channel area initialization.");
+			}
+			if(K==0){
+				loads->AddObject(new Channel(i+1,0.,i,iomodel));
+			}
+			else{
+				loads->AddObject(new Channel(i+1,channelarea[i],i,iomodel));
+			}
+			iomodel->DeleteData(1,"md.initialization.channelarea");
 		}
 	}
Index: /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp	(revision 24131)
+++ /issm/trunk-jpl/src/c/classes/Loads/Channel.cpp	(revision 24132)
@@ -39,5 +39,6 @@
 }
 /*}}}*/
-Channel::Channel(int channel_id,int i,int index,IoModel* iomodel){/*{{{*/
+Channel::Channel(int channel_id,IssmDouble channelarea,int index,IoModel* iomodel){/*{{{*/
+//Channel::Channel(int channel_id,int i,int index,IoModel* iomodel)
 
 	this->id=channel_id;
@@ -48,6 +49,8 @@
 
 	/*Set channel cross section to 0*/
-	this->S    = 0.;
-	this->Sold = 0.;
+	//this->S    = 0.;
+	//this->Sold = 0.;
+	this->S    = channelarea;
+	this->Sold = channelarea;
 
 	/*Get edge info*/
Index: /issm/trunk-jpl/src/c/classes/Loads/Channel.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Channel.h	(revision 24131)
+++ /issm/trunk-jpl/src/c/classes/Loads/Channel.h	(revision 24132)
@@ -40,5 +40,6 @@
 		/*Channel constructors,destructors {{{*/
 		Channel();
-		Channel(int numericalflux_id,int i,int index,IoModel* iomodel);
+		Channel(int numericalflux_id,IssmDouble channelarea,int index,IoModel* iomodel);
+		//Channel(int numericalflux_id,int i,int index,IoModel* iomodel);
 		~Channel();
 		/*}}}*/
Index: /issm/trunk-jpl/src/m/classes/initialization.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.m	(revision 24131)
+++ /issm/trunk-jpl/src/m/classes/initialization.m	(revision 24132)
@@ -19,4 +19,5 @@
 		watercolumn         = NaN;
 		hydraulic_potential = NaN;
+		channelarea         = NaN;
 	end
 	methods
@@ -94,4 +95,5 @@
 					md = checkfield(md,'fieldname','initialization.watercolumn','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
 					md = checkfield(md,'fieldname','initialization.hydraulic_potential','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+					md = checkfield(md,'fieldname','initialization.channelarea','NaN',1,'Inf',1,'>=',0,'size',[md.mesh.numberofedges 1]);
 				end
 			end
@@ -126,4 +128,5 @@
 			fielddisplay(self,'watercolumn','subglacial water sheet thickness (for Shreve and GlaDS) [m]');
 			fielddisplay(self,'hydraulic_potential','Hydraulic potential (for GlaDS) [Pa]');
+			fielddisplay(self,'channelarea','subglacial water channel area (for GlaDS) [m2]');
 
 		end % }}}
@@ -142,4 +145,5 @@
 			WriteData(fid,prefix,'object',self,'fieldname','epl_thickness','format','DoubleMat','mattype',1);
 			WriteData(fid,prefix,'object',self,'fieldname','watercolumn','format','DoubleMat','mattype',1);
+			WriteData(fid,prefix,'object',self,'fieldname','channelarea','format','DoubleMat','mattype',1);
 			WriteData(fid,prefix,'object',self,'fieldname','hydraulic_potential','format','DoubleMat','mattype',1);
 
@@ -172,4 +176,5 @@
 			writejs1Darray(fid,[modelname '.initialization.watercolumn'],self.watercolumn);
 			writejs1Darray(fid,[modelname '.initialization.hydraulic_potential'],self.hydraulic_potential);
+			writejs1Darray(fid,[modelname '.initialization.channel'],self.channelarea);
 
 		end % }}}