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

Data source: Toriumi & Takasao 2017, ApJ, 850, 39
https://ui.adsabs.harvard.edu/abs/2017ApJ...850...39T/abstract
Data and document: Shin Toriumi (ISAS/JAXA)
Contact: toriumi.shin@jaxa.jp
Date: April 5, 2023

IDL save files of 2D slices at z=0 from t=0 through 300
and 3D cubes above z=0 at limited cadence are provided.
A quadrupolar active region is created with a sheared PIL in between.

Units (non-dimensional) are
 Length: 170 km
 Magnetic field strength: 250 Gauss
 Velocity: 6.8 km/s

All data have been reduced to float to save the data size,
but the original simulation was done using double precision.

To read the Bz slices and plot magnetogram at final state (t=300),
IDL> restore,'2d/bz.sav'
IDL> restore,'2d/x.sav'
IDL> restore,'2d/y.sav'
IDL> contour,bz[*,*,300]>(-0.5)<0.5,x,y,/fill,nlevels=100,/iso,xs=1,ys=1

3D cubes for (Bx, By, Bz) are found at the following directory.
Note that ID of the file names is *half* of the actual time step.
Therefore, to read and plot the magnetogram (bottom boundary of Bz cube) at t=300,
IDL> restore,'3d/bz3d.0150.sav'
IDL> restore,'3d/x3d.sav'
IDL> restore,'3d/y3d.sav'
IDL> restore,'3d/z3d.sav'
IDL> contour,bz3d[*,*,0]>(-0.5)<0.5,x3d,y3d,/fill,nlevels=100,/iso,xs=1,ys=1

To plot the abosolute field strength on the vertical crosssection at t=300,
IDL> restore,'3d/bx3d.0150.sav'
IDL> restore,'3d/by3d.0150.sav'
IDL> b3d=sqrt(bx3d^2+by3d^2+bz3d^2)
IDL> contour,alog10(reform(b3d[*,256,*]>0.001<1.)),x3d,z3d $
IDL> ,xr=[-150,150],yr=[0,150],/fill,nlevels=100,/iso,xs=1,ys=1


#EOF