Changeset 21890 for issm/trunk-jpl/src

Timestamp:

07/27/17 17:04:23 (8 years ago)

Author:

seroussi

Message:

BUG: yet another problem fixed with melting

File:

• issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp

@@ -520 +520 @@
 
                 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){