def METHOD_NAME(modelHoldupNotDynamic):
    slab = modelHoldupNotDynamic.fs.slab
    nx = len(slab.ixnodes)
    nz = len(slab.iznodes)
    nt = len(slab.flowsheet().time)
    comp_dict = common_components(nt, nz, nx)
    comp_dict[pyo.Var]["int_energy_density_solid"] = nx * nz * nt
    comp_dict[pyo.Param] = {}
    comp_dict[pyo.Param]["dcedt_solid"] = nt * nx * nz
    comp_dict[pyo.Constraint]["int_energy_density_solid_eqn"] = nt * nx * nz
    soc_testing._build_test_utility(
        slab,
        comp_dict=comp_dict,
        references=[
            "temperature_z",
            "current_density",
            "length_z",
            "length_y",
            "temperature_deviation_x0",
            "heat_flux_x0",
            "temperature_deviation_x1",
            "heat_flux_x1",
        ],
    )
    assert degrees_of_freedom(slab) == 0