Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Masstransport/UpdateElementsMasstransport.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Masstransport/UpdateElementsMasstransport.cpp	(revision 15951)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Masstransport/UpdateElementsMasstransport.cpp	(revision 15952)
@@ -45,4 +45,5 @@
 	iomodel->FetchDataToInput(elements,BathymetryEnum);
 	iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum);
+	iomodel->FetchDataToInput(elements,MaskGroundediceLevelsetEnum);
 	iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);
 	iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum);
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp	(revision 15951)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp	(revision 15952)
@@ -96,4 +96,5 @@
 	iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);
 	iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum);
+	iomodel->FetchDataToInput(elements,MaskGroundediceLevelsetEnum);
 	iomodel->FetchDataToInput(elements,MaterialsRheologyBEnum);
 	iomodel->FetchDataToInput(elements,MaterialsRheologyNEnum);
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Transient/UpdateElementsTransient.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Transient/UpdateElementsTransient.cpp	(revision 15951)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Transient/UpdateElementsTransient.cpp	(revision 15952)
@@ -10,38 +10,3 @@
 
 void	UpdateElementsTransient(Elements* elements, Parameters* parameters,IoModel* iomodel,int analysis_counter,int analysis_type){
-
-	bool        isgl;
-	int         migration_style;
-	IssmDouble  rho_ice,rho_water;
-	IssmDouble *phi   = NULL;
-
-	/*Fetch data needed*/
-	iomodel->Constant(&isgl,TransientIsgroundinglineEnum);
-	iomodel->Constant(&migration_style,GroundinglineMigrationEnum);
-
-	if(isgl==1 && (migration_style==SubelementMigrationEnum || migration_style==SubelementMigration2Enum)){
-
-		/*get parameters and constants: */
-		iomodel->Constant(&rho_ice,MaterialsRhoIceEnum);
-		iomodel->Constant(&rho_water,MaterialsRhoWaterEnum);
-		iomodel->FetchData(2,ThicknessEnum,BathymetryEnum);
-
-		/*Create phi vector */
-		phi=xNew<IssmDouble>(iomodel->numberofvertices);
-		for(int i=0;i<iomodel->numberofvertices;i++){
-			if(iomodel->my_vertices[i]){
-				phi[i] = iomodel->Data(ThicknessEnum)[i] + rho_water/rho_ice * iomodel->Data(BathymetryEnum)[i];
-			}
-		}
-
-		/*Update elements: */
-		for(int i=0;i<elements->Size();i++){
-			Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
-			element->InputCreate(phi,element->Sid(),iomodel,iomodel->numberofvertices,1,1,MaskGroundediceLevelsetEnum,1);
-		}
-
-		/*Free ressources:*/
-		iomodel->DeleteData(2,ThicknessEnum,BathymetryEnum);
-		xDelete<IssmDouble>(phi);
-	}
 }