Source code for pyro.burgers_viscous.problems.tophat

"""Initialize a circular "tophat" region that will drive a shock."""

from pyro.util import msg

DEFAULT_INPUTS = "inputs.tophat"




def init_data(myd, rp):
    """ initialize the tophat burgers problem """

    if rp.get_param("driver.verbose"):
        msg.bold("initializing the tophat burgers problem...")

    u = myd.get_var("x-velocity")
    v = myd.get_var("y-velocity")

    xmin = myd.grid.xmin
    xmax = myd.grid.xmax

    ymin = myd.grid.ymin
    ymax = myd.grid.ymax

    xctr = 0.5*(xmin + xmax)
    yctr = 0.5*(ymin + ymax)

    u[:, :] = 0.0
    v[:, :] = 0.0

    R = 0.1

    inside = (myd.grid.x2d - xctr)**2 + (myd.grid.y2d - yctr)**2 < R**2

    u[inside] = 1.0
    v[inside] = 1.0





def finalize():
    """ print out any information to the user at the end of the run """