def METHOD_NAME():
    TestSection("DDF Spatial Collapse with Ghost Zones")
    DeleteAllPlots()
    OpenDatabase(silo_data_path("bigsil.silo"))
    DefineScalarExpression("mesh_x_nodal","coord(mesh)[0]")
    DefineScalarExpression("mesh_y_nodal","coord(mesh)[1]")
    DefineScalarExpression("mesh_z_nodal","coord(mesh)[2]")
    DefineScalarExpression("mesh_x_zonal","recenter(coord(mesh)[0])")
    DefineScalarExpression("mesh_y_zonal","recenter(coord(mesh)[1])")
    DefineScalarExpression("mesh_z_zonal","recenter(coord(mesh)[2])")
    AddPlot("Pseudocolor","dist")
    DrawPlots()
    atts = visit.ConstructDDFAttributes()
    atts.name = "ddf_dist_1d"
    atts.varForReductionOperator = "dist"
    atts.varnames = ("mesh_x_zonal",)
    atts.binBoundaries = (0,1)
    atts.numBins = (21,)
    ddf(atts,"dist","sum")
    atts = visit.ConstructDDFAttributes()
    atts.name = "ddf_dist_2d"
    atts.varForReductionOperator = "dist"
    atts.varnames = ("mesh_x_zonal", "mesh_y_zonal")
    atts.binBoundaries = (0,1, 0,1)
    atts.numBins = (21,21)
    ddf(atts,"dist","sum")
    # plot 1d result
    DeleteAllPlots()
    OpenDatabase("ddf_dist_1d.ult")
    AddPlot("Curve","dist_sum_1d")
    DrawPlots()
    Test("ddf_spatial_1d")
    # plot 2d result
    DeleteAllPlots()
    OpenDatabase("ddf_dist_2d.vtk")
    AddPlot("Pseudocolor","dist_sum_2d")
    DrawPlots()
    Test("ddf_spatial_2d")