Changeset 27342

Timestamp:

10/27/22 10:07:59 (2 years ago)

Author:

vverjans

Message:

BUG: fixing indexing in Armaprocess

Location:

issm/trunk-jpl/src/c

Files:

• classes/Elements/Element.cpp (modified) (2 diffs)
• modules/FrontalForcingsx/FrontalForcingsx.cpp (modified) (1 diff)

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

@@ -130 +130 @@
                         /*Compute polynomial with parameters of indperiod*/
                         if(indperiod==0) telapsed_break = tatstep;
-                        else             telapsed_break = tatstep-datebreaks_basin[indperiod];
+                        else             telapsed_break = tatstep-datebreaks_basin[indperiod-1];
                         for(int j=0;j<numparams;j++)   sumpoly[s] = sumpoly[s]+polyparams_basin[indperiod+j*numperiods]*pow(telapsed_break,j);
                 }
@@ -2767 +2767 @@
                 /*Compute intercept+trend terms with parameters of indperiod*/
       if(indperiod==0) telapsed_break = time;
-      else             telapsed_break = time-datebreaks_b[indperiod];
+      else             telapsed_break = time-datebreaks_b[indperiod-1];
       piecewiselin = intercepts_b[indperiod]+trends_b[indperiod]*telapsed_break;
         }

issm/trunk-jpl/src/c/modules/FrontalForcingsx/FrontalForcingsx.cpp

@@ -101 +101 @@
    xDelete<IssmDouble>(polyparams);
    xDelete<IssmDouble>(datebreaks);
+   xDelete<IssmDouble>(monthdatebreaks);
+   xDelete<IssmDouble>(monthintercepts);
+   xDelete<IssmDouble>(monthtrends);
 }/*}}}*/