Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 24975)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 24976)
@@ -6010,5 +6010,21 @@
 	}
 
-	if (masks->isiceonly[this->lid]){
+
+	if(masks->isoceanin[this->lid]){
+		/*From Sg, recover water sea level rise:*/
+		S=0; for(int i=0;i<NUMVERTICES;i++) S+=Sg[this->vertices[i]->Sid()]/NUMVERTICES;
+
+		/*convert to kg/m^2:*/
+		S=rho_water*S;
+
+		for(int i=0;i<gsize;i++){
+			Up[i]+=S*GU[i];
+			if(horiz){
+				North[i]+=S*GN[i];
+				East[i]+=S*GE[i];
+			}
+		}
+	}
+	else if (masks->isiceonly[this->lid]){
 
 		/*Compute ice thickness change: */
@@ -6025,19 +6041,4 @@
 				North[i]+=I*GN[i];
 				East[i]+=I*GE[i];
-			}
-		}
-	}
-	else if(masks->isoceanin[this->lid]){
-		/*From Sg, recover water sea level rise:*/
-		S=0; for(int i=0;i<NUMVERTICES;i++) S+=Sg[this->vertices[i]->Sid()]/NUMVERTICES;
-
-		/*convert to kg/m^2:*/
-		S=rho_water*S;
-
-		for(int i=0;i<gsize;i++){
-			Up[i]+=S*GU[i];
-			if(horiz){
-				North[i]+=S*GN[i];
-				East[i]+=S*GE[i];
 			}
 		}