Differences between revisions 23 and 26 (spanning 3 versions)
⇤ ← Revision 23 as of 2010-01-18 14:15:24 →
Size: 1693
Editor: AnnaFriedrichs
Comment:
← Revision 26 as of 2010-01-18 14:38:17 → ⇥
Size: 2357
Editor: AnnaFriedrichs
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
== Daten ==
Line 9: Line 11:
== Skript ==
Line 10: Line 14:
#!/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()

# 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()
Line 57: Line 94:
 * Drucklevels  * Drucklevels - ok
Line 59: Line 96:
 * MISSING Values  * MISSING Values -ok
Line 61: Line 98:
 * Koordinaten  * Koordinaten - von der Grenzschichtgruppe
Line 63: Line 100:
 * Profile ploten  * Profile plotten - formal ok
Line 65: Line 102:
 * Grenzschichtdicke in hPa bestimmen  * Grenzschichtdicke in hPa bestimmen - sind überhaupt Inversionen sichtbar?

Daten

Temperaturprofile: AIRS

Hilfreiche Links:

Documentation

Ftp

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

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 - ok
  • MISSING Values -ok
  • Koordinaten - von der Grenzschichtgruppe
  • Profile plotten - formal ok
  • Grenzschichtdicke in hPa bestimmen - sind überhaupt Inversionen sichtbar?
  • Grenzschichtdicke in Meter umrechnen (barometrische Höhenformel)

LehreWiki: OpenSource2010/Project/Project Idea2010/Grenzschichtdicke (last edited 2011-01-17 09:45:59 by anonymous)