gs2crmod 0.5.2

Sign up to get free protection for your applications and to get access to all the features.
Files changed (45) hide show
  1. data/.document +5 -0
  2. data/Gemfile +13 -0
  3. data/LICENSE.txt +20 -0
  4. data/README.md +4 -0
  5. data/README.rdoc +19 -0
  6. data/Rakefile +56 -0
  7. data/VERSION +1 -0
  8. data/ext/extconf.rb +9 -0
  9. data/ext/gs2crmod_ext.c +366 -0
  10. data/gs2crmod.gemspec +98 -0
  11. data/include/gs2crmod_ext.h +58 -0
  12. data/lib/gs2crmod/astrogk/astrogk.rb +201 -0
  13. data/lib/gs2crmod/astrogk/calculations.rb +57 -0
  14. data/lib/gs2crmod/astrogk/check_convergence.rb +7 -0
  15. data/lib/gs2crmod/astrogk/deleted_variables.rb +76 -0
  16. data/lib/gs2crmod/astrogk/graphs.rb +13 -0
  17. data/lib/gs2crmod/astrogk/gsl_data.rb +13 -0
  18. data/lib/gs2crmod/astrogk/gsl_tools.rb +182 -0
  19. data/lib/gs2crmod/astrogk/ingen.rb +18 -0
  20. data/lib/gs2crmod/astrogk/input_file_tools.rb +7 -0
  21. data/lib/gs2crmod/astrogk/namelist_tools.rb +14 -0
  22. data/lib/gs2crmod/astrogk/namelists.rb +2800 -0
  23. data/lib/gs2crmod/astrogk/properties.rb +17 -0
  24. data/lib/gs2crmod/astrogk/species_dependent_namelists.rb +228 -0
  25. data/lib/gs2crmod/astrogk/test_gs2.rb +231 -0
  26. data/lib/gs2crmod/astrogk.rb +200 -0
  27. data/lib/gs2crmod/calculations.rb +780 -0
  28. data/lib/gs2crmod/check_convergence.rb +179 -0
  29. data/lib/gs2crmod/deleted_variables.rb +916 -0
  30. data/lib/gs2crmod/graphs.rb +1899 -0
  31. data/lib/gs2crmod/graphs_rdoc.rb +556 -0
  32. data/lib/gs2crmod/gs2.rb +1143 -0
  33. data/lib/gs2crmod/gsl_data.rb +1181 -0
  34. data/lib/gs2crmod/gsl_data_3d.rb +705 -0
  35. data/lib/gs2crmod/gsl_tools.rb +187 -0
  36. data/lib/gs2crmod/ingen.rb +218 -0
  37. data/lib/gs2crmod/namelists.rb +5142 -0
  38. data/lib/gs2crmod/properties.rb +22 -0
  39. data/lib/gs2crmod/species_dependent_namelists.rb +228 -0
  40. data/lib/gs2crmod/test_gs2.rb +231 -0
  41. data/lib/gs2crmod.rb +2 -0
  42. data/lib/gs2crmod_extension.rb +1 -0
  43. data/test/helper.rb +18 -0
  44. data/test/test_gs2crmod.rb +7 -0
  45. 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