diff --git a/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.C b/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.C
index f3814b08..d9875fac 100644
--- a/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.C
+++ b/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.C
@@ -851,6 +851,27 @@ const surfaceScalarField& gerhardsRecModel::exportSurfaceScalarField(word fieldn
 }
 
 
+PtrList<volScalarField>& gerhardsRecModel::exportVolScalarFieldList(word fieldname)
+{
+    const label fieldI = getVolScalarFieldIndex(fieldname);
+
+    return volScalarFieldList_[fieldI];
+}
+
+PtrList<volVectorField>& gerhardsRecModel::exportVolVectorFieldList(word fieldname)
+{
+    const label fieldI = getVolVectorFieldIndex(fieldname);
+
+    return volVectorFieldList_[fieldI];
+}
+
+PtrList<surfaceScalarField>& gerhardsRecModel::exportSurfaceScalarFieldList(word fieldname)
+{
+    const label fieldI = getSurfaceScalarFieldIndex(fieldname);
+
+    return surfaceScalarFieldList_[fieldI];
+}
+
 
 
 SymmetricSquareMatrix<scalar>& gerhardsRecModel::recurrenceMatrix()
diff --git a/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.H b/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.H
index 81fb2ccd..af0b0e98 100644
--- a/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.H
+++ b/src/recurrence/recModel/gerhardsRecModel/gerhardsRecModel.H
@@ -100,9 +100,9 @@ protected:
 
     scalar checkTimeStep();
 
-    inline label getVolScalarFieldIndex(word, label) const;
-    inline label getVolVectorFieldIndex(word, label) const;
-    inline label getSurfaceScalarFieldIndex(word, label) const;
+    inline label getVolScalarFieldIndex(word, label=0) const;
+    inline label getVolVectorFieldIndex(word, label=0) const;
+    inline label getSurfaceScalarFieldIndex(word, label=0) const;
 
     void readFieldSeries();
 
@@ -143,7 +143,11 @@ public:
     const volVectorField& exportVolVectorField(word, label);
     const surfaceScalarField& exportSurfaceScalarField(word, label);
 
-  //  tmp<surfaceScalarField> exportAveragedSurfaceScalarField(word, scalar, label index = -1);
+    PtrList<volScalarField>& exportVolScalarFieldList(word);
+    PtrList<volVectorField>& exportVolVectorFieldList(word);
+    PtrList<surfaceScalarField>& exportSurfaceScalarFieldList(word);
+
+    //  tmp<surfaceScalarField> exportAveragedSurfaceScalarField(word, scalar, label index = -1);
     void exportAveragedVolVectorField(volVectorField&, word, scalar, label index = -1) const;
 
     SymmetricSquareMatrix<scalar>& recurrenceMatrix();