Add and use helper

CliMA · Jul 3, 2024 · fcd5725 · fcd5725
1 parent a800430
commit fcd5725
Show file tree

Hide file tree

Showing 13 changed files with 40 additions and 59 deletions.
diff --git 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
@@ -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,14 +95,14 @@ 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
 
 """ [Brown2002](@cite) """
 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
@@ -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
@@ -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
@@ -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
@@ -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
@@ -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
@@ -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
@@ -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
@@ -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
@@ -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
@@ -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
@@ -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))