def METHOD_NAME(self):
    mesh = discretize.TensorMesh([30, 30], x0=[-0.5, -1.0])
    sigma = np.random.rand(mesh.nC)
    model = np.log(sigma)
    # prob = DC.Simulation3DCellCentered(mesh, rhoMap=maps.ExpMap(mesh))
    # prob1 = DC.Simulation3DCellCentered(mesh, rhoMap=maps.ExpMap(mesh))
    rx = DC.Rx.Pole(
        utils.ndgrid([mesh.cell_centers_x, np.r_[mesh.cell_centers_y.max()]])
    )
    rx1 = DC.Rx.Pole(
        utils.ndgrid([mesh.cell_centers_x, np.r_[mesh.cell_centers_y.min()]])
    )
    src = DC.Src.Dipole(
        [rx],
        np.r_[-0.25, mesh.cell_centers_y.max()],
        np.r_[0.25, mesh.cell_centers_y.max()],
    )
    src1 = DC.Src.Dipole(
        [rx1],
        np.r_[-0.25, mesh.cell_centers_y.max()],
        np.r_[0.25, mesh.cell_centers_y.max()],
    )
    survey = DC.Survey([src])
    simulation0 = DC.simulation.Simulation3DCellCentered(
        mesh=mesh, survey=survey, rhoMap=maps.ExpMap(mesh)
    )
    survey1 = DC.Survey([src1])
    simulation1 = DC.simulation.Simulation3DCellCentered(
        mesh=mesh, survey=survey1, rhoMap=maps.ExpMap(mesh)
    )
    dobs0 = simulation0.make_synthetic_data(model)
    dobs1 = simulation1.make_synthetic_data(model)
    self.mesh = mesh
    self.model = model
    self.survey0 = survey
    self.sim0 = simulation0
    self.survey1 = survey1
    self.sim1 = simulation1
    # self.dmis0 = data_misfit.L2DataMisfit(self.survey0)
    self.dmis0 = data_misfit.L2DataMisfit(data=dobs0, simulation=simulation0)
    self.dmis1 = data_misfit.L2DataMisfit(data=dobs1, simulation=simulation1)
    self.dmiscombo = self.dmis0 + self.dmis1