尝试绘制T-LnP图

MeteoInfo

尝试绘制T-LnP图，从micaps第5类数据中读取了某个站的探空数据（不同高度的温度、露点、风向、风速）。

目前绘制的图形，温度、露点、风向杆廓线、干绝热线：


脚本程序：

1. #Read data from MICAPS 5 data file
   
2. fn = 'D:/Temp/micaps/16051108.000'
   
3. f = open(fn)
   
4. for i in range(2):
   
5.     f.readline()
   
6. line = f.readline().strip()
   
7. n = int(line.split()[4]) / 6
   
8. pres = []    #Pressure
   
9. height = []  #Height
   
10. temp = []    #Temperature
    
11. temp_d = []  #Dew point
    
12. wdir = []    #Wind direction
    
13. wspeed = []  #Wind speed
    
14. for i in range(n):
    
15.     line = f.readline().strip()
    
16.     data = line.split()
    
17.     for j in range(len(data)):
    
18.         if data[j] == '9999':
    
19.             data[j] = nan
    
20.         else:
    
21.             data[j] = float(data[j])   
    
22.     pres.append(data[0])
    
23.     height.append(data[1])
    
24.     temp.append(data[2])
    
25.     temp_d.append(data[3])
    
26.     wdir.append(data[4])
    
27.     wspeed.append(data[5])
    
28. pres = array(pres)
    
29. y = log(1000./pres)
    
30. 
    
31. #Plot
    
32. #Plot dry lapse lines
    
33. for t0 in range(-90, 200, 10):
    
34.     t = meteo.dry_lapse(pres, t0 + 273.15) - 273.15
    
35.     plot(t, y, color=[255,204,0])
    
36. #Plot temperature line
    
37. line1 = plot(temp, y, 'r', linewidth=2, label='Temperature')
    
38. #Plot dew point line
    
39. line2 = plot(temp_d, y, 'g', linewidth=2, label='Dew Point')
    
40. #Plot wind barbs
    
41. x = zeros(len(y))
    
42. x = x + 10
    
43. barbs(x, y, wdir, wspeed, isuv=False, colors='b')
    
44. yticks(y[1:], pres[1:].astype('int'))
    
45. ylabel('Pressure (hPa)')
    
46. xlabel(r'$\rm{Temperature} \ (^{\circ} \ C)$')
    
47. legend([line1, line2])
    
48. grid()
    
49. xlim(-90, 20)
    
50. ylim(y[1], y[len(y)-1])
    
51. title('T-LnP figure example')

#Read data from MICAPS 5 data file<br /> fn = 'D:/Temp/micaps/16051108.000'<br /> f = open(fn)<br /> for i in range(2):<br />     f.readline()<br /> line = f.readline().strip()<br /> n = int(line.split()[4]) / 6<br /> pres = []    #Pressure<br /> height = []  #Height<br /> temp = []    #Temperature<br /> temp_d = []  #Dew point<br /> wdir = []    #Wind direction<br /> wspeed = []  #Wind speed<br /> for i in range(n):<br />     line = f.readline().strip()<br />     data = line.split()<br />     for j in range(len(data)):<br />         if data[j] == '9999':<br />             data[j] = nan<br />         else:<br />             data[j] = float(data[j])    <br />     pres.append(data[0])<br />     height.append(data[1])<br />     temp.append(data[2])<br />     temp_d.append(data[3])<br />     wdir.append(data[4])<br />     wspeed.append(data[5])<br /> pres = array(pres)<br /> y = log(1000./pres)<br /> <br /> #Plot<br /> #Plot dry lapse lines<br /> for t0 in range(-90, 200, 10):<br />     t = meteo.dry_lapse(pres, t0 + 273.15) - 273.15<br />     plot(t, y, color=[255,204,0])<br /> #Plot temperature line<br /> line1 = plot(temp, y, 'r', linewidth=2, label='Temperature')<br /> #Plot dew point line<br /> line2 = plot(temp_d, y, 'g', linewidth=2, label='Dew Point')<br /> #Plot wind barbs<br /> x = zeros(len(y))<br /> x = x + 10<br /> barbs(x, y, wdir, wspeed, isuv=False, colors='b')<br /> yticks(y[1:], pres[1:].astype('int'))<br /> ylabel('Pressure (hPa)')<br /> xlabel(r'$\rm{Temperature} \ (^{\circ} \ C)$')<br /> legend([line1, line2])<br /> grid()<br /> xlim(-90, 20)<br /> ylim(y[1], y[len(y)-1])<br /> title('T-LnP figure example')

Iron_Lee

楼主barbs函数是自己写的吗？能分享下么？

MeteoInfo

Iron_Lee 发表于 2016-7-22 08:34
楼主barbs函数是自己写的吗？能分享下么？

MeteoInfo是开源软件

Iron_Lee

MeteoInfo 发表于 2016-7-22 08:35
MeteoInfo是开源软件

！！！！！！
之前没有用过，可以尝试学习一下，谢谢楼主。

chennizi

楼主果然是勇于创新，好样的！NCL里面应该会有绘制TLOGP的类似程序，里面会有方法吧

OomaulenOo

grads也有画t-lnp图的gs，可以参考一下~

jzd19851102

坐等metroinfo的效果超过这个图

KIMO23

谢谢您的分享，T-LnP图的样式还是如jzd19851102网友所附的图形格式违常用格式，再请您参考一下~~~

MeteoInfo

我想问的是用1楼脚本中的数据怎么计算出干绝热线？

兰溪之水

MeteoInfo 发表于 2016-7-22 12:52
我想问的是用1楼脚本中的数据怎么计算出干绝热线？

其他的线都不是额外的数据计算的，比如干绝热线可以用干绝热递减率计算~