Hack to avoid including AstroLib due to it breaking compat

eford · eford · commit 0e4564f19e24 · 2021-04-21T20:58:19.000-04:00
diff --git a/Project.toml b/Project.toml
@@ -4,7 +4,6 @@ authors = ["Eric B. Ford, Christian Gilbertson, Joe Ninan, Michael L. Palumbo II
 version = "0.2.0"
 
 [deps]
-AstroLib = "c7932e45-9af1-51e7-9da9-f004cd3a462b"
 CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
 DSP = "717857b8-e6f2-59f4-9121-6e50c889abd2"
 DataFrames = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
@@ -37,7 +36,6 @@ TemporalGPs = "e155a3c4-0841-43e1-8b83-a0e4f03cc18f"
 ThreadedIterables = "11d239b0-c0b9-11e8-1935-d5cfa53abb03"
 
 [compat]
-AstroLib = "0.4"
 CSV = "0.7, 0.8"
 DSP = "0.5, 0.6"
 DataFrames = "0.20, 0.21, 0.22"
@@ -60,12 +58,12 @@ Query = "1.0"
 RvSpectMLBase = "0.2"
 Scalpels = "0.1"
 SpecialFunctions = "0.10, 1.0, 1.1"
-StaticArrays = "0.12, 1.0"
+StaticArrays = "0.8, 0.9, 0.10, 0.11, 0.12, 1"
 StatsBase = "0.33"
 Stheno = "0.6"
-TemporalGPs = "0.3"
+TemporalGPs = "0.3, 0.4, 0.5"
 ThreadedIterables = "0.2"
-julia = "1.6"
+julia = "1.5, 1.6"
 
 [extras]
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
diff --git a/src/RvSpectML.jl b/src/RvSpectML.jl
@@ -35,7 +35,8 @@ export Scalpels
 using LinearAlgebra, Statistics
 using DataFrames, Query
 using Dates
-using AstroLib  # For solar position
+#using AstroLib  # For solar position
+include("astrolib/astrolib.jl")
 # Packages that are being used and likely can be shared
 # using Distributions, Interpolations, MultivariateStats, PDMats
 # Packages we might use soon
diff --git a/src/astrolib/astrolib.jl b/src/astrolib/astrolib.jl
@@ -0,0 +1,20 @@
+# AstroLib.jl[https://github.com/JuliaAstro/AstroLib.jl] was breaking compatability with other packages
+# So this is a hack to include the just needed functions and avoid extra dependancies
+
+module AstroLib
+
+using Dates
+const JULIANCENTURY = 36_525
+const ct2lst_c  = (280.46061837, 360.98564736629, 0.000387933, 38710000.0)
+include("jdcnv.jl")
+const J2000 = Int(jdcnv(2000, 01, 01, 12)) # 2451545
+
+include("sec2rad.jl")
+include("ct2lst.jl")
+export ct2lst
+include("hadec2altaz.jl")
+export hadec2altaz
+include("sunpos.jl")
+export sunpos
+
+end
diff --git a/src/astrolib/ct2lst.jl b/src/astrolib/ct2lst.jl
@@ -0,0 +1,111 @@
+# This file is a part of AstroLib.jl. License is MIT "Expat".
+# Copyright (C) 2016 Mosè Giordano.
+
+function ct2lst(long::T, jd::T) where {T<:AbstractFloat}
+    t0 = jd - J2000
+    t  = t0 / JULIANCENTURY
+    # Compute GST in seconds.
+    θ = ct2lst_c[1] + (ct2lst_c[2]*t0) + t*t*(ct2lst_c[3] - t / ct2lst_c[4])
+    return mod((θ + long)/15, 24)
+end
+
+"""
+    ct2lst(longitude, jd) -> local_sidereal_time
+    ct2lst(longitude, tz, date) -> local_sidereal_time
+
+### Purpose ###
+
+Convert from Local Civil Time to Local Mean Sidereal Time.
+
+### Arguments ###
+
+The function can be called in two different ways.  The only argument common to
+both methods is `longitude`:
+
+* `longitude`: the longitude in degrees (east of Greenwich) of the place for which the local
+  sidereal time is desired.  The Greenwich mean sidereal time (GMST) can be found by setting
+  longitude = `0`.
+
+The civil date to be converted to mean sidereal time can be specified either by
+providing the Julian days:
+
+* `jd`: this is number of Julian days for the date to be converted.
+
+or the time zone and the date:
+
+* `tz`: the time zone of the site in hours, positive East of the
+  Greenwich meridian (ahead of GMT).  Use this parameter to easily
+  account for Daylight Savings time (e.g. -4=EDT, -5 = EST/CDT).
+* `date`: this is the local civil time with type `DateTime`.
+
+### Output ###
+
+The local sidereal time for the date/time specified in hours.
+
+### Method ###
+
+The Julian days of the day and time is question is used to determine the number
+of days to have passed since 2000-01-01.  This is used in conjunction with the
+GST of that date to extrapolate to the current GST; this is then used to get the
+LST.  See Astronomical Algorithms by Jean Meeus, p. 84 (Eq. 11-4) for the
+constants used.
+
+### Example ###
+
+Find the Greenwich mean sidereal time (GMST) on 2008-07-30 at 15:53 in
+Baltimore, Maryland (longitude=-76.72 degrees).  The timezone is EDT or tz=-4
+
+```jldoctest
+julia> using AstroLib, Dates
+
+julia> lst = ct2lst(-76.72, -4, DateTime(2008, 7, 30, 15, 53))
+11.356505172312609
+
+julia> sixty(lst)
+3-element StaticArrays.SArray{Tuple{3},Float64,1,3} with indices SOneTo(3):
+ 11.0
+ 21.0
+ 23.418620325392112
+```
+
+Find the Greenwich mean sidereal time (GMST) on 2015-11-24 at 13:21 in
+Heidelberg, Germany (longitude=08° 43' E).  The timezone is CET or tz=1.
+Provide `ct2lst` only with the longitude of the place and the number of
+Julian days.
+
+```jldoctest
+julia> using AstroLib, Dates
+
+julia> longitude=ten(8, 43); # Convert longitude to decimals.
+
+julia> jd = jdcnv(DateTime(2015, 11, 24, 13, 21) - Dates.Hour(1));
+# Get number of Julian days. Remember to subtract the time zone in
+# order to convert local time to UTC.
+
+julia> lst = ct2lst(longitude, jd) # Calculate Greenwich Mean Sidereal Time.
+17.140685171005316
+
+julia> sixty(lst)
+3-element StaticArrays.SArray{Tuple{3},Float64,1,3} with indices SOneTo(3):
+ 17.0
+  8.0
+ 26.466615619137883
+```
+
+### Notes ###
+
+Code of this function is based on IDL Astronomy User's Library.
+"""
+ct2lst(long::Real, jd::Real) = ct2lst(promote(float(long), float(jd))...)
+
+function ct2lst(long::T, tz::T, date::DateTime) where {T<:AbstractFloat}
+    # In order to handle time zones, package "TimeZones.jl" is much better, but
+    # here we need only to add the time zone to UTC time.  All time zones I know
+    # are either integer, or ±30 or ±45 minutes, so it should be safe enough to
+    # convert hours to minutes and subtract minutes from "time".
+    date = date - Dates.Minute(round(Int, tz*60))
+    return ct2lst(long, jdcnv(date))
+end
+
+ct2lst(long::Real, tz::Real, date::DateTime) =
+    ct2lst(promote(float(long), float(tz))..., date)
diff --git a/src/astrolib/hadec2altaz.jl b/src/astrolib/hadec2altaz.jl
@@ -0,0 +1,96 @@
+# This file is a part of AstroLib.jl. License is MIT "Expat".
+# Copyright (C) 2016 Mosè Giordano.
+
+function _hadec2altaz(ha::T, dec::T, lat::T, ws::Bool) where {T<:AbstractFloat}
+    sh, ch = sincos(deg2rad(ha))
+    sd, cd = sincos(deg2rad(dec))
+    sl, cl = sincos(deg2rad(lat))
+
+    x = -ch*cd*sl + sd*cl
+    y = -sh*cd
+    z = ch*cd*cl + sd*sl
+    r = hypot(x, y)
+
+    # Now get altitude, azimuth
+    az  = rad2deg(mod2pi(atan(y, x)))
+    alt = rad2deg(atan(z, r))
+    # Convert azimuth to West from South, if desired
+    if ws
+        az = mod(az + 180, 360)
+    end
+    return alt, az
+end
+
+"""
+    hadec2altaz(ha, dec, lat[, ws=true]) -> alt, az
+
+### Purpose ###
+
+Convert Hour Angle and Declination to Horizon (Alt-Az) coordinates.
+
+### Explanation ###
+
+Can deal with the NCP singularity.  Intended mainly to be used by program
+`eq2hor`.
+
+### Arguments ###
+
+Input coordinates may be either a scalar or an array, of the same dimension.
+
+* `ha`: the local apparent hour angle, in degrees.  The hour angle is the time
+  that right ascension of 0 hours crosses the local meridian.  It is
+  unambiguously defined.
+* `dec`: the local apparent declination, in degrees.
+* `lat`: the local geodetic latitude, in degrees, scalar or array.
+* `ws` (optional boolean keyword): if true, the output azimuth is measured West
+  from South.  The default is to measure azimuth East from North.
+
+`ha` and `dec` can be given as a 2-tuple `(ha, dec)`.
+
+### Output ###
+
+2-tuple `(alt, az)`
+
+* `alt`: local apparent altitude, in degrees.
+* `az`: the local apparent azimuth, in degrees.
+
+The output coordinates are always floating points and have the same type (scalar
+or array) as the input coordinates.
+
+### Example ###
+
+Arcturus is observed at an apparent hour angle of 336.6829 and a declination of
+19.1825 while at the latitude of +43° 4' 42''.  What are the local altitude and
+azimuth of this object?
+
+```jldoctest
+julia> using AstroLib
+
+julia> alt, az = hadec2altaz(336.6829, 19.1825, ten(43, 4, 42))
+(59.08617155005685, 133.3080693440254)
+```
+
+### Notes ###
+
+`altaz2hadec` converts Horizon (Alt-Az) coordinates to Hour Angle and
+Declination.
+
+Code of this function is based on IDL Astronomy User's Library.
+"""
+hadec2altaz(ha::Real, dec::Real, lat::Real; ws::Bool=false) =
+    _hadec2altaz(promote(float(ha), float(dec), float(lat))..., ws)
+
+hadec2altaz(hadec::Tuple{Real, Real}, lat::Real; ws::Bool=false) =
+    hadec2altaz(hadec..., lat, ws=ws)
+
+function hadec2altaz(ha::AbstractArray{R}, dec::AbstractArray{<:Real},
+                     lat::AbstractArray{<:Real}; ws::Bool=false) where {R<:Real}
+    @assert length(ha) == length(dec) == length(lat)
+    typeha = float(R)
+    alt = similar(ha, typeha)
+    az  = similar(ha, typeha)
+    for i in eachindex(ha)
+        alt[i], az[i] = hadec2altaz(ha[i], dec[i], lat[i], ws=ws)
+    end
+    return alt, az
+end
diff --git a/src/astrolib/jdcnv.jl b/src/astrolib/jdcnv.jl
@@ -0,0 +1,53 @@
+# This file is a part of AstroLib.jl. License is MIT "Expat".
+# Copyright (C) 2016 Mosè Giordano.
+
+"""
+    jdcnv(date) -> julian_days
+
+### Purpose ###
+
+Convert proleptic Gregorian Calendar date in UTC standard to number of Julian
+days.
+
+### Explanation ###
+
+Takes the given proleptic Gregorian date in UTC standard and returns the number
+of Julian calendar days since epoch `-4713-11-24T12:00:00`.
+
+### Argument ###
+
+* `date`: date in proleptic Gregorian Calendar.  Each element can be either a `DateTime` or
+  anything that can be converted directly to `DateTime`.
+
+### Output ###
+
+Number of Julian days, as a floating point.
+
+### Example ###
+
+Find the Julian days number at 2016 August 23, 03:39:06.
+
+```jldoctest
+julia> using AstroLib, Dates
+
+julia> jdcnv(DateTime(2016, 08, 23, 03, 39, 06))
+2.4576236521527776e6
+
+julia> jdcnv(2016, 08, 23, 03, 39, 06)
+2.4576236521527776e6
+
+julia> jdcnv("2016-08-23T03:39:06")
+2.4576236521527776e6
+```
+
+### Notes ###
+
+This is the inverse of `daycnv`.
+
+`get_juldate` returns the number of Julian days for current time.  It is
+equivalent to `jdcnv(now(Dates.UTC))`.
+
+For the conversion of Julian date to number of Julian days, use `juldate`.
+"""
+const jdcnv = Dates.datetime2julian
+jdcnv(date...) = jdcnv(DateTime(date...))
diff --git a/src/astrolib/sec2rad.jl b/src/astrolib/sec2rad.jl
@@ -0,0 +1,32 @@
+# This file is a part of AstroLib.jl. License is MIT "Expat".
+# Copyright (C) 2016 Mosè Giordano.
+
+"""
+    sec2rad(sec) -> radians
+
+### Purpose ###
+
+Convert from seconds to radians.
+
+### Argument ###
+
+* `sec`: number of seconds.
+
+### Output ###
+
+The number of radians corresponding to `sec`.
+
+### Example ###
+
+```jldoctest
+julia> using AstroLib
+
+julia> sec2rad(3600 * 30)
+0.5235987755982988
+```
+
+### Notes ###
+
+Use `rad2sec` to convert radians to seconds.
+"""
+sec2rad(sec::Real) = deg2rad(sec / 3600)
diff --git a/src/astrolib/sunpos.jl b/src/astrolib/sunpos.jl
