Merge pull request #62 from CliMA/gb/update_interpolations

Update interpolations
CliMA · Dec 5, 2023 · f3de10f · f3de10f
2 parents f7be333 + 605515e
commit f3de10f
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/Project.toml b/Project.toml
@@ -13,5 +13,5 @@ Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
 Artifacts = "<0.0.1, 1"
 Dates = "<0.0.1, 1"
 DelimitedFiles = "<0.0.1, 1"
-Interpolations = "0.13, 0.14"
+Interpolations = "0.13, 0.14, 0.15"
 julia = "1.7, 1.8, 1.9"
diff --git a/src/Insolation.jl b/src/Insolation.jl
@@ -12,7 +12,7 @@ export orbital_params
 """
     OrbitalData
 
-The parameters vary due to Milankovitch cycles. 
+The parameters vary due to Milankovitch cycles.
 
 Orbital parameters from the Laskar 2004 paper are
 lazily downloaded from Caltech Box to the
@@ -31,11 +31,11 @@ function OrbitalData()
     x, _ = readdlm(datapath, ',', Float64, header = true)
     t_range = (x[1, 1] * 1e3):1e3:(x[end, 1] * 1e3) # array of every 1 kyr to range of years
     e_spline_etp =
-        CubicSplineInterpolation(t_range, x[:, 2], extrapolation_bc = NaN)
+        cubic_spline_interpolation(t_range, x[:, 2], extrapolation_bc = NaN)
     γ_spline_etp =
-        CubicSplineInterpolation(t_range, x[:, 3], extrapolation_bc = NaN)
+        cubic_spline_interpolation(t_range, x[:, 3], extrapolation_bc = NaN)
     ϖ_spline_etp =
-        CubicSplineInterpolation(t_range, x[:, 4], extrapolation_bc = NaN)
+        cubic_spline_interpolation(t_range, x[:, 4], extrapolation_bc = NaN)
 
     E = typeof(e_spline_etp)
     G = typeof(γ_spline_etp)