Skip to content

Magnetopause script updates#326

Closed
jreimi wants to merge 22 commits intofmihpc:devfrom
jreimi:streamline_magnetopause
Closed

Magnetopause script updates#326
jreimi wants to merge 22 commits intofmihpc:devfrom
jreimi:streamline_magnetopause

Conversation

@jreimi
Copy link
Copy Markdown
Contributor

@jreimi jreimi commented Jun 5, 2025

Updates and documentation for magnetopause finding with streamlines.


#Create streamline seeds (starting points for streamlines)
seedN = 50 #seeds per row, final seed count will be seedN*seedN !
seedN = 25 #seeds per row, final seed count will be seedN*seedN !
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This should be a (default) parameter for the function

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

e.g. def make_streamlines(vlsvFileName, seedN=25):

@@ -185,27 +148,27 @@ def make_streamlines(vlsvFileName):
streamline_seeds[k] = [20, i, j]
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

20 -> another default parameter

Copy link
Copy Markdown
Member

@alhom alhom left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

As these are aimed to be included in the analysator modules, I'd like a top-level callable function like find_magnetopause(vlsvfile, seedN=25, x0=20 Re, ...) (with inputs in SI) that wraps the existing parts of main()s (~below) that actually gather the magnetopause, and return the triangle mesh. Ditto for magnetopause2d (incl. default parameter comments).

def main():
    ## STREAMLINES
    streams = make_streamlines(fileLocation+fileN)
    ## MAGNETOPAUSE
    magnetopause = make_magnetopause(streams)

@jreimi jreimi marked this pull request as ready for review June 13, 2025 11:41
@alhom alhom changed the base branch from master to dev June 16, 2025 07:05
magnetopause2d
--------------
:doc:`magnetopause2d`
magnetopause_sw_streamline_2d
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

These are now documented by autodoc in pyCalculations, rather have the example scripts described here (or rather, similarly under a new examples page!)


------------

magnetopause2d
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Restore these to point towards your plot scripts

@alhom
Copy link
Copy Markdown
Member

alhom commented Aug 29, 2025

Included in #338.

@alhom alhom closed this Aug 29, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants