test_gas.py 文件源码

def test_total_number_large_aperture(self):
        """Test column density for aperture >> lengthscale."""
        Q = 1 / u.s
        v = 1 * u.km / u.s
        rho = 1000 * u.km
        parent = 10 * u.km
        N = Haser(Q, v, parent).total_number(rho)
        ideal = Q * parent / v
        assert np.isclose(N, ideal.decompose().value)

    # TEST FAILS
    #
    # def test_total_number_rho_NJ78(self):
    #     """Reproduce Newburn and Johnson 1978.

    #     Species, N observed, Q/v (km**-1)
    #     CN, 6.4e26, 5.8e23
    #     C3, 8.3e28, 9.0e23
    #     C2, 7.8e27, 5.9e24

    #     Cannot reproduce C3 quoted in paper.

    #     """

    #     #Nobs = [6.41756750e26, 8.63191842e+28, 7.81278300e27]
    #     #Nobs = [6.4e26, 4.2e27, 7.8e27]
    #     Nobs = [6.4e26, 8.3e28, 7.8e27]
    #     parent = [1.4e4, 0, 1.0e4] * u.km
    #     daughter = [1.7e5, 4.6e4, 7.6e4] * u.km
    #     Q = [5.8e23, 9.0e23, 5.9e24] / u.s
    #     rho = 3300 * u.km

    #     N = np.zeros(3)
    #     for i in range(3):
    #         coma = Haser(Q[i], 1 * u.km / u.s, parent[i], daughter[i])
    #         N[i] = coma.total_number(rho)

    #     assert np.allclose(N, Nobs)