logo

Reading Output Into IDL

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'.

Back