Code to construct sunspot number series from the series recorded by many observers with their different personal/instrumental biases
The main module.
A module with useful routines needed to determine the optimal sunspot number, including the reading routine for the CSV export from Slunce software for the CESLOPOL observers.
Contains the configuration directives influencing the performance of the code
A pseudocode showing how to do some trivial plots and comparisons of the results
A compiled database used in the study submitted to Solar Physics journal, it contains example observations.
WDC-SILSO 2.0 sunspot number
Revised group number
A mock-up script (pseudo-code) to demonstrate how to fill in the database Pandas DataFrame compatible with the main code.
databaze={'Pozorovatel':[], 'Pristroj':[], 'Datum':[], 'UT':[], 'Q':[], 'Obr.':[], 'g':[], 'f':[], 'R':[]}
- 'Pozorovatel': STRING: The name of the observer
- 'Pristroj': STRING: A unique identification of the observer's instrument, it may also identify the station
- 'Datum': DATETIME timestamp: The date (the day) of the observation
- 'UT': FLOAT: the time of observation in UT in hours
- 'Q': FLOAT: the quality of observations, following the empirical table of visibility of structures
- 'Obr.': FLOAT: "obraz" (image), the another index assessing the seeing
- 'g': FLOAT: the total number of groups
- 'f': FLOAT: the total number a sunspots
- 'R': FLOAT: the relative number as 10g+f
The mandatory fields are 'Pozorovatel', 'Pristroj', 'Datum', 'Q', 'g', 'f', 'R', other are optional
target={'Datum':[], 'g':[], 'f':[], 'R':[], 'Q':[5]}
- 'Datum': DATETIME timestamp: The date (the day) of the observation
- 'g': FLOAT: the total number of groups
- 'f': FLOAT: the total number a sunspots
- 'R': FLOAT: the relative number as 10g+f
- 'Q': FLOAT: the quality of observations, default is 5 in the case of reduction to optimal conditions
uncertainty={'Datum':[], 'g':[], 'f':[], 'R':[]}
- 'Datum': DATETIME timestamp: The date (the day) of the observation
- 'g': FLOAT: uncertainty of the total number of groups
- 'f': FLOAT: uncertainty of the total number a sunspots
- 'R': FLOAT: uncertainty the relative number as 10g+f
k_table={'Datum':[], 'Jmeno':[], 'k_f':[], 'k_g':[], 'Pruchod':[]}
- 'Datum': DATETIME timestamp: The date (the day) of the observation
- 'Jmeno': STRING: the name of the observer
- 'k_f': FLOAT: observer's conversion coefficient to targer series for f (total sunspots number)
- 'k_g': FLOAT: observer's conversion coefficient to targer series for g (number of sunspot groups)
- 'Pruchod': FLOAT: identification of the iteration for which these coefficients are valid. Integer values indicate the full-interation (forth and back), whereas the values ending with 0.5 indicate the half-iteration (forth only)