def METHOD_NAME():
    SetWindowLayout(2)
    fmt = "Blueprint"
    db = "shaping.root"
    for a in Argv():
        if ".mfem_root" in a:
            db = a
            fmt = "MFEM"
            print(f"Using {db}")
            break
        elif ".root" in a:
            db = a
            fmt = "Blueprint"
            print(f"Using {db}")
            break
    OpenDatabase(db)
    md = GetMetaData(db)
    # Look at the scalars to get the materials. Note that we sort the
    # material names in the end because the order of scalars in the file
    # could differ depending on the order in which the VFs were added to
    # the dataset.
    mats = []
    for i in range(md.GetNumScalars()):
        name = md.GetScalars(i).name
        if fmt == "MFEM":
            key = "vol_frac_"
            pos = name.find(key)
            if pos == 0:
                matname=name[pos+len(key):]
                #print(f"name={name}, pos={pos}, matname={matname}")
                mats.append(matname)
        else:
            key = "/vol_frac_"
            pos = name.find(key)
            if pos != -1:
                matname=name[pos+len(key):]
                #print(f"name={name}, pos={pos}, matname={matname}")
                mats.append(matname)
    mats = sorted(mats)
    SetActiveWindow(1)
    PlotVolumeFractions(mats, fmt)
    SetActiveWindow(2)
    DeleteAllPlots()
    PlotMaterials(mats, fmt)
    SetActiveWindow(1)
    swa = GetSaveWindowAttributes()
    swa.saveTiled = 1
    swa.width = 2048
    SetSaveWindowAttributes(swa)
    SaveWindow()
    sys.exit()