def METHOD_NAME():
    """
    Test loading a recent gizmo snapshot that doesn't have cooling/metallicity
    The gizmo_zeldovich file has no metallicity field on the gas particles
    but is a cosmological dataset run using GIZMO_version=2022. There are two
    versions of the file, with GIZMO_version (_wg) and without GIZMO_version
    (_ng). Check that both load as gizmo datasets and correctly pull the
    cosmological variables. This test should get simpler when the file switches
    to pytest.
    """
    for fn in [zeld_wg, zeld_ng]:
        ds = yt.load(fn)
        assert isinstance(ds, GizmoDataset)
        assert ds.cosmological_simulation
        assert ds.omega_matter == 1.0
        assert ds.omega_lambda == 0.0
        # current_time is calculated from the cosmology so this checks if that
        # was calculated correctly
        assert_allclose_units(ds.current_time, 1672.0678 * Myr)