In [1]:
from scipy.interpolate import interp2d

def TRS398_table7():
    table = dict()
    table['R50'] = [
        1.0, 1.4, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5,
        6.0, 7.0, 8.0, 10.0, 13.0, 16.0, 20.0
    ]
    table['depth/R50'] = [
        0.02, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35,
        0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75,
        0.80, 0.85, 0.90, 0.95, 1.00, 1.05, 1.10, 1.15,
        1.20
    ]
    table['contents'] = [
        [1.076, 1.060, 1.042, 1.030, 1.020, 1.012, 1.004, 0.997, 0.991, 0.986,
         0.980, 0.971, 0.963, 0.950, 0.935, 0.924, 0.914],
        [1.078, 1.061, 1.044, 1.032, 1.022, 1.014, 1.006, 1.000, 0.994, 0.988,
         0.983, 0.974, 0.965, 0.952, 0.937, 0.926, 0.916],
        [1.080, 1.064, 1.047, 1.036, 1.026, 1.018, 1.010, 1.004, 0.998, 0.992,
         0.987, 0.978, 0.970, 0.957, 0.942, 0.931, 0.920],
        [1.083, 1.067, 1.050, 1.039, 1.030, 1.022, 1.014, 1.008, 1.002, 0.997,
         0.992, 0.983, 0.975, 0.961, 0.946, 0.935, 0.924],
        [1.085, 1.070, 1.053, 1.043, 1.034, 1.026, 1.019, 1.012, 1.006, 1.001,
         0.996, 0.987, 0.979, 0.966, 0.951, 0.940, 0.929],
        [1.088, 1.073, 1.057, 1.046, 1.037, 1.030, 1.023, 1.017, 1.011, 1.006,
         1.001, 0.992, 0.984, 0.971, 0.956, 0.945, 0.933],
        [1.091, 1.076, 1.060, 1.050, 1.041, 1.034, 1.027, 1.021, 1.016, 1.010,
         1.006, 0.997, 0.989, 0.976, 0.961, 0.950, 0.938],
        [1.093, 1.079, 1.064, 1.054, 1.045, 1.038, 1.032, 1.026, 1.020, 1.015,
         1.011, 1.002, 0.995, 0.982, 0.966, 0.955, 0.943],
        [1.096, 1.082, 1.067, 1.058, 1.049, 1.042, 1.036, 1.030, 1.025, 1.020,
         1.016, 1.007, 1.000, 0.987, 0.972, 0.960, 0.948],
        [1.099, 1.085, 1.071, 1.062, 1.054, 1.047, 1.041, 1.035, 1.030, 1.025,
         1.021, 1.013, 1.006, 0.993, 0.978, 0.966, 0.953],
        [1.102, 1.089, 1.075, 1.066, 1.058, 1.051, 1.046, 1.040, 1.035, 1.031,
         1.027, 1.019, 1.012, 0.999, 0.984, 0.971, 0.959],
        [1.105, 1.092, 1.078, 1.070, 1.062, 1.056, 1.051, 1.045, 1.041, 1.036,
         1.032, 1.025, 1.018, 1.005, 0.990, 0.977, 0.964],
        [1.108, 1.095, 1.082, 1.074, 1.067, 1.061, 1.056, 1.051, 1.046, 1.042,
         1.038, 1.031, 1.024, 1.012, 0.996, 0.984, 0.970],
        [1.111, 1.099, 1.086, 1.078, 1.072, 1.066, 1.061, 1.056, 1.052, 1.048,
         1.044, 1.037, 1.030, 1.018, 1.003, 0.990, 0.976],
        [1.114, 1.102, 1.090, 1.082, 1.076, 1.071, 1.066, 1.062, 1.058, 1.054,
         1.050, 1.043, 1.037, 1.025, 1.010, 0.997, 0.983],
        [1.117, 1.105, 1.094, 1.087, 1.081, 1.076, 1.072, 1.067, 1.064, 1.060,
         1.057, 1.050, 1.044, 1.033, 1.017, 1.004, 0.989],
        [1.120, 1.109, 1.098, 1.091, 1.086, 1.081, 1.077, 1.073, 1.070, 1.066,
         1.063, 1.057, 1.051, 1.040, 1.025, 1.012, 0.996],
        [1.123, 1.112, 1.102, 1.096, 1.091, 1.087, 1.083, 1.080, 1.076, 1.073,
         1.070, 1.064, 1.059, 1.048, 1.033, 1.019, 1.004],
        [1.126, 1.116, 1.107, 1.101, 1.096, 1.092, 1.089, 1.086, 1.083, 1.080,
         1.077, 1.072, 1.067, 1.056, 1.041, 1.028, 1.011],
        [1.129, 1.120, 1.111, 1.106, 1.102, 1.098, 1.095, 1.092, 1.090, 1.087,
         1.085, 1.080, 1.075, 1.065, 1.050, 1.036, 1.019],
        [1.132, 1.124, 1.115, 1.111, 1.107, 1.104, 1.101, 1.099, 1.097, 1.095,
         1.092, 1.088, 1.083, 1.074, 1.059, 1.045, 1.028],
        [1.136, 1.127, 1.120, 1.116, 1.113, 1.110, 1.108, 1.106, 1.104, 1.102,
         1.100, 1.096, 1.092, 1.083, 1.069, 1.055, 1.037],
        [1.139, 1.131, 1.125, 1.121, 1.118, 1.116, 1.115, 1.113, 1.112, 1.110,
         1.109, 1.105, 1.102, 1.093, 1.079, 1.065, 1.046],
        [1.142, 1.135, 1.129, 1.126, 1.124, 1.123, 1.122, 1.120, 1.119, 1.118,
         1.117, 1.114, 1.111, 1.104, 1.090, 1.075, 1.056],
        [1.146, 1.139, 1.134, 1.132, 1.130, 1.129, 1.129, 1.128, 1.128, 1.127,
         1.126, 1.124, 1.121, 1.115, 1.101, 1.086, 1.066]
    ]
    
    return table

In [2]:
data = TRS398_table7()

In [3]:
stop_power_interp = interp2d(
    data['R50'], data['depth/R50'], data['contents'], kind='cubic'
)

In [4]:
stop_power_interp(2.0, 0.10)


Out[4]:
array([ 1.047])

In [ ]: