Index: /issm/trunk/src/m/classes/model.m
===================================================================
--- /issm/trunk/src/m/classes/model.m	(revision 7646)
+++ /issm/trunk/src/m/classes/model.m	(revision 7647)
@@ -123,4 +123,5 @@
 		 meltingpoint=0;
 		 beta=0;
+		 
 		 mixed_layer_capacity=0;
 		 thermal_exchange_velocity=0;
@@ -143,7 +144,9 @@
 		 elementonicesheet=NaN;
 		 elementonwater=NaN;
+		 elementonnuna=NaN;
 		 gridoniceshelf=NaN;
 		 gridonicesheet=NaN;
 		 gridonwater=NaN;
+		 gridonnuna=NaN;
 		 surface=NaN;
 		 thickness=NaN;
@@ -159,4 +162,5 @@
 		 spctemperature=NaN;
 		 spcthickness=NaN;
+		 spcwatercolumn=NaN;
 
 		 %Observations
@@ -256,5 +260,14 @@
 		 melting_rate_correction_apply=NaN;
 		 pressure=NaN;
-
+		 
+		 %hydrology
+		 watercolumn=NaN;
+		 hydro_n=NaN;
+		 hydro_CR=NaN;
+		 hydro_p=NaN;
+		 hydro_q=NaN;
+		 hydro_gamma=0;
+		 hydro_kn=0;
+		 
 		 %Parallelisation parameters
 		 cluster=NaN;
@@ -703,4 +716,13 @@
 			 md.kff=1;
 
+			 %hydrology:  from Johnson's 2002 thesis, section 3.5.4
+			 md.hydro_n=.02;
+			 md.hydro_CR=2;
+			 md.hydro_p=2;
+			 md.hydro_q=1;
+			 md.hydro_gamma=1/(md.hydro_n*md.hydro_CR^md.hydro_p*(md.rho_water*md.g)^md.hydro_q);
+			 md.hydro_kn=0;
+
+
 			 %i/o:
 			 md.io_gather=1;