; .r start
; .r r
;
; Find largest possible cube
;
nmin=min([nx,ny,nz])
l2=l1+nmin-1
m2=m1+nmin-1
n2=n1+nmin-1
nx=nmin & ny=nmin & nz=nmin
mx=l2+4 & my=m2+4 & mz=n2+4
aa_sub=aa(0:l2+nghostx,0:m2+nghosty,0:n2+nghostz,*)
uu_sub=uu(0:l2+nghostx,0:m2+nghosty,0:n2+nghostz,*)
lnrho_sub=lnrho(0:l2+nghostx,0:m2+nghosty,0:n2+nghostz)
;
; Calculate magnetic field
;
bb_sub=curl(aa_sub)
;
; Find mean values
;
bb_mean=fltarr(3)
uu_mean=fltarr(3)
bb_mean(0)=mean(bb_sub(l1:l2,m1:m2,n1:n2,0))
bb_mean(1)=mean(bb_sub(l1:l2,m1:m2,n1:n2,1))
bb_mean(2)=mean(bb_sub(l1:l2,m1:m2,n1:n2,2))
B0=(bb_mean(0)^2.+bb_mean(1)^2.+bb_mean(2)^2.)^(0.5)
uu_mean(0)=mean(uu_sub(l1:l2,m1:m2,n1:n2,0))
uu_mean(1)=mean(uu_sub(l1:l2,m1:m2,n1:n2,1))
uu_mean(2)=mean(uu_sub(l1:l2,m1:m2,n1:n2,2))
U0=(uu_mean(0)^2.+uu_mean(1)^2.+uu_mean(2)^2.)^(0.5)
;
; Find the magnetic field minus the mean
;
b=bb_sub & u=uu_sub
b(*,*,*,0)=bb_sub(*,*,*,0)-bb_mean(0)
b(*,*,*,1)=bb_sub(*,*,*,1)-bb_mean(1)
b(*,*,*,2)=bb_sub(*,*,*,2)-bb_mean(2)
u(*,*,*,0)=uu_sub(*,*,*,0)-uu_mean(0)
u(*,*,*,1)=uu_sub(*,*,*,1)-uu_mean(1)
u(*,*,*,2)=uu_sub(*,*,*,2)-uu_mean(2)
;
; Find rms values
;
brms=rms(b(l1:l2,m1:m2,n1:n2,*))
urms=rms(u(l1:l2,m1:m2,n1:n2,*))
;
; Print results
;
print,'The size of the box is',nmin,' mesh points in each direction'
print,'B0=',B0
print,'U0=',U0
print,'brms=',brms
print,'urms=',urms
END