Changeset 24279

Timestamp:

10/25/19 13:36:55 (6 years ago)

Author:

schlegel

Message:

CHG: update grain radius when merging cells

Location:

issm/trunk-jpl

Files:

• src/c/modules/SurfaceMassBalancex/Gembx.cpp (modified) (2 diffs)
• test/Archives/Archive243.arch (modified) ( previous)
• test/Archives/Archive244.arch (modified) ( previous)
• test/Archives/Archive252.arch (modified) ( previous)

issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp

@@ -1198 +1198 @@
                         dfall = min(max(1.29 - 0.17*V,0.20),1.0);
                         sfall = min(max(0.08*V + 0.38,0.5),0.9);
-                        refall = (0.1 + (1.0-dfall)*0.25 + (0.5-sfall)*0.1)/2.0;
+                        refall = max((0.1 + (1.0-dfall)*0.25 + (0.5-sfall)*0.1)/2.0,Gdntol);
 
                         // if snow depth is greater than specified min dz, new cell created
@@ -1228 +1228 @@
                                 gdn[0] = (dfall * P + gdn[0] * mInit[0])/mass;
                                 gsp[0] = (sfall * P + gsp[0] * mInit[0])/mass;
-                                refall = (0.1 + (1.0-dfall)*0.25 + (0.5-sfall)*0.1)/2.0;
+                                re[0] = max((0.1 + (1.0-gdn[0])*0.25 + (0.5-gsp[0])*0.1)/2.0,Gdntol);
                         }
                 }