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. В большинстве случаев эта функция может работать хорошо.
Интересно, как решить эту проблему. Большое спасибо!
1 ответ
Решение
Это ошибка в MetPy. Можете ли вы открыть проблему с помощью MetPy Issue Tracker?