-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmephit_g33353_2900_EQH.in
102 lines (102 loc) · 3.97 KB
/
mephit_g33353_2900_EQH.in
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
89
90
91
92
93
94
95
96
97
98
99
100
101
102
&scalars
! 0: None, 1: Err, 2: Warn, 3: Info, 4: Debug
config%log_level = 4
! 0: single iteration, 1: direct iterations, 2: preconditioned iterations
config%runmode = 2
! 0: EPS 2018 paper profile, 1: parabolic profile, 2: EQDSK profile
config%pres_prof = 2
! 0: only PS current, 1: rot B_0, 2: EQDSK profile
config%curr_prof = 2
! 0: surface integral, 1: field line integral, 2: EQDSK profile
config%q_prof = 1
! 0: geometric series around resonance, 1: Gaussian profile
config%refinement_scheme = 0
! 0: iMHD, 1: FLRE 2nd order
config%currn_model = 1
! 0: Viktor Nemov's code, 1: GPEC, 2: pre-computed Fourier transform
config%vac_src = 2
! prefactor in Biot-Savart law, e.g. 5 windings and conversion from A to statA / c_0
config%Biot_Savart_prefactor = 5.0d-1
!> file containing coil field data
config%coil_file = '/proj/plasma/CODE/lainer_p/MEPHIT/build/run/AUG_B_coils.h5'
!> file containing coil currents
config%currents_file = '/proj/plasma/DATA/BALANCE/COIL/33353/33353.2900_coil_markl.dat'
!> file containing equilibrium density profile
config%dens_file = '/proj/plasma/DATA/BALANCE/PROF/33353/33353.2900_ne_PED_MMARKL_rho_pol.dat'
!> file containing equilibrium electron temperature profile
config%temp_e_file = '/proj/plasma/DATA/BALANCE/PROF/33353/33353.2900_Te_PED_MMARKL_rho_pol.dat'
!> file containing equilibrium ion temperature profile
config%temp_i_file = '/proj/plasma/DATA/BALANCE/PROF/33353/33353.2900_Ti_PED_MMARKL_rho_pol.dat'
!> file containing equilibrium radial electric field profile
config%E_r_file = '/proj/plasma/DATA/BALANCE/PROF/33353/33353_2900_Er.dat'
! mass number of ions
config%m_i = 2.0
! toroidal mode number of perturbation
config%n = 2
! maximum number of iteration steps
config%niter = 75
! maximum relative error in fixed-point iteration
config%iter_rel_err = 1d-9
! maximum number of iterations in Arnoldi method
config%nkrylov = 100
! threshold for eigenvalues in preconditioner
config%ritz_threshold = 0.5d0
! maximim relative error for eigenvalues in preconditioner
config%ritz_rel_err = 1d-9
! maximum distance between flux surfaces along theta = 0
config%max_Delta_rad = 1.0d0
! maximum number of points per flux surface (0: no maximum)
config%pol_max = 0
! minimum number of points per flux surface (0: no minimum)
config%pol_min = 0
! ignore resonance position where q = 1, which is usually spurious
config%ignore_q1_res = .true.
! enable damping of Pfirsch-Schlueter current near resonances
config%damp = .true.
! maximal number of eigenvectors to be saved
config%max_eig_out = 10
! maximum poloidal mode number for Fourier transform of results
config%m_max = 24
! number of points in sweep over electron fluid resonance
config%resonance_sweep = 0
! offset added to radial electric field
config%offset_E_r = 0.0d0
/
&arrays
refinement(3) = 1.5d0
refinement(4) = 1.5d0
refinement(5) = 1.5d0
refinement(6) = 1.5d0
refinement(7) = 1.5d0
refinement(8) = 1.5d0
refinement(9) = 1.5d0
refinement(10) = 0.0d0
refinement(11) = 0.0d0
add_fine(3) = 4
add_fine(4) = 4
add_fine(5) = 4
add_fine(6) = 4
add_fine(7) = 4
add_fine(8) = 4
add_fine(9) = 2
add_fine(10) = 0
add_fine(11) = 0
Delta_rad_res(3) = 1.300211470239295e-02
Delta_rad_res(4) = 1.782989651072676e-02
Delta_rad_res(5) = 1.786554488832890e-02
Delta_rad_res(6) = 1.177360561550166e-02
Delta_rad_res(7) = 1.330771149330919e-02
Delta_rad_res(8) = 2.770845072558438e-02
Delta_rad_res(9) = 5.347772320733589e-02
Delta_rad_res(10) = 6.066328245075899e-02
Delta_rad_res(11) = 3.464590203985631e-02
sheet_current_factor(3) = 1d5
sheet_current_factor(4) = 1d5
sheet_current_factor(5) = 1d5
sheet_current_factor(6) = 1d5
sheet_current_factor(7) = 1d5
sheet_current_factor(8) = 1d5
sheet_current_factor(9) = 1d5
sheet_current_factor(10) = 0d0
sheet_current_factor(11) = 0d0
/