Write a script to compute admin area of each DRA locality

Create a python script called somthing like findDraLocalityAdminAreas.py

The script should read the following inputs:
1. DRA locality coverage containing polygons for each dra locality.
2. Admin Area coverage containing polygons for a given set of admin areas (e.g., Community Health Service Areas).
3. A parameter called adminAreaType which contains the official code of the set of admin areas in the coverage (e.g., "CHSA")

The script should perform a spatial overlay of each DRA locality polygon against the admin area coverage to determine the list of intersecting admin areas.

Script should output the following for each dra locality:

1. localityName
2. adminAreaType
3. adminAreaCode

If there is more than one admin area intersecting a locality, set adminAreaCode to a special character like asterisk ("*") to represent "more than one". For example, if the dra locality polygon for Victoria intersects five community health service areas, set Victoria's adminAreaCode to "*", If the dra locality polygon for Ahousat intersects a single CHSA called "NBAH", set Ahousat's adminAreaCode to "NBAH"
If there is no intersecting admin area for a given dra locality polygon, set adminAreaCode to the empty string ("")

Start by locating existing Python scripts and the repository's conventions for reading coverage inputs and writing results. Implement the named script around the three stated inputs and spatial overlay, then verify outputs for localities with multiple, single, and no intersecting admin areas.

Written by the indexing model from the issue text.