diff --git a/docs/source/devguide/index.rst b/docs/source/devguide/index.rst
index a8ea5434c0..83b741c43c 100644
--- a/docs/source/devguide/index.rst
+++ b/docs/source/devguide/index.rst
@@ -14,3 +14,4 @@ as debugging.
 
     structures
     xml-fortran
+    statepoint
diff --git a/docs/source/devguide/statepoint.rst b/docs/source/devguide/statepoint.rst
new file mode 100644
index 0000000000..cd2785fc5a
--- /dev/null
+++ b/docs/source/devguide/statepoint.rst
@@ -0,0 +1,205 @@
+.. _devguide_statepoint:
+
+======================================
+State Point Binary File Specifications
+======================================
+
+----------
+Revision 2
+----------
+
+**integer(4) REVISION_STATEPOINT**
+
+    Revision of the binary state point file. Any time a change is made in the
+    format of the state-point file, this integer is incremented.
+
+**integer(4) VERSION_MAJOR**
+
+    Major version number for OpenMC
+
+**integer(4) VERSION_MINOR**
+
+    Minor version number for OpenMC
+
+**integer(4) VERSION_RELEASE**
+
+    Release version number for OpenMC
+
+**integer(4) run_mode**
+
+    run mode used. The modes are described in constants.F90.
+
+**integer(8) n_particles**
+
+    Number of particles used per generation.
+
+**integer(4) n_batches**
+
+    Total number of batches (active + inactive).
+
+**integer(4) n_inactive**
+
+    Number of inactive batches
+
+**integer(4) gen_per_batch**
+
+    Number of generations per batch for criticality calculations
+
+**integer(4) current_batch**
+
+    The number of batches already simulated.
+
+*do i = 1, current_batch*
+
+    **real(8) k_batch(i)**
+
+        k-effective for the i-th batch
+
+    if (entropy_on)
+
+        **real(8) entropy for the i-th batch**
+
+**integer(4) N_GLOBAL_TALLIES**
+
+    Number of global tally scores
+
+*do i = 1, N_GLOBAL_TALLIES*
+
+    **real(8) global_tallies(i) % sum**
+
+        Accumulated sum for the i-th global tally
+
+*do i = 1, N_GLOBAL_TALLIES*
+
+    **real(8) global_tallies(i) % sum_sq**
+
+        Accumulated sum of squares for the i-th global tally
+
+**integer(4) n_tallies**
+
+*do i = 1, n_tallies*
+
+    **integer(4) size(tallies(i) % scores, 1)**
+
+        Total number of score bins for the i-th tally
+
+    **integer(4) size(tallies(i) % scores, 2)**
+
+        Total number of filter bins for the i-th tally
+
+*do i = 1, n_tallies*
+
+    *do k = 1, size(tallies(i) % scores, 2)*
+
+        *do j = 1, size(tallies(i) % scores, 1)*
+
+            **real(8) tallies(i) % scores(j,k) % sum**
+            
+                Accumulated sum for the j-th score and k-th filter of the i-th
+                tally
+
+            **real(8) tallies(i) % scores(j,k) % sum_sq**
+
+                Accumulated sum of squares for the j-th score and k-th filter of
+                the i-th tally
+
+----------
+Revision 1
+----------
+
+**integer(4) REVISION_STATEPOINT**
+
+    Revision of the binary state point file. Any time a change is made in the
+    format of the state-point file, this integer is incremented.
+
+**integer(4) VERSION_MAJOR**
+
+    Major version number for OpenMC
+
+**integer(4) VERSION_MINOR**
+
+    Minor version number for OpenMC
+
+**integer(4) VERSION_RELEASE**
+
+    Release version number for OpenMC
+
+**integer(4) run_mode**
+
+    run mode used. The modes are described in constants.F90.
+
+**integer(8) n_particles**
+
+    Number of particles used per generation.
+
+**integer(4) n_batches**
+
+    Total number of batches (active + inactive).
+
+**integer(4) n_inactive**
+
+    Number of inactive batches
+
+**integer(4) gen_per_batch**
+
+    Number of generations per batch for criticality calculations
+
+**integer(4) current_batch**
+
+    The number of batches already simulated.
+
+*do i = 1, current_batch*
+
+    **real(8) k_batch(i)**
+
+        k-effective for the i-th batch
+
+    if (entropy_on)
+
+        **real(8) entropy for the i-th batch**
+
+**integer(4) N_GLOBAL_TALLIES**
+
+    Number of global tally scores
+
+*do i = 1, N_GLOBAL_TALLIES*
+
+    **real(8) global_tallies(i) % sum**
+
+        Accumulated sum for the i-th global tally
+
+*do i = 1, N_GLOBAL_TALLIES*
+
+    **real(8) global_tallies(i) % sum_sq**
+
+        Accumulated sum of squares for the i-th global tally
+
+**integer(4) n_tallies**
+
+*do i = 1, n_tallies*
+
+    **integer(4) size(tallies(i) % scores, 1)**
+
+        Total number of score bins for the i-th tally
+
+    **integer(4) size(tallies(i) % scores, 2)**
+
+        Total number of filter bins for the i-th tally
+
+    *do k = 1, size(tallies(i) % scores, 2)*
+
+        *do j = 1, size(tallies(i) % scores, 1)*
+
+            **real(8) tallies(i) % scores(j,k) % sum**
+            
+                Accumulated sum for the j-th score and k-th filter of the i-th
+                tally
+
+    *do k = 1, size(tallies(i) % scores, 2)*
+
+        *do j = 1, size(tallies(i) % scores, 1)*
+
+            **real(8) tallies(i) % scores(j,k) % sum_sq**
+
+                Accumulated sum of squares for the j-th score and k-th filter of
+                the i-th tally