Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 8537)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 8538)
@@ -991,5 +991,5 @@
 	double     Ke_gg_gaussian2[numdof][numdof]  ={0.0};
 	double     Ke_gg_gaussian3[numdof][numdof]  ={0.0};
-	double     hydro_gamma;
+	double     hydro_gamma,hydro_p,hydro_q;
 	double     DL_scalar1;
 	double     DL_scalar2;
@@ -1003,7 +1003,9 @@
 	GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
 	this->parameters->FindParam(&dt,DtEnum);
-	
+
 	/*retrieve material parameters: */
 	hydro_gamma=matpar->GetGamma();
+	hydro_p    =matpar->GetHydroP();
+	hydro_q    =matpar->GetHydroQ();
 
 	/* Start  looping on the number of gaussian points: */
Index: /issm/trunk/src/c/objects/IoModel.cpp
===================================================================
--- /issm/trunk/src/c/objects/IoModel.cpp	(revision 8537)
+++ /issm/trunk/src/c/objects/IoModel.cpp	(revision 8538)
@@ -213,4 +213,6 @@
 	IoModelFetchData(&this->hydro_gamma,iomodel_handle,"hydro_gamma");
 	IoModelFetchData(&this->hydro_kn,iomodel_handle,"hydro_kn");
+	IoModelFetchData(&this->hydro_p,iomodel_handle,"hydro_p");
+	IoModelFetchData(&this->hydro_q,iomodel_handle,"hydro_q");
 	IoModelFetchData(&this->meltingpoint,iomodel_handle,"meltingpoint");
 	IoModelFetchData(&this->referencetemperature,iomodel_handle,"referencetemperature");
Index: /issm/trunk/src/c/objects/IoModel.h
===================================================================
--- /issm/trunk/src/c/objects/IoModel.h	(revision 8537)
+++ /issm/trunk/src/c/objects/IoModel.h	(revision 8538)
@@ -172,4 +172,6 @@
 		double hydro_gamma;
 		double hydro_kn;
+		double hydro_p;
+		double hydro_q;
 		double meltingpoint;
 		double referencetemperature;
Index: /issm/trunk/src/c/objects/Materials/Matpar.cpp
===================================================================
--- /issm/trunk/src/c/objects/Materials/Matpar.cpp	(revision 8537)
+++ /issm/trunk/src/c/objects/Materials/Matpar.cpp	(revision 8538)
@@ -67,4 +67,6 @@
 	this->gamma=matpar_gamma; 
 	this->kn=matpar_kn; 
+	this->hydro_p=iomodel->hydro_p; 
+	this->hydro_q=iomodel->hydro_q; 
 
 }
@@ -402,4 +404,14 @@
 }
 /*}}}1*/
+/*FUNCTION Matpar::GetHydroP {{{1*/
+double Matpar::GetHydroP(){
+	return hydro_p;
+}
+/*}}}1*/
+/*FUNCTION Matqar::GetHydroQ {{{1*/
+double Matqar::GetHydroQ(){
+	return hydro_q;
+}
+/*}}}1*/
 /*FUNCTION Matpar::TMeltingPoint {{{1*/
 double Matpar::TMeltingPoint(double pressure){
Index: /issm/trunk/src/c/objects/Materials/Matpar.h
===================================================================
--- /issm/trunk/src/c/objects/Materials/Matpar.h	(revision 8537)
+++ /issm/trunk/src/c/objects/Materials/Matpar.h	(revision 8538)
@@ -27,7 +27,10 @@
 		double  thermal_exchange_velocity;
 		double  g;
+
 		/*hydrology: */
 		double  kn; 
 		double  gamma;
+		double  hydro_p;
+		double  hydro_q;
 
 	public:
@@ -81,4 +84,6 @@
 		double GetGamma();
 		double GetKn();
+		double GetHydroP();
+		double GetHydroQ();
 		double TMeltingPoint(double pressure);
 		double PureIceEnthalpy(double pressure);