cubestats.f90: use double precision where necessary