def initialize(self):

        my_grid = grid_setup(self.rp, ng=4)

        # create the variables
        my_data = patch.CellCenterData2d(my_grid)
        bc = bc_setup(self.rp)[0]
        my_data.register_var("density", bc)
        my_data.create()

        self.cc_data = my_data

        # now set the initial conditions for the problem
        exec(self.problem_name + '.init_data(self.cc_data, self.rp)')
Initialize the grid and variables for advection and set the initial
dataset/ETHPy150Open zingale/pyro2/advection/simulation.py/Simulation.initialize