Output-profiles.md

+The output 1D profiles are written in a plain ASCII text file, with blank separated values in tabular format.
+
+The first 21 lines document the set of parameters used for the run. Line 22 contains the data column headers. All the 22 header lines start with a leading `#` character. Actual data start at line number 23.
+
+The file is written by default on unit=48 (unit=648 in JINTRAC).
+
+The number of points in the profiles are set with the [parameter](code-parameters) `nrho`, while the [parameter](code-parameters) `ipec` flags how the radial grid is spaced:
+- `ipec`=1: even spacing in $`\rho_{\rm pol}`$
+- `ipec`=2: even spacing in $`\rho_{\rm tor}`$
+
+The content of each column in the file is given here below.
+
+| Column | Label    |      Variable      | Units  | Description                     |
+|   ---: | :---     |        :---:       | :---:  | :----------------------   |
+|      1 | `rhop `  | $`\rho_{\rm pol}`$ |   -    | Normalized minor radius   |
+|      2 | `rhot`   | $`\rho_{\rm tor}`$ |   -    | Normalized minor radius   |
+|      3 | `Jphi`   | $`J_\phi`$     | MA/m$`^2`$ | Driven toroidal current density        |
+|      4 | `Jcdb`   | $`J_{\rm cd}`$ | MA/m$`^2`$ | Driven current density, JETTO definition     |
+|      5 | `dPdV`   | $`dP/dV`$      | MW/m$`^3`$ | Absorbed power density     |
+|      6 | `Icdins` | $`I_{\rm cd,ins}`$ |   MA   | Cumulative integral of driven current        |
+|      7 | `Pins`   | $`P_{\rm ins}`$    |   MW   | Cumulative integral of absorbed power |
+
+> All the values in the output files are given in COCOS=3, independently of the `icocos` value set in the [code parameters](code-parameters).
+
+## Definitions
+
+- Normalized minor radii
+    ```math
+    \rho_{\rm pol} = \sqrt{\frac{\psi - \psi_0}{\psi_a - \psi_0}}, \qquad
+    \rho_{\rm tor} = \sqrt{\frac{\Phi}{\Phi_a}}
+    ```
+
+- Flux surface average
+
+    ```math
+    \langle x \rangle = \left. \oint x \frac{dl}{B_p} \middle/
+                               \oint \frac{dl}{B_p} \right.
+    ```
+- Flux-averaged parallel current density
+
+    ```math
+    \langle J_{\parallel}\rangle = \left \langle
+                                     \frac{{\bf J} \cdot {\bf B}}{B}
+                                   \right \rangle
+       = \frac{\langle {{\bf J} \cdot {\bf B}}\rangle}
+              {{\langle B^2 \rangle/}{\langle B \rangle}}.
+    ```
+
+- Toroidal current density
+
+    ```math
+    J_{\phi} = \frac{f(\psi)}{\langle B \rangle}
+               \frac{\langle {1/R^2} \rangle}{\langle{1/R} \rangle}
+               {\langle J_\parallel \rangle }
+    ```
+- JETTO definition of the driven current density
+
+    ```math
+    J_{\rm cd} = \frac{\langle {\bf J} \cdot {\bf B} \rangle}
+                      {\langle {\bf J} \rangle}
+    ```
+- Driven current cumulative integral
+
+    ```math
+    I_{\mathrm{cd,ins}}(\rho) = \int_0^\rho J_\phi \frac{dA}{d\rho'} d\rho'
+    ```
+
+- Absorbed power cumulative integral
+    ```math
+    P_{\mathrm{ins}}(\rho) = \int_0^\rho \frac{dP}{dV} \frac{dV}{d\rho'} d\rho'
+    ```
+
+Where:
+- $`\Phi`$ and $`\psi`$ are the toroidal and poloidal magnetic flux respectively
+- $`\Phi_a`$ and $`\psi_a`$ their values at the last closed flux surface
+- $`\Phi_0=0`$ and $`\psi_0`$ their values value at the plasma magnetic axis
+- $`R`$ the major radius
+- $`\phi`$ the toroidal angle
+- $`{\bf J}`$ the EC driven current density vector
+- $`{\bf B}`$ the magnetic field
+- $`B_\phi`$ and $`B_p`$ its toroidal and poloidal components
+- $`f(\psi) = B_\phi R`$ the poloidal flux function
+- $`V(\rho)`$ the volume enclosed by the flux surface at $`\rho`$
+- $`A(\rho)`$ its cross-section area
\ No newline at end of file