PRO movie,device=device if (n_elements(device) eq 0) then device='x' if (device eq 'z') then begin set_plot,'z' device,set_resolution=[600,619] endif else begin set_plot,'x' window,xsize=600,ysize=619 device,pseudo_color=8,decompose=0 endelse help,/device print,'Windows done' restore,/verb,'intt_bbt_h5.save' s=size(intt) nt=s[3] print,'nt=',nt print,'times:',min(times,max=max),max !p.background=255 !p.thick=2 !x.thick=2 !y.thick=2 !p.charsize=1.5 !p.charthick=1.7 !x.margin=[6,2] hlim,intt,mini,maxi,cut=.99 print,'Int:',min(mini),max(maxi) inttb=bytscl(intt,min=min(mini),max=max(maxi)) plot,maxi oplot,mini hlim,bxt,minx,maxx,cut=.99 print,'Bx:',min(minx),max(maxx) plot,maxx hlim,byt,miny,maxy,cut=.99 print,'By:',min(miny),max(maxy) plot,maxy hlim,bzt,minz,maxz,cut=.99 print,'Bz:',min(minz),max(maxz) plot,maxz minb=min([min(minx),min(miny),min(minz)]) maxb=max([max(maxx),max(maxy),max(maxz)]) print,'B:',minb,maxb Bxtb=bytscl(Bxt,min=min(minb),max=max(maxb)) Bytb=bytscl(Byt,min=min(minb),max=max(maxb)) Bztb=bytscl(Bzt,min=min(minb),max=max(maxb)) for t=0,nt-1 do begin ix=0 iy=0 panel,ix,iy,nx=2,ny=2,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1 plot_image,inttb[*,*,t],[0,48],[0,48],ytit='!5Y (Mm)' iy=1 panel,ix,iy,nx=2,ny=2,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1 plot_image,bxtb[*,*,t],[0,48],[0,48],xtitle='!5X (Mm)',ytit='!5Y (Mm)' ix=1 iy=0 panel,ix,iy,nx=2,ny=2,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1 plot_image,bytb[*,*,t],[0,48],[0,48] iy=1 panel,ix,iy,nx=2,ny=2,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1 plot_image,bztb[*,*,t],[0,48],[0,48],xtitle='!5X (Mm)' xyouts,9,7,string(times[t]/.6/60.,'(f5.2)')+' Hrs',color=2,/device im=tvrd() if (t le 9) then begin write_jpeg,'intbb4.h5.00'+string(t,format='(i1)')+'.jpg',im endif else if (t le 99) then begin write_jpeg,'intbb4.h5.0'+string(t,format='(i2)')+'.jpg',im endif else if (t le 999) then begin write_jpeg,'intbb4.h5.'+string(t,format='(i3)')+'.jpg',im endif else begin write_jpeg,'intbb4.h5.'+string(t,format='(i4)')+'.jpg',im endelse wait,.1 endfor return end