after bounce ($t\equiv 0$) of our model 3D-H-1). The contours on the cross sections in the $x = 0$ (back right), $y = 0$ (back bottom), and $z = 0$ (back left) planes are, respectively, projected on the sidewalls of the graphs. For each snapshot, the length of white line is indicated at right bottom text.[]{data-label="fig1"}](f1d.eps "fig:"){width=".45\linewidth"} ![Three dimensional plots of entropy per baryon (top panel), $\tau_{\rm res}/\tau_{\rm heat}$ (bottom left panel) that is the ratio of the residency to the neutrino heating timescale (see the text for details), and the net neutrino heating rate (bottom right panel, in unit of ${\rm erg}~{\rm cm}^{-3}~{\rm s}^{-1}$) for three snapshots (top and bottom left: $t=230$ ms, and bottom right: $t=150$ ms measured after bounce ($t\equiv 0$) of our model 3D-H-1). The contours on the cross sections in the $x = 0$ (back right), $y = 0$ (back bottom), and $z = 0$ (back left) planes are, respectively, projected on the sidewalls of the graphs. For each snapshot, the length of white line is indicated at right bottom text.[]{data-label="fig1"}](f1b.eps "fig:"){width=".45\linewidth"}

![image](f2b.eps){width=".28\linewidth"} ![image](f2d.eps){width=".28\linewidth"} ![image](f2new.eps){width=".4\linewidth"}

Numerical Methods and Models
============================

Here we briefly summarize several major updates of the code that we have implemented after 