fix: Add missing quarter-mode mass function to BBKS transfer function…

… class

source/structure_formation.transfer_function.BBKS.F90

@@ -48,6 +48,7 @@ A bbks transfer function class.
      procedure :: value                 => bbksValue
      procedure :: logarithmicDerivative => bbksLogarithmicDerivative
      procedure :: halfModeMass          => bbksHalfModeMass
+     procedure :: quarterModeMass       => bbksQuarterModeMass
      procedure :: epochTime             => bbksEpochTime
   end type transferFunctionBBKS
 
@@ -267,6 +268,26 @@ Compute the mass corresponding to the wavenumber at which the transfer function
     return
   end function bbksHalfModeMass
 
+  double precision function bbksQuarterModeMass(self,status)
+    !!{
+    Compute the mass corresponding to the wavenumber at which the transfer function is suppressed by a factor of four relative
+    to a \gls{cdm} transfer function. Not supported in this implementation.
+    !!}
+    use :: Galacticus_Error, only : Galacticus_Error_Report, errorStatusFail
+    implicit none
+    class  (transferFunctionBBKS), intent(inout)           :: self
+    integer                      , intent(  out), optional :: status
+    !$GLC attributes unused :: self
+
+    bbksQuarterModeMass=0.0d0
+    if (present(status)) then
+       status=errorStatusFail
+    else
+       call Galacticus_Error_Report('not supported by this implementation'//{introspection:location})
+    end if
+    return
+  end function bbksQuarterModeMass
+
   double precision function bbksEpochTime(self)
     !!{
     Return the cosmic time at the epoch at which this transfer function is defined.