-
Notifications
You must be signed in to change notification settings - Fork 3
/
data16
211 lines (204 loc) · 6.17 KB
/
data16
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
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
&PEDBAR
NBAR = 2
NPED = 5
FCOND = 0.650
PEDFAC = 0.50
DPERP = 0.5
XNU = 2.5
ZEFF = 1.6
Z0 = .6
CFZINJ = 0.00
CFZINT = 0.058
CFZINJTB = 0.0
CFZINTTB = 0.058
QZMULTCORE = 1.0
ALPHAT1 = 0.05
A1 = 0.0
A2 = 1.0
PLASMACUR = 1.22
PBEAM = 4.89
ENBI = 80.0
PHIIN = 0.0
RADMULT = 1.
TAVMULT = 1.
RHOGAUSS = 0.45
SIGMA =.1
ZEFFC = 2.5
TAUE = 1.E-2
TAURATIO = 0.50
IOPTCONF = 0
C89 = 1.
C89N = 1.
IGRADPED = 2
CHEI = 0.96
RXLTIXLTE = 2.0
CHIREDGE = 3.65
IOPTPINCH = 0
CPINCH = 0.0
VPINCH = 0.0
CC = 2.5
fluxheatianom = 3.0e4
fluxheateanom = 5.0e4
FI = 0.13
DTBI = 0.3
HEATFRACE = 0.4
IOPTLN = 1
CIE = 2.0
CNE = 1.0
CHIRTBI = 1.0E-1
CHIRTBE = 1.0E-1
HN = 1.2
IOPTPART = 1
AMASS = 2.0
XNCTRPED =1.67
TCTRPED = 2.80
NNEUT = 40
IOPTPED = 1
ERAD = -5.E3
VPHI = 5.E3
VTHET = 5.E3
VRAD = 1.E3
XNUBAR = 2.5
XLQ = -1.E0
XLE = 1.E2
XLVTHET = 1.E2
XLVPHI = 1.E2
XNAV = 2.E20
XNRHOPED = 4.0
RSYN = .75
IPROF = 2
JJOPTPED = 9
FSHAPE = 1.0
IOPTGRAD = 1
RCOEF = 1.0
GRADTPEDIS = 0.0
IOPTPEDWIDTH = 2
CPED = 2.E-2
CPED1 = 1.50
CPED2 = 0.25
FION = 1.35
ISHAPE = 1
IDRAKE = 0
SHEARQ = 1.67
CBALL = 3.0
FISEP = 0.3
chitbli = 1.0
chitble = 1.0
SHEAR0 = 1.6
cshear = 0.6
ASS = 1.0
IOPTTAUE =89
IOPTTRAN = 0
MODPED = 3
CMHD = 3.0
ccrit = 1.0
cmfp = 1.0
/
NBAR = # POLOIDAL GYRORADII ACROSS TRANSPORT BARRIER
NPED = # POLOIDAL GYRORADII ACROSS PEDESTAL
FI = DTBI = FI*CHITBI
FCOND = CONDUCTIVE FRACTION OF HEAT THROUGH SEPARATRIX
XNU = TEMP DEPENDENCE OF CHIR--NU FACTOR
ZEFF = ZEFF
Z0 = PARAMETER IN FRICTION C2 FORMULA RELATED TO IMP Z
CFZINJ = FZCORE FOR INJECTED IMPURITY IN CORE
CFZINT = FZCORE FOR INTRINSIC IMPURITY IN CORE
CFZINJTB = FZTRANBARRIER FOR INJECTED IMPURITY
CFZINTTB = FZTRANBARRIER FOR INTRINSIC IMPURITY
QZMULTCORE = INTRINSIC IMPURITY RADIATION MULTIPLIER
PLASMACUR = PLASMA CURRENT (MA)
AMINOR = PLASMA RADIUS (m)
RMAJOR = MAJOR RADIUS (m)
ALPHAT = PARABOLIC TO ALPHAT POWER Te(r)
PBEAM = BEAM POWER (MW)
ENBI = BEAM PARTICLE ENERGY (keV)
PHIIN = BEAM ANGLE WITH RESPECT TO INCIDENT NORMAL (RAD)
RHOGAUSS = CENTER OF GAUSSIAN DISTR OF BEAM POWER
SIGMA = SIGMA IN GAUSSIAN DISTRIBUTION
NZIMP = IMPURITY ATOMIC #
FZ = IMPURITY CONCENTRATION FRACTION
dperp = perp diffusion coef (m^2/s)
alphan = parabolic to alphan power n(r)
rcoef = recycling coefficient
XNCTRPED = RATIO CENTRAL TO PEDESTAL DENSITIES
TCTRPED = RATIO CENTRAL TO PEDESTAL TEMPERATURES
IOPTPED = 0, USE INPUT GRADIENTS IN TRANSPORT BARRIER
= 1, CALCULATE GRADIENTS IN TRANSPORT BARRIER
TAUE = ENERGY CONFINEMENT TIME (s)
TAURATIO = TAUP/TAUE
HCONF = H-MODE CONFINEMENT ENHANCEMENT (ITER 89-P)
C89 = CONSTANT FOR CONFINEMENT DEGRADATION
CHEI = CHI-E/CHI-I FOR CONFINEMENT DEGRADATION CALC.
HN = PARTICLE CONFINEMENT ENHANCEMENT FACTOR
AMASS = MASS OF PLASMA ION (amu)
CHIRTB I&E = TRANS BARRIER CHIR FOR GROWTH RATE CALC (m^2/s)
ERAD = RADIAL ELECTRIC FIELD IN EDGE (V/m)
VTHET = V-THETA IN EDGE (m/s)
VRAD = V-RADIAL IN EDGE (m/s)
XNUBAR = TEMP COEF EXPONENT OF CHI IN TRANS BARRIER
XLQ = -(dQ/dr)/Q IN TRANSPORT BARRIER (m-1)
XLE = -(dE/dr)/E IN TRANSPORT BARRIER (m-1)
XLVTHET = -(dVtheta/dr)/Vtheta IN TRANSPORT BARRIER (m-1)
XLVPHI = -(dVphi/dr)/Vphi IN TRANSPORT BARRIER (m-1)
XNAV = AVERAGE ION DENSITY FOR STABILITY ANALYSIS
XNRHOPED = # POLOIDAL RHO ACROSS TRANS BARRIER
RSYN = WALL REFLECTION COEFFICIENT FOR SYNCHROTRON RADIATION
IPROF = 1, (1-RHO)^ALPHAT TEMP PROFILE
= 2, (1-RHO*2)^ALPHAT TEMP PROFILE
JJOPTPED = 1, BALLOONING LIMIT, EQIPARTION PEDESTAL MODEL,
CALCULATE CHIREDGE TO REMOVE HEAT & CHIRCORE
= 2, BALLOONING & HEAT REMOVAL LIMITS PEDESTAL MODEL
USE INPUT CHIR EDGE AND CORE
= 3, SAME AS #2 BUT SIMPLER ALGORITHMS
= 4, nped determined by ionization source in TB, Tped by heat removal
= 5, nped = f*nav, tped from heat removal, adjust nped for ballooning limit
= 6, nped = f*nav, tped from heat removal, adjust tped for ballooning limit
= 7, nped = f*nav, tped from heat removal
= 8, nped = nsep, tped = tsep, no pedestal
= 9, NPED=F*NAV,TPED FROM BALLOON STAB, CHIREDGE FROM HEAT REMOVAL
=10, CALCULATED FROM THEORY
JOPTPED2 = 1 EVALUATE PEDESTAL WITH TB AV N & T
= 2 EVALUATE PEDESTAL WITH PEDESTAL N & T
IOPTGRAD = 1 CALCULATE GRAD SCALE LENGTHS FOR MARFE DEN LIMIT
= 2 USE INPUT " " " " " " "
IOPTPEDWIDTH= 1 IGITKHANOV FIT OF PEDESTAL WIDTH
= 2 GROEBNER & OSBORNE FIT OF PEDESTAL WIDTH
= 4 ROGERS & DRAKE
= 5 IDEAL BALLOONING
CPED = SCALING CONSTANT TO USE WITH GROEBNER & OSBORNE FIT OF DELTB
**********cped = 0.02 ********** works for DIII-D
CPED1 = RATIO TIPED/TEPED
CPED2 = SCALING MULTIPLE FOR CRITICAL PRESSURE GRAD FORMULA
FION = RATIO TI/TE GRADIENT SCALE LENGTHS
CHEI = RATIO CHI-E TO CHI-I IN TRANSPORT BARRIER
RXLTIXLTE = RATIO ION-T AND ELECTRON-T GRADIENT SCALE LENGTHS IN TB
VPINCH = RADIAL PINCH VELOCITY IN TB(m/s)
DTBI = ION DIFFUSION COEFFICIENT IN TB (m2/s)
IOPTLN = 2 CALCULATE LN FROM FLUXPART USING DIFFUSION + PINCH
= 1 CALCULATE LN FROM FLUXPART USING NEUTRAL IONIZATION
CIE = LTI/LTE
CNE = LN/LTE
HEATFRACE = FRACTION OF HEAT FLUX ACROSS SEPARATRIX THAT IS ELECTRON
IOPTPINCH = 0 USE INPUT PINCH VELOCITY VPINCH (M/S)
= 1 USE INPUT CPINCH
IGRADPED = 1 CALCULATE PEDESTAL GRADIENTS AT SEPARATRIX
= 2 CALCULATE PEDESTAL GRADIENTS AT TRANBAR MIDPOINT
= 3 CALCULATE PEDESTAL GRAD INCLUDING CONV. HEAT FLUXES FROM PEDESTAL
CC = 3/2 OR 5/2 THE CONVECTIVE MULTIPLIER
ISHAPE = 1, DPCRIT SHAPE FACTOR = (1+1.2DELTA)(1+KAPPA^2)
2, = (1+KAPPA^2(1+2DELTA^2-1.3DELTA^3)
SHEARQ = SHEAR PARAMETER
CBALL = MULTIPLIER TO USE WITH IDEAL BALLOONING DELTB5
FISEP = ION FRACTION OF ENERGY CROSSING SEPARATRIX
IDRAKE = 0 PCRIT CONSTRAINT ON XLP
= 1 DELTA = 0.5DELTB ON RHS ROGERS&DRAKE
= 2 DP/DR = P/0.5DELTB ON LHS ROGERS&DRAKE
chitbli&e = L-mode transport barrier chi's for L-H critical flux
SHEAR0 = EDGE SHEAR BEFORE BOOTSTRAP CURRENT CORRECTION
IOPTTAUE = 89 ITER89P TAUE SCALING LAW
= 98 ITER98(y,2) TAUE SCALING LAW
MODPED = 0, DELTTEMP=DELTDENS=DELTPRESS
= 1, DELTDENS=NEUT MFP,DELTTEMP=DELTPRESS
= 2, DELTTEMP=DELTDENS=NEUT MFP
= 3, DELTPRESS=DELTTEMP=DELTDENS BALLOON-PEELING
CMHD = BALLOONING SECOND STABILITY ACCESS FACTOR