Added in a field modification for reading in restart field and updated CHANGELOG.md

CHANGELOG.md

@@ -22,9 +22,12 @@ v4.4.1
 - Remove the duplicate output of errorutil (!756)
 - Fix BLAS CMake dependencies (!763)
 - Fix interpolation issue with Lagrange basis functions (!768)
+- Fix issue with average fields not working with different polynomial order
+  fields (!776)
 
 **FieldConvert**:
 - Fix issue with FieldConvert when range flag used (!761)
+- Fix issue with field ordering in the interppointdatatofld module (!754)
 
 **NekMesh**:
 - Fix memory consumption issue with Gmsh output (!747, !762)
@@ -33,9 +36,6 @@ v4.4.1
 - Add manifold meshing option (!756)
 - Fix issue with older rea input files (!765)
 
-**FieldConvert:**
-- Fix issue with field ordering in the interppointdatatofld module (!754)
-
 v4.4.0
 ------
 **Library**:

library/SolverUtils/Filters/FilterFieldConvert.cpp

@@ -194,9 +194,10 @@ void FilterFieldConvert::v_Initialise(
         // Extract fields to output
         for (int j = 0; j < m_variables.size(); ++j)
         {
+            ASSERTL1(j < pFields.num_elements(), "Number of variables is larger than the size of pFields");
             for (int i = 0; i < fieldData.size(); ++i)
             {
-                pFields[0]->ExtractDataToCoeffs(
+                pFields[j]->ExtractDataToCoeffs(
                     fieldDef[i],
                     fieldData[i],
                     m_variables[j],