Change in averaging method for IPVariable over time

NonLinearSolver/nlsolver/nonLinearMechSolver.cpp

@@ -3838,6 +3838,7 @@ void nonLinearMechSolver::computePeriodAveragedIPDataOnPhysical(const double tst
                     {
                         double avgValue = 0.0;
                         int count = 0;
+                        double sum_dt = 0.0;
                         // here can check GP against std::map coordinates for correctness
@@ -3848,12 +3849,17 @@ void nonLinearMechSolver::computePeriodAveragedIPDataOnPhysical(const double tst
                             if(t >= tstart && t <= tend)
                             {
                                 count++;
+                                sum_dt += t;
                                 // compute average of IPVariable over given time period
-                                avgValue += IPVariableOverTime.at(t).at(eleNum).at(j);
+                                //avgValue += IPVariableOverTime.at(t).at(eleNum).at(j);
+
+                                // sum(w_AV * dt)
+                                avgValue += (IPVariableOverTime.at(t).at(eleNum).at(j)) * t;
                             }
                         }
-                        avgValue /= count;
+                        //avgValue /= count;
+                        avgValue /= sum_dt;
                         tempMap.emplace(j, avgValue);
                     }
                     avgIPVariable.emplace(eleNum, tempMap);

[Vinayak Gholap]

average = sum(w_AV * dt)/sum(dt)