GSF — Solar Modulation

import%20marimo%0A%0A__generated_with%20%3D%20%220.23.15%22%0Aapp%20%3D%20marimo.App(%0A%20%20%20%20width%3D%22full%22%2C%0A%20%20%20%20app_title%3D%22GSF%20%E2%80%94%20Solar%20Modulation%22%2C%0A%20%20%20%20css_file%3D%22gallery.css%22%2C%0A%20%20%20%20html_head_file%3D%22gallery_head.html%22%2C%0A)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%0A%20%20%20%20return%20(mo%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20Solar%20Modulation%20Effects%20on%20Cosmic%20Ray%20Flux%0A%0A%20%20%20%20Galactic%20cosmic%20rays%20are%20decelerated%20by%20the%20solar%20wind%20before%20reaching%20Earth.%20GSF%20uses%20the%20force-field%20approximation%2C%20parameterized%20by%20a%20single%20potential%20%24%5Cphi%24%20(GV)%3A%20higher%20%24%5Cphi%24%20%E2%86%92%20stronger%20modulation%2C%20lower%20observed%20flux.%20Typical%20values%3A%20%24%5Cphi%20%5Capprox%200.4%24%20GV%20at%20solar%20minimum%2C%20%24%5Cphi%20%5Capprox%201.0%24%20GV%20at%20solar%20maximum.%20The%20unmodulated%20Local%20Interstellar%20Spectrum%20(LIS)%20corresponds%20to%20%24%5Cphi%20%3D%200%24.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Aasync%20def%20_()%3A%0A%20%20%20%20import%20sys%0A%0A%20%20%20%20NOTEBOOK%20%3D%20%22solar_modulation%22%20%20%23%20used%20by%20the%20download%20buttons%20at%20the%20end%0A%20%20%20%20IN_WASM%20%3D%20%22pyodide%22%20in%20sys.modules%0A%20%20%20%20SITE%20%3D%20%22%22%0A%0A%20%20%20%20if%20IN_WASM%3A%0A%20%20%20%20%20%20%20%20%23%20WASM%20(docs%20gallery)%3A%20TEMPORARY%20until%20globalsplinefit%20is%20on%20PyPI%20%E2%80%94%0A%20%20%20%20%20%20%20%20%23%20install%20the%20wheel%20published%20under%20%2Fwheels%2F%20on%20the%20same%20site%20the%0A%20%20%20%20%20%20%20%20%23%20notebook%20is%20served%20from%20(also%20works%20on%20staging%20or%20local%20hosts).%0A%20%20%20%20%20%20%20%20from%20urllib.parse%20import%20urljoin%0A%0A%20%20%20%20%20%20%20%20import%20marimo%20as%20_mo%0A%20%20%20%20%20%20%20%20import%20micropip%0A%0A%20%20%20%20%20%20%20%20%23%20notebook_location()%20resolution%20depth%20differs%20between%20marimo%20export%0A%20%20%20%20%20%20%20%20%23%20layouts%20%E2%80%94%20try%20site-root%20%2Fwheels%2F%20from%20every%20plausible%20depth%2C%20and%0A%20%20%20%20%20%20%20%20%23%20keep%20the%20depth%20that%20worked%20as%20the%20site%20root%20for%20other%20assets.%0A%20%20%20%20%20%20%20%20_base%20%3D%20str(_mo.notebook_location())%20%2B%20%22%2F%22%0A%20%20%20%20%20%20%20%20_err%20%3D%20None%0A%20%20%20%20%20%20%20%20for%20_up%20in%20(%22%22%2C%20%22..%2F%22%2C%20%22..%2F..%2F%22%2C%20%22..%2F..%2F..%2F%22)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20await%20micropip.install(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20urljoin(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_base%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_up%20%2B%20%22wheels%2Fglobalsplinefit-2.0.0a1-py3-none-any.whl%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20SITE%20%3D%20urljoin(_base%2C%20_up)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20None%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20break%0A%20%20%20%20%20%20%20%20%20%20%20%20except%20Exception%20as%20_e%3A%20%20%23%20noqa%3A%20BLE001%20%E2%80%94%20404%20lands%20as%20generic%20error%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_err%20%3D%20_e%0A%20%20%20%20%20%20%20%20if%20_err%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%20_err%0A%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20from%20matplotlib.colors%20import%20Normalize%0A%0A%20%20%20%20from%20globalsplinefit%20import%20(%0A%20%20%20%20%20%20%20%20MODEL_VERSIONS%2C%0A%20%20%20%20%20%20%20%20GSFEnergy%2C%0A%20%20%20%20%20%20%20%20GSFEnergyPerNucleon%2C%0A%20%20%20%20%20%20%20%20GSFRigidity%2C%0A%20%20%20%20)%0A%0A%20%20%20%20plt.rcParams%5B%22figure.figsize%22%5D%20%3D%20(10%2C%206)%0A%0A%20%20%20%20def%20show(fig%3DNone)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Stamp%20and%20return%20a%20figure%20as%20the%20cell%20output.%0A%0A%20%20%20%20%20%20%20%20marimo%20renders%20a%20cell's%20last%20expression%3B%20%60plt.show()%60%20produces%20no%0A%20%20%20%20%20%20%20%20output%20at%20all%20in%20app%20view%2C%20so%20every%20plot%20cell%20ends%20with%20%60show(fig)%60.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20fig%20%3D%20plt.gcf()%20if%20fig%20is%20None%20else%20fig%0A%20%20%20%20%20%20%20%20%23%20TEMPORARY%20pre-publication%20marker%20(_gsf_preliminary)%3A%20diagonal%0A%20%20%20%20%20%20%20%20%23%20PRELIMINARY%20watermark%20on%20every%20figure.%20Delete%20this%20block%20at%20the%0A%20%20%20%20%20%20%20%20%23%20GSF%202026%20release%20%E2%80%94%20keep%20the%20%60return%20fig%60.%0A%20%20%20%20%20%20%20%20fig.text(%0A%20%20%20%20%20%20%20%20%20%20%20%200.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%200.5%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22PRELIMINARY%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ha%3D%22center%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20va%3D%22center%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20rotation%3D30%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D42%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontweight%3D%22bold%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22gray%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.18%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D1000%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20return%20fig%0A%0A%20%20%20%20def%20autoscale(ax%2C%20*series%2C%20pad%3D1.6%2C%20lift%3D1.0)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Fit%20the%20y-axis%20to%20the%20positive%2C%20finite%20values%20in%20%60series%60.%0A%0A%20%20%20%20%20%20%20%20%60lift%60%20raises%20the%20lower%20limit%20by%20that%20factor%20%E2%80%94%20the%20decades%20a%20steeply%0A%20%20%20%20%20%20%20%20falling%20group%20spends%20near%20the%20bottom%20carry%20no%20information%20and%20only%0A%20%20%20%20%20%20%20%20compress%20the%20part%20of%20the%20plot%20that%20does.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20_v%20%3D%20np.concatenate(%5Bnp.asarray(s%2C%20float).ravel()%20for%20s%20in%20series%5D)%0A%20%20%20%20%20%20%20%20_v%20%3D%20_v%5Bnp.isfinite(_v)%20%26%20(_v%20%3E%200)%5D%0A%20%20%20%20%20%20%20%20if%20_v.size%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20ax.set_ylim(_v.min()%20%2F%20pad%20*%20lift%2C%20_v.max()%20*%20pad)%0A%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20GSFEnergy%2C%0A%20%20%20%20%20%20%20%20GSFEnergyPerNucleon%2C%0A%20%20%20%20%20%20%20%20GSFRigidity%2C%0A%20%20%20%20%20%20%20%20MODEL_VERSIONS%2C%0A%20%20%20%20%20%20%20%20NOTEBOOK%2C%0A%20%20%20%20%20%20%20%20Normalize%2C%0A%20%20%20%20%20%20%20%20SITE%2C%0A%20%20%20%20%20%20%20%20autoscale%2C%0A%20%20%20%20%20%20%20%20np%2C%0A%20%20%20%20%20%20%20%20plt%2C%0A%20%20%20%20%20%20%20%20show%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(MODEL_VERSIONS%2C%20mo)%3A%0A%20%20%20%20version%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dlist(MODEL_VERSIONS)%2C%20value%3D%222026.0%22%2C%20label%3D%22model%20version%22%0A%20%20%20%20)%0A%20%20%20%20exponent%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20start%3D0.0%2C%0A%20%20%20%20%20%20%20%20stop%3D3.0%2C%0A%20%20%20%20%20%20%20%20step%3D0.1%2C%0A%20%20%20%20%20%20%20%20value%3D2.6%2C%0A%20%20%20%20%20%20%20%20label%3D%22flux%20scaling%20%20E%5Ek%22%2C%0A%20%20%20%20%20%20%20%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20plot_size%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20start%3D0.6%2C%0A%20%20%20%20%20%20%20%20stop%3D1.6%2C%0A%20%20%20%20%20%20%20%20step%3D0.1%2C%0A%20%20%20%20%20%20%20%20value%3D1.0%2C%0A%20%20%20%20%20%20%20%20label%3D%22plot%20size%22%2C%0A%20%20%20%20%20%20%20%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20mo.hstack(%5Bversion%2C%20exponent%2C%20plot_size%5D%2C%20justify%3D%22start%22%2C%20gap%3D2)%0A%20%20%20%20return%20exponent%2C%20plot_size%2C%20version%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo%2C%20plot_size)%3A%0A%20%20%20%20%23%20Drives%20--gsf-plot-scale%2C%20which%20gallery.css%20multiplies%20into%20every%20figure's%0A%20%20%20%20%23%20width.%20Scaling%20the%20rendered%20raster%20(rather%20than%20the%20figsize)%20is%20what%20keeps%0A%20%20%20%20%23%20it%20proportional%3A%20fonts%2C%20line%20widths%20and%20the%20aspect%20ratio%20all%20ride%20along.%0A%20%20%20%20mo.Html(f%22%3Cstyle%3E%3Aroot%7B%7B--gsf-plot-scale%3A%7Bplot_size.value%7D%7D%7D%3C%2Fstyle%3E%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(exponent)%3A%0A%20%20%20%20energy_exponent%20%3D%20exponent.value%20%20%23%20display%20scaling%2C%20driven%20by%20the%20slider%0A%20%20%20%20return%20(energy_exponent%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Three%20model%20types%2C%20three%20solar%20conditions%0A%0A%20%20%20%20Solar%20conditions%3A%20LIS%2C%20solar%20minimum%20(2009)%2C%20solar%20maximum%20(1991).%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(GSFEnergy%2C%20GSFEnergyPerNucleon%2C%20GSFRigidity%2C%20np%2C%20version)%3A%0A%20%20%20%20grid%20%3D%20np.logspace(np.log10(0.5)%2C%20np.log10(200)%2C%20200)%20%20%23%20GeV%20or%20GV%0A%20%20%20%20energy%20%3D%20rigidity%20%3D%20energy_per_nucleon%20%3D%20grid%0A%0A%20%20%20%20gsf_energy%20%3D%20GSFEnergy(version%3Dversion.value)%0A%20%20%20%20gsf_rigidity%20%3D%20GSFRigidity(version%3Dversion.value)%0A%20%20%20%20gsf_nucleon%20%3D%20GSFEnergyPerNucleon(version%3Dversion.value)%0A%0A%20%20%20%20%23%20time_interval%20ends%20are%20EXCLUSIVE%2C%20so%20a%20calendar%20year%20runs%20to%0A%20%20%20%20%23%20January%20of%20the%20next%20one%0A%20%20%20%20time_periods%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22LIS%22%3A%20%22LIS%22%2C%0A%20%20%20%20%20%20%20%20%22Solar%20Minimum%20(2009)%22%3A%20(200901%2C%20201001)%2C%0A%20%20%20%20%20%20%20%20%22Solar%20Maximum%20(1991)%22%3A%20(199101%2C%20199201)%2C%0A%20%20%20%20%7D%0A%20%20%20%20groups%20%3D%20%5B%22p%22%2C%20%22He%22%2C%20%22O*%22%2C%20%22Fe*%22%5D%0A%20%20%20%20group_colors%20%3D%20%5B%22red%22%2C%20%22orange%22%2C%20%22green%22%2C%20%22blue%22%5D%0A%20%20%20%20line_styles%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22LIS%22%3A%20%22-%22%2C%0A%20%20%20%20%20%20%20%20%22Solar%20Minimum%20(2009)%22%3A%20%22--%22%2C%0A%20%20%20%20%20%20%20%20%22Solar%20Maximum%20(1991)%22%3A%20%22%3A%22%2C%0A%20%20%20%20%7D%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20group_colors%2C%0A%20%20%20%20%20%20%20%20groups%2C%0A%20%20%20%20%20%20%20%20gsf_energy%2C%0A%20%20%20%20%20%20%20%20gsf_nucleon%2C%0A%20%20%20%20%20%20%20%20gsf_rigidity%2C%0A%20%20%20%20%20%20%20%20line_styles%2C%0A%20%20%20%20%20%20%20%20rigidity%2C%0A%20%20%20%20%20%20%20%20time_periods%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Flux%20comparisons%0A%0A%20%20%20%20Compute%20and%20plot%20per-group%20%2B%20total%20flux%20for%20each%20model%20and%20time%20period.%20Fluxes%20are%20scaled%20by%20%24E%5E%7B2.6%7D%24%20to%20flatten%20the%20spectrum.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20energy%2C%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20groups%2C%0A%20%20%20%20gsf_energy%2C%0A%20%20%20%20gsf_nucleon%2C%0A%20%20%20%20gsf_rigidity%2C%0A%20%20%20%20rigidity%2C%0A%20%20%20%20time_periods%2C%0A)%3A%0A%20%20%20%20%23%20Collect%20fluxes%20for%20each%20(model%2C%20period%2C%20group)%20combination%0A%20%20%20%20models%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22GSFEnergy%22%3A%20(gsf_energy%2C%20energy)%2C%0A%20%20%20%20%20%20%20%20%22GSFRigidity%22%3A%20(gsf_rigidity%2C%20rigidity)%2C%0A%20%20%20%20%20%20%20%20%22GSFEnergyPerNucleon%22%3A%20(gsf_nucleon%2C%20energy_per_nucleon)%2C%0A%20%20%20%20%7D%0A%20%20%20%20fluxes%20%3D%20%7B%7D%0A%20%20%20%20for%20_mname%2C%20(model%2C%20_grid_)%20in%20models.items()%3A%0A%20%20%20%20%20%20%20%20fluxes%5B_mname%5D%20%3D%20%7B%7D%0A%20%20%20%20%20%20%20%20for%20_period%2C%20ti%20in%20time_periods.items()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20fluxes%5B_mname%5D%5B_period%5D%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20g%3A%20model.flux(_grid_%2C%20g%2C%20time_interval%3Dti)%20for%20g%20in%20groups%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20fluxes%5B_mname%5D%5B_period%5D%5B%22Total%22%5D%20%3D%20model.total_flux(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_grid_%2C%20time_interval%3Dti%0A%20%20%20%20%20%20%20%20%20%20%20%20)%20%20%23%20fluxes%5Bmodel_name%5D%5Bperiod%5D%5Bgroup_or_%22Total%22%5D%0A%20%20%20%20return%20(fluxes%2C)%0A%0A%0A%40app.cell%0Adef%20_(energy_exponent%2C%20fluxes%2C%20group_colors%2C%20groups%2C%20line_styles)%3A%0A%20%20%20%20def%20plot_groups(ax%2C%20model_name%2C%20grid_%2C%20xlabel%2C%20scale_label)%3A%0A%20%20%20%20%20%20%20%20period_handles%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20i%2C%20(group%2C%20_color)%20in%20enumerate(zip(groups%2C%20group_colors%2C%20strict%3DFalse))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_period%2C%20style%20in%20line_styles.items()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_flux%20%3D%20fluxes%5Bmodel_name%5D%5B_period%5D%5Bgroup%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(line%2C)%20%3D%20ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20grid_%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_flux%20*%20grid_**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dstyle%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20i%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20period_handles.append((line%2C%20_period))%0A%20%20%20%20%20%20%20%20group_handles%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20ax.plot(%5B%5D%2C%20%5B%5D%2C%20color%3Dc%2C%20lw%3D3%2C%20label%3Dlbl)%5B0%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20c%2C%20lbl%20in%20zip(group_colors%2C%20%5B%22H*%22%2C%20%22He*%22%2C%20%22O*%22%2C%20%22Fe*%22%5D%2C%20strict%3DFalse)%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20ax.set_xlabel(xlabel)%0A%20%20%20%20%20%20%20%20ax.set_ylabel(scale_label)%0A%20%20%20%20%20%20%20%20ax.set_title(model_name)%0A%20%20%20%20%20%20%20%20ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20%20%20%20%20ax.set_xlim(grid_%5B0%5D%2C%20grid_%5B-1%5D)%0A%20%20%20%20%20%20%20%20period_legend%20%3D%20ax.legend(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5Bh%20for%20h%2C%20_%20in%20period_handles%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5Blbl%20for%20_%2C%20lbl%20in%20period_handles%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20loc%3D%22lower%20right%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20ax.legend(handles%3Dgroup_handles%2C%20loc%3D%22upper%20right%22%2C%20ncol%3D4)%0A%20%20%20%20%20%20%20%20ax.add_artist(period_legend)%0A%0A%20%20%20%20return%20(plot_groups%2C)%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20autoscale%2C%0A%20%20%20%20energy%2C%0A%20%20%20%20energy_exponent%2C%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20plot_groups%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20rigidity%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20axes%20%3D%20plt.subplots(3%2C%201%2C%20figsize%3D(9%2C%2012))%0A%20%20%20%20plot_groups(%0A%20%20%20%20%20%20%20%20axes%5B0%5D%2C%0A%20%20%20%20%20%20%20%20%22GSFEnergy%22%2C%0A%20%20%20%20%20%20%20%20energy%2C%0A%20%20%20%20%20%20%20%20%22Energy%20%5BGeV%5D%22%2C%0A%20%20%20%20%20%20%20%20f%22Flux%20%C3%97%20%24E%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22%2C%0A%20%20%20%20)%0A%20%20%20%20plot_groups(%0A%20%20%20%20%20%20%20%20axes%5B1%5D%2C%0A%20%20%20%20%20%20%20%20%22GSFRigidity%22%2C%0A%20%20%20%20%20%20%20%20rigidity%2C%0A%20%20%20%20%20%20%20%20%22Rigidity%20%5BGV%5D%22%2C%0A%20%20%20%20%20%20%20%20f%22Flux%20%C3%97%20%24R%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22%2C%0A%20%20%20%20)%0A%20%20%20%20plot_groups(%0A%20%20%20%20%20%20%20%20axes%5B2%5D%2C%0A%20%20%20%20%20%20%20%20%22GSFEnergyPerNucleon%22%2C%0A%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%22Energy%2Fnucleon%20%5BGeV%2Fnucleon%5D%22%2C%0A%20%20%20%20%20%20%20%20f%22Nucleon%20flux%20%C3%97%20%24E%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22%2C%0A%20%20%20%20)%0A%20%20%20%20for%20_ax%20in%20axes%3A%0A%20%20%20%20%20%20%20%20autoscale(_ax%2C%20*%5B_l.get_ydata()%20for%20_l%20in%20_ax.get_lines()%5D%2C%20lift%3D10)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20Total%20flux%20across%20model%20types%0A%0A%20%20%20%20Overlay%20the%20total%20flux%20from%20all%20three%20models%20on%20a%20single%20axis.%20Note%20the%20x-axis%20units%20differ%20by%20model%20(GeV%20vs%20GV%20vs%20GeV%2Fnucleon).%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20autoscale%2C%0A%20%20%20%20energy%2C%0A%20%20%20%20energy_exponent%2C%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20fluxes%2C%0A%20%20%20%20line_styles%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20rigidity%2C%0A%20%20%20%20show%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3D(10%2C%206))%0A%20%20%20%20plot_models%20%3D%20%5B%0A%20%20%20%20%20%20%20%20(%22GSFEnergy%22%2C%20energy%2C%20%22black%22)%2C%0A%20%20%20%20%20%20%20%20(%22GSFRigidity%22%2C%20rigidity%2C%20%22purple%22)%2C%0A%20%20%20%20%20%20%20%20(%22GSFEnergyPerNucleon%22%2C%20energy_per_nucleon%2C%20%22brown%22)%2C%0A%20%20%20%20%5D%0A%20%20%20%20for%20_mname%2C%20_grid_%2C%20_color%20in%20plot_models%3A%0A%20%20%20%20%20%20%20%20for%20_period%2C%20style%20in%20line_styles.items()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_flux%20%3D%20fluxes%5B_mname%5D%5B_period%5D%5B%22Total%22%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20label%20%3D%20_mname%20if%20_period%20%3D%3D%20%22LIS%22%20else%20None%0A%20%20%20%20%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_grid_%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_flux%20*%20_grid_**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ls%3Dstyle%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3Dlabel%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.8%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20%2F%20Rigidity%20%2F%20Energy%20per%20nucleon%22)%0A%20%20%20%20_ax.set_ylabel(f%22Total%20flux%20%C3%97%20%24X%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%20(units%20vary%20by%20model)%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.legend(loc%3D%22lower%20right%22)%0A%20%20%20%20autoscale(_ax%2C%20*%5B_l.get_ydata()%20for%20_l%20in%20_ax.get_lines()%5D%2C%20lift%3D10)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Solar%20cycle%20evolution%0A%0A%20%20%20%20Sample%2017%20monthly%20time%20points%20across%202013%E2%80%932024%20and%20compute%20the%20total%20nucleon%20flux%20at%20each.%20The%20color%20scale%20maps%20to%20%24%5Cphi%24%20(warmer%20%3D%20stronger%20modulation).%20The%20second%20figure%20shows%20%24%5Cphi%24%20and%20the%20flux-suppression%20ratio%20(modulated%20%2F%20LIS)%20at%201.5%2C%2010%2C%20and%20100%20GeV%2Fnucleon.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(gsf_energy%2C%20np)%3A%0A%20%20%20%20%23%2017%20points%20spanning%202013-01%20to%202024-12.%20The%20end%20of%20a%20time_interval%20is%0A%20%20%20%20%23%20exclusive%2C%20so%20a%20December%20point%20ends%20in%20January%20of%20the%20next%20year.%0A%20%20%20%20time_points%20%3D%20%5B%5D%0A%20%20%20%20phi_values%20%3D%20%5B%5D%0A%20%20%20%20dates%20%3D%20%5B%5D%0A%20%20%20%20for%20_idx%20in%20np.linspace(0%2C%20(2024%20-%202013%20%2B%201)%20*%2012%20-%201%2C%2017)%3A%0A%20%20%20%20%20%20%20%20month_idx%20%3D%20int(round(_idx))%0A%20%20%20%20%20%20%20%20year%20%3D%202013%20%2B%20month_idx%20%2F%2F%2012%0A%20%20%20%20%20%20%20%20month%20%3D%20(month_idx%20%25%2012)%20%2B%201%0A%20%20%20%20%20%20%20%20if%20year%20in%20gsf_energy.phi%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20end_year%2C%20end_month%20%3D%20(year%20%2B%201%2C%201)%20if%20month%20%3D%3D%2012%20else%20(year%2C%20month%20%2B%201)%0A%20%20%20%20%20%20%20%20%20%20%20%20time_points.append((year%20*%20100%20%2B%20month%2C%20end_year%20*%20100%20%2B%20end_month))%0A%20%20%20%20%20%20%20%20%20%20%20%20phi_values.append(float(gsf_energy.phi%5Byear%5D%5Bmonth%20-%201%5D))%0A%20%20%20%20%20%20%20%20%20%20%20%20dates.append(f%22%7Byear%7D-%7Bmonth%3A02d%7D%22)%0A%0A%20%20%20%20print(%0A%20%20%20%20%20%20%20%20f%22%7Blen(time_points)%7D%20time%20points%2C%20%CF%86%20range%3A%20%7Bmin(phi_values)%3A.3f%7D%E2%80%93%7Bmax(phi_values)%3A.3f%7D%20GV%22%0A%20%20%20%20)%0A%20%20%20%20return%20dates%2C%20phi_values%2C%20time_points%0A%0A%0A%40app.cell%0Adef%20_(energy_per_nucleon%2C%20gsf_nucleon%2C%20np%2C%20time_points)%3A%0A%20%20%20%20time_series_fluxes%20%3D%20np.array(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20gsf_nucleon.total_flux(energy_per_nucleon%2C%20time_interval%3Dti)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20ti%20in%20time_points%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20flux_lis%20%3D%20gsf_nucleon.total_flux(energy_per_nucleon%2C%20time_interval%3D%22LIS%22)%0A%20%20%20%20return%20flux_lis%2C%20time_series_fluxes%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20Normalize%2C%0A%20%20%20%20autoscale%2C%0A%20%20%20%20dates%2C%0A%20%20%20%20energy_exponent%2C%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20flux_lis%2C%0A%20%20%20%20phi_values%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A%20%20%20%20time_series_fluxes%2C%0A)%3A%0A%20%20%20%20cmap%20%3D%20plt.get_cmap(%22coolwarm%22)%0A%20%20%20%20norm%20%3D%20Normalize(vmin%3Dmin(phi_values)%2C%20vmax%3Dmax(phi_values))%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3D(8%2C%205))%0A%20%20%20%20for%20_flux%2C%20phi%2C%20_date%20in%20zip(time_series_fluxes%2C%20phi_values%2C%20dates%2C%20strict%3DFalse)%3A%0A%20%20%20%20%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_flux%20*%20energy_per_nucleon**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3Dcmap(norm(phi))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.7%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D1.5%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax.loglog(%0A%20%20%20%20%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20%20%20%20%20flux_lis%20*%20energy_per_nucleon**energy_exponent%2C%0A%20%20%20%20%20%20%20%20%22k-%22%2C%0A%20%20%20%20%20%20%20%20lw%3D3%2C%0A%20%20%20%20%20%20%20%20label%3D%22LIS%20(%CF%86%3D0)%22%2C%0A%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Energy%20%5BGeV%2Fnucleon%5D%22)%0A%20%20%20%20_ax.set_ylabel(f%22Nucleon%20flux%20%C3%97%20%24E%5E%7B%7B%7Benergy_exponent%7D%7D%7D%24%22)%0A%20%20%20%20_ax.set_title(%22Nucleon%20flux%20through%20the%20solar%20cycle%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax.set_xlim(energy_per_nucleon%5B0%5D%2C%20energy_per_nucleon%5B-1%5D)%0A%20%20%20%20autoscale(_ax%2C%20flux_lis%20*%20energy_per_nucleon**energy_exponent)%0A%20%20%20%20_ax.legend(loc%3D%22lower%20right%22)%0A%20%20%20%20sm%20%3D%20plt.cm.ScalarMappable(cmap%3Dcmap%2C%20norm%3Dnorm)%0A%20%20%20%20sm.set_array(%5B%5D)%0A%20%20%20%20plt.colorbar(sm%2C%20ax%3D_ax%2C%20label%3D%22%CF%86%20%5BGV%5D%22)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20dates%2C%0A%20%20%20%20energy_per_nucleon%2C%0A%20%20%20%20flux_lis%2C%0A%20%20%20%20np%2C%0A%20%20%20%20phi_values%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20show%2C%0A%20%20%20%20time_series_fluxes%2C%0A)%3A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3D(9%2C%205))%0A%20%20%20%20_ax2%20%3D%20_ax.twinx()%0A%20%20%20%20years%20%3D%20%5Bint(d%5B%3A4%5D)%20%2B%20(int(d%5B5%3A7%5D)%20-%201)%20%2F%2012%20for%20d%20in%20dates%5D%0A%20%20%20%20_ax.plot(years%2C%20phi_values%2C%20%22ro-%22%2C%20lw%3D2%2C%20ms%3D6%2C%20label%3D%22%CF%86%22)%0A%20%20%20%20_ax.set_xlabel(%22Year%22)%0A%20%20%20%20_ax.set_ylabel(%22%CF%86%20%5BGV%5D%22%2C%20color%3D%22red%22)%0A%20%20%20%20_ax.tick_params(axis%3D%22y%22%2C%20labelcolor%3D%22red%22)%0A%20%20%20%20_ax.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20for%20energy_val%2C%20_color%20in%20zip(%0A%20%20%20%20%20%20%20%20%5B1.5%2C%2010.0%2C%20100.0%5D%2C%20%5B%22blue%22%2C%20%22green%22%2C%20%22orange%22%5D%2C%20strict%3DFalse%0A%20%20%20%20)%3A%0A%20%20%20%20%20%20%20%20idx%20%3D%20np.argmin(np.abs(energy_per_nucleon%20-%20energy_val))%0A%20%20%20%20%20%20%20%20ratio%20%3D%20time_series_fluxes%5B%3A%2C%20idx%5D%20%2F%20flux_lis%5Bidx%5D%0A%20%20%20%20%20%20%20%20_ax2.plot(%0A%20%20%20%20%20%20%20%20%20%20%20%20years%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ratio%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22o-%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D_color%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ms%3D4%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3Df%22%7Benergy_per_nucleon%5Bidx%5D%3A.1f%7D%20GeV%2Fnuc%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax2.axhline(1.0%2C%20color%3D%22gray%22%2C%20ls%3D%22--%22%2C%20alpha%3D0.5)%0A%20%20%20%20_ax2.set_ylabel(%22Flux%20ratio%20(modulated%20%2F%20LIS)%22)%0A%20%20%20%20_ax2.set_ylim(0%2C%201.2)%0A%20%20%20%20_ax2.legend(loc%3D%22lower%20right%22)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Solar-cycle%20averaging%3A%20how%20many%20%24%5Cphi%24%20bins%3F%0A%0A%20%20%20%20%60time_interval%3DNone%60%20returns%20a%20Solar%20Cycle%2024%20average.%20The%20monthly%20%24%5Cphi%24%20values%0A%20%20%20%20of%20the%20interval%20are%20binned%20into%20%60solar_cycle_average_bins%60%20weighted%0A%20%20%20%20representatives%20(each%20bin's%20mean%20%24%5Cphi%24%2C%20weighted%20by%20its%20month%20count)%20and%20the%0A%20%20%20%20*modulated%20flux*%20is%20averaged%20over%20them%20%E2%80%94%20the%20same%20period%20average%20the%20fitter%0A%20%20%20%20uses.%20Because%20flux%20is%20nonlinear%20in%20%24%5Cphi%24%2C%20this%20%22mean%20of%20modulated%22%20is%20not%0A%20%20%20%20%24%5Cmathrm%7BFlux%7D(%5Clangle%5Cphi%5Crangle)%24.%0A%0A%20%20%20%20The%20weighted%20mean%20%24%5Cphi%24%20is%20exact%20at%20every%20setting%2C%20so%20the%20first-order%0A%20%20%20%20suppression%20is%20always%20right%20and%20only%20the%20curvature%20is%20approximated.%20At%201%20GeV%2C%0A%20%20%20%201%20bin%20is%20off%20by%205.2%25%20from%20the%20full%20monthly%20average%2C%206%20bins%20by%200.5%25%2C%20and%20the%20default%2012%20by%200.1%25%3B%20all%20are%20below%200.3%25%20above%0A%20%20%20%2010%20GeV.%20%60None%60%20averages%20every%20month%20explicitly.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(GSFEnergy%2C%20energy_exponent%2C%20np%2C%20plt%2C%20show)%3A%0A%20%20%20%20gsf_one%20%3D%20GSFEnergy(solar_cycle_average_bins%3D1)%0A%20%20%20%20gsf_def%20%3D%20GSFEnergy(solar_cycle_average_bins%3D12)%20%20%23%20the%20default%0A%20%20%20%20gsf_exact%20%3D%20GSFEnergy(solar_cycle_average_bins%3DNone)%0A%20%20%20%20energy_subset%20%3D%20np.logspace(0%2C%201%2C%20500)%0A%0A%20%20%20%20flux_one%20%3D%20gsf_one.total_flux(energy_subset)%0A%20%20%20%20flux_def%20%3D%20gsf_def.total_flux(energy_subset)%0A%20%20%20%20flux_exact%20%3D%20gsf_exact.total_flux(energy_subset)%0A%20%20%20%20rel_one%20%3D%20100%20*%20np.abs(flux_one%20-%20flux_exact)%20%2F%20flux_exact%0A%20%20%20%20rel_def%20%3D%20100%20*%20np.abs(flux_def%20-%20flux_exact)%20%2F%20flux_exact%0A%0A%20%20%20%20_fig%2C%20(ax1%2C%20_ax2)%20%3D%20plt.subplots(2%2C%201%2C%20sharex%3DTrue%2C%20figsize%3D(7%2C%205))%0A%20%20%20%20ax1.loglog(%0A%20%20%20%20%20%20%20%20energy_subset%2C%0A%20%20%20%20%20%20%20%20energy_subset**energy_exponent%20*%20flux_one%2C%0A%20%20%20%20%20%20%20%20%22r-%22%2C%0A%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20label%3D%221%20bin%3A%20Flux%24(%5C%5Clangle%5C%5Cphi%5C%5Crangle)%24%22%2C%0A%20%20%20%20)%0A%20%20%20%20ax1.loglog(%0A%20%20%20%20%20%20%20%20energy_subset%2C%0A%20%20%20%20%20%20%20%20energy_subset**energy_exponent%20*%20flux_def%2C%0A%20%20%20%20%20%20%20%20%22-%22%2C%0A%20%20%20%20%20%20%20%20color%3D%22C2%22%2C%0A%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20label%3D%2212%20bins%20(default)%22%2C%0A%20%20%20%20)%0A%20%20%20%20ax1.loglog(%0A%20%20%20%20%20%20%20%20energy_subset%2C%0A%20%20%20%20%20%20%20%20energy_subset**energy_exponent%20*%20flux_exact%2C%0A%20%20%20%20%20%20%20%20%22b--%22%2C%0A%20%20%20%20%20%20%20%20lw%3D2%2C%0A%20%20%20%20%20%20%20%20label%3D%22all%20months%3A%20%24%5C%5Clangle%24Flux%24(%5C%5Cphi)%5C%5Crangle%24%22%2C%0A%20%20%20%20)%0A%20%20%20%20ax1.set_ylabel(%22Total%20flux%20%C3%97%20%24E%5E%7B2.6%7D%24%22)%0A%20%20%20%20ax1.legend()%0A%20%20%20%20ax1.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax2.semilogx(energy_subset%2C%20rel_one%2C%20%22r-%22%2C%20lw%3D2%2C%20label%3D%221%20bin%22)%0A%20%20%20%20_ax2.semilogx(energy_subset%2C%20rel_def%2C%20%22-%22%2C%20color%3D%22C2%22%2C%20lw%3D2%2C%20label%3D%2212%20bins%22)%0A%20%20%20%20_ax2.set_xlabel(%22Energy%20%5BGeV%5D%22)%0A%20%20%20%20_ax2.set_ylabel(%22Deviation%20from%20full%20average%20%5B%25%5D%22)%0A%20%20%20%20_ax2.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20_ax2.legend(fontsize%3D8%2C%20loc%3D%22upper%20right%22)%0A%20%20%20%20_ax2.set_ylim(0%2C%20rel_one.max()%20*%201.1)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20show()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(NOTEBOOK%2C%20SITE%2C%20mo)%3A%0A%20%20%20%20if%20SITE%3A%20%20%23%20running%20in%20the%20browser%3A%20link%20the%20sources%20published%20next%20to%20it%0A%20%20%20%20%20%20%20%20_py%20%3D%20f%22%7BSITE%7Dgallery%2F%7BNOTEBOOK%7D%2F%7BNOTEBOOK%7D.py%22%0A%20%20%20%20%20%20%20%20_ipynb%20%3D%20f%22%7BSITE%7Dgallery%2F%7BNOTEBOOK%7D%2F%7BNOTEBOOK%7D.ipynb%22%0A%20%20%20%20else%3A%20%20%23%20running%20locally%3A%20hand%20over%20the%20file%20on%20disk%0A%20%20%20%20%20%20%20%20_py%20%3D%20(mo.notebook_dir()%20%2F%20f%22%7BNOTEBOOK%7D.py%22).read_bytes()%0A%20%20%20%20%20%20%20%20_ipynb%20%3D%20None%0A%0A%20%20%20%20_buttons%20%3D%20%5B%0A%20%20%20%20%20%20%20%20mo.download(%0A%20%20%20%20%20%20%20%20%20%20%20%20_py%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20filename%3Df%22%7BNOTEBOOK%7D.py%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22marimo%20notebook%20(.py)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mimetype%3D%22text%2Fx-python%22%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%5D%0A%20%20%20%20if%20_ipynb%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20_buttons.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.download(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_ipynb%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20filename%3Df%22%7BNOTEBOOK%7D.ipynb%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22Jupyter%20notebook%20(.ipynb)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mimetype%3D%22application%2Fx-ipynb%2Bjson%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(r%22%22%22%23%23%23%20Take%20this%20notebook%20with%20you%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(_buttons%2C%20justify%3D%22start%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A