def METHOD_NAME(self, *args, **options):
    """
    Find neighborhoods in requested municipality store them as a dataset.
    Note: storing as a dataset in this context means that a we derive
    area type names for the neighborhoods in the municipality under
    consideration. We store a copy of all relevant neighborhoods (both the
    Dutch "wijk" and "buurt" types of neighborhood) associated with the
    newly created AreaTypes.
    """
    cbs_wijk_dataset = REQUIRED_DATASETS['CBS_WIJK_DATASET']
    cbs_buurt_dataset = REQUIRED_DATASETS['CBS_BUURT_DATASET']
    municipal_code = options['municipal_code']
    if not self._cbs_data_present():
        raise CommandError('Required datasets are not available, exiting!')
    if not self._municipal_code_is_correct(municipal_code):
        raise CommandError('Provided municipal code is not present in CBS data, exiting!')
    self.stdout.write('Required CBS data are available.')
    # Datatype for our newly loaded data, note we drop the old copy.
    wijk_dataset_code, buurt_dataset_code = self._derive_dataset_codes(municipal_code)
    wijk_area_type, _ = AreaType.objects.get_or_create(code=wijk_dataset_code, name=wijk_dataset_code)
    buurt_area_type, _ = AreaType.objects.get_or_create(code=buurt_dataset_code, name=buurt_dataset_code)
    with transaction.atomic():
        to_remove_w = Area.objects.filter(_type=wijk_area_type)
        to_remove_b = Area.objects.filter(_type=buurt_area_type)
        self.stdout.write(f'\nThere are {to_remove_w.count()} old (wijk) area instances to remove')
        self.stdout.write(f'There are {to_remove_b.count()} old (buurt) area instances to remove')
        to_remove_w.delete()
        to_remove_b.delete()
        # We duplicate the neighborhoods in the municipality we are interested in, we select the geometry for
        # the municipality and use it to find the neighborhoods ("wijk" and "buurt") we want to copy.
        municipality_area_type = AreaType.objects.get(code=REQUIRED_DATASETS['CBS_MUNICIPAL_BORDERS_DATASET'])
        municipality = Area.objects.get(_type=municipality_area_type, code=municipal_code)
        self.stdout.write(f'\nNeighborhoods for the municipality of "{municipality.name}" will be copied')
        to_copy_w = Area.objects.filter(_type__code=cbs_wijk_dataset, geometry__intersects=municipality.geometry)
        to_copy_b = Area.objects.filter(_type__code=cbs_buurt_dataset, geometry__intersects=municipality.geometry)
        self.stdout.write(f'There are {to_copy_w.count()} new "wijk" area instances to copy.')
        self.stdout.write(f'"Wijk" area instances will be copied to dataset/AreaType {wijk_area_type.name}.')
        self.stdout.write(f'There are {to_copy_b.count()} new "buurt area instances to copy.')
        self.stdout.write(f'"Buurt" area instances will be copied to dataset/AreaType {buurt_area_type.name}.')
        self._copy_areas_to_area_type(to_copy_w, wijk_area_type)
        self._copy_areas_to_area_type(to_copy_b, buurt_area_type)
        self.stdout.write('Done.')