leifdenby / tephigram_python

Python package for plotting tephigrams
BSD 3-Clause "New" or "Revised" License
9 stars 5 forks source link

Plotting U of Wyoming sounding data #5

Open winash12 opened 5 years ago

winash12 commented 5 years ago

I have data from University of Wyoming soundings - [http://weather.uwyo.edu/upperair/sounding.html]. Can I use your package to plot this data and create a tephigram ?

leifdenby commented 5 years ago

Yes, that was exactly the data I used when developing this code :) Although when you say "Plotting U" do you mean plotting velocity? I haven't added any wind profile plots yet unfortunately

winash12 commented 5 years ago

I am getting this error

Traceback (most recent call last): File "test_plot.py", line 70, in test_plot_test_parcel() File "test_plot.py", line 62, in test_plot_test_parcel parcel_info = tephigram.plot_test_parcel(z=z, P=P, T=T, RH=RH) File "/usr/local/lib/python3.6/dist-packages/tephigram-0.2.0-py3.6.egg/tephigram/plotter.py", line 396, in plot_test_parcel F_parcel, parcel_info = parcel_integrator.find_EL(dqv0=dqv0, dT0=dT0, raise_exception=False, include_info=True, z0=z[0]) File "/usr/local/lib/python3.6/dist-packages/tephigram-0.2.0-py3.6.egg/tephigram/profile_integration.py", line 457, in find_EL dp=dp) File "/usr/local/lib/python3.6/dist-packages/tephigram-0.2.0-py3.6.egg/tephigram/profile_integration.py", line 223, in _integrate while not stopping_criterion(F_current=F[-1]): File "/usr/local/lib/python3.6/dist-packages/tephigram-0.2.0-py3.6.egg/tephigram/profile_integration.py", line 420, in call p=F_current[Var.p], T=self.thermodynamic_parcel_integration.profile.temp(z=z), File "/usr/local/lib/python3.6/dist-packages/tephigram-0.2.0-py3.6.egg/tephigram/profile_integration.py", line 105, in return lambda z: scipy.interpolate.interp1d(x=self.z, y=y_discrete)(z) File "/usr/local/lib/python3.6/dist-packages/scipy-1.1.0-py3.6-linux-x86_64.egg/scipy/interpolate/polyint.py", line 79, in call y = self._evaluate(x) File "/usr/local/lib/python3.6/dist-packages/scipy-1.1.0-py3.6-linux-x86_64.egg/scipy/interpolate/interpolate.py", line 664, in _evaluate below_bounds, above_bounds = self._check_bounds(x_new) File "/usr/local/lib/python3.6/dist-packages/scipy-1.1.0-py3.6-linux-x86_64.egg/scipy/interpolate/interpolate.py", line 696, in _check_bounds raise ValueError("A value in x_new is above the interpolation " ValueError: A value in x_new is above the interpolation range.

for this data file

10410 EDZE Essen Observations at 00Z 20 Mar 2019

-----------------------------------------------------------------------------

PRES HGHT TEMP DWPT RELH MIXR DRCT SKNT THTA THTE THTV

hPa m C C % g/kg deg knot K K K

-----------------------------------------------------------------------------

1014.0 147 5.2 0.6 72 3.96 170 2 277.2 288.3 277.9 1013.0 155 5.0 0.6 73 3.96 172 2 277.1 288.2 277.8 1003.0 233 7.6 -0.4 57 3.72 194 8 280.5 291.1 281.1 1000.0 257 7.4 0.4 61 3.96 200 10 280.6 291.8 281.2 959.0 600 4.5 -0.2 71 3.94 245 10 281.0 292.2 281.7 925.0 895 2.0 -0.8 82 3.92 250 8 281.4 292.5 282.0 890.0 1204 -0.7 -1.8 92 3.78 255 9 281.7 292.5 282.3 858.0 1496 -2.5 -6.3 75 2.79 260 10 282.8 290.9 283.3 850.0 1570 -2.9 -7.5 71 2.57 270 12 283.1 290.6 283.5 833.0 1730 -3.9 -9.9 63 2.17 278 11 283.7 290.1 284.1 817.0 1882 -5.2 -10.2 68 2.16 285 10 283.9 290.3 284.2 797.0 2077 -6.9 -10.6 75 2.15 308 10 284.1 290.5 284.4 781.0 2235 -6.1 -24.1 23 0.70 327 11 286.6 288.9 286.7 745.0 2604 -5.7 -27.5 16 0.54 10 12 290.9 292.7 291.0 743.0 2626 -5.7 -27.7 16 0.53 9 12 291.1 292.9 291.2 710.0 2980 -7.1 -24.5 24 0.75 350 14 293.4 295.9 293.6 700.0 3091 -7.5 -23.5 27 0.83 355 16 294.1 296.9 294.3 677.0 3347 -8.6 -26.1 23 0.67 10 19 295.8 298.0 295.9 565.0 4736 -14.5 -40.5 9 0.20 2 22 304.5 305.2 304.5 553.0 4898 -15.7 -32.7 22 0.44 1 23 304.9 306.5 305.0 550.0 4939 -16.1 -24.1 50 1.00 1 23 304.9 308.3 305.1 546.0 4994 -16.1 -44.1 7 0.14 1 23 305.6 306.1 305.6 539.0 5091 -16.9 -21.8 66 1.25 0 23 305.7 309.9 306.0 536.0 5132 -17.3 -21.6 69 1.28 0 23 305.8 310.0 306.0 532.0 5189 -17.7 -24.7 54 0.98 1 24 305.9 309.2 306.1 530.0 5217 -18.1 -20.8 79 1.38 1 25 305.8 310.4 306.0 500.0 5650 -20.5 -24.4 71 1.07 5 31 308.0 311.6 308.2 471.0 6086 -23.5 -27.9 67 0.82 15 37 309.6 312.5 309.8 443.0 6533 -26.5 -31.5 63 0.62 13 40 311.2 313.5 311.4 400.0 7260 -32.3 -36.9 64 0.41 10 45 312.9 314.4 313.0 392.0 7402 -33.5 -36.9 71 0.41 10 47 313.2 314.7 313.2 367.0 7861 -37.9 -42.5 62 0.25 12 55 313.3 314.2 313.3 363.0 7937 -38.1 -48.1 34 0.14 12 56 314.0 314.5 314.0 350.0 8187 -39.9 -51.9 27 0.09 13 60 314.8 315.2 314.9 312.0 8957 -46.8 -55.8 35 0.06 15 74 315.8 316.0 315.8 300.0 9220 -49.1 -57.1 39 0.06 10 74 316.0 316.3 316.1 292.0 9397 -50.7 -57.7 43 0.05 11 74 316.2 316.4 316.2 286.0 9533 -50.3 -67.3 12 0.02 11 74 318.7 318.7 318.7 260.0 10150 -55.7 -62.7 41 0.03 14 74 319.5 319.7 319.5 255.0 10274 -56.7 -64.7 36 0.03 15 74 319.8 319.9 319.9 250.0 10400 -57.7 -66.7 31 0.02 10 70 320.2 320.2 320.2 247.0 10476 -57.9 -68.9 23 0.01 10 69 321.0 321.0 321.0 237.0 10735 -56.9 -69.9 18 0.01 12 66 326.3 326.4 326.3 220.0 11198 -60.2 -70.5 25 0.01 15 60 328.1 328.2 328.1 213.0 11400 -61.7 -70.7 29 0.01 5 56 328.9 329.0 328.9 206.0 11606 -60.9 -73.3 18 0.01 355 62 333.4 333.4 333.4 203.0 11696 -60.5 -74.5 14 0.01 360 55 335.4 335.4 335.4 201.0 11758 -55.3 -77.3 5 0.01 3 51 344.5 344.6 344.5 200.0 11790 -54.7 -79.7 3 0.00 5 49 346.0 346.0 346.0 198.0 11855 -54.3 -80.8 2 0.00 10 45 347.6 347.6 347.6 196.0 11920 -53.9 -81.9 2 0.00 8 43 349.3 349.3 349.3 188.0 12191 -53.3 -82.3 2 0.00 0 37 354.4 354.4 354.4 167.0 12961 -51.7 -83.4 1 0.00 355 45 369.3 369.3 369.3 159.0 13280 -51.0 -83.8 1 0.00 350 41 375.7 375.7 375.7 158.0 13321 -50.9 -83.9 1 0.00 352 40 376.5 376.5 376.5 150.0 13660 -51.9 -83.9 1 0.00 5 33 380.4 380.4 380.4 143.0 13963 -52.9 -84.4 1 0.00 340 27 383.9 384.0 383.9 133.0 14422 -54.4 -85.2 1 0.00 345 43 389.3 389.3 389.3 118.0 15181 -56.9 -86.5 1 0.00 355 27 398.3 398.3 398.3 100.0 16230 -60.3 -88.3 1 0.00 15 35 410.9 411.0 410.9

Station information and sounding indices

Station identifier: EDZE

Station number: 10410

Observation time: 190320/0000

Station latitude: 51.40

Station longitude: 6.97

Station elevation: 147.0

Showalter index: 15.24

Lifted index: 14.61

LIFT computed using virtual temperature: 14.72

SWEAT index: 54.39

K index: -5.90

Cross totals index: 13.00

Vertical totals index: 17.60

Totals totals index: 30.60

Convective Available Potential Energy: 0.00

CAPE using virtual temperature: 0.00

Convective Inhibition: 0.00

CINS using virtual temperature: 0.00

Bulk Richardson Number: 0.00

Bulk Richardson Number using CAPV: 0.00

Temp [K] of the Lifted Condensation Level: 271.96

Pres [hPa] of the Lifted Condensation Level: 898.74

Mean mixed layer potential temperature: 280.40

Mean mixed layer mixing ratio: 3.92

1000 hPa to 500 hPa thickness: 5393.00

Precipitable water [mm] for entire sounding: 10.40

winash12 commented 5 years ago

Any update on this?