This folder contains the original data used to create the "Real-time ionospheric tomography using GNSS-TEC" figure. The three main types of data are as follows 1. Ionospheric electron number density distribution (1320.tomo, 1800.tomo, 2040.tomo) 2. Computation time (log_SPSOLVE.txt, log_RHRCIT_zeroX0.txt, log_RHRCIT.txt) 3. Eigenvalue distributions (all_eigs_PCG_M.txt, all_eigs_CG_M.txt) 1. About Ionospheric Electron Number Density Distribution files *************************************************************************************************************** All files are saved in the same format as .txt, although the extension is saved as .tomo. There is no problem to replace the extension directly. The file description format is as follows -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- # # UTC : Time the file was recorded. # # Number of Boxel (All area) # Number of Boxes # Latitude : Number of boxes in latitude direction # Longitude : Number of boxes in longitude direction # Height : Number of boxes in altitude direction # # *** Plain List *** # Describe the position of the grid (voxel-to-voxel boundaries) in order. # Latitude # +00.00 # ... Describe the positions of the grid in order of latitude (degrees). # +75.00 # Longitude # +100.00 # ... Describe the positions of the grid in order of longitude (degrees). # +175.00 # Height # +070.00 # ... Describe the positions of the grid in order of altitude (km). # +1010.00 # ... # END OF HEADER The end of the header part. The following is the data part. The data is described over multiple lines. All values are in [TECU/km]. The data in one line describes the electron number density of voxels at the same altitude and latitude in decreasing order of longitude (from west to east). Each new line increases the latitude by one (from south to north). A blank line is inserted when the altitude is increased by one. Example) (electron number density at altitude 0~1 latitude 0~1 longitude 0~1) (electron number density at altitude 0~1 latitude 0~1 longitude 1~2) ... (electron number density at altitude 0~1 latitude 0~1 longitude L-1~L) (electron number density at altitude 0~1 latitude 1~2 longitude 0~1) (electron number density at altitude 0~1 latitude 1~2 longitude 1~2) ... (electron number density at altitude 0~1 latitude 1~2 longitude L-1~L) ... ... (electron number density at altitude 0~1 latitude B-1~B longitude 0~1) (electron number density at altitude 0~1 latitude B-1~B longitude 1-2) ... (electron number density at altitude 0~1 latitude B-1~B longitude L-1~L) (electron number density at altitude 1~2 latitude 0~1 longitude 0~1) (electron number density at altitude 1~2 latitude 0~1 longitude 1~2) ... (electron number density at altitude 1~2 latitude 0~1 longitude L-1~L) (electron number density at altitude 1~2 latitude 1~2 longitude 0~1) (electron number density at altitude 1~2 latitude 1~2 longitude 1~2) ... (electron number density at altitude 1~2 latitude 1~2 longitude L-1~L) ... ... (electron number density at altitude 1~2 latitude B-1~B longitude 0~1) (electron number density at altitude 1~2 latitude B-1~B longitude 1~2) ... (electron number density at altitude 1~2 latitude B-1~B, longitude L-1~L) ... ... -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 2. About the computation time files This is the original computation log file from which the computation time for each method shown in Figure 2 in the paper was derived. For each method, the program start time, the problem setting phase (Setting Process End), and the end time of the seeking phase (Solving Process End) are described in chronological order. The name of the file indicates which method each file records the results of. Figure 2 in the paper is obtained by taking the time differences of the records in this file. The format of the files is as follows -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- ProgramStart : 2024-03-12 12:35:33.225393 # Start time of program SettingProcessEnd : jp 2023 196 CIT41_TID_lam=3e8_alpha=12e-1_SPSOLVE 1204 : 2024-03-12 12:36:38.423390 # End time of the problem-setting phase of the 1200 epoch of 2023/196 SettingProcessEnd : jp 2023 196 CIT41_TID_lam=3e8_alpha=12e-1_SPSOLVE 1200 : 2024-03-12 12:36:38.540633 # End time of the problem-setting phase of the 1204 epoch of 2023/196 ... ... SettingProcessEnd : jp 2023 196 CIT41_TID_lam=3e8_alpha=12e-1_SPSOLVE 2158 : 2024-03-12 15:08:19.389898 # End time of the problem-setting phase of the 2168 epoch of 2023/196 SolveProcessEnd : jp 2023 196 CIT41_TID_lam=3e8_alpha=12e-1_SPSOLVE 1200 : 2024-03-12 15:10:27.278615 # End time of the seeking phase of the 1200 epoch of 2023/196 SolveProcessEnd : jp 2023 196 CIT41_TID_lam=3e8_alpha=12e-1_SPSOLVE 1202 : 2024-03-12 15:12:43.140098 # End time of the seeking phase of the 1202 epoch of 2023/196 ... ... SolveProcessEnd : jp 2023 196 CIT41_TID_lam=3e8_alpha=12e-1_SPSOLVE 2158 : 2024-03-13 09:00:28.187528 # End time of the seeking phase of the 2158 epoch of 2023/196 -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 3. About the Eigenvalue Distribution Files This file records the values of the eigenvalue distribution of the coefficient matrix used in Figure 5 in the paper. These values were computed using scipy.sparse.linalg.eigs. Figure 5 is obtained by sorting the numbers in this file in ascending order. The file names represent the following information CG : eigenvalue distribution of the coefficient matrix without preprocessing PCG : eigenvalue distribution of the coefficient matrix with preprocessing The format of the files is as follows -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 1364993903848816.2 # List of Eigenvalues 1340874318177390.0 1340874318124925.0 ... ... 46698439.19192739 46695707.610381745 46694991.7337396 --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------