Index: /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 21889)
+++ /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 21890)
@@ -520,6 +520,13 @@
 
 		if(migration_style==SubelementMigrationEnum){
-			if (phi<0.00000001) mb=gmb;
-			else mb=(1-phi)*fmb+phi*gmb; // phi is the fraction of grounded ice so (1-phi) is floating
+			if (phi>0.999999999) mb=gmb;
+			else{
+				mb=(1-phi)*fmb+phi*gmb; // phi is the fraction of grounded ice so (1-phi) is floating
+				printf("mg %g \n",mb);
+				printf("phi %g \n",phi);
+				printf("fmb%g \n",fmb);
+				printf("gmb%g \n",gmb);
+				printf("\n");
+			}
 		}
 		else if(migration_style==SubelementMigration2Enum){