Hallo Bernhard
Super!
Noch eine Anregung, was haelst Du davon, auf der rechten Seite des Plots eine Hoehenskala anzufuegen (natuerlich berechnet aus P und T)?
Gruss Willi
Werbung
Umfrage Skew-t LopP Plot Payerne
- Bernhard Oker
- Moderator
- Beiträge: 6238
- Registriert: Do 16. Aug 2001, 11:02
- Geschlecht: männlich
- Wohnort: 8902 Urdorf
- Hat sich bedankt: 1582 Mal
- Danksagung erhalten: 2575 Mal
Umfrage Skew-t LopP Plot Payerne
Wäre nicht schwer die Höhen zu den schon eingezeichneten Druckwerten (1000 - 100 mB) auf der rechten Seite hinzuzufügen.
Die Formel zur Berechnung der Höhe aus P und T habe ich aber noch nicht gefunden. In den Temps aus dem Internet stehen die Höhen jeweils schon drin. Hast Du die Formel gerade im Kopf?
Gruss Bernhard
Die Formel zur Berechnung der Höhe aus P und T habe ich aber noch nicht gefunden. In den Temps aus dem Internet stehen die Höhen jeweils schon drin. Hast Du die Formel gerade im Kopf?
Gruss Bernhard
Bernhard Oker - Urdorf (ZH/CH) - Meine Webseiten "Never Stop Chasing!"
Umfrage Skew-t LopP Plot Payerne
Hoi Bernhard,
suchst Du etwa die barometrische Höhenformel?
Sie lautet:
z = [ (R * T(m) ) / g ] * ln [ p / p(0) ]
wobei z die Höhe, R die Gaskonstante für Luft (287 J/(Kg*K)), T(m) die Mitteltemperatur zwischen den Druckniveaus p(0) und p, p(0) der Luftdruck am Boden, p der Luftdruck in der Höhe z und g die Schwerebeschleunigung der Erde sind.
suchst Du etwa die barometrische Höhenformel?
Sie lautet:
z = [ (R * T(m) ) / g ] * ln [ p / p(0) ]
wobei z die Höhe, R die Gaskonstante für Luft (287 J/(Kg*K)), T(m) die Mitteltemperatur zwischen den Druckniveaus p(0) und p, p(0) der Luftdruck am Boden, p der Luftdruck in der Höhe z und g die Schwerebeschleunigung der Erde sind.
- Willi
- Administrator
- Beiträge: 8881
- Registriert: Fr 10. Aug 2001, 16:16
- Wohnort: 8143 Sellenbüren
- Hat sich bedankt: 3994 Mal
- Danksagung erhalten: 3892 Mal
- Kontaktdaten:
Umfrage Skew-t LopP Plot Payerne
Hoi Bernhard
Du hast ja die Höhen aus den Standardniveaus,
am besten machst Du dazwischen eine Interpolation mit der barometrischen Höhenformel.
Gruss Willi
Du hast ja die Höhen aus den Standardniveaus,
am besten machst Du dazwischen eine Interpolation mit der barometrischen Höhenformel.
Gruss Willi
Gruss Willi
Immer da wenn's wettert
Immer da wenn's wettert
- Bernhard Oker
- Moderator
- Beiträge: 6238
- Registriert: Do 16. Aug 2001, 11:02
- Geschlecht: männlich
- Wohnort: 8902 Urdorf
- Hat sich bedankt: 1582 Mal
- Danksagung erhalten: 2575 Mal
Umfrage Skew-t LopP Plot Payerne
Bin inzwischen fündig geworden bei den Fortran Formeln von wo ich auch alle anderen her habe. Die So berechneten Werte stimmen auf den Meter genau mit denen überein welche im Sounding aus dem Internet eingetragen sind.
function calcHeight(p1 : real;
p2 : real;
t1 : real;
t2 : real;
td1 : real;
td2 : real) : real;
// g.s. stipanuk 1973 original version.
// reference stipanuk paper entitled:
// "algorithms for generating a skew-t, log p
// diagram and computing selected meteorological
// quantities."
// atmospheric sciences laboratory
// u.s. army electronics command
// white sands missile range, new mexico 88002
// 33 pages
// baker, schlatter 17-may-1982
// this function returns the thickness of a layer bounded by pressure
// p(1) at the bottom and pressure pt at the top.
// on input:
// p = pressure (mb). note that p1 > p2.
// t = temperature (celsius)
// td = dew point (celsius)
// on output:
// z = geometric thickness of the layer (m)
// the algorithm involves numerical integration of the hydrostatic
// equation from p1 to p2. it is described on p.15 of stipanuk
// (1973).
const
// C1 = 0.001 * (1.0 / eps - 1.0) where eps = 0.62197 is the ratio of the
// molecular weight of water to that of
// dry air. the factor 1000.0 converts the
// mixing ratio w from g/kg to a dimensionless ratio.
// C2 = r / (2.0 * g) where r is the gas constant for dry air
// (287 kg/joule/deg k) and g is the acceleration
// due to the earth's gravity (9.8 m/s**2). the
// factor of 2 is used in averaging two virtual temperatures.
C1 = 0.0006078;
C2 = 14.64285;
var
a1, a2, z : real;
begin
t1 := t1 + 273.15;
t2 := t2 + 273.15;
a1 := t2 * (1.0 + C1 * w(td2, p2));
a2 := t1 * (1.0 + C1 * w(td1, p1));
z := z + C2 * (a1 + a2) * ln(p1 / p2);
Result := z;
end;
function w(t, p : real) : real;
var
x : real;
// This function returns the mixing ratio [g/kg] given
// the dew point [°C] and pressure [mB].
// If the temperture is input instead of the
// dew point, then saturation mixing ratio [g/kg] is returned.
// The formula is found in most meteorological texts.
begin
x := esat(t);
Result := 622 * x / (p - x);
end;
function esat(t : real) : real;
var
p1, p2, c1 : real;
// This function returns the saturation vapor pressure over
// water [mB] given the temperature [°C].
// The algorithm is due to nordquist, w.s.,1973: "numerical approxima-
// tions of selected meteorological parameters for cloud physics prob-
// lems," ecom-5475, atmospheric sciences laboratory, u.s. army
// electronics command, white sands missile range, new mexico 88002.
begin
t := t + 273.15;
p1 := 11.344 - (0.0303998 * t);
p2 := 3.49149 - (1302.8844 / t);
c1 := 23.832241 - (5.02808 * log10(t));
Result := power(10, (c1 - 1.3816E-7 * power(10, p1) + 8.1328E-3 * power(10, p2) - 2949.076 / t));
end;
Hoffentlich ist es einigermassen verständlich
Gruss Bernhard
function calcHeight(p1 : real;
p2 : real;
t1 : real;
t2 : real;
td1 : real;
td2 : real) : real;
// g.s. stipanuk 1973 original version.
// reference stipanuk paper entitled:
// "algorithms for generating a skew-t, log p
// diagram and computing selected meteorological
// quantities."
// atmospheric sciences laboratory
// u.s. army electronics command
// white sands missile range, new mexico 88002
// 33 pages
// baker, schlatter 17-may-1982
// this function returns the thickness of a layer bounded by pressure
// p(1) at the bottom and pressure pt at the top.
// on input:
// p = pressure (mb). note that p1 > p2.
// t = temperature (celsius)
// td = dew point (celsius)
// on output:
// z = geometric thickness of the layer (m)
// the algorithm involves numerical integration of the hydrostatic
// equation from p1 to p2. it is described on p.15 of stipanuk
// (1973).
const
// C1 = 0.001 * (1.0 / eps - 1.0) where eps = 0.62197 is the ratio of the
// molecular weight of water to that of
// dry air. the factor 1000.0 converts the
// mixing ratio w from g/kg to a dimensionless ratio.
// C2 = r / (2.0 * g) where r is the gas constant for dry air
// (287 kg/joule/deg k) and g is the acceleration
// due to the earth's gravity (9.8 m/s**2). the
// factor of 2 is used in averaging two virtual temperatures.
C1 = 0.0006078;
C2 = 14.64285;
var
a1, a2, z : real;
begin
t1 := t1 + 273.15;
t2 := t2 + 273.15;
a1 := t2 * (1.0 + C1 * w(td2, p2));
a2 := t1 * (1.0 + C1 * w(td1, p1));
z := z + C2 * (a1 + a2) * ln(p1 / p2);
Result := z;
end;
function w(t, p : real) : real;
var
x : real;
// This function returns the mixing ratio [g/kg] given
// the dew point [°C] and pressure [mB].
// If the temperture is input instead of the
// dew point, then saturation mixing ratio [g/kg] is returned.
// The formula is found in most meteorological texts.
begin
x := esat(t);
Result := 622 * x / (p - x);
end;
function esat(t : real) : real;
var
p1, p2, c1 : real;
// This function returns the saturation vapor pressure over
// water [mB] given the temperature [°C].
// The algorithm is due to nordquist, w.s.,1973: "numerical approxima-
// tions of selected meteorological parameters for cloud physics prob-
// lems," ecom-5475, atmospheric sciences laboratory, u.s. army
// electronics command, white sands missile range, new mexico 88002.
begin
t := t + 273.15;
p1 := 11.344 - (0.0303998 * t);
p2 := 3.49149 - (1302.8844 / t);
c1 := 23.832241 - (5.02808 * log10(t));
Result := power(10, (c1 - 1.3816E-7 * power(10, p1) + 8.1328E-3 * power(10, p2) - 2949.076 / t));
end;
Hoffentlich ist es einigermassen verständlich
Gruss Bernhard
Bernhard Oker - Urdorf (ZH/CH) - Meine Webseiten "Never Stop Chasing!"
- Bernhard Oker
- Moderator
- Beiträge: 6238
- Registriert: Do 16. Aug 2001, 11:02
- Geschlecht: männlich
- Wohnort: 8902 Urdorf
- Hat sich bedankt: 1582 Mal
- Danksagung erhalten: 2575 Mal
Umfrage Skew-t LopP Plot Payerne
Neu gestaltete Temps mit Höhenangaben sind jetzt Online.
http://www.geocities.com/okerbernhard/c ... lookch.htm
Gruss Bernhard
http://www.geocities.com/okerbernhard/c ... lookch.htm
Gruss Bernhard
Bernhard Oker - Urdorf (ZH/CH) - Meine Webseiten "Never Stop Chasing!"
-
- Beiträge: 2348
- Registriert: So 19. Aug 2001, 20:08
- Geschlecht: männlich
- Wohnort: 9056 Gais
- Hat sich bedankt: 774 Mal
- Danksagung erhalten: 2333 Mal
- Kontaktdaten:
Umfrage Skew-t LopP Plot Payerne
Super, Bernhard, danke für die Arbeit!
Auf dass die dunkelgrünen Flächen noch
ein paar Mal schön gross werden...
Gruss
Markus
Auf dass die dunkelgrünen Flächen noch
ein paar Mal schön gross werden...
Gruss
Markus
- Bernhard Oker
- Moderator
- Beiträge: 6238
- Registriert: Do 16. Aug 2001, 11:02
- Geschlecht: männlich
- Wohnort: 8902 Urdorf
- Hat sich bedankt: 1582 Mal
- Danksagung erhalten: 2575 Mal
Umfrage Skew-t LopP Plot Payerne
Demo Version des Programs kann unter folgendem Link heruntergeladen werden:
http://www.geocities.com/okerbernhard/s ... nalyse.zip
MfG Bernhard
http://www.geocities.com/okerbernhard/s ... nalyse.zip
MfG Bernhard
Bernhard Oker - Urdorf (ZH/CH) - Meine Webseiten "Never Stop Chasing!"