refrakcje

 Śnieżnik -> Alpy 



Radiosondy

University of Wyoming Atmospheric Science Radiosonde Archive


https://weather.uwyo.edu/upperair/sounding.shtml


Wykorzystałem skrypt napisany dla Stoja

27.12.2025 Prostejov


27.12.2025 Wiedeń





Symulacja refrakcja 0.45
spłaszczyłem nieco panoramę zmnieniając współczynnik wzrostu na 0.5

Zdjęcie z DO z nałożonymi azymutami
próba nazwania szczytów




Referencyjny 1

Ve Vrsich 590
49.586310444, 16.627129256
Odległość: 70.84 km
Azymut: 192.92°


Referencyjny 2
Hluboke Dvory
49.387976875, 16.501679423
Odległość: 94.43 km
Azymut: 195.34°

Maly Chlum
49.451275951, 16.585407257
Odległość: 86.15 km
Azymut: 192.68°

Radostice
49.130456080, 16.472290710
Odległość: 122.76 km
Azymut: 192.83°

Hochwechsel
47.530354880, 15.913831085
Odległość: 305.40 km
Azymut: 193.27°

Stuchleck
47.574765472, 15.789550790
Odległość: 302.77 km
Azymut: 195.19°

Heukuppe
47.689290428, 15.689361647
Odległość: 292.48 km
Azymut: 197.24°

Waxriegel (Schneeberg)
47.760299656, 15.830312201
Odległość: 282.04 km
Azymut: 195.63°

Hoher Hengst
47.756595502, 15.865440022
Odległość: 281.78 km
Azymut: 195.10°


Klosterwappen
47.767278230, 15.804367216
Odległość: 281.79 km
Azymut: 196.06°

Hoher Eisernes Tor
48.007209162, 16.134916257
Odległość: 250.09 km
Azymut: 192.23°

Niederwechsel MAX
47.518301721, 15.951186797
Odległość: 306.11 km
Azymut: 192.70°

Frauenalpen
47.559942080, 15.964364169
Odległość: 301.38 km
Azymut: 192.70°

Kampstein
47.561526381, 15.986778345
Odległość: 300.86 km


Hocheck
47.996773703, 15.951770049
Odległość: 254.31 km









Opis działania skryptu – Refrakcja i warstwy inwersji

Skrypt analizuje dane z sondowań atmosferycznych w celu wykrywania warstw inwersji temperatury oraz obliczania współczynnika refrakcji atmosferycznej. Dodatkowo klasyfikuje typ mirażu, który może wystąpić.

1. Dane wejściowe

Skrypt przyjmuje dane z sondowania atmosferycznego zawierające:

  • Ciśnienie \( P \) [hPa]
  • Wysokość \( h \) [m]
  • Temperaturę \( T \) [°C]
  • Punkt rosy \( T_d \) [°C]

2. Wykrywanie warstw inwersji

Warstwa inwersji jest obecna, gdy: \[ T_{i} > T_{i-1} \] co oznacza wzrost temperatury z wysokością – zjawisko odwrotne do normalnego spadku temperatury.

3. Obliczanie współczynnika refrakcji

Indeks refrakcji atmosfery \( n \)

Zgodnie z meteorologicznym modelem indeksu refrakcji: \[ n = 1 + 10^{-6} \left( \frac{77.6 \cdot P}{T_K} + \frac{373000 \cdot e}{T_K^2} \right) \] gdzie:

  • \( P \) – ciśnienie w hPa
  • \( T_K = T + 273.15 \) – temperatura w kelwinach
  • \( e \) – ciśnienie pary wodnej (patrz niżej)

Ciśnienie pary wodnej \( e \)

Obliczane na podstawie punktu rosy: \[ e = 6.112 \cdot \exp \left( \frac{17.67 \cdot T_d}{T_d + 243.5} \right) \]

Gradient współczynnika refrakcji

Obliczamy zmianę indeksu refrakcji w danej warstwie: \[ \frac{dn}{dh} = \frac{n_2 - n_1}{h_2 - h_1} \]

Współczynnik refrakcji \( k \)

Wzór ogólny na obliczenie współczynnika refrakcji: \[ k = - \frac{R}{H} \cdot \frac{dn}{dh} \] gdzie:

  • \( R = 6\,371\,000 \, \text{m} \) – promień Ziemi
  • \( H = 8000 \, \text{m} \) – typowa skala wysokości atmosfery

Ostatecznie wartość \( k \) jest mnożona przez 1000 i zaokrąglana do 3 miejsc po przecinku.

4. Klasyfikacja mirażu

Na podstawie wartości współczynnika refrakcji:

  • \( k > 0.1 \) – miraż górny
  • \( k < -0.1 \) – miraż dolny
  • Inaczej – brak mirażu

5. Wyniki

Skrypt wypisuje dane o każdej wykrytej warstwie:

  • Zakres wysokości
  • Zmiana temperatury
  • Wartość \( k \)
  • Typ mirażu (jeśli występuje)

6. Eksport i wykres

Użytkownik może zapisać dane do pliku CSV. Dodatkowo wartości \( k \) są przedstawiane na wykresie słupkowym (biblioteka Chart.js).


2025. Analiza optycznych efektów atmosferycznych na podstawie danych radiosond.



Komentarze