Skip to content

WIP: Updates for v3 NRL data#73

Draft
mathomp4 wants to merge 5 commits intodevelopfrom
feature/72-nrl-v3-data
Draft

WIP: Updates for v3 NRL data#73
mathomp4 wants to merge 5 commits intodevelopfrom
feature/72-nrl-v3-data

Conversation

@mathomp4
Copy link
Copy Markdown
Member

@mathomp4 mathomp4 commented Aug 7, 2025

Closes #72

NRL released v03r00 of the NRLSSI2 ancillary text files, which replace the legacy sunspot_darkening_function/facular_brightening_function columns with Bolfac and Bolspot. The v02r01 files are no longer updated past 2023.

Changes

  • Mg_SB_from_daily_file.py: Updated to detect v03r00 files and convert to legacy Mg/SB format via an empirical least-squares linear mapping learned from the v02/v03 overlap period (1882–2023)
  • TSI_Mg_SB_merged_from_daily.py: Updated plotting to shade the post-2024 v03-derived region; restore ops DATADIR/OUTDIR as default
  • README.md: Substantially expanded with v03r00 download instructions, ordering constraints, and full documentation of the conversion approach including fitted coefficients and caveats

Conversion approach

Since there is no official first-principles formula to convert v03 quantities to legacy units, an empirical linear fit is used over the overlapping historical period (1882–2023):

[Mg_SB] Fitted v03->legacy:
  Mg:  Mg_v02 ≈ 0.007506 * Bolfac + 0.150331  (R^2=0.9049, n=51864)
  SB:  SB_v02 ≈ 1941.389784 * Bolspot + 42.906433  (R^2=0.9940, n=51864)

The SB mapping is very tight (R²=0.994); the Mg mapping has more scatter (R²=0.905) and post-2023 Mg values should be understood as approximations.

TSInMgSB_plot_all_final

@mathomp4 mathomp4 self-assigned this Aug 7, 2025
@mathomp4 mathomp4 added 0 diff The changes in this pull request have verified to be zero-diff with the target branch. Non 0-diff The changes in this pull request are non-zero-diff labels Aug 7, 2025
- Expand README with v03r00 data source instructions, ordering
  constraints, and full section I documenting the empirical
  Bolfac/Bolspot -> legacy Mg/SB conversion approach
- Fix ancilliary -> ancillary typo and other minor typos
- Restore ops DATADIR/OUTDIR as default in TSI_Mg_SB_merged_from_daily.py;
  keep personal paths as commented reference
Use headers, code blocks, tables, blockquotes, and lists
for better rendering on GitHub.
@mathomp4
Copy link
Copy Markdown
Member Author

I've updated this PR with the latest information and better readme, etc.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

0 diff The changes in this pull request have verified to be zero-diff with the target branch. Non 0-diff The changes in this pull request are non-zero-diff

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant