## page was renamed from OpenSource2010/Project/Project Idea/Grenzschichtdicke
== Daten ==

Temperaturprofile: [[http://disc.gsfc.nasa.gov/AIRS/ | AIRS]]

Hilfreiche Links: 

[[http://disc.sci.gsfc.nasa.gov/AIRS/documentation/v5_docs/AIRS_V5_Release_User_Docs/V5_L3_QuickStart.pdf | Documentation]]

[[ftp://acdisc.gsfc.nasa.gov/ftp/data/s4pa/Aqua_AIRS_Level3/AIRX3STD.005 | Ftp]]

== Skript ==

{{{#!python
#!/bin/tcsh

from pyhdf import SD
from pylab import *
from mpl_toolkits.basemap import Basemap

fid=SD.SD('AIRS.2003.03.30.L3.RetStd001.v5.0.14.0.G07215021620.hdf')
T=fid.select('Temperature_A')
Temp=T.get()


#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]


# Fehlende Profile finden - Wert 1 
testnan=Temp[0,:,:]
testnan=testnan<0

# Temperaturprofile
figure(1)
subplot(1,2,1)
plot(Temp[:,16,165],Lv)
ylim((1000.0,1.0))
title('over ice, lat=-15,lon=74')
show()
subplot(1,2,2)
plot(Temp[:,16,192],Lv)
ylim((1000.0,1.0))
title('over open water, lat=12,lon=74')
show()
savefig('Soundings.png',dpi=75)

# Kartendarstellung der vorhandenen Profile
figure(2)
m = Basemap(width=2400000,height=1600000,projection='stere',lat_ts=77.0,lon_0=15.0,lat_0=77.0,resolution='l')
m.drawcoastlines()

for i in range(0,30):
    for j in range(160,200):
        if testnan[i,j]==0:
            x,y=m(j-180.,90.-i)
            m.plot(x,y,'r.')

x,y=m(-3.,77.)
m.plot(x,y,'b.')
x,y=m(10.,77.)
m.plot(x,y,'b.')
title('grid points with profiles(red), spectrum (blue)')

show()
savefig('Map.png',dpi=75)
}}}

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

 * (./) Karte mit vorhandenen Profilen - gibt es bei Wolken nie Daten? ja 
 * (./) Daten zum testen:

   [[http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Hornsund.2009073.terra.250m | 14.03.09]]

   [[http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Hornsund.2009077.terra.250m | 18.03.09]]

   [[http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Hornsund.2009089.terra.250m | 30.03.09]]

   [[http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Hornsund.2009090.terra.250m | 31.03.09]]

   [[http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Hornsund.2009091.terra.250m | 01.04.09]]

 * (./) Koordinaten - von der Grenzschichtgruppe: Breite 77°-79°; Länge 2° bis 10°


 * (./) Profile plotten

 * (./) Grenzschichtdicke in hPa bestimmen: Werte: jeweils 14 Profile
 
 * (./) sind überhaupt Inversionen sichtbar? nein - Bodendruck herausfinden

{{attachment:Rtavn00120030330.png}}

 * (./) Grenzschichtdicke in Meter umrechnen barometrische Höhenformel: {{{#!latex
$\Delta z = \frac{R_L T}{g}ln\left(\frac{p}{p_s}\right)$}}}
R=287,1J/kgK; Annahme: T=const.; p=Druckniveau der Wolkenobergrenze, ps=Bodendruck aus Karte (s.o.)