real(dp) function cepa_aqcc(run) integer, intent(in) :: run cepa_aqcc = diagsft(run) * (1.0_dp - aqcc_factor) end function cepa_aqcc