Changeset 26199

Timestamp:

04/16/21 14:20:43 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: hack for Nicole

File:

• issm/trunk-jpl/src/c/classes/Inputs/TransientInput.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/classes/Inputs/TransientInput.cpp

@@ -430 +430 @@
                 //int num_intervals = floor((time-time0)/deltat); //Cannot do that because of AD!
                 int num_intervals = reCast<int,IssmDouble>(fabs(time-time0)/deltat);
-                if(time<time0) num_intervals = -num_intervals-1;
-
-                /*Now offset time so that we do the right interpolation below*/
-                time = time - num_intervals*deltat;
+                //if(time<time0) num_intervals = -num_intervals-1;
+                
+                if(fabs(time-time0)/deltat == reCast<IssmDouble>(num_intervals)){
+                        /*Hack to make sure we always cover the last value of the series (discussion with Nicole)*/
+                        time = time1;
+                }
+                else{
+                        /*Now offset time so that we do the right interpolation below*/
+                        time = time - num_intervals*deltat;
+                }
         }