diff --git a/CHANGES.md b/CHANGES.md
index 7c6b05f31d..19e9f97746 100644
--- a/CHANGES.md
+++ b/CHANGES.md
@@ -1,5 +1,8 @@
 # 23.11
 
+  * The ca_rad_source hook in Fortran has been removed. The existing
+    problem_rad_source() hook in C++ can be used instead. (#2626)
+
   * The compile option USE_AUX_UPDATE has been removed. If you want to
     manually update the auxiliary parameters, you can use an external
     source term or you can use the problem post-timestep hook. (#2614)
diff --git a/Source/radiation/Make.package b/Source/radiation/Make.package
index 5c1866acf5..977a82fee5 100644
--- a/Source/radiation/Make.package
+++ b/Source/radiation/Make.package
@@ -47,5 +47,3 @@ CEXE_sources += RadDerive.cpp
 CEXE_headers += RadDerive.H
 CEXE_headers += rad_util.H
 CEXE_headers += blackbody.H
-
-ca_F90EXE_sources += rad_source.F90
diff --git a/Source/radiation/RAD_F.H b/Source/radiation/RAD_F.H
index 55e0154615..079dedba62 100644
--- a/Source/radiation/RAD_F.H
+++ b/Source/radiation/RAD_F.H
@@ -65,16 +65,6 @@ extern "C" {
 }
 #endif
 
-extern "C"
-{
-  void ca_rad_source
-    (const int* lo, const int* hi, 
-     BL_FORT_FAB_ARG_3D(rhs),
-     const amrex::Real* dx, amrex::Real dt,
-     amrex::Real time, int igroup);
-
-}
-
 #ifdef __cplusplus
 extern "C" {
 #endif
diff --git a/Source/radiation/RadSolve.cpp b/Source/radiation/RadSolve.cpp
index 56157ec59a..e671cf66f1 100644
--- a/Source/radiation/RadSolve.cpp
+++ b/Source/radiation/RadSolve.cpp
@@ -1102,11 +1102,6 @@ void RadSolve::levelRhs(int level, MultiFab& rhs, const MultiFab& jg,
 
           problem_rad_source(i, j, k, rhs_arr, geomdata, time, delta_t, igroup);
       });
-
-      ca_rad_source(AMREX_INT_ANYD(bx.loVect()), AMREX_INT_ANYD(bx.hiVect()),
-                    BL_TO_FORTRAN_ANYD(rhs[ri]),
-                    AMREX_REAL_ANYD(dx), delta_t, time, igroup);
-
   }
 }
 
diff --git a/Source/radiation/rad_source.F90 b/Source/radiation/rad_source.F90
deleted file mode 100644
index bb39ab1043..0000000000
--- a/Source/radiation/rad_source.F90
+++ /dev/null
@@ -1,39 +0,0 @@
-module rad_source_module
-
-  use amrex_fort_module, only: rt => amrex_real
-
-  implicit none
-
-contains
-
-  subroutine ca_rad_source(lo, hi, &
-                           rhs, rhs_lo, rhs_hi, &
-                           dx, dt, time, igroup) &
-                           bind(C, name="ca_rad_source")
-
-    implicit none
-
-    integer,  intent(in   ) :: lo(3), hi(3) 
-    integer,  intent(in   ) :: rhs_lo(3), rhs_hi(3)
-    real(rt), intent(inout) :: rhs(rhs_lo(1):rhs_hi(1),rhs_lo(2):rhs_hi(2),rhs_lo(3):rhs_hi(3))
-    real(rt), intent(in   ) :: dx(3)
-    real(rt), intent(in   ), value :: dt, time
-    integer,  intent(in   ), value :: igroup
-
-    integer :: i, j, k
-
-    do k = lo(3), hi(3)
-       do j = lo(2), hi(2)
-          do i = lo(1), hi(1)
-
-             ! This is a dummy routine that does nothing, but could be implemented
-             ! by a problem setup to inject a radiation source to group igroup.
-             ! This routine must add to rhs.
-
-          end do
-       end do
-    end do
-
-  end subroutine ca_rad_source
-
-end module rad_source_module