[go: nahoru, domu]

Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

increase MPOL of VMEC input by one so that DESC eq has requisite reso… #1005

Draft
wants to merge 12 commits into
base: master
Choose a base branch
from

Conversation

dpanici
Copy link
Collaborator
@dpanici dpanici commented Apr 28, 2024

…lution to capture shafranov shift

without this increase, the auto continuation solve from the command line yields an equilibrium with barely any shafranov shift

Things to check

  • If the continuation method steps are weird or have changed
  • if the non continuation method returns the same or not
  • if it is indeed the MPOL thing (MPOL to MPOL-1 upon loading VMEC was changed about 2 years ago, but 6 months ago we lowered the resolution from 12 to 8... so worth checking the code at that point)
  • if ANSI vs FRINGE matters here

Copy link
codecov bot commented Apr 28, 2024

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 94.97%. Comparing base (4745460) to head (8cc0c75).

Additional details and impacted files
@@            Coverage Diff             @@
##           master    #1005      +/-   ##
==========================================
- Coverage   94.98%   94.97%   -0.01%     
==========================================
  Files          87       87              
  Lines       21743    21743              
==========================================
- Hits        20652    20651       -1     
- Misses       1091     1092       +1     

see 1 file with indirect coverage changes

Copy link
Collaborator
@ddudt ddudt left a comment

Choose a reason for hiding this comment

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

I think the following files also have to be updated to stay self-consistent with this change:

  • input.HELIOTRON_desc_no_continuation
  • input.HELIOTRON_desc_no_continuation_output.h5
  • basic_equilibrium.ipynb
  • advanced_equilibrium_continuation.ipynb

Copy link

Check out this pull request on  ReviewNB

See visual diffs & provide feedback on Jupyter Notebooks.


Powered by ReviewNB

@dpanici
Copy link
Collaborator Author
dpanici commented Apr 29, 2024

I think the following files also have to be updated to stay self-consistent with this change:

  • input.HELIOTRON_desc_no_continuation
  • input.HELIOTRON_desc_no_continuation_output.h5
  • basic_equilibrium.ipynb
  • advanced_equilibrium_continuation.ipynb

the files that need updating have ben updated. The no continuation input does not need updating as the resolution was not changed for that one.

@dpanici dpanici requested a review from ddudt April 29, 2024 18:49
@f0uriest
Copy link
Member

Feels really weird that going from m=8 to m=9 makes such a big difference?

@ddudt
Copy link
Collaborator
ddudt commented Apr 30, 2024

Feels really weird that going from m=8 to m=9 makes such a big difference?

@dpanici are you sure this PR is necessary? The notebooks on master look like they already have the correct Shafranov shift

@dpanici
Copy link
Collaborator Author
dpanici commented Apr 30, 2024

Feels really weird that going from m=8 to m=9 makes such a big difference?

@dpanici are you sure this PR is necessary? The notebooks on master look like they already have the correct Shafranov shift

Yea I am sure. If you rerun the examples on master, you completely lose the shafranov shift. Try it for yourself and see

@dpanici
Copy link
Collaborator Author
dpanici commented Apr 30, 2024

Feels really weird that going from m=8 to m=9 makes such a big difference?

I also agree but this is the result I see, if the point of the notebook is to just use the automatic contiunation that by default is used when running on a vmec input file, it currently results in an equilibrium with barely any shafranov shift on master after running python -m desc input.HELIOTRON

This is fixed if we up the resolution by 1. The current saved input.HELIOTRON_output.h5 on master has L=M=8 resolution, so possibly what happened is it is from before we made the change so that when we load a VMEC input file, we set MPOL to be one less than what VMEC's MPOL is

image

@f0uriest
Copy link
Member
f0uriest commented May 1, 2024

ok, im still a bit worried this is indicative of some larger issue. Can you check to see when it broke? Was it just because we changed the default m from vmec?

@dpanici dpanici marked this pull request as draft May 6, 2024 15:27
@dpanici
Copy link
Collaborator Author
dpanici commented May 6, 2024

Note to self, must make a new environment locally to load in older DESC versions

@dpanici
Copy link
Collaborator Author
dpanici commented Jul 3, 2024

make as an issue then merge this in @dpanici

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.

None yet

4 participants