Skip to content

Add solovev_free_bdy test case input files. #461

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.

Sign up for GitHub

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

Draft
jons-pf wants to merge 1 commit into
base: main
Choose a base branch
from
Draft

Add solovev_free_bdy test case input files. #461

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
38 changes: 38 additions & 0 deletions src/vmecpp/cpp/vmecpp/test_data/coils.solovev
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,38 @@
! Coils for a free-boundary equilibrium similar to the Solov'ev equilibrium in Hirshman & Whitson (1983).
! Also used in the free-boundary SPEC verification paper by Hudson et al (2020).

&MGRID_NLI
MGRID_EXT = 'solovev'
MGRID_MODE = 'R'
LSTELL_SYM = .false.
IR = 201
JZ = 201
KP = 1
Rmin = 2.0
Rmax = 6.0
Zmin = -2.0
Zmax = 2.0
/

** coils_dot_starts_below **
periods 1
begin filament
mirror NIL
0.0 0.0 -1.0e6 1
0.0 0.0 1.0e6 1
1.0e6 0.0 1.0e6 1
1.0e6 0.0 -1.0e6 1
0.0 0.0 -1.0e6 0 1 CS
0.867845 0.0 1.86006 1 2 PF_01
4.49628 0.0 3.69583 1 3 PF_02
8.63810 0.0 -2.46839 1 4 PF_03
8.63847 0.0 2.46836 1 5 PF_03
4.49620 0.0 -3.69589 1 6 PF_02
0.867842 0.0 -1.86005 1 7 PF_01
1.19984 0.0 3.34106 1 8 PF_01
2.62980 0.0 5.24924 1 9 PF_02
9.65488 0.0 -3.29570 1 10 PF_03
9.65501 0.0 3.29541 1 11 PF_03
2.62966 0.0 -5.24925 1 12 PF_02
1.19984 0.0 -3.34102 1 13 PF_01
end
54 changes: 54 additions & 0 deletions src/vmecpp/cpp/vmecpp/test_data/input.solovev_free_bdy
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,54 @@
&INDATA
! SOLOV'EV Equilibrium from Hirshman Whitson 1983 BETA~0.1
LFREEB = T

mgrid_file = 'mgrid_solovev.nc'
EXTCUR( 1) = 3.884526409876309E+06 ! CS

EXTCUR( 2) = -2.935577123737952E+05 ! PF_01
EXTCUR( 7) = -2.935531536308510E+05 ! PF_01

EXTCUR( 8) = -3.560639108717275E+05 ! PF_01
EXTCUR(13) = -3.560589388468855E+05 ! PF_01

EXTCUR( 3) = -1.734851853677043E+04 ! PF_02
EXTCUR( 6) = -1.734993103183817E+04 ! PF_02

EXTCUR( 9) = -6.588434719283084E+04 ! PF_02
EXTCUR(12) = -6.588300858364606E+04 ! PF_02

EXTCUR( 4) = 6.002137016973160E+04 ! PF_03
EXTCUR( 5) = 6.002540940490887E+04 ! PF_03

EXTCUR(10) = -1.154387774712987E+04 ! PF_03
EXTCUR(11) = -1.153546510755219E+04 ! PF_03

NVACSKIP = 6

DELT = 9.00E-01
LASYM = F
NFP = 1
NCURR = 0
MPOL = 6
NTOR = 0
NZETA = 1
NS_ARRAY = 16 32
NITER_ARRAY = 10000 20000
FTOL_ARRAY = 1E-10 1E-14
NSTEP = 100
GAMMA = 0.000000E+00
PHIEDGE = 1.00000000000000E+00
BLOAT = 1.000000E+00
CURTOR = 0.00000000000000E+00
SPRES_PED = 1.00000000000000E+00
PMASS_TYPE = "power_series"
AM = 0.125 -0.125
AI = 0.49 -0.3

RAXIS_CC = 3.999
ZAXIS_CS = 0.000

RBC( 0,0) = 3.999 ZBS( 0,0) = 0.000
RBC( 0,1) = 1.026 ZBS( 0,1) = 1.580
RBC( 0,2) = 0.068 ZBS( 0,2) = -0.010
/