Index: /issm/trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 17269)
+++ /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 17270)
@@ -123,22 +123,4 @@
 		}
 		
-		if(isgroundingline){
-			if(VerboseSolution()) _printf0_("   computing new grounding line position\n");
-			GroundinglineMigrationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
-			if(save_results){
-				int outputs[3] = {SurfaceEnum,BedEnum,MaskGroundediceLevelsetEnum};
-				femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3);
-			}
-		}
-		if(isgia){
-			if(VerboseSolution()) _printf0_("   computing glacial isostatic adjustment\n");
-			#ifdef _HAVE_GIA_
-			gia_core(femmodel);
-			#else
-			_error_("ISSM was not compiled with gia capabilities. Exiting");
-			#endif
-
-		}
-
 		if(islevelset){
 			if(VerboseSolution()) _printf0_("   computing movement of ice boundaries\n");
@@ -173,4 +155,21 @@
 			femmodel->UpdateVertexPositionsx();
 		}
+		
+		if(isgroundingline){
+			if(VerboseSolution()) _printf0_("   computing new grounding line position\n");
+			GroundinglineMigrationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
+			if(save_results){
+				int outputs[3] = {SurfaceEnum,BedEnum,MaskGroundediceLevelsetEnum};
+				femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3);
+			}
+		}
+		if(isgia){
+			if(VerboseSolution()) _printf0_("   computing glacial isostatic adjustment\n");
+			#ifdef _HAVE_GIA_
+			gia_core(femmodel);
+			#else
+			_error_("ISSM was not compiled with gia capabilities. Exiting");
+			#endif
+		}
 
 		/*unload results*/