Python вызывает metpy.calc.lfc (), поднимает IndexError: индекс 0 выходит за пределы оси 0 с размером 0

Моя программа Python повысила IndexError: index 0 is out of bounds for axis 0 with size 0 когда я вызываю metpy.calc.lft().

Мой исходный код:

lfc_p, lfc_T = mpcalc.lfc(P,T,Td)

Полная информация:

Traceback (most recent call last):   
   File "cal_sounding_data.py", line 289, in <module>
     newset_dic[k] = specStationSounding(k,dataset_dic)
   File "cal_sounding_data.py", line 213, in specStationSounding
     result_dic = soundingCalculation(num,p,t,td,ws,wd)   
   File "cal_sounding_data.py", line 142, in soundingCalculation
     lfc_p, lfc_T = mpcalc.lfc(P,T,Td)   
   File "/anaconda2/lib/python2.7/site-packages/metpy/xarray.py", line 138, in wrapper
     return func(*args, **kwargs)   
   File "/anaconda2/lib/python2.7/site-packages/metpy/calc/thermo.py", line 378, in lfc
     return x[0], y[0]   
   File "/anaconda2/lib/python2.7/site-packages/pint/quantity.py", line 1281, in __getitem__
     value = self._magnitude[key] 
IndexError: index 0 is out of bounds for axis 0 with size 0

Я отлаживал его с помощью pdb и вставлял следующий код, чтобы поймать ошибку:

try:
    lfc_p, lfc_T = mpcalc.lfc(P,T,Td)
except IndexError:
    pdb.set_trace()

Я запустил эту программу, а затем проверил переменные P, T и Td, когда она столкнулась с этой ошибкой.

Это показывает, что:

(Pdb) P
<Quantity([ 1004.  1000.   943.   928.   925.   850.   839.   749.   700.   699.
   603.   500.   404.   400.   363.   306.   300.   250.   213.   200.
   176.   150.   120.   105.   100.    70.    66.    58.    50.    40.
    33.    30.    23.    20.    16.], 'hectopascal')>
(Pdb) T
<Quantity([ 24.2  24.   20.2  21.6  21.4  20.4  20.2  14.4  13.2  13.    6.8  -3.3
 -13.1 -13.7 -17.9 -25.5 -26.9 -37.9 -46.7 -48.7 -52.1 -58.9 -67.3 -66.5
 -66.7 -65.1 -66.1 -60.9 -60.5 -57.7 -50.1 -50.3 -50.1 -47.9 -43.1], 'degC')>
(Pdb) Td
<Quantity([  2.19000000e+01   2.21000000e+01   1.92000000e+01   2.05000000e+01
   2.04000000e+01   1.84000000e+01   1.74000000e+01   8.40000000e+00
  -2.80000000e+00  -3.00000000e+00  -1.52000000e+01  -2.03000000e+01
  -2.91000000e+01  -2.77000000e+01  -2.49000000e+01  -3.95000000e+01
  -4.19000000e+01  -5.19000000e+01  -6.07000000e+01  -6.27000000e+01
  -6.51000000e+01  -7.19000000e+01   9.99999000e+05   9.99999000e+05
   9.99999000e+05   9.99999000e+05   9.99999000e+05   9.99999000e+05
   9.99999000e+05   9.99999000e+05   9.99999000e+05   9.99999000e+05
   9.99999000e+05   9.99999000e+05   9.99999000e+05], 'degC')>

Таким образом, кажется, что P, T и Td верны. В конце концов, они на самом деле не оси с размером 0. Большую часть времени эта функция может работать хорошо.

Интересно, как решить эту проблему. Большое спасибо!