Index: /issm/trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 27313)
+++ /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 27314)
@@ -159,11 +159,9 @@
 	femmodel->parameters->FindParam(&isstochasticforcing,StochasticForcingIsStochasticForcingEnum);
 
-#if defined(_HAVE_OCEAN_ )
+#if defined(_HAVE_OCEAN_)
 	if(isoceancoupling) OceanExchangeDatax(femmodel,false);
 #endif
 
-	if(isstochasticforcing){
-		StochasticForcingx(femmodel);
-	}
+	if(isstochasticforcing) StochasticForcingx(femmodel);
 
 	if(isthermal && domaintype==Domain3DEnum){
@@ -190,27 +188,15 @@
 
 	/* Using Hydrology dc  coupled we need to compute smb in the hydrology inner time loop*/
-	if(issmb) {
-		smb_core(femmodel);
-	}
-
-	if(ishydrology){
-		hydrology_core(femmodel);
-	}
-
-	if(isstressbalance && (step%sb_coupling_frequency==0 || step==1) ) {
-		stressbalance_core(femmodel);
-	}
-
-	if(isdamageevolution) {
-		damage_core(femmodel);
-	}
-
-	if(ismovingfront)	{
-		movingfront_core(femmodel);
-	}
-
-	if(isdebris){
-		debris_core(femmodel);
-	}
+	if(issmb) smb_core(femmodel);
+
+	if(ishydrology) hydrology_core(femmodel);
+
+	if(isstressbalance && (step%sb_coupling_frequency==0 || step==1)) stressbalance_core(femmodel);
+
+	if(isdamageevolution) damage_core(femmodel);
+
+	if(ismovingfront)	movingfront_core(femmodel);
+
+	if(isdebris) debris_core(femmodel);
 
 	/* from here on, prepare geometry for next time step*/
@@ -219,16 +205,13 @@
 		bmb_core(femmodel);
 		masstransport_core(femmodel);
-		femmodel->UpdateVertexPositionsx();
-	}
-
-	if(isoceantransport){
-		oceantransport_core(femmodel);
-	}
-
-	if(isgroundingline){
-		groundingline_core(femmodel);
-	}
-
-	/*esa: */
+	}
+
+	if(isoceantransport) oceantransport_core(femmodel);
+
+	if(isgroundingline) groundingline_core(femmodel);
+
+	/*Update mesh vertices now that we have changed the geometry*/
+	if(ismasstransport || isgroundingline) femmodel->UpdateVertexPositionsx();
+
 	if(isesa) esa_core(femmodel);
 
@@ -243,8 +226,5 @@
 
 	/*Sampling: */
-	if(issampling){
-		if(VerboseSolution()) _printf0_("   computing Gaussian random field\n");
-		sampling_core(femmodel);
-	}
+	if(issampling) sampling_core(femmodel);
 
 	/*Any requested output that needs to be saved?*/