def METHOD_NAME():
    TestSection("Save image plus depth")
    db = silo_data_path("multi_rect3d.silo")
    OpenDatabase(db)
    AddPlot("Contour", "u")
    DrawPlots()
    v = GetView3D()
    v.viewNormal = (-0.483301, 0.412135, 0.772376)
    v.focus = (0.5, 0.5, 0.5)
    v.viewUp = (0.218992, 0.911123, -0.349139)
    v.viewAngle = 30
    v.parallelScale = 0.866025
    v.nearPlane = -1.73205
    v.farPlane = 1.73205
    v.imagePan = (0, 0)
    v.imageZoom = 1
    v.perspective = 1
    v.eyeAngle = 2
    v.centerOfRotationSet = 0
    v.centerOfRotation = (0.5, 0.5, 0.5)
    v.axis3DScaleFlag = 0
    v.axis3DScales = (1, 1, 1)
    v.shear = (0, 0, 1)
    v.windowValid = 1
    SetView3D(v)
    # Save against a solid background.
    Test("pixeldata_1_00")
    SavePixelData("pdz0", 0, 1, 0, 0)
    # Save against a solid background adding axes
    a = GetAnnotationAttributes()
    a.axes3D.bboxFlag = 1
    SetAnnotationAttributes(a)
    Test("pixeldata_1_01")
    SavePixelData("pdz1", 0, 1, 0, 0)
    # Save against a gradient background 
    a.backgroundMode = a.Gradient
    a.gradientColor1 = (180, 100, 180, 255)
    a.gradientColor2 = (0, 0, 0, 255)
    a.axes3D.bboxFlag = 0
    a.foregroundColor = (255,255,255,255)
    SetAnnotationAttributes(a)
    Test("pixeldata_1_02")
    SavePixelData("pdz2", 0, 1, 0, 0)
    # Save against a gradient background adding axes
    a.backgroundMode = a.Gradient
    a.gradientColor1 = (180, 100, 180, 255)
    a.gradientColor2 = (0, 0, 0, 255)
    a.axes3D.bboxFlag = 1
    a.foregroundColor = (255,255,255,255)
    SetAnnotationAttributes(a)
    Test("pixeldata_1_03")
    SavePixelData("pdz3", 0, 1, 0, 0)
    DeleteAllPlots()
    CloseDatabase(db)
    # Open the files we made.
    idx = 4
    for db in ("pdz0.png", "pdz1.png", "pdz2.png", "pdz3.png"):
        OpenDatabase(db)
        AddPlot("Pseudocolor", "depth")
        pc = GetPlotOptions()
        pc.max = 220.
        pc.maxFlag = 1
        SetPlotOptions(pc)
        DrawPlots()
        v = GetView2D()
        v.viewportCoords = (0, 1, 0, 1)
        SetView2D(v)
        Test("pixeldata_1_%02d" % idx)
        DeleteAllPlots()
        CloseDatabase(db)
        idx = idx + 1