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

Timestamp:

02/16/18 14:00:20 (7 years ago)

Author:

schlegel

Message:

CHG: decreasing gemb thermo dt to 1/11

File:

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

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

@@ -655 +655 @@
 
         // determine minimum acceptable delta t (diffusion number > 1/2) [s]
-        // NS: 2.16.18 divided dt by 9 for stability (changed 3 to 27)
+        // NS: 2.16.18 divided dt by 11 for stability (changed 3 to 33)
         dt=1e12; 
-        for(int i=0;i<m;i++)dt = fmin(dt,CI * pow(dz[i],2) * d[i]  / (27 * K[i]));
+        for(int i=0;i<m;i++)dt = fmin(dt,CI * pow(dz[i],2) * d[i]  / (33 * K[i]));
 
         // smallest possible even integer of 60 min where diffusion number > 1/2