SPARC, MW RC

Author: islent

Project.toml

@@ -14,6 +14,7 @@ Dates = "ade2ca70-3891-5945-98fb-dc099432e06a"
 Dierckx = "39dd38d3-220a-591b-8e3c-4c3a8c710a94"
 Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
+Measurements = "eff96d63-e80a-5855-80a2-b1b0885c5ab7"
 Optim = "429524aa-4258-5aef-a3af-852621145aeb"
 PhysicalConstants = "5ad8b20f-a522-5ce9-bfc9-ddf1d5bda6ab"
 PhysicalParticles = "78724250-4eb2-11e9-0f34-13a7aab864c5"
@@ -36,6 +37,7 @@ DataFrames = "1"
 Dierckx = "0.5"
 Distributions = "0.25"
 DocStringExtensions = "0.8, 0.9"
+Measurements = "2.14.1"
 Optim = "1"
 PhysicalConstants = "0.2"
 PhysicalParticles = "1"

src/AstroIC.jl

@@ -17,6 +17,7 @@ using QuadGK
 using Optim
 using CSV
 using UnitfulAngles
+using Measurements
 
 @reexport using PhysicalParticles
 @reexport using AstroIO
@@ -75,6 +76,13 @@ include("solarsystem.jl")
 include("data/MW_satellites.jl")
 export load_data_MW_satellites
 
+include("data/SPARC.jl")
+export load_SPARC_LTGs_RC, load_SPARC_LTGs_data, load_li2018_SPARC
+export load_SPARC_Xray_ETGs_data, load_SPARC_rotating_ETGs_data, load_SPARC_rotating_ETGs_RC, load_SPARC_rotating_ETGs_rotmod
+
+include("data/MilkyWay.jl")
+export load_MW_RC_Eilers2019, load_MW_RC_Mroz2019, load_MW_RC_stddev_W21, load_MW_RC_DS_W21
+
 include("precompile.jl")
 
 """

src/data/MilkyWay.jl

@@ -0,0 +1,35 @@
+"""
+$(TYPEDSIGNATURES)
+"""
+function load_MW_RC_Eilers2019(filename = joinpath(@__DIR__, "MilkyWay/MW_RC_Eilers2019.dat"))
+    dfEilers2019 = DataFrame(CSV.File(filename; header = false, delim=" ", ignorerepeated = true))
+    rename!(dfEilers2019, [:r, :v, :σ_low, :σ_high])
+    return dfEilers2019
+end
+
+"""
+$(TYPEDSIGNATURES)
+"""
+function load_MW_RC_Mroz2019(filename = joinpath(@__DIR__, "MilkyWay/MW_RC_Mroz2019.dat"))
+    dfMroz2019 = DataFrame(CSV.File(filename; header = false, delim=" ", ignorerepeated = true))
+    rename!(dfMroz2019, [:r, :v, :σ_low, :σ_high])
+    return dfMroz2019
+end
+
+"""
+$(TYPEDSIGNATURES)
+"""
+function load_MW_RC_stddev_W21(filename = joinpath(@__DIR__, "MilkyWay/MW_RC_stddev_W21.txt"))
+    dfstddev_W21 = DataFrame(CSV.File(filename; header = false, delim=" ", ignorerepeated = true))
+    rename!(dfstddev_W21, [:id, :r, :v, :σ_v, :σ_r])
+    return dfstddev_W21
+end
+
+"""
+$(TYPEDSIGNATURES)
+"""
+function load_MW_RC_DS_W21(filename = joinpath(@__DIR__, "MilkyWay/MW_dsvel-W21-RC.txt"))
+    dfDS_W21 = DataFrame(CSV.File(filename; header = false, skipto=3, delim="\t", ignorerepeated = true))
+    rename!(dfDS_W21, [:r, :v_b, :v_CDM, :v_QUMOND, :v_MOG])
+    return dfDS_W21
+end

src/data/MilkyWay/MW_RC_Eilers2019.dat

@@ -0,0 +1,38 @@
+5.27 226.83 1.91 1.90
+5.74 230.80 1.43 1.35
+6.23 231.20 1.70 1.10
+6.73 229.88 1.44 1.32
+7.22 229.61 1.37 1.11
+7.82 229.91 0.92 0.88
+8.19 228.86 0.80 0.67
+8.78 226.50 1.07 0.95
+9.27 226.20 0.72 0.62
+9.76 225.94 0.42 0.52
+10.26 225.68 0.44 0.40
+10.75 224.73 0.38 0.41
+11.25 224.02 0.33 0.54
+11.75 223.86 0.40 0.39
+12.25 222.23 0.51 0.37
+12.74 220.77 0.54 0.46
+13.23 220.92 0.57 0.40
+13.74 217.47 0.64 0.51
+14.24 217.31 0.77 0.66
+14.74 217.60 0.65 0.68
+15.22 217.07 1.06 0.80
+15.74 217.38 0.84 1.07
+16.24 216.14 1.20 1.48
+16.74 212.52 1.39 1.43
+17.25 216.41 1.44 1.85
+17.75 213.70 2.22 1.65
+18.24 207.89 1.76 1.88
+18.74 209.60 2.31 2.77
+19.22 206.45 2.54 2.36
+19.71 201.91 2.99 2.26
+20.27 199.84 3.15 2.89
+20.78 198.14 3.33 3.37
+21.24 195.30 5.99 6.50
+21.80 213.67 15.38 12.18
+22.14 176.97 28.58 18.57
+22.73 193.11 27.64 19.05
+23.66 176.63 18.67 16.74
+24.82 198.42 6.50 6.12