|
Size: 731
Comment:
|
Size: 2430
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| == Daten == |
|
| Line 9: | Line 11: |
| {{{!#python | == Skript == {{{#!python #!/bin/tcsh |
| Line 12: | Line 18: |
| from mpl_toolkits.basemap import Basemap | #from mpl_toolkits.basemap import Basemap |
| Line 18: | Line 24: |
| #Bild imshow(Temp[0,0:90,130:180],vmin=170,vmax=300) |
#Drucklevel Lv=[1000.0,925.0,850.0,700.0,600.0,500.0,400.0,300.0,250.0,200.0,150.0,100.0, 70.0, 50.0,30.0, 20.0, 15.0, 10.0, 7.0,5.0, 3.0,2.0,1.5,1.0] # Ausschnitt definieren sued=30 nord=0 west=130 ost=180 # Kartendarstellung # global figure(1) imshow(Temp[0,:,:],vmin=170,vmax=300) |
| Line 21: | Line 40: |
| show() | |
| Line 22: | Line 42: |
| #Basemap m = Basemap(projection='ortho',lon_0=0.0,lat_0=75.0,resolution='l') m.bluemarble() |
# Ausschnitt figure(2) imshow(Temp[0,nord:sued,west:ost],vmin=170,vmax=300) colorbar() show() |
| Line 26: | Line 48: |
| #PLOTBEFEHL FEHLT NOCH - Drucklevels auslesen}}} | # Fehlende Profile finden - Wert 1 testnan=Temp[0,nord:sued,west:ost] testnan=testnan<0 # Temperaturprofil figure(3) subplot(1,2,1) plot(Temp[:,30,150],Lv) ylim((1000.0,1.0)) show() subplot(1,2,2) plot(Temp[:,25,150],Lv) ylim((1000.0,1.0)) show() }}} Drucklevel: ||Index of T||TempPresLvls|| || Levels || in mb of T Levels || || 1 || 1000.0|| || 2 || 925.0|| || 3 || 850.0|| || 4 || 700.0|| || 5 || 600.0|| || 6 || 500.0|| || 7 || 400.0|| || 8 || 300.0|| || 9 || 250.0|| || 10 || 200.0|| || 11 || 150.0|| || 12 || 100.0|| || 13 || 70.0|| || 14 || 50.0|| || 15 || 30.0|| || 16 || 20.0|| || 17 || 15.0|| || 18 || 10.0|| || 19 || 7.0|| || 20 || 5.0|| || 21 || 3.0|| || 22 || 2.0|| || 23 || 1.5|| || 24 || 1|| == Notizen == * Drucklevels - (./) * MISSING Values - (./) * Koordinaten - von der Grenzschichtgruppe * Profile plotten - formal ok * Grenzschichtdicke in hPa bestimmen - sind überhaupt Inversionen sichtbar? - Bodendruck herausfinden * Grenzschichtdicke in Meter umrechnen (barometrische Höhenformel) {{attachment: Rtavn00120030330.png}} |
Daten
Temperaturprofile: AIRS
Hilfreiche Links:
Skript
1 #!/bin/tcsh
2
3 from pyhdf import SD
4 from pylab import *
5 #from mpl_toolkits.basemap import Basemap
6
7 fid=SD.SD('AIRS.2003.03.30.L3.RetStd001.v5.0.14.0.G07215021620.hdf')
8 T=fid.select('Temperature_A')
9 Temp=T.get()
10
11
12 #Drucklevel
13 Lv=[1000.0,925.0,850.0,700.0,600.0,500.0,400.0,300.0,250.0,200.0,150.0,100.0,
14 70.0, 50.0,30.0, 20.0, 15.0, 10.0, 7.0,5.0, 3.0,2.0,1.5,1.0]
15
16 # Ausschnitt definieren
17 sued=30
18 nord=0
19 west=130
20 ost=180
21
22 # Kartendarstellung
23 # global
24 figure(1)
25 imshow(Temp[0,:,:],vmin=170,vmax=300)
26 colorbar()
27 show()
28
29 # Ausschnitt
30 figure(2)
31 imshow(Temp[0,nord:sued,west:ost],vmin=170,vmax=300)
32 colorbar()
33 show()
34
35 # Fehlende Profile finden - Wert 1
36 testnan=Temp[0,nord:sued,west:ost]
37 testnan=testnan<0
38
39 # Temperaturprofil
40 figure(3)
41 subplot(1,2,1)
42 plot(Temp[:,30,150],Lv)
43 ylim((1000.0,1.0))
44 show()
45 subplot(1,2,2)
46 plot(Temp[:,25,150],Lv)
47 ylim((1000.0,1.0))
48 show()
Drucklevel:
Index of T |
|
Levels |
in mb of T Levels |
1 |
1000.0 |
2 |
925.0 |
3 |
850.0 |
4 |
700.0 |
5 |
600.0 |
6 |
500.0 |
7 |
400.0 |
8 |
300.0 |
9 |
250.0 |
10 |
200.0 |
11 |
150.0 |
12 |
100.0 |
13 |
70.0 |
14 |
50.0 |
15 |
30.0 |
16 |
20.0 |
17 |
15.0 |
18 |
10.0 |
19 |
7.0 |
20 |
5.0 |
21 |
3.0 |
22 |
2.0 |
23 |
1.5 |
24 |
1 |
Notizen
Drucklevels -
MISSING Values -
- Koordinaten - von der Grenzschichtgruppe
- Profile plotten - formal ok
- Grenzschichtdicke in hPa bestimmen - sind überhaupt Inversionen sichtbar? - Bodendruck herausfinden
- Grenzschichtdicke in Meter umrechnen (barometrische Höhenformel)
