def METHOD_NAME():
    setup_plot()
    AddOperator("DataBinning")
    datts = DataBinningAttributes()
    datts.numDimensions = datts.Two
    datts.dim1BinBasedOn = datts.Variable
    datts.dim1Var = "mesh_x_zonal"
    datts.dim1SpecifyRange = 0
    datts.dim1NumBins = 10
    datts.dim2BinBasedOn = datts.Variable
    datts.dim2Var = "mesh_y_zonal"
    datts.dim2SpecifyRange = 0
    datts.dim2NumBins = 10
    datts.outOfBoundsBehavior = datts.Clamp
    datts.reductionOperator = datts.Sum
    datts.varForReduction = "mass"
    datts.emptyVal = 0
    datts.outputType = datts.OutputOnBins
    SetOperatorOptions(datts)
    DrawPlots()
    # we have to export b/c we can't query the 
    # result of the operated created expr ...
    ofname  = "dbin_mass_sum_using_coords_exprs"
    eatts = ExportDBAttributes()
    eatts.db_type = "VTK"
    eatts.filename = ofname
    ExportDatabase(eatts)
    DeleteAllPlots()    
    dbin_varname = "%s_%s_%dd" % ("mass","sum",2)
    ofname += ".vtk"
    orig_vtk_var = "SCALARS %s float" % "operators/DataBinning"
    ddf_vtk_var  = "SCALARS %s float" % dbin_varname
    data = open(ofname).read()
    f = open(ofname, "w")
    data = data.replace(orig_vtk_var,ddf_vtk_var)
    f.write(data)
    f.close()
    OpenDatabase(ofname)
    AddPlot("Pseudocolor","mass_sum_2d")
    DrawPlots()
    Test("ddf_vs_dbinning_dbin_coords_exprs_result")
    res = query("Variable Sum")
    DeleteAllPlots()
    CloseDatabase(ofname)
    return res