diff --git a/Source/hydro/outline.txt b/Source/hydro/outline.txt
deleted file mode 100644
index 2a7f80919e..0000000000
--- a/Source/hydro/outline.txt
+++ /dev/null
@@ -1,55 +0,0 @@
-
-   1. MFIter loop:
-
-      * convert U_avg to U_cc
-      * construct q_bar from U_avg (using ctoprim)
-      * construct q_cc from U_cc (using ctoprim)
-
-   2. MFIter loop
-
-      * construct q_avg from q_bar and q_cc (using ctoprim)
-
-   3. in C++, loop over dimensions -- to save memory, we do each separately
-
-      a. main Fortran routine (accepts construct q_avg and q_bar)
-
-         * compute flattening coefficients from q_bar
-         * find face-average interface states from q_avg  incl/ flattening
-         * solve Riemann problem and return face-averaged interface states
-
-      b. separate MFIter:
-
-         * construct face-centered interface states
-         * compute face-centered and face-average fluxes
-
-      c. separate MFIter:
-
-         * compute final fluxes
-
-
-Extra MFs:
-
-  1. q_bar -- needed throughout the entire process
-  2. q_int_fc and q_int_avg
-  3. F_fc to hold the face-centered fluxes
-
-
-
-Also needed:
-
-  -- for the Riemann solver, we return the interface state of only a subset
-     of variable NGDNV, not all NQ -- this means that we will need a separate
-     Riemann solver here to get all interface states
-
-  -- how do we compute the pdivu term to 4th order -- this is needed
-     for the internal energy evolution
-
-  -- how do we do source terms to 4th order (like gravity)
-
-  -- we need to make sure that at the end of evolution Source_Type has
-     the proper old and new sources, to allow for multilevel support
-
-  -- we need to fix pradial
-
-  -- for non-Cartesian, we need to have our Laplacian (and transverse Lap) have
-     the proper area factors