function cc_triples_factor() result(factor) real(dp) :: factor ! for now just set that to 1 factor = 1.0_dp ! essentially we would need to have the influence of the ! quadrupels and the quintuples.. but those numbers are hard to ! estimate.. ! todo end function cc_triples_factor