From 5e825ffdf48511c9e5574a3e136aa1a0374453b7 Mon Sep 17 00:00:00 2001
From: Bryan Herman <bherman@mit.edu>
Date: Mon, 7 Nov 2011 22:19:30 -0500
Subject: [PATCH] added in cmfd execute information for testing flux extraction
 from tally object

---
 src/DEPENDENCIES     |  1 +
 src/cmfd_execute.F90 | 17 +++++++++++++++++
 src/main.f90         |  4 ++++
 3 files changed, 22 insertions(+)

diff --git a/src/DEPENDENCIES b/src/DEPENDENCIES
index 32b3406f08..6a0fc05dba 100644
--- a/src/DEPENDENCIES
+++ b/src/DEPENDENCIES
@@ -118,6 +118,7 @@ interpolation.o: search.o
 logging.o: constants.o
 logging.o: global.o
 
+main.o: cmfd_execute.o
 main.o: constants.o
 main.o: global.o
 main.o: initialize.o
diff --git a/src/cmfd_execute.F90 b/src/cmfd_execute.F90
index 80bcded5d2..3cd608714f 100644
--- a/src/cmfd_execute.F90
+++ b/src/cmfd_execute.F90
@@ -9,6 +9,23 @@ module cmfd_execute
 
 contains
 
+!===============================================================================
+! EXECUTE_CMFD is the highest level routine that controls CMFD calculation 
+!===============================================================================
+
+  subroutine execute_cmfd()
+
+    ! allocate cmfd object
+    call allocate_cmfd()
+
+    ! calculate all cross sections based on reaction rates from last batch
+    call compute_xs()
+
+    ! print out flux for debugginb
+    write(11,*) cmfd % flux
+
+  end subroutine
+
 !===============================================================================
 ! ALLOCATE_CMFD allocates all of the space for the cmfd object based on tallies
 !===============================================================================
diff --git a/src/main.f90 b/src/main.f90
index d60083f862..5b90a527a1 100644
--- a/src/main.f90
+++ b/src/main.f90
@@ -1,5 +1,6 @@
 program main
 
+  use cmfd_execute,    only: execute_cmfd
   use constants
   use global
   use initialize,      only: initialize_run
@@ -132,6 +133,9 @@ contains
     ! Calculate statistics for tallies
     call tally_statistics()
 
+    ! call cmfd calculation
+    call execute_cmfd()
+
     if (master) call header("SIMULATION FINISHED", 1)
 
   end subroutine run_problem