function cc_quads_factor() result(factor) real(dp) :: factor ! see above factor = 1.0_dp end function cc_quads_factor