Kalibnracija magnitude glede na standardne magnitude B,V

Uporabili bomo Vezier in iskali naše zvezde po koordinatah, za referenčni katalog izberemo:

    Montgomery K.A., Marschall L.A., Janes K.A.
    (1993) Astron. J. 106, 181
    CCD photometry of the old open cluster M67
    1993AJ....106..181M 

To je katalog 54 iz http://www.univie.ac.at/webda/ Dostop do magnitud pa je pri Webdi je malo bolj neroden
ker moramo zvezde iskat v dveh tabelah, pa tudi nima vgrajenega iskalnika za koordinate. 

Zato raje uporabimo http://vizier.u-strasbg.fr/viz-bin/VizieR-3 
za katalog vpišemo "Montgomery K.A., Marschall L.A., Janes K.A."  ali pa kar "J/AJ/106/181" in izberemo table3.  

Zvezde lahko iščemo posamič, ali pa izberemo zavihek "List Of Tagret" tu pa označimo 
"single table" in "add your input as first column" ter pripravimo datoteko, 
iz LO_M67VB.wcs izberemo 30 ne najsvetlejših zvezd, naprimer od zaporecne številke 20 do 50. 
V datoteki za upload morajo biti samo koordinate (brez magnitud).

Iz zadtkov sestavimo tabelo v takšni obliki:

#  ----------------izmerjeno------------   ---katalog----  -----izracun----
#    rek          dek        v      b         V      B-V       v-V     b-B
132.89311342  11.852570599  14.14  15.15    10.940  0.220  
 132.8100906   11.74987854  14.21  15.37    11.078  0.430  
132.88588056  11.813980279  14.29  15.26    11.063  0.190  
 132.8211742  11.845753505  14.35  16.08    11.200  1.080  
 132.8367605  11.890405852  14.46  15.60    11.251  0.415  
132.92669177  11.856178819  14.48  16.20    11.330  1.070  
....

Iz tabele izračunamo, recimo v preglednici, aditivni konstanti za V in B
                V            B
'POVPREČJE  3,179279      3,899375
'STD.DEV.   0,038882      0,029192


Če kakšna zvezda močno izstopa jo izbrišemo iz tabele. Lako tudi preverimo ker so nekatere zvezde 
znane kot spremenjive, teh seveda ne vzamemo za kalibracijo. 

Zgaj v gnuplotu lahko narišem HR diagram:

gnuplot> plot [0:2][18:10] "LO_M67VB.wcs"  u ($4-3.89-$3+3.179):($3-3.179)