Changeset 27862 for issm/trunk-jpl/src/c/classes/Elements/Element.cpp

Timestamp:

07/29/23 06:15:00 (21 months ago)

Author:

vverjans

Message:

CHG: ARMA process should start from starttime and not from 0

File:

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

issm/trunk-jpl/src/c/classes/Elements/Element.cpp

@@ -160 +160 @@
                         for(int j=0;j<numparams;j++)   sumpoly[s] = sumpoly[s]+polyparams_basin[indperiod+j*numperiods]*pow(telapsed_break,j);
                 }
-                else for(int j=0;j<numparams;j++) sumpoly[s] = sumpoly[s]+polyparams_basin[j*numperiods]*pow(tatstep,j);
+                else for(int j=0;j<numparams;j++) sumpoly[s] = sumpoly[s]+polyparams_basin[j*numperiods]*pow(tatstep-starttime,j);
         }
 
@@ -2462 +2462 @@
                 if(lapserates_b[ii]!=0) isadjustsmb=true;
         }
-
+        
         /*Adjust SMB if any lapse rate value is non-zero*/
         if(isadjustsmb){