From 364199eb2847d67e15d8bc3544a9038b5b3ce177 Mon Sep 17 00:00:00 2001
From: Bryan Herman <bherman@mit.edu>
Date: Wed, 17 Sep 2014 17:25:06 -0400
Subject: [PATCH] added input description of write_initial_source and allowed
 binary extension for name

---
 docs/source/usersguide/input.rst | 42 ++++++++++++++++----------------
 src/source.F90                   |  4 +++
 2 files changed, 25 insertions(+), 21 deletions(-)

diff --git a/docs/source/usersguide/input.rst b/docs/source/usersguide/input.rst
index 5a775be8e9..4fc49b35f5 100644
--- a/docs/source/usersguide/input.rst
+++ b/docs/source/usersguide/input.rst
@@ -433,6 +433,13 @@ attributes/sub-elements:
 
       *Default*: 0.988 2.249
 
+  :write_initial:
+    An element specifying whether to write out the intial source bank used
+    at the beginning of the first batch. The output file is named 
+    "initial_source.binary(h5)"
+
+      *Default*: false
+
 ``<state_point>`` Element
 -------------------------
 
@@ -1303,15 +1310,6 @@ Currently, it allows users to accelerate fission source convergence during
 inactive neutron batches. To run CMFD, the ``<run_cmfd>`` element in
 ``settings.xml`` should be set to "true".
 
-``<atoli>`` Element
---------------------
-
-The ``<atoli>`` element specifies the absolute inner tolerance on the Gauss-Seidel
-when performing CMFD calculations. It is only used in the standalone CMFD solver
-and not when PETSc is active. 
-
-  *Default*: 1.e-10
-
 ``<begin>`` Element
 -------------------
 
@@ -1338,9 +1336,9 @@ The ``<display>`` element sets one additional CMFD output column. Options are:
 ``<dhat_reset>`` Element
 ------------------------
 
-The ``<dhat_reset>`` element controls whether dhat nonlinear CMFD parameters
-should be reset to zero before solving CMFD eigenproblem. It can be turned on with 
-"true" and off with "false".
+The ``<dhat_reset>`` element controls whether :math:`\widehat{D}` nonlinear 
+CMFD parameters should be reset to zero before solving CMFD eigenproblem. 
+It can be turned on with "true" and off with "false".
 
   *Default*: false
 
@@ -1362,6 +1360,17 @@ It can be turned on with "true" and off with "false".
 
   *Default*: false
 
+``<gauss_seidel_tolerance>`` Element
+--------------------
+
+The ``<gauss_seidel_tolerance>`` element specifies two parameters. The first is
+the absolute inner tolerance for Gauss-Seidel iterations when performing CMFD
+and the second is the relative inner tolerance for Gauss-Seidel iterations
+for CMFD calculations. It is only used in the standalone CMFD power iteration
+solver and not when PETSc is active. 
+
+  *Default*: 1.e-10 1.e-5
+
 ``<ksp_monitor>`` Element
 -------------------------
 
@@ -1455,15 +1464,6 @@ This option can be turned on with "true" and turned off with "false".
 
   *Default*: false
 
-``<rtoli>`` Element
---------------------
-
-The ``<rtoli>`` element specifies the relative inner tolerance on the Gauss-Seidel
-when performing CMFD calculations. It is only used in the standalone CMFD solver
-and not when PETSc is active. 
-
-  *Default*: 1.e-5
-
 ``<run_adjoint>`` Element
 -------------------------
 
diff --git a/src/source.F90 b/src/source.F90
index f6de6042b1..fd889abc42 100644
--- a/src/source.F90
+++ b/src/source.F90
@@ -81,7 +81,11 @@ contains
     if (write_initial_source) then
       message = 'Writing out initial source guess...'
       call write_message(1)
+#ifdef HDF5
       filename = 'initial_source.h5'
+#else
+      filename = 'initial_source.binary'
+#endif
       call sp % file_create(filename, serial = .false.)
       call sp % write_source_bank()
       call sp % file_close()