Umfrage Skew-t LopP Plot Payerne

[Willi]

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

[Bernhard Oker]

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

[Dani Murer]

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.

[Willi]

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

[Bernhard Oker]

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

[Bernhard Oker]

Neu gestaltete Temps mit Höhenangaben sind jetzt Online.
http://www.geocities.com/okerbernhard/c ... lookch.htm

Gruss Bernhard

[Markus Pfister]

Super, Bernhard, danke für die Arbeit!
Auf dass die dunkelgrünen Flächen noch
ein paar Mal schön gross werden...

Gruss

Markus

[Bernhard Oker]

Demo Version des Programs kann unter folgendem Link heruntergeladen werden:
http://www.geocities.com/okerbernhard/s ... nalyse.zip

MfG Bernhard