diff --git a/analysis/runAnalysis.cpp b/analysis/runAnalysis.cpp
index a98a6448..a2b3801d 100644
--- a/analysis/runAnalysis.cpp
+++ b/analysis/runAnalysis.cpp
@@ -412,6 +412,7 @@ public:
 
     	PROFILE_START("Compute subphase",1);
     	Averages.Full();
+    	Averages.Write(timestep);
     	PROFILE_STOP("Compute subphase",1);
     }
 private:
@@ -598,6 +599,7 @@ void runAnalysis::finish( )
     d_wait_blobID.reset();
     d_wait_analysis.reset();
     d_wait_vis.reset();
+    d_wait_subphase.reset();
     d_wait_restart.reset();
     // Syncronize
     MPI_Barrier( d_comm );
diff --git a/analysis/runAnalysis.h b/analysis/runAnalysis.h
index d9c96e6b..df73fafd 100644
--- a/analysis/runAnalysis.h
+++ b/analysis/runAnalysis.h
@@ -15,7 +15,7 @@ typedef std::shared_ptr<std::vector<BlobIDType>> BlobIDList;
 
 // Types of analysis
 enum class AnalysisType : uint64_t { AnalyzeNone=0, IdentifyBlobs=0x01, CopyPhaseIndicator=0x02, 
-    CopySimState=0x04, ComputeAverages=0x08, CreateRestart=0x10, WriteVis=0x20 };
+    CopySimState=0x04, ComputeAverages=0x08, CreateRestart=0x10, WriteVis=0x20, ComputeSubphase=0x40 };
 
 
 //! Class to run the analysis in multiple threads