Changeset 26197

Timestamp:

04/16/21 11:33:10 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added cycling inputs

File:

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

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

@@ -420 +420 @@
         this->parameters->FindParam(&cycle,TimesteppingCycleForcingEnum);
 
-        if(cycle) _error_("not implemented yet");
+        /*Change input time if we cycle through the forcing*/
+        IssmDouble time0 = this->timesteps[0];
+        IssmDouble time1 = this->timesteps[this->numtimesteps - 1];
+        if(cycle && (time<time0 || time>time1)){
+
+                /*See by how many intervals we have to offset time*/
+                IssmDouble deltat = time1-time0;
+
+                //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;
+        }
 
         /*Figure step out*/
@@ -480 +494 @@
                 this->current_input->AXPY(input2,alpha2);
         }
-
 }/*}}}*/
 void TransientInput::SetAverageAsCurrentTimeInput(IssmDouble start_time,IssmDouble end_time, int averaging_method){/*{{{*/