def METHOD_NAME():
    ds_mesh, _ = _get_mesh()
    gf = GeometricFeatures()
    fc_mask = gf.read(componentName='ocean', objectType='region',
                      featureNames=['North Atlantic Ocean'])
    pool = _get_pool()
    ds_masks = compute_mpas_region_masks(
        ds_mesh, fc_mask, maskTypes=('cell', 'edge', 'vertex'), pool=pool)
    for dim in ['nCells', 'nEdges', 'nVertices', 'nRegions']:
        assert dim in ds_masks.dims
    for var in ['regionCellMasks', 'regionEdgeMasks', 'regionVertexMasks',
                'regionNames']:
        assert var in ds_masks.data_vars