Gyrooutput - Theory Group

GYRO Output

NEO Basics

NEO Input

NEO Output

TGYRO Basics

TGYRO Input

TGYRO Output

1 GYRO Data Output

GYRO Data Output

This section describes the contents of the major GYRO output files.

Poloidal output points

For $k = 0, \ldots, N_\theta -1$ :

$\theta_k = \begin{cases} 0 , & \mbox{if } N_\theta = 1 \\ -\pi + 2 \pi (k-1)/N_\theta , & \mbox{if } N_\theta > 1 \end{cases}$

Radial output points

Radial indices are stored in the vector r(n_x) (see profile_vugyro.out).

Array Dimensions

Reading the output files requires that the following dimensions are known to the reading program:

n_time: Number of discrete output times (one must read t.out to obtain this).
n_field: Number of potential fields (profile_vugyro.out)
n_spec: Number of species (profile_vugyro.out)
n_kinetic: Number of kinetic species (profile_vugyro.out)
n_x: Number of radial output points ( $N_r\,\!$ ; profile_vugyro.out)
n_theta_plot: Number of poloidal output points ( $N_\theta\,\!$ ; profile_vugyro.out)
n_n: Number of toroidal modes ( $N_n\,\!$ ; profile_vugyro.out)
n_energy: Number of energy gridpoints (profile_vugyro.out)
n_lambda = n_pass + n_trap: Number of pitch-angle gridpoints (profile_vugyro.out)

Output file structure

profile_vugyro.out

Format: Mixed scalars and vectors in one monolithic column of ASCII data.
Description: Equilibrium profile data, intended for use by vugyro.

n_x n_theta_section n_pass n_trap n_energy n_theta_plot n0 n_n d_n n_explicit_damp nonlinear_flag electron_method n_field n_ion n_kinetic n_spec field_r0_flag field_r0_grid n_grid_exp boundary_method r(n_x) q(n_x) r_s(n_x) q_s(n_x) dlntdr_s(n_spec,n_x) dlnndr_s(n_spec,n_x) tem_s(n_spec,n_x) den_s(n_spec,n_x) phi_doppler(n_x) aspect_s(n_x) delta_s(n_x) kappa_s(n_x) shift_s(n_x) shat_s(n_x) s_delta_s(n_x) s_kappa_s(n_x) beta_unit_s(n_x) pgamma_s(n_spec,n_x) b_unit_s(n_x) dr_eodr(n_x) grad_r_s(n_x) surf_hat_s(n_x) z_eff_s(n_x) nu_s(n_x) gamma_eb_s(n_x) er_exp_s(n_x) chi_gb_norm(n_x) chi_i_exp(n_x) chi_e_exp(n_x) diff_to_flow_e1(n_x) diff_to_flow_e2(n_x) eta_i_tot_exp(n_x) diff_to_flow_mi(n_x) aolvi_exp(n_x) diff_ne_exp(n_x) diff_to_flow_ne aolne_exp(n_x) diff_to_flow_heating(n_x) lambda(n_lambda) energy(n_energy) lambda_tp kt_rho(n_n) rho_s zcharge(n_spec) n_fine

geometry_arrays.out

Format: Rectangular array (8=i,n_fine,n_r)
Description: Radial and poloidal structure of geometry arrays; i=0 $:\nu\,\!$; i=1 $:\mbox{gsin}\,\!$; i=2 $:\mbox{gcos1}\,\!$; i=3 $:B\,\!$; i=4 $:J_\psi B\,\!$; i=5 $:|\nabla r|\,\!$; i=6 $:G_q\,\!$; i=7 $:\Theta\,\!$

t.out

Format: Rectangular array (2=i,n_time)
Description: Simulation time index and absolute time; i=0 $:\quad t/\Delta t\,\!$ (an integer); i=1 $:\quad ({\bar c}_s/a)t\,\!$

freq.out

Format: Rectangular array (4=i,n_n,n_time)
Description: Normalized mode frequency and growth rate; i=0 $:\quad (a/{\bar c}_s)\omega_{R,n}\,\!$; i=1 $:\quad (a/{\bar c}_s)\gamma_n\,\!$; i=2 $: \quad \mbox{error in} \; (a/{\bar c}_s)\omega_{R,n}\,\!$; i=3 $: \quad \mbox{error in} \; (a/{\bar c}_s)\gamma_n\,\!$

diff.out

Format: Rectangular array (n_kinetic,n_field,2=i,n_time)
Description: Gyrobohm-normalized diffusivities averaged over radius and summed over mode number.; i=0 $:\quad D_\sigma/\chi_{GB}\,\!$ (particle diffusivity); i=1 $:\quad \chi_\sigma/\chi_{GB}\,\!$ (energy diffusivity)

diff_i.out

Format: Rectangular array (n_kinetic,n_field,2=i,n_x,n_time)
Description: Gyrobohm-normalized diffusivities as a function of radius, summed over mode number.; i=0 $:\quad D_\sigma(r)/\chi_{GB}\,\!$ (particle diffusivity); i=1 $:\quad \chi_\sigma(r)/\chi_{GB}\,\!$ (energy diffusivity)

diff_n.out

Format: Rectangular array (n_kinetic,n_field,2=i,n_n,n_time)
Description: Gyrobohm-normalized diffusivities averaged over radius for each mode number.; i=0 $:\quad D_{\sigma,n}/\chi_{GB}\,\!$ (particle diffusivity); i=1 $:\quad \chi_{\sigma,n}/\chi_{GB}\,\!$ (energy diffusivity)

u.out

Format: Rectangular array (2,n_theta_plot,n_r,n_field=i_field,n_n,n_time).
Description: Potential expansion coefficients.; i_field=0 $:\quad \frac{e\delta \phi_n}{{\bar T}_e}\,\!$ (electrostatic potential); i_field=1 $:\quad \frac{{\bar c}_s}{c}\frac{e\delta A_{\lVert n}}{{\bar T}_e}\,\!$ (electromagnetic potential)

moment_n.out

Format: Rectangular array (2,n_theta_plot,n_r,n_kinetic,n_n,n_time).
Description: Density moment expansion coefficients $\frac{\delta n_{\sigma,n}}{{\bar n}_e}\,\!$ .

moment_e.out

Format: Rectangular array (2,n_theta_plot,n_r,n_kinetic,n_n,n_time).
Description: Energy moment expansion coefficients $\frac{\delta E_{\sigma,n}}{{\bar n}_e {\bar T}_e}\,\!$ .

flux_velocity.out

Format: Rectangular array (n_energy,n_lambda,n_kinetic,n_field=i_field,2=i,n_n,n_time).
Description: Velocity-space flux densities : $\quad \Gamma = \int d\varepsilon \int d\lambda\, \Gamma(\varepsilon,\lambda)\; , \quad Q = \int d\varepsilon \int d\lambda \, Q(\varepsilon,\lambda)$; i=0 $: \quad \Gamma_{\sigma,n}(\varepsilon,\lambda) \,\!$ (particle flux); i=1 $: \quad Q_{\sigma,n}(\varepsilon,\lambda) \,\!$ (energy flux); i_field=0 $: \quad~$ electrostatic component; i_field=1 $: \quad~$ electromagnetic component