read the data with the idl routine read_slice.pro in the directory. Example: IDL>file='slice_0100.dat' IDL>read_slice,file,nph,nth,time,ph,lat,br,bth,bph,vr,vth,vph the data read in are: nph: dimension in phi nth: dimension in latitude time: time in sec ph(nph): phi grid (in rad) lat(nth): latitude grid (in rad) br(nph,nth): Br in G bth(nph,nth): Bth in G (positive southward) bph(nph,nth): Bph in G vr(nph,nth): Vr in cm/s vth(nph,nth): Vth in cm/s (positive southward) vph(nph,nth): Vph in cm/s I have written 101 snapshots at 3 hour intervals, tracking an area centered on the emerging region (tracked rotation rate is 477 nhz), with data (Br, Bth, Bph, Vr, Vth, Vph) written on a phi-latitude grid of about 34.4 degree in total widths (400Mm). The slice is at a raidus of r= 6.6605524E10 cm, about 30 Mm below the photosphere. A movie of the slices is movie_sm_50to150_10framespersec.mov in the current directory Note: the time (sec) is relative to some chosen time which is about the time the central main emerging region begins to emerge.