Differences
This shows you the differences between two versions of the page.
Both sides previous revision Previous revision | Next revision Both sides next revision | ||
data:solar_satire [2016/08/08 15:25] jypeter Added warning about the time axis |
data:solar_satire [2016/08/09 14:16] jypeter Added more details and python code to deal with the time axis |
||
---|---|---|---|
Line 26: | Line 26: | ||
* 1st array: **wavelength array** in [nm], listing the center of each wavelength bin\\ '' | * 1st array: **wavelength array** in [nm], listing the center of each wavelength bin\\ '' | ||
* 2nd array: **wavelength bin** in [nm], listing the bin width of each wavelength bin\\ '' | * 2nd array: **wavelength bin** in [nm], listing the bin width of each wavelength bin\\ '' | ||
- | * 3rd array: **time** in [year] (floating numbers)\\ For 14C: '' | + | * 3rd array: **time** in [year] (floating numbers) |
+ | * 69235 time steps for 14C: '' | ||
+ | * 61595 time steps for 10Be: '' | ||
+ | 2015.995, 2015.997]'' | ||
* 4th array: **SSI reconstruction** in [W m-2 nm-1]. SSI is average SSI in corresponding bin. | * 4th array: **SSI reconstruction** in [W m-2 nm-1]. SSI is average SSI in corresponding bin. | ||
Line 52: | Line 55: | ||
</ | </ | ||
- | * The following python code shows how to deal with the original compressed text data. You can also check the {{: | + | * The following python code shows how to deal with the original compressed text data. You can also check the {{: |
+ | * <code python># Get directly the data from the bz2 compressed file | ||
file_in = bz2.BZ2File(input_full_path) | file_in = bz2.BZ2File(input_full_path) | ||
Line 86: | Line 90: | ||
tsi = np.dot(ssi, wl_bin) | tsi = np.dot(ssi, wl_bin) | ||
</ | </ | ||
+ | * The following python example shows how to determine the indices of specific years, before Jan 1850 1st and after. Once you know the indices of Jan 1st and Dec 31st, and if a year is a leap year or not | ||
+ | * <code python> # 1 value per day, AFTER (and including) Jan 1st 1850 | ||
+ | >>> | ||
+ | >>> | ||
+ | >>> | ||
+ | (8605, 8970) | ||
+ | >>> | ||
+ | 365 | ||
+ | >>> | ||
+ | >>> | ||
+ | (1850.0, 1850.9973, 1851.0) | ||
+ | |||
+ | # 1 value per year, strictly BEFORE Jan 1st 1850 | ||
+ | >>> | ||
+ | >>> | ||
+ | (6805, 50.5)</ | ||
+ | |||
===== References ===== | ===== References ===== | ||