;
; This is a program for calculating the filling factor defined as
; f=<b^2>^2/<b^4>
;
PRO filling_decay,bb,time,f=f,nx=nx
;
default,nx,256
;
b2=dot2(bb(3:nx+2,3:nx+2,3:nx+2,*))
b4=b2*b2

ngrid=(double(nx))^3
b2_mean=total(b2)/ngrid
b4_mean=total(b4)/ngrid

f=b2_mean*b2_mean/b4_mean
print,'f=',f

file='filling_decay_t='+string(time,FORMAT='(F7.2)')+'.sav'
save,file=file,f,time

END