def _calc_gamma(self, V_a):
        V0, gamma0, gamma_inf, beta = self.eos_mod.get_param_values(
            param_names=['V0','gamma0','gamma_inf','beta'])

        gamma_a = gamma_inf + (gamma0-gamma_inf)*(V_a/V0)**beta
        return gamma_a