diff --git a/src/Utility/HDF5Writer.C b/src/Utility/HDF5Writer.C
index 99fdb19f..9c84b947 100644
--- a/src/Utility/HDF5Writer.C
+++ b/src/Utility/HDF5Writer.C
@@ -338,6 +338,8 @@ void HDF5Writer::writeSIM (int level, const DataEntry& entry,
     int loc = sim->getLocalPatchIndex(i+1);
     if (loc > 0) // we own the patch
     {
+      sim->extractPatchSolution(Vectors(1,*sol), loc-1);
+      // order is important - do this last to ensure correct vector is loaded
       size_t ndof1 = sim->extractPatchSolution(*sol,loc-1);
       Vector& psol = const_cast<IntegrandBase*>(prob)->getSolution();
       if (results & DataExporter::PRIMARY) {