Using
  using [1,1,1] block setup for computation on domain [0:31:1,0:31:1,0:31:1]
  solving eeq
  time increments from [0, 1.79769e+308], cfl 0.5
  starting at t = 0, i = 1
  cell physical/total domain [0:30:1,0:30:1,0:30:1], [-2:32:1,-2:32:1,-2:32:1]
  face  physical/total domain [0:30:1,0:30:1,0:30:1], [-2:32:1,-2:32:1,-2:32:1]
  periodic boundaries in X Y Z
i = 1	 t = 0.00428322	 dt = 0.00428322i = 2	 t = 0.0081095	 dt = 0.00382628i = 3	 t = 0.011952	 dt = 0.00384251i = 4	 t = 0.0157823	 dt = 0.00383029Correctness:
	sum(rh) difference = 3.04226e-10
	sum(vx) = 0.0540387
	sum(vy) = 0.0505179
	sum(vz) = 0.0485836
	sum(rh*T) difference = -18.5667
exit 0