geq Module

Name of the specific implementation

Has this instance been initialised

Does this implementation have the full or half theta grid? If only half then we assume up-down symmetry and use this to set gradients etc.

Typically the number of radial and theta grid points.

Typically the number of radial and theta grid points.

2D map of co-ordinates to R, Z or B

2D map of co-ordinates to R, Z or B

2D map of co-ordinates to R, Z or B

2D maps of theta and eqpsi

2D maps of theta and eqpsi

Minor radius gradient in cylindrical coordinates $(\hat{e}_R, \hat{e}_Z, \hat{e}_\zeta)$

Theta gradient in cylindrical coordinates $(\hat{e}_R, \hat{e}_Z, \hat{e}_\zeta)$

B gradient in cylindrical coordinates $(\hat{e}_R, \hat{e}_Z, \hat{e}_\zeta)$

Minor radius gradient in Bishop coordinates $(\hat{e}_{r_\psi}, \hat{e}_l, \hat{e}_\zeta)$

Theta gradient in Bishop coordinates $(\hat{e}_{r_\psi}, \hat{e}_l, \hat{e}_\zeta)$

B gradient in Bishop coordinates $(\hat{e}_{r_\psi}, \hat{e}_l, \hat{e}_\zeta)$

Equilibrium psi (not normalised) and pressure grids

Equilibrium psi (not normalised) and pressure grids

Normalised psi, fp, q, beta, diameter and R_centre values on grid

Normalised psi, fp, q, beta, diameter and R_centre values on grid

Normalised psi, fp, q, beta, diameter and R_centre values on grid

Normalised psi, fp, q, beta, diameter and R_centre values on grid

Normalised psi, fp, q, beta, diameter and R_centre values on grid

Normalised psi, fp, q, beta, diameter and R_centre values on grid

Calculate the derivative of rp w.r.t. R and Z and return the modulus sqrt(drp/dR ^ 2 + drp/dZ^2). I.e. return |grad rp|. Parameters are: - rgrid: Submodule radial grid (r_circ for EFIT, rp otherwise) - theta: grid of theta - gradf(:,1): |grad psi| - gradf(:,2): 0.0 - char: if char = 'R', return |grad pressure| instead - char: if char = 'T', then return theta gradient in bishop form - gradf(:,1): (dtheta/dZ d rp/dR - dtheta/dR drp/dZ)/ |grad rp| - gradf(:,2): (dtheta/dR d rp/dR + dtheta/dZ drp/dZ)/ |grad rp| - rp: value of rp on the flux surface where we want the grad - nth: number of theta points - ntgrid: lower index of theta array is -ntgrid