Index: /issm/trunk-jpl/src/m/classes/hydrologydc.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 23663)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 23664)
@@ -148,5 +148,5 @@
 
 	def defaultoutputs(self,md): # {{{
-		list = ['SedimentHeadHydrostep','SedimentHeadResidual','EffectivePressureHydrostep']
+		list = ['SedimentHeadHydrostep','SedimentHeadResidual','EffectivePressureHydrostep','HydrologydcMaskThawedNode','HydrologydcMaskThawedElt']
 		if self.isefficientlayer==1:
 			list.extend(['EplHeadHydrostep','HydrologydcMaskEplactiveNode','HydrologydcMaskEplactiveElt','EplHeadSlopeX','EplHeadSlopeY','HydrologydcEplThicknessHydrostep'])
@@ -204,5 +204,5 @@
 			md = checkfield(md,'fieldname','hydrology.epl_max_thickness','numel',[1],'>',0.)
 			md = checkfield(md,'fieldname','hydrology.epl_initial_thickness','numel',[1],'>',0.)
-			md = checkfield(md,'fieldname','hydrology.epl_colapse_thickness','numel',[1],'>',0.)
+			md = checkfield(md,'fieldname','hydrology.epl_colapse_thickness','numel',[1],'>',0.,'<',self.epl_initial_thickness)
 			md = checkfield(md,'fieldname','hydrology.epl_thick_comp','numel',[1],'values',[0,1])
 			md = checkfield(md,'fieldname','hydrology.eplflip_lock','>=',0.,'numel',[1])
Index: /issm/trunk-jpl/src/m/classes/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.py	(revision 23663)
+++ /issm/trunk-jpl/src/m/classes/model.py	(revision 23664)
@@ -783,12 +783,14 @@
 
 		# Hydrologydc variables
-		if type(md.hydrology) is 'hydrologydc':
-			md.hydrology.spcsediment_head=project2d(md,md.hydrology.spcsediment_head,1)
-			md.hydrology.sediment_transmitivity=project2d(md,md.hydrology.sediment_transmitivity,1)
-			md.hydrology.basal_moulin_input=project2d(md,md.hydrology.basal_moulin_input,1)
-			md.hydrology.mask_thawed_node=project2d(md,md.hydrology.mask_thawed_node,1)
-			if md.hydrology.isefficientlayer == 1:
-				md.hydrology.mask_eplactive_node=project2d(md,md.hydrology.mask_eplactive_node,1)
-				md.hydrology.spcepl_head=project2d(md,md.hydrology.spcepl_head,1)
+		hydrofields=md.hydrology.__dict__.keys()
+		for field in hydrofields:
+			try:
+				isvector=np.logical_or(np.shape(md.hydrology.__dict__[field])[0]==md.mesh.numberofelements,
+															 np.shape(md.hydrology.__dict__[field])[0]==md.mesh.numberofvertices)
+			except IndexError:
+				isvector=False
+			#we colpase only fields that are vertices or element based
+			if isvector:
+				md.hydrology.__dict__[field]=project2d(md,md.hydrology.__dict__[field],1)
 
 		#boundary conditions