-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathgpec.in
More file actions
88 lines (77 loc) · 6.43 KB
/
gpec.in
File metadata and controls
88 lines (77 loc) · 6.43 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
&GPEC_INPUT
dcon_dir = '' ! Path to DCON outputs. Individual file paths are appended to this
idconfile = 'euler.bin' ! Eigenfunctions output by DCON run
ieqfile = 'psi_in.bin' ! Flux surfaces output by DCON run
rdconfile = 'globalsol.bin' ! Galerkin functions output by resistive DCON run
gal_flag = t ! Use the Galerkin eigenfunctions from resistive DCON instead of the ideal DCON shooting method solutions
jac_in="hamada" ! Coordinate system of external perturbation (affects data_flag and harmonic_flag). Options include: pest,boozer,hamada
jsurf_in=1 ! True(1) if external perturbation is area weighted (i.e. flux)
tmag_in =0 ! True(1) if external perturbation has toroidal angle defined by jac_in. False(0) if machine angle.
mthsurf =0 ! Number of poloidal angle gridpoints. Minimum of twice the Nyquist condition from DCON m limits is enforced.
coil_flag=t ! Calculate external perturbation from coils in coil.in
data_flag=f ! Apply perturbation from thrid party code output
data_type="surfmn" ! Name of said code. Choose: surfmn,??
infile="path/to/file" ! Path to said output
nmin=0 ! Extra parameters needed to read surfmn files.
nmax=64 ! Extra parameters needed to read surfmn files.
mmin=-64 ! Extra parameters needed to read surfmn files.
mmax=64 ! Extra parameters needed to read surfmn files.
harmonic_flag=f ! Apply self-defined perturbations
cosmn(3)=1e-4 ! Real part of m=3 component (n defined by DCON run).
sinmn(6)=1e-4 ! Imaginary part of m=6 component. etc.
displacement_flag=f ! Perturbations are displacement (m). Default is field (T).
fixed_boundary_flag=f ! Total perturbation = external perturbation
mode_flag = .F. ! Return only a single DCON eigenfunction defined by mode (ignores external pert.)
mode = 0 ! DCON eigen function index (1 is least stable)
filter_types = 'w' ! Isolate the filter_modes of energy (w), reluctance (r), permiability (p) or singular coupling (s) eigenvectors/singular-vectors
filter_modes = 0 ! Number of modes isolated, negative values isolate backwards from the last mode
/
&GPEC_CONTROL
resp_index=1 ! Integer 0-4, determining method of inductance matrix formation. 0 uses energy identity, and is recommended for ideal MHD. Use 1 when using resistive DCON.
sing_spot=3.38e-3 ! Minimum distance to singular surface in (m-nq) space
sing_npsi=100 ! Number of points between rational surfaces used for b-normal splines across sin_spot
reg_flag=f ! Regularize the solution near resonant surfaces
reg_spot=5e-2 ! Regularize by (m-nq)^2/((m-nq)^2+reg_spot^2)
chebyshev_flag=f ! Fit Chebyshev polynomials to the displacement.
nche=20 ! Number of Chebyshev polynomials used
use_classic_splines = f ! Use a classical cubic spline instead of tri-diagonal solution for splines with extrapolation boundary conditions
/
&GPEC_OUTPUT
jac_out="hamada" ! Coordinate system of outputs
jsurf_out=0 ! True(1) for area weighted outputs (only appends standard gpec_control output)
tmag_out=1 ! True(1) for magnetic toroidal angle outputs (affects singcoup and singfld, appends gpec_control)
mlim_out=64 ! If positive, Fourier decomposed outputs will span -mlim_out<=m<=mlim_out. Default (<=0) outputs on internal m range.
resp_flag=t ! Output energy, reluctance, inductance, and permeability eigenvalues, eigenvectors, and matrices
filter_flag=t ! Outputs energy, reluctance, pereability, and (optionally) the singular-coupling eigenmodes on the control surface
singcoup_flag=t ! Calculate coupling of each m to resonant surfaces
singfld_flag=t ! Output of resonant surface quantities (flux,current,etc.)
vsingfld_flag=t ! Output of vacuum resonant surface quantities
singthresh_flag=f ! Calculate Callen resonant field penetration threshold (requires kinetic profiles from pentrc.in)
pmodb_flag=f ! Outputs delta-B_Lagrangian
xbnormal_flag=t ! Outputs normal displacement and field in plasma
vbnormal_flag=f ! Outputs vacuum normal displacement and field in plasma
xclebsch_flag=t ! Outputs clebsch coordinate displacement for PENT
dw_flag=f ! Outputs self-consistent energy and torque profiles (when kin_flag=t from DCON)
eqbrzphi_flag=f ! Outputs equilibrium field on regular (r,z) grid
brzphi_flag=t ! Outputs perturbed field on regular (r,z) grid
xrzphi_flag=f ! Outputs displacement on regular (r,z) grid
vbrzphi_flag=f ! Outputs field on regular (r,z) grid due to total boundary surface current (not a real field)
vvbrzphi_flag=f ! Outputs field on regular (r,z) grid due to exeternal perturbation surface current (not a real field)
vsbrzphi_flag=f ! Outputs the field on regular (r,z) grid due to singular surface currents
ss_flag(3)=f ! Outputs the field on regular (r,z) grid due to the 3rd singular surface current, etc.
xbrzphifun_flag=f ! Outputs r,z, and perturbation vectors on psi,theta grid
arzphifun_flag=f ! Outputs r,z, and vector potential on psi,theta grid
nr = 64 ! Number of radial grid points in rzphi outputs
nz = 124 ! Number of vertical grid points in rzphi outputs
fun_flag=t ! Write outputs to (r,z) in addition to default (psi,m)
flux_flag=f ! Write outputs to (psi,theta) in addition to default (psi,m)
bin_flag=f ! Write binary outputs for use with xdraw
bin_2d_flag=f ! Write 2D binary outputs for use with xdraw
out_ahg2msc=f ! Output deprecated ahg2msc.out files. This is used to communicate with vacuum, but is now done through memory.
max_linesout=0 ! Maximum length of ascii data tables enforced by increasing radial step size (iff .gt. 0).
verbose=t ! Print run log to terminal
/
&GPEC_DIAGNOSE
timeit=t ! Print timer splits for major subroutines
radvar_flag=f ! Map various radial variables (rho,psi_tor) on psi_n grid
/