Zobrazit diskusi v plném vzhledu
Dotazy ohledně radarů, družic a celkově o interpelaci meteorologických dat.
Odeslat odpověď

Kimatická simulace s jednou hladinou na Icosahedrální síti

20. 05. 2015 22:34

První dobře vypadacící výsledek - klimatické modelování 2D s využitím šestiúhelníkové sítě, stupeň rozlišení 5 - 40 962 bodů. Výpočet trval poměrně dlouho - přibližně 60 dní. Nejvíce výpočetního času sebrala interpolace (z 40962 bodů na 2000x1000, jeden obrázek trval okolo 23 s, celkem 32777 map pro 6 veličin - 196 662 map).

Pro exoplanetu Tau Cet e (ještě nebyla existence zcela prokázána) zvoleny parametry: Tlak 2 bar, tepelná kapacita 1007 J/K/kg, rovnovážná teplota 289 K (Země 255 K), Sklon rotační osy 37° (Země 23,45°), délka siderického roku 14525568 s/ 168.12 našich dní, perioda rotace 158 400 s, délka dne 160 127 s, eccentricita 0,05, rozdíl perihelium-rovnodenost 15°.

Teplota
" http://uloz.to/xHpgWTPK/qx-000000040962 ... d-temp-avi "

Forcing teploty
" http://uloz.to/xGAhCtrQ/qx-000000040962 ... rcingt-avi "

Změna teploty
" http://uloz.to/xk2jqKEA/qx-000000040962 ... rtialt-avi "

Insolace
" http://uloz.to/xdrVAtK1/qx-000000040962 ... ated-s-avi "

Radiační bilance
" http://uloz.to/xMEpUhxC/qx-000000040962 ... ated-q-avi "

K - Sinus zenitového úhlu
" http://uloz.to/xA3LH2Fm/qx-000000040962 ... ated-k-avi "

Jiný běh modelu - s atmosférou jen 0.2 bar a stejnými ostatními hodnotami.
Mapy jsou bez interpolace, neboť interpolace trvá doluho. Video vytvořeno z 182 345 map.
Teplota " http://uloz.to/xvRQGhuS/map-classic-mola-temp-test-avi "
Ro rychlejší download " https://mega.co.nz/#!VxJgFQqT!LyKmgxe6j ... KGVN6sEnik "

Parametry programu :
SigmaST=5.670373E-8 ! St.-B. konstanta
Mplanet=2.565E+25 Hmotnost planety
a=11582000.0 ! Rovníkový poloměr
asp=21000.0 ! Rozdíl poloměru rovník a póly
g=(6.6732E-11)*Mplanet/(a*a) ! povrchová gravitace, hmotnost 2.565E25 kg/4.25 Me, poloměr 11 582 km, hustota 3941.4
Efektivní tíhové zrychlení:
!gm=((6.6732E-11)*Mplanet/((a-(asp*cos(souradnice(ID,2))))**2)-((2.0*Pi*(1.0/DelkaDne))**2)*((a-(asp*cos(souradnice(ID,2))))*cos(souradnice(ID,2)))-(U0*U0+V0*V0+W0*W0)/((a-(asp*cos(souradnice(ID,2))))**2)-(KorF*U0*(-sin(souradnice(ID,3)))))
Sklon=0.64577182323790194346176558434079 ! Sklon rotační osy, 37 deg
e=0.05 ! ekcentricita
DelkaRoku=14525568.0 ! Sidercká, 168.12 Earth days
DelkaDne=158400.0 ! Siderická délka dne - perioda rotace, 44 hodin
DelkaDneSyn=160127.34 ! Sinodická délka dne, delší o 1727.34 s
EquinoxShift=0.26179939 ! Úhlový rozdíl v dráze mezi periheliem a rovnodenností, perihelium 15°před rovnodenností
HlAxis=82579200000.0 ! 8.25792E+10 m or 0.552 AU Hlavní polosa
KonstT=44.682256021 ! Konstanta teplota vs. výška H/T
Tsl=5344.0 ! Povrchová teplota hvězdy
aSL=551928000.0 ! POloměr hězdy, 0.793 Slunečního
Time=0.0 ! Čas start
dT=4.0 ! Časový krok
TimeW=14525568.0 ! První zápis hodnot (po 1 roce)
TimeWG=14525568.0 ! Čas pro první zápis globálních průměrů
TimeWSU=14525568.0 ! Čas pro start výpoču dlouhodobých průměrů
TimeWO=1329.5 ! Čas mezi zápisem hodnot pro vytvoření mapy
TimeWOG=41.55 !Čas mezi zápisi globálních průměrů
TimeMax=58102272.0 ! Délka běhu modelu - 4 místní roky
CasSU=0 ! Proměnná času pro dlouhodobé průměry, 0 na začátku, integer
PHIpsu=159686.38 ! průměrný geopotenciál na začátku, tlak 2 bar, hustota 1.252455, g 12.76 ms-2
Hsu=PHIpsu/g ! Průměrná výška na začaátku (tloušťka atmosféry)
EPSm=0.55 ! ckoeficien emisivitz vyzařování (může být 0-1)
Albedo=0.303 ! Albedo - průměr. Pro albedo 303 rovnovážná teplota 280.08 K (rovnovážná subsellární teplota 396.09 K) - hustota atmosféry 1.3954 kg/m3.
AlbedoAtm=0.434 ! Albedo atmosféry (pro nekonečnou vzduchovou hmotu, pro nulovou vzuchovou hmotu je to 0.303)
OpAtm=0.785 ! Propustnost atmosféry (vzduchová hmota 1, v zenitu)
Rho=1.252455 ! Hustota atmosféry homogenní, odpovídající 1 bar a 280.08 K, plynová konstanta R 288.05 pro suchý vzduch (R = 288.844 for mixing ratio 10 g/kg)
KpMerna=1007.0 ! Měrná tepelná kapacita cp - J/(kg.K).

Re: Kimatická simulace s jednou hladinou na Icosahedrální sí

25. 07. 2015 16:27

Nový kodek nalezen a soubory jsou nyní podstatně menší. Jedná se o h265, kodování nejpomalejší, quality 17. Rozlišení 2000x1000 30 fps. Velikost je nesrovnatelně nižší, kodování videa je výpočetně extrémně náročné. Staré velké soubory budou smazány.

Teplota (K)
http://uloz.to/xfCFumK1/qx-000000040962 ... d-temp-mp4

Forcing teploty dT/dt (K/s)
http://uloz.to/xc11MtP6/qx-000000040962 ... rcingt-mp4

Insolace S (W/m2)
http://uloz.to/xnd3crkh/qx-000000040962 ... ated-s-mp4

Radiační bilance Q (W/m2)
http://uloz.to/xGpp1gny/qx-000000040962 ... ated-q-mp4

Sinus výšky Slunce nad obzorem K
http://uloz.to/xwvvYUo4/qx-000000040962 ... ated-k-mp4
Odeslat odpověď