To read TRANSP output from the netCDF file into an IDL session one can use the routine ~dkeel/IDL/CDF_read/read_transp_jet.pro
written by Dave Keeling
if one is on the Jet Analysis Cluster.
If the TRANSP run has been archived so that the netCDF file is stored in the data warehouse then this routine can be called as follows:
data = read_transp_jet('runID','signal')
Where runID
is the run ID e.g. 84543K01 and signal
is the signal name e.g. ne for the electron density. This call will return a data structure with the following
fields:
NAME | STRING | 'NE' |
INFO | STRING | 'Data mapped against X' |
DATA | FLOAT | Array[50, 401] |
X | FLOAT | Array[50, 401] |
T | FLOAT | Array[401] |
UNITS | STRING | 'N/CM**3 ' |
RADTYPE | STRING | 'X' |
IFLAG | INT | -1 |
DIMENSIONS | LONG | 2 |
NT | LONG | 401 |
NX | LONG | 50 |
Most of these variables are self explanatory, RADTYPE
indicates the radial coordinate used for this quantity, some quantities, for instance profiles, are stored at the grid centre
points in which case this will have the value 'X'
whereas others, for instance fluxes, will be stored at the edges of the grid cells in which case this variable will have the
value 'XB'
.