diff --git a/Source/driver/Castro.cpp b/Source/driver/Castro.cpp
index 39e836b3f0..7721244b93 100644
--- a/Source/driver/Castro.cpp
+++ b/Source/driver/Castro.cpp
@@ -24,7 +24,6 @@
 
 #ifdef RADIATION
 #include <Radiation.H>
-#include <RAD_F.H>
 #endif
 
 #ifdef AMREX_PARTICLES
diff --git a/Source/hydro/Castro_ctu_rad.cpp b/Source/hydro/Castro_ctu_rad.cpp
index 00ee52923e..464f52226a 100644
--- a/Source/hydro/Castro_ctu_rad.cpp
+++ b/Source/hydro/Castro_ctu_rad.cpp
@@ -3,7 +3,6 @@
 
 #include "Radiation.H"
 #include "RadHydro.H"
-#include "RAD_F.H"
 #include "fluxlimiter.H"
 
 using namespace amrex;
diff --git a/Source/radiation/MGFLD.cpp b/Source/radiation/MGFLD.cpp
index 8709ad344e..9be9efb67c 100644
--- a/Source/radiation/MGFLD.cpp
+++ b/Source/radiation/MGFLD.cpp
@@ -1,5 +1,4 @@
 #include <Radiation.H>
-#include <RAD_F.H>
 #include <rad_util.H>
 #include <filter.H>
 #include <blackbody.H>
diff --git a/Source/radiation/MGFLDRadSolver.cpp b/Source/radiation/MGFLDRadSolver.cpp
index c622affdf7..250de9b1ae 100644
--- a/Source/radiation/MGFLDRadSolver.cpp
+++ b/Source/radiation/MGFLDRadSolver.cpp
@@ -5,9 +5,6 @@
 #include <Radiation.H>
 #include <RadSolve.H>
 
-
-#include <RAD_F.H>
-
 #include <iostream>
 #include <iomanip>
 
diff --git a/Source/radiation/MGRadBndry.cpp b/Source/radiation/MGRadBndry.cpp
index 84d08e3273..f69f4ac773 100644
--- a/Source/radiation/MGRadBndry.cpp
+++ b/Source/radiation/MGRadBndry.cpp
@@ -10,8 +10,6 @@
 #include <omp.h>
 #endif
 
-#include <RAD_F.H>
-
 using namespace amrex;
 
 int         MGRadBndry::ngroups = 1;
diff --git a/Source/radiation/Make.package b/Source/radiation/Make.package
index 93c5f935d9..5a0023f4da 100644
--- a/Source/radiation/Make.package
+++ b/Source/radiation/Make.package
@@ -29,8 +29,6 @@ CEXE_headers += HABEC.H
 CEXE_headers += filter.H
 CEXE_headers += filt_prim.H
 
-FEXE_headers += RAD_F.H
-
 CEXE_sources += trace_ppm_rad.cpp
 
 CEXE_headers += fluxlimiter.H
diff --git a/Source/radiation/RAD_F.H b/Source/radiation/RAD_F.H
deleted file mode 100644
index a1768d1aa1..0000000000
--- a/Source/radiation/RAD_F.H
+++ /dev/null
@@ -1,85 +0,0 @@
-#ifndef CASTRO_RAD_F_H
-#define CASTRO_RAD_F_H
-
-#include <AMReX_ArrayLim.H>
-
-#ifdef BL_FORT_USE_UNDERSCORE
-#  define   FORT_RADBNDRY2 rbndr2_
-
-#elif defined(BL_FORT_USE_UPPERCASE)
-
-#  define   FORT_RADBNDRY2 RBNDR2
-
-#endif
-
-#include <AMReX_BLFort.H>
-
-#ifdef __cplusplus
-extern "C" {
-#endif
-  void ca_filt_prim
-    (const int lo[], const int hi[],
-     BL_FORT_FAB_ARG(Stmp),
-     BL_FORT_FAB_ARG(Snew),
-     BL_FORT_FAB_ARG(mask),
-     const int* filter_prim_T, const int* filter_prim_S,
-     const int* domlo, const int* domhi,
-     const amrex::Real* delta, const amrex::Real* xlo, const amrex::Real* problo,
-     const amrex::Real* time, const int* level);
-#ifdef __cplusplus
-}
-#endif
-
-#ifdef __cplusplus
-extern "C"
-{
-#endif
-
-  void ca_flux_face2center
-   (const int lo[], const int hi[],
-    BL_FORT_FAB_ARG(test), 
-    const BL_FORT_FAB_ARG(flx),
-    const amrex::Real* x, const int* xlo, const int* xhi, 
-    const int* ntest, const int* idim, const int* iflx);
-
-#ifdef __cplusplus
-}
-#endif
-// </ MGFLD>
-
-#ifdef __cplusplus
-extern "C" {
-#endif
-  void fkpn(const int* lo, const int* hi,
-            BL_FORT_FAB_ARG_3D(fkp),       
-            amrex::Real con, amrex::Real em, amrex::Real en,
-            amrex::Real ep, amrex::Real xnu, amrex::Real tf,
-            BL_FORT_FAB_ARG_3D(temp), 
-            BL_FORT_FAB_ARG_3D(state),
-            amrex::Real temp_offset);
-
-  void rosse1(const int* lo, const int* hi,
-              amrex::Real con, amrex::Real em, amrex::Real en,
-              amrex::Real ep, amrex::Real nu,
-              amrex::Real tf, amrex::Real kfloor,
-              BL_FORT_FAB_ARG_3D(state),
-              BL_FORT_FAB_ARG_3D(kappar));
-
-  void rosse1s(const int* lo, const int* hi,
-               amrex::Real con, amrex::Real em, amrex::Real en,
-               amrex::Real ep,  amrex::Real scon,
-               amrex::Real sem, amrex::Real sen,
-               amrex::Real sep, amrex::Real nu,
-               amrex::Real tf, amrex::Real kfloor,
-               BL_FORT_FAB_ARG_3D(state),
-               BL_FORT_FAB_ARG_3D(kappar));
-  
-  void FORT_RADBNDRY2(amrex::Real* bf, ARLIM_P(blo), ARLIM_P(bhi), 
-                      int* tfab, ARLIM_P(dlo), ARLIM_P(dhi),
-                      const amrex::Real* dx, const amrex::Real* xlo, const amrex::Real& time);
-
-#ifdef __cplusplus
-}
-#endif
-
-#endif
diff --git a/Source/radiation/RadBndry.cpp b/Source/radiation/RadBndry.cpp
index 47af9e86f8..aa742febef 100644
--- a/Source/radiation/RadBndry.cpp
+++ b/Source/radiation/RadBndry.cpp
@@ -9,8 +9,6 @@
 #include <omp.h>
 #endif
 
-#include <RAD_F.H>
-
 using namespace amrex;
 
 int         RadBndry::first = 1;
diff --git a/Source/radiation/RadMultiGroup.cpp b/Source/radiation/RadMultiGroup.cpp
index 6d10bced88..83a1811ffc 100644
--- a/Source/radiation/RadMultiGroup.cpp
+++ b/Source/radiation/RadMultiGroup.cpp
@@ -4,8 +4,6 @@
 
 #include <Radiation.H>
 
-#include <RAD_F.H>
-
 #include <AMReX_ParmParse.H>
 
 #include <iostream>
diff --git a/Source/radiation/RadPlotvar.cpp b/Source/radiation/RadPlotvar.cpp
index 15ed08f304..33b90393a7 100644
--- a/Source/radiation/RadPlotvar.cpp
+++ b/Source/radiation/RadPlotvar.cpp
@@ -2,8 +2,6 @@
 #include <Radiation.H>
 #include <fluxlimiter.H>
 
-#include <RAD_F.H>
-
 using namespace amrex;
 
 void Radiation::save_lambda_in_plotvar(int level, const Array<MultiFab,AMREX_SPACEDIM>& lambda)
diff --git a/Source/radiation/RadSolve.cpp b/Source/radiation/RadSolve.cpp
index ea8e78d0df..372384d415 100644
--- a/Source/radiation/RadSolve.cpp
+++ b/Source/radiation/RadSolve.cpp
@@ -7,7 +7,6 @@
 #include <Radiation.H>  // for access to static physical constants only
 #include <rad_util.H>
 #include <problem_rad_source.H>
-#include <RAD_F.H>
 
 #include <iostream>
 
diff --git a/Source/radiation/Radiation.cpp b/Source/radiation/Radiation.cpp
index f1932f2a74..81a69d64bd 100644
--- a/Source/radiation/Radiation.cpp
+++ b/Source/radiation/Radiation.cpp
@@ -5,7 +5,6 @@
 #include <rad_util.H>
 #include <filt_prim.H>
 
-#include <RAD_F.H>
 #include <AMReX_PROB_AMR_F.H>
 
 #include <opacity.H>
diff --git a/Source/radiation/SGRadSolver.cpp b/Source/radiation/SGRadSolver.cpp
index 7ad805a38c..eca2c5db37 100644
--- a/Source/radiation/SGRadSolver.cpp
+++ b/Source/radiation/SGRadSolver.cpp
@@ -5,9 +5,6 @@
 #include <Radiation.H>
 #include <RadSolve.H>
 
-
-#include <RAD_F.H>
-
 #include <iostream>
 
 #ifdef _OPENMP