Index: /issm/trunk-jpl/src/c/classes/Inputs/TransientInput.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/TransientInput.cpp	(revision 26198)
+++ /issm/trunk-jpl/src/c/classes/Inputs/TransientInput.cpp	(revision 26199)
@@ -430,8 +430,14 @@
 		//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;
+		}
 	}