Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 23372)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 23373)
@@ -607,9 +607,9 @@
 				//thickness[i] = old_thickness[i]*(1.0+opening-closing)/(1.0-opening+closing);
 
-				if(element->nodes[i]->Sid()==2299){
-					printf("for node %i \n",element->nodes[i]->Sid());
-					printf(" old thickness is  %g \n",old_thickness[i]);
-					printf(" new thickness is  %g \n",thickness[i]);
-				}
+				//if(element->nodes[i]->Sid()==2299){
+				//	printf("for node %i \n",element->nodes[i]->Sid());
+				//	printf(" old thickness is  %g \n",old_thickness[i]);
+				//	printf(" new thickness is  %g \n",thickness[i]);
+				//}
 				/*Take care of otherthikening*/
 				if(thickness[i]>max_thick){
@@ -710,6 +710,6 @@
 			/* If epl thickness gets under colapse thickness, close the layer */
 			if(epl_thickness[i]<colapse_thick){
-				if(element->nodes[i]->Sid()==2299){
-					printf("masking node %i with thickness %f\n",element->nodes[i]->Sid(),epl_thickness[i]);}
+				//if(element->nodes[i]->Sid()==2299){
+					//printf("masking node %i with thickness %f\n",element->nodes[i]->Sid(),epl_thickness[i]);}
 				vec_mask->SetValue(basalelement->nodes[i]->Sid(),0.,INS_VAL);
 				recurence->SetValue(basalelement->nodes[i]->Sid(),1.,INS_VAL);
Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 23372)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 23373)
@@ -10,5 +10,4 @@
 	return 1;
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/
 
@@ -29,6 +28,6 @@
 
 	/*retrieve some parameters: */
-	/* bool   issmb; */
-	/* iomodel->FindConstant(&issmb,"md.transient.issmb"); */
+	 bool   issmb;
+	iomodel->FindConstant(&issmb,"md.transient.issmb");
 	iomodel->FindConstant(&hydrology_model,"md.hydrology.model");
 
@@ -65,7 +64,7 @@
 		parameters->AddObject(new DoubleParam(HydrologydcSedimentlimitEnum,sedimentlimit));
 	}
-	/* if(!issmb){ */
-	/* 	parameters->AddObject(iomodel->CopyConstantObject("md.smb.model",SmbEnum)); */
-	/* } */
+	if(!issmb){
+		parameters->AddObject(iomodel->CopyConstantObject("md.smb.model",SmbEnum));
+	}
 
   /*Requested outputs*/
@@ -75,5 +74,4 @@
   iomodel->DeleteData(&requestedoutputs,numoutputs,"md.hydrology.requested_outputs");
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::UpdateElements(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type){/*{{{*/
 
@@ -83,6 +81,4 @@
 	/*Fetch data needed: */
 	iomodel->FindConstant(&hydrology_model,"md.hydrology.model");
-	bool   issmb;
-	iomodel->FindConstant(&issmb,"md.transient.issmb");
 
 	/*Now, do we really want DC?*/
@@ -117,10 +113,6 @@
 		iomodel->FetchDataToInput(elements,"md.initialization.epl_head",EplHeadHydrostepEnum);
 	}
-	if(!issmb){
-		iomodel->FetchDataToInput(elements,"md.smb.model",SmbEnum);
-}
-
-}/*}}}*/
-
+
+}/*}}}*/
 void HydrologyDCInefficientAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/
 
@@ -138,5 +130,4 @@
 	iomodel->DeleteData(2,"md.mesh.vertexonbase","md.mesh.vertexonsurface");
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/
 
@@ -148,5 +139,4 @@
 	IoModelToConstraintsx(constraints,iomodel,"md.hydrology.spcsediment_head",HydrologyDCInefficientAnalysisEnum,P1Enum);
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/
 
@@ -183,14 +173,11 @@
 	_error_("not implemented");
 }/*}}}*/
-
 ElementVector* HydrologyDCInefficientAnalysis::CreateDVector(Element* element){/*{{{*/
 	/*Default, return NULL*/
 	return NULL;
 }/*}}}*/
-
 ElementMatrix* HydrologyDCInefficientAnalysis::CreateJacobianMatrix(Element* element){/*{{{*/
 _error_("Not implemented");
 }/*}}}*/
-
 ElementMatrix* HydrologyDCInefficientAnalysis::CreateKMatrix(Element* element){/*{{{*/
 
@@ -318,5 +305,4 @@
 	return Ke;
 }/*}}}*/
-
 ElementVector* HydrologyDCInefficientAnalysis::CreatePVector(Element* element){/*{{{*/
 
@@ -464,5 +450,4 @@
 	return pe;
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::GetB(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss){/*{{{*/
 	/*Compute B  matrix. B=[B1 B2 B3] where Bi is of size 3*NDOF2.
@@ -492,13 +477,10 @@
 	xDelete<IssmDouble>(dbasis);
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element){/*{{{*/
 	element->GetSolutionFromInputsOneDof(solution,SedimentHeadHydrostepEnum);
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::GradientJ(Vector<IssmDouble>* gradient,Element* element,int control_type,int control_index){/*{{{*/
 	_error_("Not implemented yet");
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/
 
@@ -595,10 +577,8 @@
 	}
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/
 	/*Default, do nothing*/
 	return;
 }/*}}}*/
-
 IssmDouble HydrologyDCInefficientAnalysis::SedimentStoring(Element* element,Gauss* gauss,Input* sed_head_input, Input* base_input){/*{{{*/
 	int unconf_scheme;
@@ -640,5 +620,4 @@
 	return sediment_storing;
 }/*}}}*/
-
 IssmDouble HydrologyDCInefficientAnalysis::SedimentTransmitivity(Element* element,Gauss* gauss,Input* sed_head_input, Input* base_input,Input* SedTrans_input){/*{{{*/
 	int unconf_scheme;
@@ -675,5 +654,4 @@
 	return sediment_transmitivity;
 }/*}}}*/
-
 void  HydrologyDCInefficientAnalysis::GetHydrologyDCInefficientHmax(IssmDouble* ph_max,Element* element, Node* innode){/*{{{*/
 
@@ -711,5 +689,4 @@
 }
 /*}}}*/
-
 IssmDouble HydrologyDCInefficientAnalysis::GetHydrologyKMatrixTransfer(Element* element){/*{{{*/
 
@@ -733,5 +710,4 @@
 	return transfer;
 }/*}}}*/
-
 IssmDouble HydrologyDCInefficientAnalysis::GetHydrologyPVectorTransfer(Element* element, Gauss* gauss, Input* epl_head_input){/*{{{*/
 
@@ -758,5 +734,4 @@
 	return transfer;
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::ElementizeEplMask(FemModel* femmodel){/*{{{*/
 
@@ -777,6 +752,5 @@
 	}
 }/*}}}*/
-
-void  HydrologyDCInefficientAnalysis::HydrologyIDSGetMask(Vector<IssmDouble>* vec_mask, Element* element){
+void  HydrologyDCInefficientAnalysis::HydrologyIDSGetMask(Vector<IssmDouble>* vec_mask, Element* element){/*{{{*/
 	bool        active_element;
 	int         domaintype;
@@ -819,5 +793,5 @@
 	xDelete<IssmDouble>(meltingrate);
 	xDelete<IssmDouble>(groundedice);
-}
+}/*}}}*/
 void HydrologyDCInefficientAnalysis::ElementizeIdsMask(FemModel* femmodel){/*{{{*/
 
@@ -838,5 +812,4 @@
 	}
 }/*}}}*/
-
 void HydrologyDCInefficientAnalysis::HydrologyIdsGetActive(Vector<IssmDouble>* active_vec, Element* element){/*{{{*/
 	/*Constants*/