diff --git a/process/fw.py b/process/fw.py
index 8f150b7280..084d575915 100644
--- a/process/fw.py
+++ b/process/fw.py
@@ -144,23 +144,7 @@ def fw_temp(
 
         # Model B is given for comparison
         # Model C is used
-        # Model A is not used.
-
-        # Model A: LeClaire formula for circular pipes (NOT USED) !
-        # !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
-        # This gives a very small temperature differential, as the wall thickness is
-        # the same all the way round, but the surface area is bigger than that of a flat wall.
-
-        # Calculate peak temperature - occurs at (r,theta) = (bfw,0)
-        # bfw = thickness/2.0d0
-
-        # po.cosine_term(afw, bfw, 0.0D0, bfw, qpp, hcoeff, tkfw, tmthet)
-
-        # deltat_solid = bfw/tkfw * (qpp/pi + qppp*bfw/2.0D0) * log(bfw/afw) &
-        #             - qppp/4.0D0/tkfw*(bfw**2-afw**2)
-
-        # tpeakfw = fwoutlet + deltat_solid + deltat_coolant + tmthet  # in K
+        # Model A LeClaire formula for circular pipes (removed GitHub #389).
 
         # Model B: Simple 1-dimensional calculation !
         # !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
diff --git a/source/fortran/pulse_variables.f90 b/source/fortran/pulse_variables.f90
index 791fdec141..eebf0b241b 100644
--- a/source/fortran/pulse_variables.f90
+++ b/source/fortran/pulse_variables.f90
@@ -19,9 +19,6 @@ module pulse_variables
   real(dp) :: bctmp
   !! first wall bulk coolant temperature (C)
 
-  real(dp) :: bfw
-  !! outer radius of each first wall structural tube (m) (0.5 * average of fwith and fwoth)
-
   real(dp) :: dtstor
   !! maximum allowable temperature change in stainless steel thermal storage block (K) (`istore=3`)
 
@@ -52,7 +49,6 @@ subroutine init_pulse_variables
     implicit none
 
     bctmp = 320.0D0
-    bfw = 0.0D0
     dtstor = 300.0D0
     istore = 1
     itcycl = 1