Franta (Pondělí 20. prosince 2010)

No já mám vzorce dva + tabulku od Ministra, která je prý z hvězdářské ročenky a odpovídá prý Valouchovým tabulkám, které nemám a ještě jsem se nedostal k tomu je sehnat. Zatím ovšem se zábývám jen armosferickou refrakcí. A s horizontem mám problém, to už jsem se svěřil dříve.

No, jo. Typicky jsem zase zapomněl, že to chce vysvětlivky. Je to funkce z astronomické knihovny, která se jmenuje libnova. Takže takhle v C vypadá zdrojový kód funkce.

/*! fn double ln_get_refraction_adj (double altitude, double atm_pres, double temp)
* param altitude The altitude of the object above the horizon in degrees
* param atm_pres Atmospheric pressure in milibars
* param temp Temperature in degrees C.
* return Adjustment in objects altitude in degrees.
*
* Calculate the adjustment in altitude of a body due to atmosphric
* refraction. This value varies over altitude, pressure and temperature.
*
* Note: Default values for pressure and teperature are 1010 mBar and 10C
* respectively.
*/
/*
double ln_get_refraction_adj (double altitude, double atm_pres, double temp)
{
long double R;

/* equ 16.3 */ /*
R = 1.0 / tan (ln_deg_to_rad (altitude + (7.31 / (altitude + 4.4))));
R -= 0.06 * sin (ln_deg_to_rad (14.7 * (R / 60.0) + 13.0));

/* take into account of atm press and temp */ /*
R *= ((atm_pres / 1010) * (283 / (273 + temp)));

/* convert from arcminutes to degrees */ /*
R /= 60.0;

return R;
}

Výsledky podle Ministra (oranžově) a podle libnova (žlutě) jsou v
pefr.net/refrakce.xls