ISEE workshop "Magnetic Flux Rope on the Sun"
led by K.D. Leka and B. Kliem

Data source: Toriumi & Hotta 2019, ApJL, 886, L21
https://ui.adsabs.harvard.edu/abs/2019ApJ...886L..21T/abstract
Data and document: Shin Toriumi (ISAS/JAXA)
Contact: toriumi.shin@jaxa.jp
Date: April 5, 2023

IDL save files of 2D slices at tau=1 and 3D cubes
at two time steps (35.0 hr and 42.3 hr) are provided.
The developed spots with a sheared PIL and
flux-rope-like structure are found in the t=42.3 hr data.

Units are
 Length: Mm
 Magnetic field strength: Gauss
 Velocity: cm/s

To read the Bz slice for t=42.3 hr and plot magnetogram,
IDL> restore,'t42.3/tau1/bz.sav'
IDL> restore,'t42.3/tau1/x.sav'
IDL> restore,'t42.3/tau1/y.sav'
IDL> contour,bz>(-500)<500,x,y,/fill,nlevels=100,/iso,xs=1,ys=1

Along with the (Bx, By, Bz) and (Vx, Vy, Vz) vectors, intensity map
is also available.
IDL> restore,'t42.3/tau1/rd.sav'
IDL> rd0=mean(rd[0:199,0:199])
IDL> contour,(rd/rd0)>0.5<1.5,x,y,/fill,nlevels=100,/iso,xs=1,ys=1

3D cubes for (Bx, By, Bz) are found at,
IDL> restore,'t42.3/3d/bz3d.sav'
IDL> restore,'t42.3/3d/x.sav'
IDL> restore,'t42.3/3d/y.sav'
IDL> restore,'t42.3/3d/z.sav'

The z=0 Mm level is at (slightly below) 241th grid, so
IDL> k=241
IDL> print,z[k]
IDL> contour,bz3d[*,*,k]>(-500)<500,x,y,/fill,nlevels=100,/iso,xs=1,ys=1

#EOF