From fd0194301578e14f1e359097ccd984abed678ed6 Mon Sep 17 00:00:00 2001
From: Gabriele Bozzola <sbozzolator@gmail.com>
Date: Tue, 2 Jul 2024 19:35:56 -0700
Subject: [PATCH 1/4] Remove bors

---
 README.md |  4 ----
 bors.toml | 10 ----------
 2 files changed, 14 deletions(-)
 delete mode 100644 bors.toml

diff --git a/README.md b/README.md
index 25ea0a5..6670796 100644
--- a/README.md
+++ b/README.md
@@ -7,7 +7,6 @@ A library of atmospheric profiles
 | **Docs Build**       | [![docs build][docs-bld-img]][docs-bld-url]   |
 | **Documentation**    | [![dev][docs-dev-img]][docs-dev-url]          |
 | **GHA CI**           | [![gha ci][gha-ci-img]][gha-ci-url]           |
-| **Bors enabled**     | [![bors][bors-img]][bors-url]                 |
 
 [docs-bld-img]: https://github.com/CliMA/AtmosphericProfilesLibrary.jl/actions/workflows/docs.yml/badge.svg
 [docs-bld-url]: https://github.com/CliMA/AtmosphericProfilesLibrary.jl/actions/workflows/docs.yml
@@ -18,8 +17,5 @@ A library of atmospheric profiles
 [gha-ci-img]: https://github.com/CliMA/AtmosphericProfilesLibrary.jl/actions/workflows/ci.yml/badge.svg
 [gha-ci-url]: https://github.com/CliMA/AtmosphericProfilesLibrary.jl/actions/workflows/ci.yml
 
-[bors-img]: https://bors.tech/images/badge_small.svg
-[bors-url]: https://app.bors.tech/repositories/41911
-
 See our [docs](https://CliMA.github.io/AtmosphericProfilesLibrary.jl/dev/) for usage and how to contribute.
 
diff --git a/bors.toml b/bors.toml
deleted file mode 100644
index d088ad1..0000000
--- a/bors.toml
+++ /dev/null
@@ -1,10 +0,0 @@
-status = [
-  "docbuild",
-  "ci 1.7.0 - ubuntu-latest",
-  "ci 1.7.0 - windows-latest",
-  "ci 1.7.0 - macOS-latest",
-]
-delete_merged_branches = true
-timeout_sec = 10800
-block_labels = [ "do-not-merge-yet" ]
-cut_body_after = "<!--"

From bf6112bdc50d53bd8d35a6f1522b8b6d5cfc66f5 Mon Sep 17 00:00:00 2001
From: Gabriele Bozzola <sbozzolator@gmail.com>
Date: Tue, 2 Jul 2024 19:45:57 -0700
Subject: [PATCH 2/4] Remove Dierckx, use Interpolations instead

---
 Project.toml                      |  4 ++--
 docs/Project.toml                 |  4 ++--
 docs/src/define_save_plots.jl     |  1 -
 docs/src/index.md                 |  2 +-
 src/AtmosphericProfilesLibrary.jl |  2 +-
 src/profiles/ARM_SGP.jl           | 16 ++++++++--------
 src/profiles/Bomex.jl             |  2 +-
 src/profiles/DryBubble.jl         | 10 +++++-----
 src/profiles/Dycoms_RF01.jl       |  2 +-
 src/profiles/GABLS.jl             |  2 +-
 src/profiles/GATE_III.jl          | 12 ++++++------
 src/profiles/LifeCycleTan2018.jl  |  2 +-
 src/profiles/Nieuwstadt.jl        |  2 +-
 src/profiles/Rico.jl              |  2 +-
 src/profiles/Soares.jl            |  2 +-
 src/profiles/TRMM_LBA.jl          | 14 +++++++-------
 16 files changed, 39 insertions(+), 40 deletions(-)

diff --git a/Project.toml b/Project.toml
index c7a7459..186dc73 100644
--- a/Project.toml
+++ b/Project.toml
@@ -4,9 +4,9 @@ authors = ["Climate Modeling Alliance"]
 version = "0.1.5"
 
 [deps]
-Dierckx = "39dd38d3-220a-591b-8e3c-4c3a8c710a94"
+Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 
 [compat]
-Dierckx = "0.5"
+Interpolations = "0.14"
 julia = "1"
diff --git a/docs/Project.toml b/docs/Project.toml
index 3c2fcdf..f060804 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -1,10 +1,10 @@
 [deps]
 AtmosphericProfilesLibrary = "86bc3604-9858-485a-bdbe-831ec50de11d"
-Dierckx = "39dd38d3-220a-591b-8e3c-4c3a8c710a94"
 Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4"
 DocumenterCitations = "daee34ce-89f3-4625-b898-19384cb65244"
 Glob = "c27321d9-0574-5035-807b-f59d2c89b15c"
+Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 
 [compat]
-Dierckx = "0.5"
+Interpolations = "0.14"
diff --git a/docs/src/define_save_plots.jl b/docs/src/define_save_plots.jl
index 733060e..63f473b 100644
--- a/docs/src/define_save_plots.jl
+++ b/docs/src/define_save_plots.jl
@@ -1,5 +1,4 @@
 import Plots
-import Dierckx
 
 # https://github.com/jheinen/GR.jl/issues/278#issuecomment-587090846
 ENV["GKSwstype"] = "nul"
diff --git a/docs/src/index.md b/docs/src/index.md
index f4e029f..e897fa0 100644
--- a/docs/src/index.md
+++ b/docs/src/index.md
@@ -14,7 +14,7 @@ This repository is meant to serve as a library of static functions in efforts to
    - `prof = z->2*z` callable by `prof(z)`
    - `prof = t->2*t` callable by `prof(t)`
    - `prof = (t,z) -> z*t` callable by `prof(t, z)`
-   - `prof = Dierckx.Spline1D(FT[0, 1, 2], FT[0, 2, 4]; k = 1)` callable by `prof(z)`
+   - `prof = Interpolations.interpolate(FT[0, 1, 2], FT[0, 2, 4])` callable by `prof(z)`
 
 ## How to contribute
 
diff --git a/src/AtmosphericProfilesLibrary.jl b/src/AtmosphericProfilesLibrary.jl
index 507e98c..61dac70 100644
--- a/src/AtmosphericProfilesLibrary.jl
+++ b/src/AtmosphericProfilesLibrary.jl
@@ -1,6 +1,6 @@
 module AtmosphericProfilesLibrary
 
-import Dierckx
+import Interpolations as Intp
 
 abstract type AbstractProfile end
 struct TimeProfile{P} <: AbstractProfile
diff --git a/src/profiles/ARM_SGP.jl b/src/profiles/ARM_SGP.jl
index 251ec26..183484e 100644
--- a/src/profiles/ARM_SGP.jl
+++ b/src/profiles/ARM_SGP.jl
@@ -7,7 +7,7 @@ end
 function ARM_SGP_θ_liq_ice(::Type{FT}) where {FT}
     z_in = ARM_SGP_z(FT)
     θ_liq_ice_in = FT[299.0, 301.5, 302.5, 303.53, 303.7, 307.13, 314.0, 343.2] # K
-    profile = Dierckx.Spline1D(z_in, θ_liq_ice_in; k = 1)
+    profile = Intp.interpolate((z_in, ), θ_liq_ice_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 """ [Brown2002](@cite) """
@@ -18,7 +18,7 @@ function ARM_SGP_q_tot(::Type{FT}) where {FT}
     z_in = ARM_SGP_z(FT)
     r_in = FT[15.2,15.17,14.98,14.8,14.7,13.5,3.0,3.0] ./ 1000 # qt should be in kg/kg
     q_tot_in = r_in ./ (1 .+ r_in)
-    profile = Dierckx.Spline1D(z_in, q_tot_in; k = 1)
+    profile = Intp.interpolate((z_in, ), q_tot_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 """ [Brown2002](@cite) """
@@ -49,7 +49,7 @@ function ARM_SGP_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
 
@@ -66,8 +66,8 @@ function ARM_SGP_dTdt(::Type{FT}) where {FT}
     AT_in = FT[0.0, 0.0, 0.0, -0.08, -0.016, -0.016] ./ 3600
     # Radiative forcing for theta [K/h] converted to [K/sec]
     RT_in = FT[-0.125, 0.0, 0.0, 0.0, 0.0, -0.1] ./ 3600
-    dTdt_A = Dierckx.Spline1D(t_in, AT_in; k = 1)
-    dTdt_R = Dierckx.Spline1D(t_in, RT_in; k = 1)
+    dTdt_A = Intp.interpolate((t_in, ), AT_in, Intp.Gridded(Intp.Linear()))
+    dTdt_R = Intp.interpolate((t_in, ), RT_in, Intp.Gridded(Intp.Linear()))
     prof = (t, z) -> if z <= 1000.0
         FT(dTdt_A(t)+dTdt_R(t))
     elseif z > 1000.0 && z <= 2000.0
@@ -83,7 +83,7 @@ function ARM_SGP_dqtdt(::Type{FT}) where {FT}
     t_in = ARM_SGP_time(FT)
     # Radiative forcing for qt converted to [kg/kg/sec]
     Rqt_in = FT[0.08, 0.02, 0.04, -0.1, -0.16, -0.3] ./ 1000 ./ 3600
-    dqtdt = Dierckx.Spline1D(t_in, Rqt_in; k = 1)
+    dqtdt = Intp.interpolate((t_in, ), Rqt_in, Intp.Gridded(Intp.Linear()))
     prof = (Π, t, z) -> if z <= 1000.0
         FT(dqtdt(t) * Π)
     elseif z > 1000.0 && z <= 2000.0
@@ -98,7 +98,7 @@ end
 function ARM_SGP_shf(::Type{FT}) where {FT}
     t_Sur_in = FT[0.0, 4.0, 6.5, 7.5, 10.0, 12.5, 14.5] .* 3600 #LES time is in sec
     shf = FT[-30.0, 90.0, 140.0, 140.0, 100.0, -10, -10] # W/m^2
-    profile = Dierckx.Spline1D(t_Sur_in, shf; k = 1)
+    profile = Intp.interpolate((t_Sur_in, ), shf, Intp.Gridded(Intp.Linear()))
     return TimeProfile(profile)
 end
 
@@ -106,6 +106,6 @@ end
 function ARM_SGP_lhf(::Type{FT}) where {FT}
     t_Sur_in = FT[0.0, 4.0, 6.5, 7.5, 10.0, 12.5, 14.5] .* 3600 #LES time is in sec
     lhf = FT[5.0, 250.0, 450.0, 500.0, 420.0, 180.0, 0.0] # W/m^2
-    profile = Dierckx.Spline1D(t_Sur_in, lhf; k = 1)
+    profile = Intp.interpolate((t_Sur_in, ), lhf, Intp.Gridded(Intp.Linear()))
     return TimeProfile(profile)
 end
diff --git a/src/profiles/Bomex.jl b/src/profiles/Bomex.jl
index 1a9e06b..cbdf900 100644
--- a/src/profiles/Bomex.jl
+++ b/src/profiles/Bomex.jl
@@ -45,7 +45,7 @@ function Bomex_tke_prescribed(::Type{FT}) where {FT}
                 0.1468, 0.1470, 0.1458, 0.1423, 0.1351, 0.1224, 0.1033, 0.0793, 0.0542, 0.0322,
                 0.0162, 0.0068, 0.0024, 0.0007, 0.0001, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
 
diff --git a/src/profiles/DryBubble.jl b/src/profiles/DryBubble.jl
index 445db6b..3aaf111 100644
--- a/src/profiles/DryBubble.jl
+++ b/src/profiles/DryBubble.jl
@@ -61,7 +61,7 @@ function DryBubble_θ_liq_ice(::Type{FT}) where {FT}
         299.9837, 299.9837, 299.9837, 299.9837, 299.9837, 299.9837,
         299.9837, 299.9837
     ]
-    profile = Dierckx.Spline1D(z_in, θ_liq_ice_in; k = 1)
+    profile = Intp.interpolate((z_in, ), θ_liq_ice_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -97,7 +97,7 @@ function DryBubble_updrafts_θ_liq_ice(::Type{FT}) where {FT}
         300.5212, 300.482 , 300.4272, 300.3875, 300.3354, 300.2968,
         300.2587, 300.2216, 300.1782, 300.1452, 300.1143, 300.0859,
         300.0603, 300.0408, 300.0211, 300.0067, 299.9963, 299.9884]
-    profile = Dierckx.Spline1D(z_in, θ_liq_in; k = 1)
+    profile = Intp.interpolate((z_in, ), θ_liq_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -114,7 +114,7 @@ function DryBubble_updrafts_area(::Type{FT}) where {FT}
         0.15 , 0.145, 0.145, 0.145, 0.14 , 0.14 , 0.14 , 0.135, 0.135,
         0.13 , 0.13 , 0.125, 0.12 , 0.115, 0.115, 0.11 , 0.105, 0.1  ,
         0.095, 0.085, 0.08 , 0.07 , 0.055, 0.04]
-    profile = Dierckx.Spline1D(z_in, Area_in; k = 1)
+    profile = Intp.interpolate((z_in, ), Area_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -132,7 +132,7 @@ function DryBubble_updrafts_w(::Type{FT}) where {FT}
         0.1601, 0.1559, 0.1589, 0.1543, 0.1496, 0.1514, 0.1464, 0.1475,
         0.1422, 0.1425, 0.1424, 0.1419, 0.1361, 0.135 , 0.1335, 0.1316,
         0.1294, 0.1302, 0.1271, 0.1264, 0.1269, 0.1256]
-    profile = Dierckx.Spline1D(z_in, W_in; k = 1)
+    profile = Intp.interpolate((z_in, ), W_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -153,6 +153,6 @@ function DryBubble_updrafts_T(::Type{FT}) where {FT}
         270.4339, 269.9094, 269.3711, 268.8465, 268.311 , 267.7877,
         267.2649, 266.7432, 266.2159, 265.698 , 265.1821, 264.6685,
         264.1574, 263.6518, 263.1461, 262.6451, 262.1476, 261.6524]
-    profile = Dierckx.Spline1D(z_in, T_in; k = 1)
+    profile = Intp.interpolate((z_in, ), T_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
diff --git a/src/profiles/Dycoms_RF01.jl b/src/profiles/Dycoms_RF01.jl
index bfa2536..4df60f7 100644
--- a/src/profiles/Dycoms_RF01.jl
+++ b/src/profiles/Dycoms_RF01.jl
@@ -30,7 +30,7 @@ function Dycoms_RF01_tke_prescribed(::Type{FT}) where {FT}
               775.0, 825.0, 875.0, 925.0, 975.0, 1025.0, 1075.0, 1125.0, 1175.0, 1225.0, 1275.0, 1325.0, 1375.0, 1425.0, 1475.0]
     tke_in = FT[0.2726, 0.5479, 0.6597, 0.7079, 0.7285, 0.7343, 0.7319, 0.7252, 0.7166, 0.7064, 0.6887, 0.6317,
                 0.6362, 0.6266, 0.5832, 0.4633, 0.0504, 0.0001, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
 
diff --git a/src/profiles/GABLS.jl b/src/profiles/GABLS.jl
index 4829936..7958728 100644
--- a/src/profiles/GABLS.jl
+++ b/src/profiles/GABLS.jl
@@ -26,7 +26,7 @@ GABLS_tke(::Type{FT}) where {FT} =
 function GABLS_tke_prescribed(::Type{FT}) where {FT}
     z_in = FT[25.0, 75.0, 125.0, 175.0, 225.0, 275.0, 325.0, 375.0]
     tke_in = FT[0.4662, 0.3873, 0.2777, 0.0277, 0.0003, 5.89e-8, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
 
diff --git a/src/profiles/GATE_III.jl b/src/profiles/GATE_III.jl
index 853fd4b..3ab028f 100644
--- a/src/profiles/GATE_III.jl
+++ b/src/profiles/GATE_III.jl
@@ -15,7 +15,7 @@ function GATE_III_q_tot(::Type{FT}) where {FT}
                0.9, 0.5, 0.25, 0.125, 0.065, 0.003, 0.0015, 0.0007,  0.0003,  0.0001,  0.0001,  0.0001,  0.0001,
                0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001, 0.0001] ./ 1000 # mixing ratio should be in kg/kg
     q_tot = r_in ./ (1 .+ r_in) # convert mixing ratio to specific humidity
-    profile = Dierckx.Spline1D(z_in, q_tot; k = 1)
+    profile = Intp.interpolate((z_in, ), q_tot, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -26,7 +26,7 @@ function GATE_III_u(::Type{FT}) where {FT}
     U_in = FT[  -1, -1.75, -2.5, -3.6, -6.0, -8.75, -11.75, -13.0, -13.1, -12.1, -11.0, -8.5, -5.0, -2.6, 0.0,
                0.5, 0.4,  0.3,   0.0,  -1.0, -2.5,   -3.5,   -4.5, -4.8, -5.0, -3.5, -2.0, -1.0, -1.0, -1.0,
                -1.5, -2.0, -2.5, -2.6, -2.7, -3.0, -3.0, -3.0] # [m/s]
-    profile = Dierckx.Spline1D(z_in, U_in; k = 1)
+    profile = Intp.interpolate((z_in, ), U_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -36,7 +36,7 @@ function GATE_III_T(::Type{FT}) where {FT}
     z_in = GATE_III_z(FT)
     T_in = FT[299.184, 294.836, 294.261, 288.773, 276.698, 265.004, 253.930, 243.662, 227.674, 214.266, 207.757, 201.973, 198.278, 197.414, 198.110, 198.110]
     z_T_in = FT[0.0, 0.492, 0.700, 1.698, 3.928, 6.039, 7.795, 9.137, 11.055, 12.645, 13.521, 14.486, 15.448, 16.436, 17.293, 22.0] .* 1000 # for km
-    profile = Dierckx.Spline1D(z_T_in, T_in; k = 1)
+    profile = Intp.interpolate((z_T_in, ), T_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -77,7 +77,7 @@ function GATE_III_dqtdt(::Type{FT}) where {FT}
                    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0] ./ (24 * 3600) ./ 1000
 
     Qtend_in = r_tend_in ./ (1 .+ r_tend_in) # convert mixing ratio to specific humidity
-    profile = Dierckx.Spline1D(z_in, Qtend_in; k = 1)
+    profile = Intp.interpolate((z_in, ), Qtend_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(profile)
 end
 
@@ -97,7 +97,7 @@ function GATE_III_dTdt(::Type{FT}) where {FT}
                   -3.7, -3.5, -3.25, -3.0, -2.8,  -2.5, -2.1,  -1.7,  -1.3,   -1.0,   -0.7,   -0.5, -0.4,
                   -0.3, -0.2, -0.1,-0.05,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0] ./ (24*3600)
 
-    profile_T = Dierckx.Spline1D(z_in, Ttend_in; k = 1)
-    profile_R = Dierckx.Spline1D(z_in, RAD_in; k = 1)
+    profile_T = Intp.interpolate((z_in, ), Ttend_in, Intp.Gridded(Intp.Linear()))
+    profile_R = Intp.interpolate((z_in, ), RAD_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(z -> profile_T(z) + profile_R(z))
 end
diff --git a/src/profiles/LifeCycleTan2018.jl b/src/profiles/LifeCycleTan2018.jl
index 71305ae..bf44330 100644
--- a/src/profiles/LifeCycleTan2018.jl
+++ b/src/profiles/LifeCycleTan2018.jl
@@ -57,7 +57,7 @@ function LifeCycleTan2018_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
 
diff --git a/src/profiles/Nieuwstadt.jl b/src/profiles/Nieuwstadt.jl
index 7a74c62..bacde1e 100644
--- a/src/profiles/Nieuwstadt.jl
+++ b/src/profiles/Nieuwstadt.jl
@@ -51,6 +51,6 @@ function Nieuwstadt_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
diff --git a/src/profiles/Rico.jl b/src/profiles/Rico.jl
index 4d472ae..65bf492 100644
--- a/src/profiles/Rico.jl
+++ b/src/profiles/Rico.jl
@@ -78,6 +78,6 @@ function Rico_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
diff --git a/src/profiles/Soares.jl b/src/profiles/Soares.jl
index 276a986..9048b0f 100644
--- a/src/profiles/Soares.jl
+++ b/src/profiles/Soares.jl
@@ -46,6 +46,6 @@ function Soares_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
diff --git a/src/profiles/TRMM_LBA.jl b/src/profiles/TRMM_LBA.jl
index a4872de..d9447b2 100644
--- a/src/profiles/TRMM_LBA.jl
+++ b/src/profiles/TRMM_LBA.jl
@@ -21,7 +21,7 @@ function TRMM_LBA_p(::Type{FT}) where {FT}
               182.3, 167.9, 154.9, 143.0, 131.1, 119.7, 108.9,
               100.1,  92.1,  84.6,  77.5,  71.4,  65.9,  60.7,
                55.9,  51.3,  47.2,  43.3,  10.3] .* 100 # LES pres is in pasc
-    prof = Dierckx.Spline1D(z_in, p_in; k = 1)
+    prof = Intp.interpolate((z_in, ), p_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(prof)
 end
 """ [Grabowski2006](@cite) """
@@ -34,7 +34,7 @@ function TRMM_LBA_T(::Type{FT}) where {FT}
              -59.16, -63.60, -67.68, -70.77, -74.41, -77.51, -80.64,
              -80.69, -80.00, -81.38, -81.17, -78.32, -74.77, -74.52,
              -72.62, -70.87, -69.19, -66.90, -66.90] .+ FT(273.15) # LES T is in deg K
-    prof = Dierckx.Spline1D(z_in, T_in; k = 1)
+    prof = Intp.interpolate((z_in, ), T_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(prof)
 end
 """ [Grabowski2006](@cite) """
@@ -47,7 +47,7 @@ function TRMM_LBA_RH(::Type{FT}) where {FT}
                45.33,  39.78,  33.78,  28.78,  24.67,  20.67,  17.67,
                17.11,  16.22,  14.22,  13.00,  13.00,  12.22,   9.56,
                 7.78,   5.89,   4.33,   3.00,   3.00]
-    prof = Dierckx.Spline1D(z_in, RH_in; k = 1)
+    prof = Intp.interpolate((z_in, ), RH_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(prof)
 end
 """ [Grabowski2006](@cite) """
@@ -60,7 +60,7 @@ function TRMM_LBA_u(::Type{FT}) where {FT}
              -9.00,  -7.77,  -5.37,  -3.88,  -1.15,  -2.36,  -9.20,
              -8.01,  -5.68,  -8.83, -14.51, -15.55, -15.36, -17.67,
             -17.82, -18.94, -15.92, -15.32, -15.32]
-    prof = Dierckx.Spline1D(z_in, u_in; k = 1)
+    prof = Intp.interpolate((z_in, ), u_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(prof)
 end
 """ [Grabowski2006](@cite) """
@@ -73,7 +73,7 @@ function TRMM_LBA_v(::Type{FT}) where {FT}
                3.14,   3.93,   7.57,   2.58,   2.50,   6.44,   6.84,
                0.19,  -2.20,  -3.60,   0.56,   6.68,   9.41,   7.03,
                5.32,   1.14,  -0.65,   5.27,   5.27]
-    prof = Dierckx.Spline1D(z_in, v_in; k = 1)
+    prof = Intp.interpolate((z_in, ), v_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(prof)
 end
 
@@ -110,7 +110,7 @@ function TRMM_LBA_tke_prescribed(::Type{FT}) where {FT}
        1.7938 , 1.56451, 1.37531, 1.17515, 0.96797, 0.61262, 0.26423,
        0.14929, 0.07465, 0.00635, 0.     , 0.     , 0.     , 0.     ,
        0.     , 0.     , 0.     , 0.     , 0.     ]
-    not_type_stable_spline = Dierckx.Spline1D(z_in, tke_in; k = 1)
+    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
     return ZProfile(x -> FT(not_type_stable_spline(x)))
 end
 
@@ -234,6 +234,6 @@ function TRMM_LBA_radiation(::Type{FT}) where {FT}
     rad_in = reduce(hcat, rad_in)::Matrix{FT}
     rad_in .= rad_in ./ 86400
     rad_in = (rad_in')::LinearAlgebra.Adjoint{FT, Matrix{FT}}
-    profile = Dierckx.Spline2D(rad_time, z_in, rad_in; kx = 1, ky = 1)
+    profile = Intp.interpolate((rad_time, z_in), rad_in, (Intp.Gridded(Intp.Linear()), Intp.Gridded(Intp.Linear())))
     return TimeZProfile(profile)
 end

From a800430167145b59a0ccfc8ec81f8c1a8e93bc85 Mon Sep 17 00:00:00 2001
From: Gabriele Bozzola <sbozzolator@gmail.com>
Date: Tue, 2 Jul 2024 19:47:13 -0700
Subject: [PATCH 3/4] Bump to 0.1.6

---
 Project.toml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/Project.toml b/Project.toml
index 186dc73..e46a4bd 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "AtmosphericProfilesLibrary"
 uuid = "86bc3604-9858-485a-bdbe-831ec50de11d"
 authors = ["Climate Modeling Alliance"]
-version = "0.1.5"
+version = "0.1.6"
 
 [deps]
 Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"

From fcd57253979ca43d75f6e980d495fa1117f6ec9f Mon Sep 17 00:00:00 2001
From: Charles Kawczynski <kawczynski.charles@gmail.com>
Date: Wed, 3 Jul 2024 12:59:50 -0400
Subject: [PATCH 4/4] Add and use helper

---
 src/AtmosphericProfilesLibrary.jl |  5 +++++
 src/profiles/ARM_SGP.jl           | 19 ++++++++-----------
 src/profiles/Bomex.jl             |  3 +--
 src/profiles/DryBubble.jl         | 15 +++++----------
 src/profiles/Dycoms_RF01.jl       |  3 +--
 src/profiles/GABLS.jl             |  3 +--
 src/profiles/GATE_III.jl          | 16 ++++++----------
 src/profiles/LifeCycleTan2018.jl  |  3 +--
 src/profiles/Nieuwstadt.jl        |  3 +--
 src/profiles/Rico.jl              |  3 +--
 src/profiles/Soares.jl            |  3 +--
 src/profiles/TRMM_LBA.jl          | 21 ++++++++-------------
 test/runtests.jl                  |  2 +-
 13 files changed, 40 insertions(+), 59 deletions(-)

diff --git a/src/AtmosphericProfilesLibrary.jl b/src/AtmosphericProfilesLibrary.jl
index 61dac70..372bf90 100644
--- a/src/AtmosphericProfilesLibrary.jl
+++ b/src/AtmosphericProfilesLibrary.jl
@@ -28,6 +28,11 @@ struct ΠZProfile{P} <: AbstractProfile
 end
 @inline (prof::ΠZProfile)(Π, z) = prof.prof(Π, z)
 
+@inline function linear_interp(z, x)
+  _interp = Intp.interpolate((z, ), x, Intp.Gridded(Intp.Linear()))
+  return Intp.extrapolate(_interp, Intp.Flat())
+end
+
 # Large data-based profiles
 include("profiles/Soares.jl")
 include("profiles/Nieuwstadt.jl")
diff --git a/src/profiles/ARM_SGP.jl b/src/profiles/ARM_SGP.jl
index 183484e..22841be 100644
--- a/src/profiles/ARM_SGP.jl
+++ b/src/profiles/ARM_SGP.jl
@@ -7,8 +7,7 @@ end
 function ARM_SGP_θ_liq_ice(::Type{FT}) where {FT}
     z_in = ARM_SGP_z(FT)
     θ_liq_ice_in = FT[299.0, 301.5, 302.5, 303.53, 303.7, 307.13, 314.0, 343.2] # K
-    profile = Intp.interpolate((z_in, ), θ_liq_ice_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, θ_liq_ice_in))
 end
 """ [Brown2002](@cite) """
 ARM_SGP_u(::Type{FT}) where {FT} = ZProfile(z -> FT(10))
@@ -18,8 +17,7 @@ function ARM_SGP_q_tot(::Type{FT}) where {FT}
     z_in = ARM_SGP_z(FT)
     r_in = FT[15.2,15.17,14.98,14.8,14.7,13.5,3.0,3.0] ./ 1000 # qt should be in kg/kg
     q_tot_in = r_in ./ (1 .+ r_in)
-    profile = Intp.interpolate((z_in, ), q_tot_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, q_tot_in))
 end
 """ [Brown2002](@cite) """
 function ARM_SGP_tke(::Type{FT}) where {FT}
@@ -49,8 +47,7 @@ function ARM_SGP_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
 
 """ [Brown2002](@cite) """
@@ -66,8 +63,8 @@ function ARM_SGP_dTdt(::Type{FT}) where {FT}
     AT_in = FT[0.0, 0.0, 0.0, -0.08, -0.016, -0.016] ./ 3600
     # Radiative forcing for theta [K/h] converted to [K/sec]
     RT_in = FT[-0.125, 0.0, 0.0, 0.0, 0.0, -0.1] ./ 3600
-    dTdt_A = Intp.interpolate((t_in, ), AT_in, Intp.Gridded(Intp.Linear()))
-    dTdt_R = Intp.interpolate((t_in, ), RT_in, Intp.Gridded(Intp.Linear()))
+    dTdt_A = linear_interp(t_in, AT_in)
+    dTdt_R = linear_interp(t_in, RT_in)
     prof = (t, z) -> if z <= 1000.0
         FT(dTdt_A(t)+dTdt_R(t))
     elseif z > 1000.0 && z <= 2000.0
@@ -83,7 +80,7 @@ function ARM_SGP_dqtdt(::Type{FT}) where {FT}
     t_in = ARM_SGP_time(FT)
     # Radiative forcing for qt converted to [kg/kg/sec]
     Rqt_in = FT[0.08, 0.02, 0.04, -0.1, -0.16, -0.3] ./ 1000 ./ 3600
-    dqtdt = Intp.interpolate((t_in, ), Rqt_in, Intp.Gridded(Intp.Linear()))
+    dqtdt = linear_interp(t_in, Rqt_in)
     prof = (Π, t, z) -> if z <= 1000.0
         FT(dqtdt(t) * Π)
     elseif z > 1000.0 && z <= 2000.0
@@ -98,7 +95,7 @@ end
 function ARM_SGP_shf(::Type{FT}) where {FT}
     t_Sur_in = FT[0.0, 4.0, 6.5, 7.5, 10.0, 12.5, 14.5] .* 3600 #LES time is in sec
     shf = FT[-30.0, 90.0, 140.0, 140.0, 100.0, -10, -10] # W/m^2
-    profile = Intp.interpolate((t_Sur_in, ), shf, Intp.Gridded(Intp.Linear()))
+    profile = linear_interp(t_Sur_in, shf)
     return TimeProfile(profile)
 end
 
@@ -106,6 +103,6 @@ end
 function ARM_SGP_lhf(::Type{FT}) where {FT}
     t_Sur_in = FT[0.0, 4.0, 6.5, 7.5, 10.0, 12.5, 14.5] .* 3600 #LES time is in sec
     lhf = FT[5.0, 250.0, 450.0, 500.0, 420.0, 180.0, 0.0] # W/m^2
-    profile = Intp.interpolate((t_Sur_in, ), lhf, Intp.Gridded(Intp.Linear()))
+    profile = linear_interp(t_Sur_in, lhf)
     return TimeProfile(profile)
 end
diff --git a/src/profiles/Bomex.jl b/src/profiles/Bomex.jl
index cbdf900..eb290ac 100644
--- a/src/profiles/Bomex.jl
+++ b/src/profiles/Bomex.jl
@@ -45,8 +45,7 @@ function Bomex_tke_prescribed(::Type{FT}) where {FT}
                 0.1468, 0.1470, 0.1458, 0.1423, 0.1351, 0.1224, 0.1033, 0.0793, 0.0542, 0.0322,
                 0.0162, 0.0068, 0.0024, 0.0007, 0.0001, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
 
 # Geostrophic velocity profiles
diff --git a/src/profiles/DryBubble.jl b/src/profiles/DryBubble.jl
index 3aaf111..5f62c12 100644
--- a/src/profiles/DryBubble.jl
+++ b/src/profiles/DryBubble.jl
@@ -61,8 +61,7 @@ function DryBubble_θ_liq_ice(::Type{FT}) where {FT}
         299.9837, 299.9837, 299.9837, 299.9837, 299.9837, 299.9837,
         299.9837, 299.9837
     ]
-    profile = Intp.interpolate((z_in, ), θ_liq_ice_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, θ_liq_ice_in))
 end
 
 """ :( """
@@ -97,8 +96,7 @@ function DryBubble_updrafts_θ_liq_ice(::Type{FT}) where {FT}
         300.5212, 300.482 , 300.4272, 300.3875, 300.3354, 300.2968,
         300.2587, 300.2216, 300.1782, 300.1452, 300.1143, 300.0859,
         300.0603, 300.0408, 300.0211, 300.0067, 299.9963, 299.9884]
-    profile = Intp.interpolate((z_in, ), θ_liq_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, θ_liq_in))
 end
 
 """ :( """
@@ -114,8 +112,7 @@ function DryBubble_updrafts_area(::Type{FT}) where {FT}
         0.15 , 0.145, 0.145, 0.145, 0.14 , 0.14 , 0.14 , 0.135, 0.135,
         0.13 , 0.13 , 0.125, 0.12 , 0.115, 0.115, 0.11 , 0.105, 0.1  ,
         0.095, 0.085, 0.08 , 0.07 , 0.055, 0.04]
-    profile = Intp.interpolate((z_in, ), Area_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, Area_in))
 end
 
 """ :( """
@@ -132,8 +129,7 @@ function DryBubble_updrafts_w(::Type{FT}) where {FT}
         0.1601, 0.1559, 0.1589, 0.1543, 0.1496, 0.1514, 0.1464, 0.1475,
         0.1422, 0.1425, 0.1424, 0.1419, 0.1361, 0.135 , 0.1335, 0.1316,
         0.1294, 0.1302, 0.1271, 0.1264, 0.1269, 0.1256]
-    profile = Intp.interpolate((z_in, ), W_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, W_in))
 end
 
 """ :( """
@@ -153,6 +149,5 @@ function DryBubble_updrafts_T(::Type{FT}) where {FT}
         270.4339, 269.9094, 269.3711, 268.8465, 268.311 , 267.7877,
         267.2649, 266.7432, 266.2159, 265.698 , 265.1821, 264.6685,
         264.1574, 263.6518, 263.1461, 262.6451, 262.1476, 261.6524]
-    profile = Intp.interpolate((z_in, ), T_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, T_in))
 end
diff --git a/src/profiles/Dycoms_RF01.jl b/src/profiles/Dycoms_RF01.jl
index 4df60f7..9425194 100644
--- a/src/profiles/Dycoms_RF01.jl
+++ b/src/profiles/Dycoms_RF01.jl
@@ -30,7 +30,6 @@ function Dycoms_RF01_tke_prescribed(::Type{FT}) where {FT}
               775.0, 825.0, 875.0, 925.0, 975.0, 1025.0, 1075.0, 1125.0, 1175.0, 1225.0, 1275.0, 1325.0, 1375.0, 1425.0, 1475.0]
     tke_in = FT[0.2726, 0.5479, 0.6597, 0.7079, 0.7285, 0.7343, 0.7319, 0.7252, 0.7166, 0.7064, 0.6887, 0.6317,
                 0.6362, 0.6266, 0.5832, 0.4633, 0.0504, 0.0001, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
 
diff --git a/src/profiles/GABLS.jl b/src/profiles/GABLS.jl
index 7958728..ce95ac2 100644
--- a/src/profiles/GABLS.jl
+++ b/src/profiles/GABLS.jl
@@ -26,8 +26,7 @@ GABLS_tke(::Type{FT}) where {FT} =
 function GABLS_tke_prescribed(::Type{FT}) where {FT}
     z_in = FT[25.0, 75.0, 125.0, 175.0, 225.0, 275.0, 325.0, 375.0]
     tke_in = FT[0.4662, 0.3873, 0.2777, 0.0277, 0.0003, 5.89e-8, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
 
 """ :( """
diff --git a/src/profiles/GATE_III.jl b/src/profiles/GATE_III.jl
index 3ab028f..3b31d0b 100644
--- a/src/profiles/GATE_III.jl
+++ b/src/profiles/GATE_III.jl
@@ -15,8 +15,7 @@ function GATE_III_q_tot(::Type{FT}) where {FT}
                0.9, 0.5, 0.25, 0.125, 0.065, 0.003, 0.0015, 0.0007,  0.0003,  0.0001,  0.0001,  0.0001,  0.0001,
                0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001,  0.0001, 0.0001] ./ 1000 # mixing ratio should be in kg/kg
     q_tot = r_in ./ (1 .+ r_in) # convert mixing ratio to specific humidity
-    profile = Intp.interpolate((z_in, ), q_tot, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, q_tot))
 end
 
 
@@ -26,8 +25,7 @@ function GATE_III_u(::Type{FT}) where {FT}
     U_in = FT[  -1, -1.75, -2.5, -3.6, -6.0, -8.75, -11.75, -13.0, -13.1, -12.1, -11.0, -8.5, -5.0, -2.6, 0.0,
                0.5, 0.4,  0.3,   0.0,  -1.0, -2.5,   -3.5,   -4.5, -4.8, -5.0, -3.5, -2.0, -1.0, -1.0, -1.0,
                -1.5, -2.0, -2.5, -2.6, -2.7, -3.0, -3.0, -3.0] # [m/s]
-    profile = Intp.interpolate((z_in, ), U_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, U_in))
 end
 
 """ [Khairoutdinov2009](@cite) """
@@ -36,8 +34,7 @@ function GATE_III_T(::Type{FT}) where {FT}
     z_in = GATE_III_z(FT)
     T_in = FT[299.184, 294.836, 294.261, 288.773, 276.698, 265.004, 253.930, 243.662, 227.674, 214.266, 207.757, 201.973, 198.278, 197.414, 198.110, 198.110]
     z_T_in = FT[0.0, 0.492, 0.700, 1.698, 3.928, 6.039, 7.795, 9.137, 11.055, 12.645, 13.521, 14.486, 15.448, 16.436, 17.293, 22.0] .* 1000 # for km
-    profile = Intp.interpolate((z_T_in, ), T_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_T_in, T_in))
 end
 
 """ [Khairoutdinov2009](@cite) """
@@ -77,8 +74,7 @@ function GATE_III_dqtdt(::Type{FT}) where {FT}
                    0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0] ./ (24 * 3600) ./ 1000
 
     Qtend_in = r_tend_in ./ (1 .+ r_tend_in) # convert mixing ratio to specific humidity
-    profile = Intp.interpolate((z_in, ), Qtend_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(profile)
+    return ZProfile(linear_interp(z_in, Qtend_in))
 end
 
 """ [Khairoutdinov2009](@cite) """
@@ -97,7 +93,7 @@ function GATE_III_dTdt(::Type{FT}) where {FT}
                   -3.7, -3.5, -3.25, -3.0, -2.8,  -2.5, -2.1,  -1.7,  -1.3,   -1.0,   -0.7,   -0.5, -0.4,
                   -0.3, -0.2, -0.1,-0.05,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0,  0.0] ./ (24*3600)
 
-    profile_T = Intp.interpolate((z_in, ), Ttend_in, Intp.Gridded(Intp.Linear()))
-    profile_R = Intp.interpolate((z_in, ), RAD_in, Intp.Gridded(Intp.Linear()))
+    profile_T = linear_interp(z_in, Ttend_in)
+    profile_R = linear_interp(z_in, RAD_in)
     return ZProfile(z -> profile_T(z) + profile_R(z))
 end
diff --git a/src/profiles/LifeCycleTan2018.jl b/src/profiles/LifeCycleTan2018.jl
index bf44330..e3d9b89 100644
--- a/src/profiles/LifeCycleTan2018.jl
+++ b/src/profiles/LifeCycleTan2018.jl
@@ -57,8 +57,7 @@ function LifeCycleTan2018_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
 
 # Large-scale cooling
diff --git a/src/profiles/Nieuwstadt.jl b/src/profiles/Nieuwstadt.jl
index bacde1e..41dfdf0 100644
--- a/src/profiles/Nieuwstadt.jl
+++ b/src/profiles/Nieuwstadt.jl
@@ -51,6 +51,5 @@ function Nieuwstadt_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
diff --git a/src/profiles/Rico.jl b/src/profiles/Rico.jl
index 65bf492..d84a186 100644
--- a/src/profiles/Rico.jl
+++ b/src/profiles/Rico.jl
@@ -78,6 +78,5 @@ function Rico_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
diff --git a/src/profiles/Soares.jl b/src/profiles/Soares.jl
index 9048b0f..486bb53 100644
--- a/src/profiles/Soares.jl
+++ b/src/profiles/Soares.jl
@@ -46,6 +46,5 @@ function Soares_tke_prescribed(::Type{FT}) where {FT}
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
                 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
diff --git a/src/profiles/TRMM_LBA.jl b/src/profiles/TRMM_LBA.jl
index d9447b2..67a6003 100644
--- a/src/profiles/TRMM_LBA.jl
+++ b/src/profiles/TRMM_LBA.jl
@@ -21,8 +21,7 @@ function TRMM_LBA_p(::Type{FT}) where {FT}
               182.3, 167.9, 154.9, 143.0, 131.1, 119.7, 108.9,
               100.1,  92.1,  84.6,  77.5,  71.4,  65.9,  60.7,
                55.9,  51.3,  47.2,  43.3,  10.3] .* 100 # LES pres is in pasc
-    prof = Intp.interpolate((z_in, ), p_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(prof)
+    return ZProfile(linear_interp(z_in, p_in))
 end
 """ [Grabowski2006](@cite) """
 function TRMM_LBA_T(::Type{FT}) where {FT}
@@ -34,8 +33,7 @@ function TRMM_LBA_T(::Type{FT}) where {FT}
              -59.16, -63.60, -67.68, -70.77, -74.41, -77.51, -80.64,
              -80.69, -80.00, -81.38, -81.17, -78.32, -74.77, -74.52,
              -72.62, -70.87, -69.19, -66.90, -66.90] .+ FT(273.15) # LES T is in deg K
-    prof = Intp.interpolate((z_in, ), T_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(prof)
+    return ZProfile(linear_interp(z_in, T_in))
 end
 """ [Grabowski2006](@cite) """
 function TRMM_LBA_RH(::Type{FT}) where {FT}
@@ -47,8 +45,7 @@ function TRMM_LBA_RH(::Type{FT}) where {FT}
                45.33,  39.78,  33.78,  28.78,  24.67,  20.67,  17.67,
                17.11,  16.22,  14.22,  13.00,  13.00,  12.22,   9.56,
                 7.78,   5.89,   4.33,   3.00,   3.00]
-    prof = Intp.interpolate((z_in, ), RH_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(prof)
+    return ZProfile(linear_interp(z_in, RH_in))
 end
 """ [Grabowski2006](@cite) """
 function TRMM_LBA_u(::Type{FT}) where {FT}
@@ -60,8 +57,7 @@ function TRMM_LBA_u(::Type{FT}) where {FT}
              -9.00,  -7.77,  -5.37,  -3.88,  -1.15,  -2.36,  -9.20,
              -8.01,  -5.68,  -8.83, -14.51, -15.55, -15.36, -17.67,
             -17.82, -18.94, -15.92, -15.32, -15.32]
-    prof = Intp.interpolate((z_in, ), u_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(prof)
+    return ZProfile(linear_interp(z_in, u_in))
 end
 """ [Grabowski2006](@cite) """
 function TRMM_LBA_v(::Type{FT}) where {FT}
@@ -73,8 +69,7 @@ function TRMM_LBA_v(::Type{FT}) where {FT}
                3.14,   3.93,   7.57,   2.58,   2.50,   6.44,   6.84,
                0.19,  -2.20,  -3.60,   0.56,   6.68,   9.41,   7.03,
                5.32,   1.14,  -0.65,   5.27,   5.27]
-    prof = Intp.interpolate((z_in, ), v_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(prof)
+    return ZProfile(linear_interp(z_in, v_in))
 end
 
 """ [Grabowski2006](@cite) """
@@ -110,8 +105,7 @@ function TRMM_LBA_tke_prescribed(::Type{FT}) where {FT}
        1.7938 , 1.56451, 1.37531, 1.17515, 0.96797, 0.61262, 0.26423,
        0.14929, 0.07465, 0.00635, 0.     , 0.     , 0.     , 0.     ,
        0.     , 0.     , 0.     , 0.     , 0.     ]
-    not_type_stable_spline = Intp.interpolate((z_in, ), tke_in, Intp.Gridded(Intp.Linear()))
-    return ZProfile(x -> FT(not_type_stable_spline(x)))
+    return ZProfile(linear_interp(z_in, tke_in))
 end
 
 
@@ -234,6 +228,7 @@ function TRMM_LBA_radiation(::Type{FT}) where {FT}
     rad_in = reduce(hcat, rad_in)::Matrix{FT}
     rad_in .= rad_in ./ 86400
     rad_in = (rad_in')::LinearAlgebra.Adjoint{FT, Matrix{FT}}
-    profile = Intp.interpolate((rad_time, z_in), rad_in, (Intp.Gridded(Intp.Linear()), Intp.Gridded(Intp.Linear())))
+    interp = Intp.interpolate((rad_time, z_in), rad_in, (Intp.Gridded(Intp.Linear()), Intp.Gridded(Intp.Linear())))
+    profile = Intp.extrapolate(interp, Intp.Flat())
     return TimeZProfile(profile)
 end
diff --git a/test/runtests.jl b/test/runtests.jl
index 699ef21..78a39d7 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -5,7 +5,7 @@ import AtmosphericProfilesLibrary as APL
 iscallable(f) = !isempty(methods(f))
 
 function test_profile(::Type{FT}, prof, name) where {FT}
-    @debug "Testing $name"
+    # println("Testing $name")
     if prof isa APL.TimeProfile
         prof(1)
         @inferred prof(Float32(1))