gs2crmod 0.5.2
Sign up to get free protection for your applications and to get access to all the features.
- data/.document +5 -0
- data/Gemfile +13 -0
- data/LICENSE.txt +20 -0
- data/README.md +4 -0
- data/README.rdoc +19 -0
- data/Rakefile +56 -0
- data/VERSION +1 -0
- data/ext/extconf.rb +9 -0
- data/ext/gs2crmod_ext.c +366 -0
- data/gs2crmod.gemspec +98 -0
- data/include/gs2crmod_ext.h +58 -0
- data/lib/gs2crmod/astrogk/astrogk.rb +201 -0
- data/lib/gs2crmod/astrogk/calculations.rb +57 -0
- data/lib/gs2crmod/astrogk/check_convergence.rb +7 -0
- data/lib/gs2crmod/astrogk/deleted_variables.rb +76 -0
- data/lib/gs2crmod/astrogk/graphs.rb +13 -0
- data/lib/gs2crmod/astrogk/gsl_data.rb +13 -0
- data/lib/gs2crmod/astrogk/gsl_tools.rb +182 -0
- data/lib/gs2crmod/astrogk/ingen.rb +18 -0
- data/lib/gs2crmod/astrogk/input_file_tools.rb +7 -0
- data/lib/gs2crmod/astrogk/namelist_tools.rb +14 -0
- data/lib/gs2crmod/astrogk/namelists.rb +2800 -0
- data/lib/gs2crmod/astrogk/properties.rb +17 -0
- data/lib/gs2crmod/astrogk/species_dependent_namelists.rb +228 -0
- data/lib/gs2crmod/astrogk/test_gs2.rb +231 -0
- data/lib/gs2crmod/astrogk.rb +200 -0
- data/lib/gs2crmod/calculations.rb +780 -0
- data/lib/gs2crmod/check_convergence.rb +179 -0
- data/lib/gs2crmod/deleted_variables.rb +916 -0
- data/lib/gs2crmod/graphs.rb +1899 -0
- data/lib/gs2crmod/graphs_rdoc.rb +556 -0
- data/lib/gs2crmod/gs2.rb +1143 -0
- data/lib/gs2crmod/gsl_data.rb +1181 -0
- data/lib/gs2crmod/gsl_data_3d.rb +705 -0
- data/lib/gs2crmod/gsl_tools.rb +187 -0
- data/lib/gs2crmod/ingen.rb +218 -0
- data/lib/gs2crmod/namelists.rb +5142 -0
- data/lib/gs2crmod/properties.rb +22 -0
- data/lib/gs2crmod/species_dependent_namelists.rb +228 -0
- data/lib/gs2crmod/test_gs2.rb +231 -0
- data/lib/gs2crmod.rb +2 -0
- data/lib/gs2crmod_extension.rb +1 -0
- data/test/helper.rb +18 -0
- data/test/test_gs2crmod.rb +7 -0
- metadata +176 -0
@@ -0,0 +1,556 @@
|
|
1
|
+
class CodeRunner::Gs2::GraphKits
|
2
|
+
# Graph of apar^2 vs time integrated over all space. No options
|
3
|
+
def apar2_vs_time_graphkit
|
4
|
+
end
|
5
|
+
# Plot the eigenfunction along the extended domain. Options mag, norm, z can be specified by using a short hand in the name of the graph, eg. efnmagnormz, efnmag, efnnorm etc. If the range is set to 0, it plots the whole eigenfunction. Otherwise it plot a small bit of it. Only specify kx or kx_index if magnetic shear is 0.
|
6
|
+
# Options:
|
7
|
+
#
|
8
|
+
# mag: Plot the magnitude, e.g. mag: true
|
9
|
+
#
|
10
|
+
# norm: Normalise the graph so that its maximum is 1, e.g. norm: true
|
11
|
+
#
|
12
|
+
# z: Plot quantities vs z = theta/shat rather than theta. See Beer, Cowley Hammet 1996, eg. z: true
|
13
|
+
#
|
14
|
+
# flip: Flip the y axis, e.g. flip: true
|
15
|
+
#
|
16
|
+
# range:
|
17
|
+
#
|
18
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
19
|
+
#
|
20
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
21
|
+
#
|
22
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
23
|
+
#
|
24
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
25
|
+
#
|
26
|
+
# strongest_non_zonal_mode: Plot the graph requested for the mode with the highest value of phi^2. Overrides ky, kx, ky_index, kx_index. Can be set true or false; e.g. strongest_non_zonal_mode: true
|
27
|
+
def efn_graphkit
|
28
|
+
end
|
29
|
+
# Plot the eigenfunction along the extended domain. Options mag, norm, z can be specified by using a short hand in the name of the graph, eg. efnmagnormz, efnmag, efnnorm etc. If the range is set to 0, it plots the whole eigenfunction. Otherwise it plot a small bit of it. Only specify kx or kx_index if magnetic shear is 0.
|
30
|
+
# Options:
|
31
|
+
#
|
32
|
+
# mag: Plot the magnitude, e.g. mag: true
|
33
|
+
#
|
34
|
+
# norm: Normalise the graph so that its maximum is 1, e.g. norm: true
|
35
|
+
#
|
36
|
+
# z: Plot quantities vs z = theta/shat rather than theta. See Beer, Cowley Hammet 1996, eg. z: true
|
37
|
+
#
|
38
|
+
# flip: Flip the y axis, e.g. flip: true
|
39
|
+
#
|
40
|
+
# range:
|
41
|
+
#
|
42
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
43
|
+
#
|
44
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
45
|
+
#
|
46
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
47
|
+
#
|
48
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
49
|
+
#
|
50
|
+
# strongest_non_zonal_mode: Plot the graph requested for the mode with the highest value of phi^2. Overrides ky, kx, ky_index, kx_index. Can be set true or false; e.g. strongest_non_zonal_mode: true
|
51
|
+
def efnim_graphkit
|
52
|
+
end
|
53
|
+
# Plot the eigenfunction along the extended domain. Options mag, norm, z can be specified by using a short hand in the name of the graph, eg. efnmagnormz, efnmag, efnnorm etc. If the range is set to 0, it plots the whole eigenfunction. Otherwise it plot a small bit of it. Only specify kx or kx_index if magnetic shear is 0.
|
54
|
+
# Options:
|
55
|
+
#
|
56
|
+
# mag: Plot the magnitude, e.g. mag: true
|
57
|
+
#
|
58
|
+
# norm: Normalise the graph so that its maximum is 1, e.g. norm: true
|
59
|
+
#
|
60
|
+
# z: Plot quantities vs z = theta/shat rather than theta. See Beer, Cowley Hammet 1996, eg. z: true
|
61
|
+
#
|
62
|
+
# flip: Flip the y axis, e.g. flip: true
|
63
|
+
#
|
64
|
+
# range:
|
65
|
+
#
|
66
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
67
|
+
#
|
68
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
69
|
+
#
|
70
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
71
|
+
#
|
72
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
73
|
+
#
|
74
|
+
# strongest_non_zonal_mode: Plot the graph requested for the mode with the highest value of phi^2. Overrides ky, kx, ky_index, kx_index. Can be set true or false; e.g. strongest_non_zonal_mode: true
|
75
|
+
def efnmag_graphkit
|
76
|
+
end
|
77
|
+
# Plot the eigenfunction along the extended domain. Options mag, norm, z can be specified by using a short hand in the name of the graph, eg. efnmagnormz, efnmag, efnnorm etc. If the range is set to 0, it plots the whole eigenfunction. Otherwise it plot a small bit of it. Only specify kx or kx_index if magnetic shear is 0.
|
78
|
+
# Options:
|
79
|
+
#
|
80
|
+
# mag: Plot the magnitude, e.g. mag: true
|
81
|
+
#
|
82
|
+
# norm: Normalise the graph so that its maximum is 1, e.g. norm: true
|
83
|
+
#
|
84
|
+
# z: Plot quantities vs z = theta/shat rather than theta. See Beer, Cowley Hammet 1996, eg. z: true
|
85
|
+
#
|
86
|
+
# flip: Flip the y axis, e.g. flip: true
|
87
|
+
#
|
88
|
+
# range:
|
89
|
+
#
|
90
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
91
|
+
#
|
92
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
93
|
+
#
|
94
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
95
|
+
#
|
96
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
97
|
+
#
|
98
|
+
# strongest_non_zonal_mode: Plot the graph requested for the mode with the highest value of phi^2. Overrides ky, kx, ky_index, kx_index. Can be set true or false; e.g. strongest_non_zonal_mode: true
|
99
|
+
def eigenfunction_graphkit
|
100
|
+
end
|
101
|
+
# 'es_heat_by_ky_vs_time' or 'es_heat_by_kx_vs_time': Electrostatic Heat Flux vs Time for a given kx or ky, integrated over the other direction
|
102
|
+
# Options:
|
103
|
+
#
|
104
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
105
|
+
#
|
106
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
107
|
+
#
|
108
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
109
|
+
#
|
110
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
111
|
+
#
|
112
|
+
# species_index: Which GS2 species to plot the graph for (1-based).
|
113
|
+
def es_heat_by_kx_vs_time_graphkit
|
114
|
+
end
|
115
|
+
# 'es_heat_by_ky_vs_time' or 'es_heat_by_kx_vs_time': Electrostatic Heat Flux vs Time for a given kx or ky, integrated over the other direction
|
116
|
+
# Options:
|
117
|
+
#
|
118
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
119
|
+
#
|
120
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
121
|
+
#
|
122
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
123
|
+
#
|
124
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
125
|
+
#
|
126
|
+
# species_index: Which GS2 species to plot the graph for (1-based).
|
127
|
+
def es_heat_by_kxy_or_mode_vs_time_graphkit
|
128
|
+
end
|
129
|
+
# 'es_heat_by_ky_vs_time' or 'es_heat_by_kx_vs_time': Electrostatic Heat Flux vs Time for a given kx or ky, integrated over the other direction
|
130
|
+
# Options:
|
131
|
+
#
|
132
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
133
|
+
#
|
134
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
135
|
+
#
|
136
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
137
|
+
#
|
138
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
139
|
+
#
|
140
|
+
# species_index: Which GS2 species to plot the graph for (1-based).
|
141
|
+
def es_heat_by_ky_vs_time_graphkit
|
142
|
+
end
|
143
|
+
# 'es_heat_by_ky_vs_time' or 'es_heat_by_kx_vs_time': Electrostatic Heat Flux vs Time for a given kx or ky, integrated over the other direction
|
144
|
+
# Options:
|
145
|
+
#
|
146
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
147
|
+
#
|
148
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
149
|
+
#
|
150
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
151
|
+
#
|
152
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
153
|
+
#
|
154
|
+
# species_index: Which GS2 species to plot the graph for (1-based).
|
155
|
+
def es_heat_by_mode_vs_time_graphkit
|
156
|
+
end
|
157
|
+
# Graph of electrostatic contribution to heat flux at a given time vs kx and ky
|
158
|
+
# Options:
|
159
|
+
#
|
160
|
+
# with: Gnuplot Option (may not apply when using other packages), e.g. with: 'lp' or with 'pm3d palette'
|
161
|
+
def es_heat_flux_vs_ky_vs_kx_graphkit
|
162
|
+
end
|
163
|
+
# Heat flux vs time for each species.
|
164
|
+
# Options:
|
165
|
+
#
|
166
|
+
# t_index_window: [begin, end], window of time indices to plot (e.g. t_index_window: [0,10])
|
167
|
+
#
|
168
|
+
# species_index: Which GS2 species to plot the graph for (1-based).
|
169
|
+
def es_heat_flux_vs_time_graphkit
|
170
|
+
end
|
171
|
+
# Momentum flux vs time for each species.
|
172
|
+
# Options:
|
173
|
+
#
|
174
|
+
# t_index_window: [begin, end], window of time indices to plot (e.g. t_index_window: [0,10])
|
175
|
+
#
|
176
|
+
# species_index: Which GS2 species to plot the graph for (1-based).
|
177
|
+
def es_mom_flux_vs_time_graphkit
|
178
|
+
end
|
179
|
+
# 'growth_rate_by_ky_vs_time' or 'growth_rate_by_kx_vs_time': Growth rate vs time for a given kx or ky, integrated over the other direction
|
180
|
+
# Options:
|
181
|
+
#
|
182
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
183
|
+
#
|
184
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
185
|
+
#
|
186
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
187
|
+
#
|
188
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
189
|
+
def growth_rate_by_kx_vs_time_graphkit
|
190
|
+
end
|
191
|
+
# 'growth_rate_by_ky_vs_time' or 'growth_rate_by_kx_vs_time': Growth rate vs time for a given kx or ky, integrated over the other direction
|
192
|
+
# Options:
|
193
|
+
#
|
194
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
195
|
+
#
|
196
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
197
|
+
#
|
198
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
199
|
+
#
|
200
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
201
|
+
def growth_rate_by_kxy_or_mode_vs_time_graphkit
|
202
|
+
end
|
203
|
+
# 'growth_rate_by_ky_vs_time' or 'growth_rate_by_kx_vs_time': Growth rate vs time for a given kx or ky, integrated over the other direction
|
204
|
+
# Options:
|
205
|
+
#
|
206
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
207
|
+
#
|
208
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
209
|
+
#
|
210
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
211
|
+
#
|
212
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
213
|
+
def growth_rate_by_ky_vs_time_graphkit
|
214
|
+
end
|
215
|
+
# 'growth_rate_by_ky_vs_time' or 'growth_rate_by_kx_vs_time': Growth rate vs time for a given kx or ky, integrated over the other direction
|
216
|
+
# Options:
|
217
|
+
#
|
218
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
219
|
+
#
|
220
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
221
|
+
#
|
222
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
223
|
+
#
|
224
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
225
|
+
def growth_rate_by_mode_vs_time_graphkit
|
226
|
+
end
|
227
|
+
# growth_rate_vs_ky or growth_rate_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /growth_rate_vs_kx_vs_ky/.
|
228
|
+
def growth_rate_vs_kx_graphkit
|
229
|
+
end
|
230
|
+
# 3D plot of growth rates vs ky and kx for phi^2
|
231
|
+
# Options:
|
232
|
+
#
|
233
|
+
# with: Gnuplot Option (may not apply when using other packages), e.g. with: 'lp' or with 'pm3d palette'
|
234
|
+
def growth_rate_vs_kx_vs_ky_graphkit
|
235
|
+
end
|
236
|
+
# growth_rate_vs_ky or growth_rate_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /growth_rate_vs_kx_vs_ky/.
|
237
|
+
def growth_rate_vs_kxy_graphkit
|
238
|
+
end
|
239
|
+
# growth_rate_vs_ky or growth_rate_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /growth_rate_vs_kx_vs_ky/.
|
240
|
+
def growth_rate_vs_ky_graphkit
|
241
|
+
end
|
242
|
+
# Graph of total heat flux vs time. No options
|
243
|
+
def hflux_tot_vs_time_graphkit
|
244
|
+
end
|
245
|
+
# Graph of the k_parallel at a given kx and ky
|
246
|
+
# Options:
|
247
|
+
#
|
248
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
249
|
+
#
|
250
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
251
|
+
#
|
252
|
+
# strongest_non_zonal_mode: Plot the graph requested for the mode with the highest value of phi^2. Overrides ky, kx, ky_index, kx_index. Can be set true or false; e.g. strongest_non_zonal_mode: true
|
253
|
+
def kpar_spectrum_graphkit
|
254
|
+
end
|
255
|
+
# ky_spectrum or kx_spectrum: Graph of phi^2 vs kx or ky
|
256
|
+
# Options:
|
257
|
+
#
|
258
|
+
# t: float, value of time at which to plot (e.g. t: 2.45)
|
259
|
+
#
|
260
|
+
# t_index: The (1-based) time index
|
261
|
+
def kx_spectrum_graphkit
|
262
|
+
end
|
263
|
+
# ky_spectrum or kx_spectrum: Graph of phi^2 vs kx or ky
|
264
|
+
# Options:
|
265
|
+
#
|
266
|
+
# t: float, value of time at which to plot (e.g. t: 2.45)
|
267
|
+
#
|
268
|
+
# t_index: The (1-based) time index
|
269
|
+
def kxy_spectrum_graphkit
|
270
|
+
end
|
271
|
+
# ky_spectrum or kx_spectrum: Graph of phi^2 vs kx or ky
|
272
|
+
# Options:
|
273
|
+
#
|
274
|
+
# t: float, value of time at which to plot (e.g. t: 2.45)
|
275
|
+
#
|
276
|
+
# t_index: The (1-based) time index
|
277
|
+
def ky_spectrum_graphkit
|
278
|
+
end
|
279
|
+
# A graph of the evolution of a single Lagrangian kx vs Eulerian kx and ky. Principally for debugging purposes
|
280
|
+
def lagrangian_kx_graphkit
|
281
|
+
end
|
282
|
+
# The potential at the outboard midplane
|
283
|
+
# Options:
|
284
|
+
#
|
285
|
+
# rgbformulae: Gnuplot Option (may not apply when using other packages), sets colour mapping. See gnuplot help set rgbformulae
|
286
|
+
#
|
287
|
+
# limit: Limit the range of quantity begin plotted - any values of the quantity outside the limits will be set to the limit: eg. limit: [0,80]
|
288
|
+
def phi0_vs_x_vs_y_graphkit
|
289
|
+
end
|
290
|
+
# 'phi2_by_ky_vs_time' or 'phi2_by_kx_vs_time': Phi^2 over time for a given kx or ky, integrated over the other direction
|
291
|
+
# Options:
|
292
|
+
#
|
293
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
294
|
+
#
|
295
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
296
|
+
#
|
297
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
298
|
+
#
|
299
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
300
|
+
def phi2_by_kx_vs_time_graphkit
|
301
|
+
end
|
302
|
+
# 'phi2_by_ky_vs_time' or 'phi2_by_kx_vs_time': Phi^2 over time for a given kx or ky, integrated over the other direction
|
303
|
+
# Options:
|
304
|
+
#
|
305
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
306
|
+
#
|
307
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
308
|
+
#
|
309
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
310
|
+
#
|
311
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
312
|
+
def phi2_by_kxy_or_mode_vs_time_graphkit
|
313
|
+
end
|
314
|
+
# 'phi2_by_ky_vs_time' or 'phi2_by_kx_vs_time': Phi^2 over time for a given kx or ky, integrated over the other direction
|
315
|
+
# Options:
|
316
|
+
#
|
317
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
318
|
+
#
|
319
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
320
|
+
#
|
321
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
322
|
+
#
|
323
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
324
|
+
def phi2_by_ky_vs_time_graphkit
|
325
|
+
end
|
326
|
+
# 'phi2_by_ky_vs_time' or 'phi2_by_kx_vs_time': Phi^2 over time for a given kx or ky, integrated over the other direction
|
327
|
+
# Options:
|
328
|
+
#
|
329
|
+
# ky: float, value of ky at which to plot (e.g. ky: 0.1)
|
330
|
+
#
|
331
|
+
# ky_index: integer, index of ky at which to plot (e.g. ky_index: 20)
|
332
|
+
#
|
333
|
+
# kx: float, value of kx at which to plot (e.g. kx: 0.1)
|
334
|
+
#
|
335
|
+
# kx_index: integer, index of kx at which to plot (e.g. kx_index: 20)
|
336
|
+
def phi2_by_mode_vs_time_graphkit
|
337
|
+
end
|
338
|
+
# Graph of phi^2 vs time integrated over all space. No options
|
339
|
+
def phi2tot_vs_time_graphkit
|
340
|
+
end
|
341
|
+
# The potential as a function of cartesian coordinates, on one specified side of the flux tube (specified using the options :coordinate (:x, :y, :theta) and :side (:min, :max))
|
342
|
+
# Options:
|
343
|
+
#
|
344
|
+
# Rgeo:
|
345
|
+
#
|
346
|
+
# n0: The toroidal mode number of the longest y mode. In effect it is the number of periodic copies of the flux tube that will fit in the torus. Periodicity requires that n0 q is also an integer. If you specify :n0 where this is not the case, q will automatically be adjusted until it is
|
347
|
+
#
|
348
|
+
# rho_star: The ratio of the reference Lamour radius to the GS2 normalising length a. Cannot be specified at the same time as n0. If specified, both n0 and q will be adjusted to ensure periodicity
|
349
|
+
#
|
350
|
+
# t_index: The (1-based) time index
|
351
|
+
#
|
352
|
+
# nakx: The number of radial wave numbers to include in the plot. In effect, it is a low pass filter which reduces the resolution in the radial direction without changing the shape of the final surface. Minimum value is 4
|
353
|
+
#
|
354
|
+
# naky: The number of kys to include in the plot. In effect, it is a low pass filter which reduces the resolution in the y direction without changing the shape of the final surface. Minimum value is 4
|
355
|
+
#
|
356
|
+
# gs2_coordinate_factor: When set to 1, plot the graph in GS2 coordinates. When set to 0 plot the graph in real space. Can be set at any value between 0 and 1: the graph will smoothly distort between the two limits
|
357
|
+
#
|
358
|
+
# xmax: The (0-based) index of the maximum value of x to include in the plot
|
359
|
+
#
|
360
|
+
# xmin: The (0-based) index of the minimum value of x to include in the plot
|
361
|
+
#
|
362
|
+
# ymax: The (0-based) index of the maximum value of y to include in the plot
|
363
|
+
#
|
364
|
+
# ymin: The (0-based) index of the minimum value of y to include in the plot
|
365
|
+
#
|
366
|
+
# thetamax: The (0-based) index of the maximum value of theta to include in the plot
|
367
|
+
#
|
368
|
+
# thetamin: The (0-based) index of the minimum value of theta to include in the plot
|
369
|
+
#
|
370
|
+
# ncopies: The number of periodic copies of the flux tube to include
|
371
|
+
#
|
372
|
+
# side:
|
373
|
+
#
|
374
|
+
# coordinate:
|
375
|
+
#
|
376
|
+
# torphi_values: An array of two values of the toroidal angle. The graph will be plotted in between those two values with poloidal cross sections at either end
|
377
|
+
def phi_flux_tube_boundary_surface_graphkit
|
378
|
+
end
|
379
|
+
# The potential as a function of y, x and theta
|
380
|
+
# Options:
|
381
|
+
#
|
382
|
+
# rgbformulae: Gnuplot Option (may not apply when using other packages), sets colour mapping. See gnuplot help set rgbformulae
|
383
|
+
#
|
384
|
+
# limit: Limit the range of quantity begin plotted - any values of the quantity outside the limits will be set to the limit: eg. limit: [0,80]
|
385
|
+
#
|
386
|
+
# t_index: The (1-based) time index
|
387
|
+
def phi_gs2_space_graphkit
|
388
|
+
end
|
389
|
+
# The potential as a function of cartesian coordinates, plus a close up section.
|
390
|
+
# Options:
|
391
|
+
#
|
392
|
+
# Rgeo:
|
393
|
+
#
|
394
|
+
# n0: The toroidal mode number of the longest y mode. In effect it is the number of periodic copies of the flux tube that will fit in the torus. Periodicity requires that n0 q is also an integer. If you specify :n0 where this is not the case, q will automatically be adjusted until it is
|
395
|
+
#
|
396
|
+
# rho_star: The ratio of the reference Lamour radius to the GS2 normalising length a. Cannot be specified at the same time as n0. If specified, both n0 and q will be adjusted to ensure periodicity
|
397
|
+
#
|
398
|
+
# t_index: The (1-based) time index
|
399
|
+
#
|
400
|
+
# nakx: The number of radial wave numbers to include in the plot. In effect, it is a low pass filter which reduces the resolution in the radial direction without changing the shape of the final surface. Minimum value is 4
|
401
|
+
#
|
402
|
+
# naky: The number of kys to include in the plot. In effect, it is a low pass filter which reduces the resolution in the y direction without changing the shape of the final surface. Minimum value is 4
|
403
|
+
#
|
404
|
+
# xmax: The (0-based) index of the maximum value of x to include in the plot
|
405
|
+
#
|
406
|
+
# xmin: The (0-based) index of the minimum value of x to include in the plot
|
407
|
+
#
|
408
|
+
# thetamax: The (0-based) index of the maximum value of theta to include in the plot
|
409
|
+
#
|
410
|
+
# thetamin: The (0-based) index of the minimum value of theta to include in the plot
|
411
|
+
#
|
412
|
+
# magnify: The magnification factor of the small section. It can take any value greater than or equal to 1
|
413
|
+
def phi_magnifying_glass_graphkit
|
414
|
+
end
|
415
|
+
# The potential as a function of cartesian coordinates
|
416
|
+
# Options:
|
417
|
+
#
|
418
|
+
# rgbformulae: Gnuplot Option (may not apply when using other packages), sets colour mapping. See gnuplot help set rgbformulae
|
419
|
+
#
|
420
|
+
# limit: Limit the range of quantity begin plotted - any values of the quantity outside the limits will be set to the limit: eg. limit: [0,80]
|
421
|
+
#
|
422
|
+
# t_index: The (1-based) time index
|
423
|
+
def phi_real_space_graphkit
|
424
|
+
end
|
425
|
+
# The potential as a function of cartesian coordinates showing a cut at one toroidal angle, with multiple periodic copies of the flux tube used to fill the whole circle..
|
426
|
+
# Options:
|
427
|
+
#
|
428
|
+
# Rgeo:
|
429
|
+
#
|
430
|
+
# n0: The toroidal mode number of the longest y mode. In effect it is the number of periodic copies of the flux tube that will fit in the torus. Periodicity requires that n0 q is also an integer. If you specify :n0 where this is not the case, q will automatically be adjusted until it is
|
431
|
+
#
|
432
|
+
# rho_star: The ratio of the reference Lamour radius to the GS2 normalising length a. Cannot be specified at the same time as n0. If specified, both n0 and q will be adjusted to ensure periodicity
|
433
|
+
#
|
434
|
+
# t_index: The (1-based) time index
|
435
|
+
#
|
436
|
+
# nakx: The number of radial wave numbers to include in the plot. In effect, it is a low pass filter which reduces the resolution in the radial direction without changing the shape of the final surface. Minimum value is 4
|
437
|
+
#
|
438
|
+
# naky: The number of kys to include in the plot. In effect, it is a low pass filter which reduces the resolution in the y direction without changing the shape of the final surface. Minimum value is 4
|
439
|
+
#
|
440
|
+
# xmax: The (0-based) index of the maximum value of x to include in the plot
|
441
|
+
#
|
442
|
+
# xmin: The (0-based) index of the minimum value of x to include in the plot
|
443
|
+
#
|
444
|
+
# thetamax: The (0-based) index of the maximum value of theta to include in the plot
|
445
|
+
#
|
446
|
+
# thetamin: The (0-based) index of the minimum value of theta to include in the plot
|
447
|
+
#
|
448
|
+
# torphi:
|
449
|
+
def phi_real_space_poloidal_plane_graphkit
|
450
|
+
end
|
451
|
+
# The potential as a function of cartesian coordinates showing showing the standard way of representing the turbulence, with two poloidal cuts and the inner and outer radial surfaces. Multiple copies of the flux tube are used to fill the space.
|
452
|
+
# Options:
|
453
|
+
#
|
454
|
+
# Rgeo:
|
455
|
+
#
|
456
|
+
# n0: The toroidal mode number of the longest y mode. In effect it is the number of periodic copies of the flux tube that will fit in the torus. Periodicity requires that n0 q is also an integer. If you specify :n0 where this is not the case, q will automatically be adjusted until it is
|
457
|
+
#
|
458
|
+
# rho_star: The ratio of the reference Lamour radius to the GS2 normalising length a. Cannot be specified at the same time as n0. If specified, both n0 and q will be adjusted to ensure periodicity
|
459
|
+
#
|
460
|
+
# t_index: The (1-based) time index
|
461
|
+
#
|
462
|
+
# nakx: The number of radial wave numbers to include in the plot. In effect, it is a low pass filter which reduces the resolution in the radial direction without changing the shape of the final surface. Minimum value is 4
|
463
|
+
#
|
464
|
+
# naky: The number of kys to include in the plot. In effect, it is a low pass filter which reduces the resolution in the y direction without changing the shape of the final surface. Minimum value is 4
|
465
|
+
#
|
466
|
+
# xmax: The (0-based) index of the maximum value of x to include in the plot
|
467
|
+
#
|
468
|
+
# xmin: The (0-based) index of the minimum value of x to include in the plot
|
469
|
+
#
|
470
|
+
# thetamax: The (0-based) index of the maximum value of theta to include in the plot
|
471
|
+
#
|
472
|
+
# thetamin: The (0-based) index of the minimum value of theta to include in the plot
|
473
|
+
#
|
474
|
+
# torphi_values: An array of two values of the toroidal angle. The graph will be plotted in between those two values with poloidal cross sections at either end
|
475
|
+
def phi_real_space_standard_representation_graphkit
|
476
|
+
end
|
477
|
+
# The potential as a function of cartesian coordinates, plotted on the six outer surfaces of constant x, y and theta.
|
478
|
+
# Options:
|
479
|
+
#
|
480
|
+
# Rgeo:
|
481
|
+
#
|
482
|
+
# n0: The toroidal mode number of the longest y mode. In effect it is the number of periodic copies of the flux tube that will fit in the torus. Periodicity requires that n0 q is also an integer. If you specify :n0 where this is not the case, q will automatically be adjusted until it is
|
483
|
+
#
|
484
|
+
# rho_star: The ratio of the reference Lamour radius to the GS2 normalising length a. Cannot be specified at the same time as n0. If specified, both n0 and q will be adjusted to ensure periodicity
|
485
|
+
#
|
486
|
+
# t_index: The (1-based) time index
|
487
|
+
#
|
488
|
+
# nakx: The number of radial wave numbers to include in the plot. In effect, it is a low pass filter which reduces the resolution in the radial direction without changing the shape of the final surface. Minimum value is 4
|
489
|
+
#
|
490
|
+
# naky: The number of kys to include in the plot. In effect, it is a low pass filter which reduces the resolution in the y direction without changing the shape of the final surface. Minimum value is 4
|
491
|
+
#
|
492
|
+
# gs2_coordinate_factor: When set to 1, plot the graph in GS2 coordinates. When set to 0 plot the graph in real space. Can be set at any value between 0 and 1: the graph will smoothly distort between the two limits
|
493
|
+
#
|
494
|
+
# xmax: The (0-based) index of the maximum value of x to include in the plot
|
495
|
+
#
|
496
|
+
# xmin: The (0-based) index of the minimum value of x to include in the plot
|
497
|
+
#
|
498
|
+
# ymax: The (0-based) index of the maximum value of y to include in the plot
|
499
|
+
#
|
500
|
+
# ymin: The (0-based) index of the minimum value of y to include in the plot
|
501
|
+
#
|
502
|
+
# thetamax: The (0-based) index of the maximum value of theta to include in the plot
|
503
|
+
#
|
504
|
+
# thetamin: The (0-based) index of the minimum value of theta to include in the plot
|
505
|
+
#
|
506
|
+
# ncopies: The number of periodic copies of the flux tube to include
|
507
|
+
def phi_real_space_surface_graphkit
|
508
|
+
end
|
509
|
+
# Graph of phi^2 at a given time vs kx and ky
|
510
|
+
# Options:
|
511
|
+
#
|
512
|
+
# with: Gnuplot Option (may not apply when using other packages), e.g. with: 'lp' or with 'pm3d palette'
|
513
|
+
def spectrum_graphkit
|
514
|
+
end
|
515
|
+
# Graph of phi^2 * ky^2 at a given time vs kpar and ky
|
516
|
+
# Options:
|
517
|
+
#
|
518
|
+
# with: Gnuplot Option (may not apply when using other packages), e.g. with: 'lp' or with 'pm3d palette'
|
519
|
+
#
|
520
|
+
# log: Plot the log of a given quantity (exact meaning varies). boolean
|
521
|
+
#
|
522
|
+
# no_zonal: Don't plot the ky=0 part (boolean, e.g. no_zonal: true)
|
523
|
+
#
|
524
|
+
# no_kpar0: Don't plot the kpar=0 part (boolean, e.g. no_kpar0: true)
|
525
|
+
def spectrum_vs_kpar_vs_ky_graphkit
|
526
|
+
end
|
527
|
+
# transient_amplification_vs_ky or transient_amplification_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /transient_amplification_vs_kx_vs_ky/.
|
528
|
+
def transient_amplification_vs_kx_graphkit
|
529
|
+
end
|
530
|
+
# transient_amplification_vs_ky or transient_amplification_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /transient_amplification_vs_kx_vs_ky/.
|
531
|
+
def transient_amplification_vs_kxy_graphkit
|
532
|
+
end
|
533
|
+
# transient_amplification_vs_ky or transient_amplification_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /transient_amplification_vs_kx_vs_ky/.
|
534
|
+
def transient_amplification_vs_ky_graphkit
|
535
|
+
end
|
536
|
+
# transient_es_heat_flux_amplification_vs_ky or transient_es_heat_flux_amplification_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /transient_es_heat_flux_amplification_vs_kx_vs_ky/.
|
537
|
+
def transient_es_heat_flux_amplification_vs_kx_graphkit
|
538
|
+
end
|
539
|
+
# transient_es_heat_flux_amplification_vs_ky or transient_es_heat_flux_amplification_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /transient_es_heat_flux_amplification_vs_kx_vs_ky/.
|
540
|
+
def transient_es_heat_flux_amplification_vs_kxy_graphkit
|
541
|
+
end
|
542
|
+
# transient_es_heat_flux_amplification_vs_ky or transient_es_heat_flux_amplification_vs_kx. Growth rates vs either ky or kx for phi^2 integrated over the other direction. For growth rates at a specific kx AND ky, see /transient_es_heat_flux_amplification_vs_kx_vs_ky/.
|
543
|
+
def transient_es_heat_flux_amplification_vs_ky_graphkit
|
544
|
+
end
|
545
|
+
# Plots vspace diagnostics. All lines shouldn't stray much above 0.1 - otherwise large amounts of the distribution function is in the higher k velocity space and velocity space is probably unresolved. (NB This graph is here temporarily (ha ha) until I add the vspace diagnostics to the NetCDF file (or the apocalypse, whichever is sooner) EGH)
|
546
|
+
def vspace_diagnostics_graphkit
|
547
|
+
end
|
548
|
+
# zonal_spectrum: Graph of kx^4 phi^2 vs kx for ky=0
|
549
|
+
# Options:
|
550
|
+
#
|
551
|
+
# t: float, value of time at which to plot (e.g. t: 2.45)
|
552
|
+
#
|
553
|
+
# t_index: The (1-based) time index
|
554
|
+
def zonal_spectrum_graphkit
|
555
|
+
end
|
556
|
+
end
|