;------------------------------------------------------------------------------
; Example illustrating how to read MAS simulation data into IDL.
;
; *** NOTE: the ps_read_mas_run is not memory efficient for big runs.
;     - you may need something like 15Gb of RAM free while it does some
;       temporary operations. edit ps_read_mas_run if that is a problem.
;------------------------------------------------------------------------------

;--- define the path to where the simulation data is located or unzipped to
MhdDir = './corona'

;--- This is the Sequence number
Seq = '002'

;--- read in the mas snapshot (creates a structure of the full simulation data)
ps_read_mas_run, MhdDir, mas, Seq=Seq

;--- examine the structure
help, mas

;--- here is the annotated output of the help command
** Structure <3289a08>, 14 tags, length=5746962776, data length=5746962776, refs=1:
;   R               DOUBLE    Array[287]             [Rs,  1-30]  (radius coordinate)
;   T               DOUBLE    Array[326]             [Rad, 0-Pi]  (theta coordinate)
;   P               DOUBLE    Array[698]             [Rad, 0-2Pi] (phi coordinate)
;   RHO             DOUBLE    Array[287, 326, 698]   [cm-3]     (electron number density, ne)
;   TEMP            DOUBLE    Array[287, 326, 698]   [MK]       (electron temperature, Te)
;   VR              DOUBLE    Array[287, 326, 698]   [km/s]     (r component of flow)
;   VT              DOUBLE    Array[287, 326, 698]   [km/s]     (t component of flow)
;   VP              DOUBLE    Array[287, 326, 698]   [km/s]     (p component of flow)
;   BR              DOUBLE    Array[287, 326, 698]   [Gauss]    (r component of field)
;   BT              DOUBLE    Array[287, 326, 698]   [Gauss]    (t component of field)
;   BP              DOUBLE    Array[287, 326, 698]   [Gauss]    (p component of field)
;   JR              DOUBLE    Array[287, 326, 698]   [A/m2]     (r component of current)
;   JT              DOUBLE    Array[287, 326, 698]   [A/m2]     (t component of current)
;   JP              DOUBLE    Array[287, 326, 698]   [A/m2]     (p component of current)