EM GeoSci Apps
Open MT_LayeredEarth on Binder:
https://mybinder.org/v2/gh/geoscixyz/em-apps/main?filepath=notebooks/em/MT_LayeredEarth.ipynb
EM GeoSci Apps
Open MT_LayeredEarth on Binder:
https://mybinder.org/v2/gh/geoscixyz/em-apps/main?filepath=notebooks/em/MT_LayeredEarth.ipynb
spherical harmonic
See visualisations at https://book.magneticearth.org/practical-field-modelling/01a_making-a-main-field-model
See eoxmagmod documentation (WIP): https://esa-vires.github.io/MagneticModel/
Space Physics book companion... has nice potential fields tool: https://spacephysics.ucla.edu/PotentialField/RotatablePotentialFields.html
Provide instructions for sharing notebooks, e.g. via https://notebooksharing.space/
Beginner: For more educative introductions, see:
Change the F10.7 source to get from GFZ instead(?) https://kp.gfz-potsdam.de/en/data
See: - https://wdc.bgs.ac.uk/ - https://wdcapi.bgs.ac.uk/docs
To do: - more explanation of AUX_OBS vs INTERMAGNET vs WDC vs SuperMAG - compare access to AUX_OBS over VirES and WDC/INTERMAGNET HAPI - Uusing WDC API for metadata (e.g. acknowledgements)
related repositories.
More to add: - https://github.com/heliophysicsPy/summer-school - https://gitlab.com/KNMI-OSS/spaceweather/notebooks
More resources at: https://swarmdisc.org/lab/
Jupyter
CodeRefinery demo: https://www.youtube.com/watch?v=cSHjTA_O8_w
Intermediate: For creating importable modules and packages, testing, documenting, logging, optimising…:
https://merely-useful.tech/py-rse/
https://alan-turing-institute.github.io/rse-course/
https://www.pyopensci.org/python-package-guide/index.html
Reference http://intermagnet.org/new_data_download.html for other access mechanisms
Spherical harmonic models
See https://github.com/klaundal/ppigrf for a better implementation
ds.plot.scatter(x="Timestamp", y="Vixv", ax=axes[0], s=0.1) ds.plot.scatter(x="Timestamp", y="Vixh", ax=axes[1], s=0.1) ds.plot.scatter(x="Timestamp", y="Viy", ax=axes[2], s=0.1) ds.plot.scatter(x="Timestamp", y="Viz", ax=axes[3], s=0.1, label="Velocities")
Change to:
ds.plot.scatter(x="Timestamp", y="Vixv", ax=axes[0], s=1, linewidths=0)
ds.plot.scatter(x="Timestamp", y="Vixh", ax=axes[1], s=1, linewidths=0)
ds.plot.scatter(x="Timestamp", y="Viy", ax=axes[2], s=1, linewidths=0)
ds.plot.scatter(x="Timestamp", y="Viz", ax=axes[3], s=1, linewidths=0, label="Velocities")
Geomagnetic Datacubes Development
Some slides to go with this: https://docs.google.com/presentation/d/1w7unlYWDDqQlJeMxn6v9MRxrOxVveOyP/edit?usp=sharing&ouid=111125506726710589088&rtpof=true&sd=true
Algorithms: see https://doi.org/10.5047/eps.2013.09.006
Evolutions: see https://doi.org/10.1029/2018JA026439
Documentation
See also https://doi.org/10.1029/2021JA030183
(link to dashboard TBD)
Prototype at http://panel.swarmdisc.org/
Swarm Measurements¶ A short demo of taking magnetic measurements from a single Swarm spacecraft and deriving a main field model directly from one month of data!
See also: https://github.com/klaundal/notebooks/tree/master/geomodel
viresclient API¶ Author: Luca Mariani Abstract: Describe the main classes and methods defined in the viresclient package
recommend referring directly to https://viresclient.readthedocs.io here
To-do:
link to https://tutorial.xarray.dev and https://foundations.projectpythia.org/core/xarray.html
Data 8 - UC Berkeley Foundations of Data Science course
TODO: replace link with https://inferentialthinking.com ?
Project Pythia - A Community Educational Resource for Geoscientists
TODO: update link to https://cookbooks.projectpythia.org/
Project Python Foundations
TODO: fix type Python->Pythia
more information
See also: https://nbviewer.org/github/pacesm/jupyter_notebooks/tree/master/EFI_TII/
anywhere on these pages
It works like this :)
TODO: process/migrate notes from old domain: https://hypothes.is/search?q=magneticearth.org
See also:
An alternative way to derive FACs: https://swarm-fac-exploration.magneticearth.org/notebooks/dual_sat_fac
Abstract: Access to the field aligned currents evaluated by the dual satellite method (level 2 product). We also show an orbit-by-orbit plot using a periodic axis to display (centred over both poles) an overview of the FAC estimates over two weeks.
Health warning
Note that the data quality varies from October 2019 onward due to the orbit changes of Swarm A & C:
These changes are yet to be made in the data products as of writing (April 2022)
TODO: add backlinks to Swarm Notebooks
Implement automated testing and generation of notebooks containing output cells using Treebeard - to be viewable at swarm-vre.readthedocs.io
out-of-date text...
Seminar series
Select the analysis interval
Different UI when using Notebook or Lab... not visible in Book. How to reconcile this?
See also, new gradient dataset: https://doi.org/10.11583/DTU.14695590.v2
(link to dashboard TBD)
http://176.58.119.29:5006/02_GVO-Geomagnetic-Virtual-Observatories
Swarm_notebooks/dashboards/04_Conjunctions.ipynb (link TBD)
Beginner
Project Pythia is already getting very good for an introduction to core packages: https://foundations.projectpythia.org/landing-page.html
Lorenz’s (1963) model: forward modelling¶
Binder link for this: https://mybinder.org/v2/gh/MagneticEarth/book.magneticearth.org/simpler_env?urlpath=lab/tree/data_assimilation/01a_lorenz.ipynb
(this has a simpler temporary environment, but won't work with the notebooks under "Geomagnetic Data & Models")
the file oamjumpsapplied.dat in the external data directory.
Sorry I didn't update this! Here it is: https://raw.githubusercontent.com/MagneticEarth/IAGA_SummerSchool2019/master/data/external/oamjumpsapplied.dat
Geomagnetic pulsations¶ .. see https://iaga-school-2019.netlify.app/notebooks/1-visualising-geomagnetic-observatory-data#geomagnetic-pulsations Hartland (HAD) second data are not available. Fetching CLF instead:
Probably just remove this section and create a dedicated page that goes into it
# Using only the radial component, perform the inversion: G = G_radius d = B_radius m = np.linalg.inv(G.T @ G) @ (G.T @ d)
Only using one component to keep it simpler. The radial component is the least contaminated so damping is less important
Estimated core GVO magnetic field, spherical polar (r,theta,phi) vector components [i.e (-C,-N,E) components in NEC frame]
NB: The source (CDF) files provide the measurements in the RTP coordinates rather than NEC. In VirES, we provide them in NEC coordinates for consistency with other Swarm products.
This has led to the misleading variable descriptions in data retrieved from VirES, stating that they are given in "spherical polar (r,theta,phi) vector components [i.e (-C,-N,E) components in NEC frame]". We have not fixed this yet
Satellite missions
Some more: https://imgur.com/bhgHscY
1964-70 Kosmos https://essd.copernicus.org/articles/12/555/2020/
Geomagnetic data
extend with paleomag databases https://imgur.com/a/hDto0s5
See also:<br> https://helioml.org/Introduction/title.html<br> https://earth-env-data-science.github.io/intro
and other examples at:<br> https://executablebooks.org/en/latest/gallery.html
k_inds
is there a nice library for accessing these? are they published data?
ESK_2003_PATH = os.path.abspath('../data/external/ESK_2003/') ESK_HOURLY_PATH = os.path.abspath('../data/external/ESK_hourly/')
will be able to replace with data loading live from vires
# Import notebook dependencies import os import sys import datetime as dt
can investigate hiding code to make the book more friendly to browse https://jupyterbook.org/interactive/hiding.html
>> USER INPUT HERE
I'll have a look for simple "date/time" input boxes through ipywidgets or similar
By David Kerridge et al.
Change to e.g. "IAGA and friends" This is displayed identically on all pages. Perhaps we can add an "authors" element to each notebook (at top or bottom)
prototype quicklook routine
The intention is that in the future we will be able to access such routine from something like:
from viresclient import SwarmQuicklooks
fig = SwarmQuicklooks.AEJxLPS.plot()
documentationn,
I found a typo!
Virtual Research Environment
Too many "virtuals"