; IDL Version 8.1, Mac OS X (darwin x86_64 m64)
; Journal File for Bob@steinr.pa.msu.edu
; Working directory: /Volumes/d2/MHDV/mhd12v/Slices
; Date: Thu Jun  9 10:54:31 2011
 
journ
restore,/ver,'int_mhd12_6km_226.save'
help,int
t=9
dx=12./2016.
print,dx
;  0.00595238
print,5/dx
;     840.000
int6=reform(int[2015-840:*,*,t])
window,xs=1800,ys=800
tvscl,congrid(int1,841/2,2016/2)
print,841*dx,2016*dx
;     5.00595      12.0000
restore,/ver,'int_mhd12_224.save'
help,int
;INT             FLOAT     = Array[1008, 1008, 20]
t=19
int12=reform(int[1008-841/2:*,*,t])
tvscl,congrid(int12,841/2,2016/2)
restore,/ver,'../../../mhdh/mhd48-1/Slices/int_48h1_1949.save'
 help,int
;    0.0238095
dx24=48./2016.
print,5/dx24
;     210.000
print,12/dx24
;     504.000
int24=reform(int[0:210,0:504,26])
help,int24
; INT24           FLOAT     = Array[211, 505]
restore,/ver,'../../../mhdh/mhd48-1/Slices/int_48h1_1843.save'
dx48=48./1008.
print,dx48
;    0.0476190
print,5/dx48
;      105.000
print,12/dx48
;      252.000
help,int
int48=reform(int[0:105,0:252,40])
ploghist,int48,/prob,n=500,yr=[1.e-7,.1],oplt=4
tvscl,congrid(int48,841/2,2016/2),2

print,min(int6,max=max),max
;    0.121109     0.558864
print,min(int24,max=max),max
;    0.130748     0.551583
print,min(int12,max=max),max
;    0.106737     0.407235
int6b=bytscl(int6,min=.13,max=.54)
int12b=bytscl(int12,min=.13,max=.54)
int24b=bytscl(int24,min=.13,max=.54)
int48b=bytscl(int48,min=.13,max=.54)
int96b=bytscl(int96,min=.13,max=.54)
int6bb=congrid(int6b,211,505)
int12bb=congrid(int12b,211,505)
int24bb=int24b
int48bb=congrid(int48b,211,505)
int96bb=congrid(int96b,211,505)
window,xs=1660,ys=800
iy=0
ix=0
panel,ix,iy,nx=5,ny=1,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1
plot_image,int6bb,[0,5],[0,12],xtitle='!5X (Mm)',ytit='!5Y (Mm)'
ix=1
panel,ix,iy,nx=5,ny=1,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1
plot_image,int12bb,[0,5],[0,12],xtitle='!5X (Mm)'
ix=2
panel,ix,iy,nx=5,ny=1,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1
plot_image,int24bb,[0,5],[0,12],xtitle='!5X (Mm)'
ix=3
panel,ix,iy,nx=5,ny=1,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1
plot_image,int48bb,[0,5],[0,12],xtitle='!5X (Mm)'
ix=4
panel,ix,iy,nx=5,ny=1,xdist=1,ydist=1,xleft=6,ybottom=4,xright=2,ytop=1
plot_image,int96bb,[0,5],[0,12],xtitle='!5X (Mm)'