Changeset 22747

Timestamp:

05/07/18 13:56:48 (7 years ago)

Author:

schlegel

Message:

CHG: check mass as done in other GEMB routines

File:

• issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp (modified) (1 diff)

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

@@ -1763 +1763 @@
         /*only in forward mode! avoid round in AD mode as it is not differentiable: */
         #ifndef _HAVE_ADOLC_
-        dm = round((mSum0 - mSum1 + mAdd)*1000);
+        dm = round((mSum0 - mSum1 + mAdd)*100.0)/100.0;
         dE = round(sumE0 - sumE1 - sumER +  addE);
         if (dm !=0  || dE !=0) _error_("mass or energy are not conserved in melt equations\n"