Changeset 22432

Timestamp:

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

Author:

schlegel

Message:

CHG: decreasing gemb thermo dt to 1/11

Location:

issm/trunk-jpl

Files:

• src/c/modules/SurfaceMassBalancex/Gembx.cpp (modified) (1 diff)
• test/Archives/Archive243.arch (modified) ( previous)
• test/Archives/Archive244.arch (modified) ( previous)
• test/NightlyRun/test243.m (modified) (1 diff)
• test/NightlyRun/test243.py (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

issm/trunk-jpl/test/NightlyRun/test243.m

@@ -45 +45 @@
 %Fields and tolerances to track changes
 field_names      ={'SmbDz','SmbT' ,'SmbD' ,'SmbRe','SmbGdn','SmbGsp','SmbA' ,'SmbEC','SmbMassBalance'};
-field_tolerances ={1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-12};
+field_tolerances ={1e-11,1e-12,1e-11,1e-11,1e-11,1e-11,1e-12,1e-12,1e-12};
 
 field_values={...

issm/trunk-jpl/test/NightlyRun/test243.py

@@ -58 +58 @@
 #Fields and tolerances to track changes
 field_names      = ['SmbDz','SmbT' ,'SmbD' ,'SmbRe','SmbGdn','SmbGsp','SmbA' ,'SmbEC','SmbMassBalance']
-field_tolerances = [5e-4,5e-5,0.0006,0.0002,1e-5,0.0003,2e-5,2e-7,1e-7]
+field_tolerances = [1e-11,1e-12,1e-11,1e-11,1e-11,1e-11,1e-12,1e-12,1e-12]
 #shape is different in python solution (fixed using reshape) which can cause test failure:
 field_values = [