From 938307853368158002ea0e59386e9efe5a7b078f Mon Sep 17 00:00:00 2001
From: Jamie Taylor <jamie.taylor@sheffield.ac.uk>
Date: Wed, 22 Mar 2023 17:48:26 +0000
Subject: [PATCH] 0.12.0

- Add support for reverse-geocoding to NUTS regions
- Use pathlib.Path throughout when working with files
- Fix bugs in gmaps.py
---
 .gitignore                                    |   1 +
 Docs/build/html/_modules/geocode/geocode.html |  70 +++++----
 README.md                                     | 116 +++++++--------
 Tests/test_geocode.py                         |  41 +++++-
 docs/build/doctrees/environment.pickle        | Bin 44958 -> 46674 bytes
 docs/build/doctrees/modules.doctree           | Bin 53378 -> 67043 bytes
 docs/build/html/.buildinfo                    |   2 +-
 docs/build/html/_modules/index.html           |   2 +-
 .../html/_static/documentation_options.js     |   2 +-
 docs/build/html/genindex.html                 |   8 +-
 docs/build/html/index.html                    |   2 +-
 docs/build/html/modules.html                  |  23 ++-
 docs/build/html/objects.inv                   | Bin 400 -> 406 bytes
 docs/build/html/search.html                   |   2 +-
 docs/build/html/searchindex.js                |   2 +-
 docs/source/conf.py                           |   2 +-
 example.py                                    |  42 +++---
 geocode/bng2latlon.py                         |  11 +-
 geocode/cache_manager.py                      |  21 +--
 geocode/cpo.py                                |  18 +--
 geocode/eurostat.py                           | 136 ++++++++++++++++++
 geocode/geocode.py                            |  66 +++++----
 geocode/gmaps.py                              |  28 ++--
 geocode/latlons2gsp.py                        |  11 +-
 geocode/latlons2llsoa.py                      |  11 +-
 geocode/ngeso.py                              |  17 +--
 geocode/ons_nrs.py                            |  59 +++-----
 geocode/postcodes2latlon.py                   |  11 +-
 geocode/utilities.py                          |   7 +-
 geocode/version.py                            |   2 +-
 setup.py                                      |   4 +-
 31 files changed, 468 insertions(+), 249 deletions(-)
 create mode 100644 geocode/eurostat.py
diff --git a/.gitignore b/.gitignore
index f168604..5245a68 100644
--- a/.gitignore
+++ b/.gitignore
@@ -106,3 +106,4 @@ venv.bak/
 # Project specific
 *.p
 geocode/google_maps/key.txt
+geocode/cache/key.txt
diff --git a/Docs/build/html/_modules/geocode/geocode.html b/Docs/build/html/_modules/geocode/geocode.html
index a7c30b9..6533c70 100644
--- a/Docs/build/html/_modules/geocode/geocode.html
+++ b/Docs/build/html/_modules/geocode/geocode.html
@@ -3,7 +3,7 @@
 <head>
   <meta charset="utf-8" />
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
-  <title>geocode.geocode &mdash; Geocode 0.11.1 documentation</title>
+  <title>geocode.geocode &mdash; Geocode 0.12.0 documentation</title>
       <link rel="stylesheet" href="../../_static/pygments.css" type="text/css" />
       <link rel="stylesheet" href="../../_static/css/theme.css" type="text/css" />
   <!--[if lt IE 9]>
@@ -80,42 +80,33 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
 <span class="kn">import</span> <span class="nn">os</span>
 <span class="kn">import</span> <span class="nn">sys</span>
 <span class="kn">import</span> <span class="nn">logging</span>
+<span class="kn">from</span> <span class="nn">pathlib</span> <span class="kn">import</span> <span class="n">Path</span>
 <span class="kn">import</span> <span class="nn">time</span> <span class="k">as</span> <span class="nn">TIME</span>
 <span class="kn">import</span> <span class="nn">argparse</span>
-<span class="kn">from</span> <span class="nn">typing</span> <span class="kn">import</span> <span class="n">Optional</span>
 <span class="kn">from</span> <span class="nn">shutil</span> <span class="kn">import</span> <span class="n">copyfile</span>
+<span class="kn">from</span> <span class="nn">typing</span> <span class="kn">import</span> <span class="n">Optional</span><span class="p">,</span> <span class="n">Iterable</span><span class="p">,</span> <span class="n">Tuple</span><span class="p">,</span> <span class="n">Union</span><span class="p">,</span> <span class="n">List</span><span class="p">,</span> <span class="n">Dict</span><span class="p">,</span> <span class="n">Literal</span>
 
-<span class="kn">from</span> <span class="nn">typing</span> <span class="kn">import</span> <span class="n">Optional</span><span class="p">,</span> <span class="n">Iterable</span><span class="p">,</span> <span class="n">Tuple</span><span class="p">,</span> <span class="n">Union</span><span class="p">,</span> <span class="n">List</span><span class="p">,</span> <span class="n">Dict</span>
 <span class="kn">import</span> <span class="nn">pyproj</span>
 
-<span class="k">try</span><span class="p">:</span>
-    <span class="kn">from</span> <span class="nn">shapely.geometry</span> <span class="kn">import</span> <span class="n">shape</span><span class="p">,</span> <span class="n">Point</span>
-    <span class="kn">from</span> <span class="nn">shapely.ops</span> <span class="kn">import</span> <span class="n">unary_union</span>
-    <span class="n">SHAPELY_AVAILABLE</span> <span class="o">=</span> <span class="kc">True</span>
-<span class="k">except</span> <span class="ne">ImportError</span><span class="p">:</span>
-    <span class="n">logging</span><span class="o">.</span><span class="n">warning</span><span class="p">(</span><span class="s2">&quot;Failed to import Shapely library - you will not be able to reverse-geocode! &quot;</span>
-                    <span class="s2">&quot;See notes in the README about installing Shapely on Windows machines.&quot;</span><span class="p">)</span>
-    <span class="n">SHAPELY_AVAILABLE</span> <span class="o">=</span> <span class="kc">False</span>
-
-<span class="n">SCRIPT_DIR</span> <span class="o">=</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">dirname</span><span class="p">(</span><span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">realpath</span><span class="p">(</span><span class="vm">__file__</span><span class="p">))</span>
-<span class="c1"># sys.path.append(SCRIPT_DIR)</span>
-
 <span class="kn">from</span> <span class="nn">.</span> <span class="n">utilities</span> <span class="kn">import</span> <span class="nn">GenericException</span>
 <span class="kn">from</span> <span class="nn">.</span> <span class="n">cpo</span> <span class="kn">import</span> <span class="nn">CodePointOpen</span>
 <span class="kn">from</span> <span class="nn">.</span> <span class="n">ngeso</span> <span class="kn">import</span> <span class="nn">NationalGrid</span>
 <span class="kn">from</span> <span class="nn">.</span> <span class="n">ons_nrs</span> <span class="kn">import</span> <span class="nn">ONS_NRS</span>
+<span class="kn">from</span> <span class="nn">.</span> <span class="n">eurostat</span> <span class="kn">import</span> <span class="nn">Eurostat</span>
 <span class="kn">from</span> <span class="nn">.</span> <span class="n">gmaps</span> <span class="kn">import</span> <span class="nn">GMaps</span>
 <span class="kn">from</span> <span class="nn">.</span> <span class="n">cache_manager</span> <span class="kn">import</span> <span class="nn">CacheManager</span>
 <span class="kn">from</span> <span class="nn">.</span> <span class="n">version</span> <span class="kn">import</span> <span class="nn">__version__</span>
 
+<span class="n">SCRIPT_DIR</span> <span class="o">=</span> <span class="n">Path</span><span class="p">(</span><span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">dirname</span><span class="p">(</span><span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">realpath</span><span class="p">(</span><span class="vm">__file__</span><span class="p">)))</span>
+
 <div class="viewcode-block" id="Geocoder"><a class="viewcode-back" href="../../modules.html#geocode.Geocoder">[docs]</a><span class="k">class</span> <span class="nc">Geocoder</span><span class="p">:</span>
 <span class="w">    </span><span class="sd">&quot;&quot;&quot;</span>
 <span class="sd">    Geocode addresses, postcodes, LLSOAs or Constituencies or reverse-geocode latitudes and</span>
 <span class="sd">    longitudes.</span>
 <span class="sd">    &quot;&quot;&quot;</span>
     <span class="k">def</span> <span class="fm">__init__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span>
-                 <span class="n">cache_dir</span><span class="p">:</span> <span class="n">Optional</span><span class="p">[</span><span class="nb">str</span><span class="p">]</span> <span class="o">=</span> <span class="kc">None</span><span class="p">,</span>
-                 <span class="n">gmaps_key_file</span><span class="p">:</span> <span class="n">Optional</span><span class="p">[</span><span class="nb">str</span><span class="p">]</span> <span class="o">=</span> <span class="kc">None</span><span class="p">)</span> <span class="o">-&gt;</span> <span class="kc">None</span><span class="p">:</span>
+                 <span class="n">cache_dir</span><span class="p">:</span> <span class="n">Optional</span><span class="p">[</span><span class="n">Path</span><span class="p">]</span> <span class="o">=</span> <span class="kc">None</span><span class="p">,</span>
+                 <span class="n">gmaps_key_file</span><span class="p">:</span> <span class="n">Optional</span><span class="p">[</span><span class="n">Path</span><span class="p">]</span> <span class="o">=</span> <span class="kc">None</span><span class="p">)</span> <span class="o">-&gt;</span> <span class="kc">None</span><span class="p">:</span>
 <span class="w">        </span><span class="sd">&quot;&quot;&quot;</span>
 <span class="sd">        Geocode addresses, postcodes, LLSOAs or Constituencies or reverse-geocode latitudes and</span>
 <span class="sd">        longitudes.</span>
@@ -132,6 +123,7 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
         <span class="bp">self</span><span class="o">.</span><span class="n">cpo</span> <span class="o">=</span> <span class="n">CodePointOpen</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">cache_manager</span><span class="p">)</span>
         <span class="bp">self</span><span class="o">.</span><span class="n">ngeso</span> <span class="o">=</span> <span class="n">NationalGrid</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">cache_manager</span><span class="p">)</span>
         <span class="bp">self</span><span class="o">.</span><span class="n">ons_nrs</span> <span class="o">=</span> <span class="n">ONS_NRS</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">cache_manager</span><span class="p">)</span>
+        <span class="bp">self</span><span class="o">.</span><span class="n">eurostat</span> <span class="o">=</span> <span class="n">Eurostat</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">cache_manager</span><span class="p">)</span>
         <span class="bp">self</span><span class="o">.</span><span class="n">gmaps</span> <span class="o">=</span> <span class="n">GMaps</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">cache_manager</span><span class="p">,</span> <span class="n">gmaps_key_file</span><span class="p">)</span>
         <span class="bp">self</span><span class="o">.</span><span class="n">status_codes</span> <span class="o">=</span> <span class="p">{</span>
             <span class="mi">0</span><span class="p">:</span> <span class="s2">&quot;Failed&quot;</span><span class="p">,</span>
@@ -147,16 +139,18 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
 
     <span class="k">def</span> <span class="fm">__exit__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="nb">type</span><span class="p">,</span> <span class="n">value</span><span class="p">,</span> <span class="n">traceback</span><span class="p">):</span>
 <span class="w">        </span><span class="sd">&quot;&quot;&quot;Context manager.&quot;&quot;&quot;</span>
-        <span class="k">pass</span>
+        <span class="bp">self</span><span class="o">.</span><span class="n">gmaps</span><span class="o">.</span><span class="n">flush_cache</span><span class="p">()</span>
 
-<div class="viewcode-block" id="Geocoder.force_setup"><a class="viewcode-back" href="../../modules.html#geocode.Geocoder.force_setup">[docs]</a>    <span class="k">def</span> <span class="nf">force_setup</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">ngeso_setup</span><span class="o">=</span><span class="kc">True</span><span class="p">,</span> <span class="n">cpo_setup</span><span class="o">=</span><span class="kc">True</span><span class="p">,</span> <span class="n">ons_setup</span><span class="o">=</span><span class="kc">True</span><span class="p">):</span>
+<div class="viewcode-block" id="Geocoder.force_setup"><a class="viewcode-back" href="../../modules.html#geocode.Geocoder.force_setup">[docs]</a>    <span class="k">def</span> <span class="nf">force_setup</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">ngeso_setup</span><span class="o">=</span><span class="kc">True</span><span class="p">,</span> <span class="n">cpo_setup</span><span class="o">=</span><span class="kc">True</span><span class="p">,</span> <span class="n">ons_setup</span><span class="o">=</span><span class="kc">True</span><span class="p">,</span> <span class="n">eurostat_setup</span><span class="o">=</span><span class="kc">True</span><span class="p">):</span>
 <span class="w">        </span><span class="sd">&quot;&quot;&quot;Download all data and setup caches.&quot;&quot;&quot;</span>
         <span class="k">if</span> <span class="n">ngeso_setup</span><span class="p">:</span>
             <span class="bp">self</span><span class="o">.</span><span class="n">ngeso</span><span class="o">.</span><span class="n">force_setup</span><span class="p">()</span>
         <span class="k">if</span> <span class="n">cpo_setup</span><span class="p">:</span>
             <span class="bp">self</span><span class="o">.</span><span class="n">cpo</span><span class="o">.</span><span class="n">force_setup</span><span class="p">()</span>
         <span class="k">if</span> <span class="n">ons_setup</span><span class="p">:</span>
-            <span class="bp">self</span><span class="o">.</span><span class="n">ons_nrs</span><span class="o">.</span><span class="n">force_setup</span><span class="p">()</span></div>
+            <span class="bp">self</span><span class="o">.</span><span class="n">ons_nrs</span><span class="o">.</span><span class="n">force_setup</span><span class="p">()</span>
+        <span class="k">if</span> <span class="n">eurostat_setup</span><span class="p">:</span>
+            <span class="bp">self</span><span class="o">.</span><span class="n">eurostat</span><span class="o">.</span><span class="n">force_setup</span><span class="p">()</span></div>
 
 <div class="viewcode-block" id="Geocoder.get_dno_regions"><a class="viewcode-back" href="../../modules.html#geocode.Geocoder.get_dno_regions">[docs]</a>    <span class="k">def</span> <span class="nf">get_dno_regions</span><span class="p">(</span><span class="bp">self</span><span class="p">):</span>
 <span class="w">        </span><span class="sd">&quot;&quot;&quot;</span>
@@ -206,12 +200,30 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
 <span class="sd">        ----------</span>
 <span class="sd">        `latlons` : iterable of strings</span>
 <span class="sd">            Specific latlons to geocode to llsoa boundaries.</span>
-<span class="sd">        </span>
 <span class="sd">        `dz` : Boolean</span>
 <span class="sd">            Indication whether to consider datazones</span>
 <span class="sd">        &quot;&quot;&quot;</span>
         <span class="k">return</span> <span class="bp">self</span><span class="o">.</span><span class="n">reverse_geocode</span><span class="p">(</span><span class="n">latlons</span><span class="p">,</span> <span class="s2">&quot;llsoa&quot;</span><span class="p">,</span> <span class="n">datazones</span><span class="o">=</span><span class="n">dz</span><span class="p">)</span></div>
 
+<div class="viewcode-block" id="Geocoder.reverse_geocode_nuts"><a class="viewcode-back" href="../../modules.html#geocode.Geocoder.reverse_geocode_nuts">[docs]</a>    <span class="k">def</span> <span class="nf">reverse_geocode_nuts</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">latlons</span><span class="p">:</span> <span class="n">List</span><span class="p">[</span><span class="n">Tuple</span><span class="p">[</span><span class="nb">float</span><span class="p">,</span> <span class="nb">float</span><span class="p">]],</span>
+                             <span class="n">level</span><span class="p">:</span> <span class="n">Literal</span><span class="p">[</span><span class="mi">0</span><span class="p">,</span> <span class="mi">1</span><span class="p">,</span> <span class="mi">2</span><span class="p">,</span> <span class="mi">3</span><span class="p">],</span>
+                             <span class="n">year</span><span class="p">:</span> <span class="n">Literal</span><span class="p">[</span><span class="mi">2003</span><span class="p">,</span> <span class="mi">2006</span><span class="p">,</span> <span class="mi">2010</span><span class="p">,</span> <span class="mi">2013</span><span class="p">,</span> <span class="mi">2016</span><span class="p">,</span> <span class="mi">2021</span><span class="p">]</span> <span class="o">=</span> <span class="mi">2021</span>
+                            <span class="p">)</span> <span class="o">-&gt;</span> <span class="n">List</span><span class="p">[</span><span class="nb">str</span><span class="p">]:</span>
+<span class="w">        </span><span class="sd">&quot;&quot;&quot;</span>
+<span class="sd">        Function to reverse geocode a collection of latlons into NUTS boundaries.</span>
+<span class="sd">        </span>
+<span class="sd">        Parameters</span>
+<span class="sd">        ----------</span>
+<span class="sd">        `latlons` : iterable of strings</span>
+<span class="sd">            Specific latlons to geocode to llsoa boundaries.</span>
+<span class="sd">        `level` : int</span>
+<span class="sd">            Specify the NUTS level, must be one of [0,1,2,3].</span>
+<span class="sd">        `year` : int</span>
+<span class="sd">            Specify the year of NUTS regulation, must be one of [2003,2006,2010,2013,2016,2021],</span>
+<span class="sd">            defaults to 2021.</span>
+<span class="sd">        &quot;&quot;&quot;</span>
+        <span class="k">return</span> <span class="bp">self</span><span class="o">.</span><span class="n">reverse_geocode</span><span class="p">(</span><span class="n">latlons</span><span class="p">,</span> <span class="s2">&quot;nuts&quot;</span><span class="p">,</span> <span class="n">level</span><span class="o">=</span><span class="n">level</span><span class="p">,</span> <span class="n">year</span><span class="o">=</span><span class="n">year</span><span class="p">)</span></div>
+
 <div class="viewcode-block" id="Geocoder.geocode_constituency"><a class="viewcode-back" href="../../modules.html#geocode.Geocoder.geocode_constituency">[docs]</a>    <span class="k">def</span> <span class="nf">geocode_constituency</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">constituencies</span><span class="p">):</span>
 <span class="w">        </span><span class="sd">&quot;&quot;&quot;</span>
 <span class="sd">        Function to geocode a collection of constituencies into latlons.</span>
@@ -311,6 +323,10 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
         <span class="k">elif</span> <span class="n">entity</span> <span class="o">==</span> <span class="s2">&quot;llsoa&quot;</span><span class="p">:</span>
             <span class="n">datazones</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s2">&quot;datazones&quot;</span><span class="p">,</span> <span class="kc">False</span><span class="p">)</span>
             <span class="k">return</span> <span class="bp">self</span><span class="o">.</span><span class="n">ons_nrs</span><span class="o">.</span><span class="n">reverse_geocode_llsoa</span><span class="p">(</span><span class="n">latlons</span><span class="o">=</span><span class="n">latlons</span><span class="p">,</span> <span class="n">datazones</span><span class="o">=</span><span class="n">datazones</span><span class="p">)</span>
+        <span class="k">elif</span> <span class="n">entity</span> <span class="o">==</span> <span class="s2">&quot;nuts&quot;</span><span class="p">:</span>
+            <span class="n">level</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s2">&quot;level&quot;</span><span class="p">)</span>
+            <span class="n">year</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s2">&quot;year&quot;</span><span class="p">,</span> <span class="mi">2021</span><span class="p">)</span>
+            <span class="k">return</span> <span class="bp">self</span><span class="o">.</span><span class="n">eurostat</span><span class="o">.</span><span class="n">reverse_geocode_nuts</span><span class="p">(</span><span class="n">latlons</span><span class="o">=</span><span class="n">latlons</span><span class="p">,</span> <span class="n">level</span><span class="o">=</span><span class="n">level</span><span class="p">,</span> <span class="n">year</span><span class="o">=</span><span class="n">year</span><span class="p">)</span>
         <span class="k">else</span><span class="p">:</span>
             <span class="k">raise</span> <span class="n">GenericException</span><span class="p">(</span><span class="sa">f</span><span class="s2">&quot;Entity &#39;</span><span class="si">{</span><span class="n">entity</span><span class="si">}</span><span class="s2">&#39; is not supported.&quot;</span><span class="p">)</span></div>
 
@@ -382,7 +398,7 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
                         <span class="n">required</span><span class="o">=</span><span class="kc">False</span><span class="p">,</span> <span class="n">help</span><span class="o">=</span><span class="s2">&quot;Force download all datasets to local cache (useful &quot;</span>
                                              <span class="s2">&quot;if running inside a Docker container i.e. run this &quot;</span>
                                              <span class="s2">&quot;as part of image build). Possible values are &quot;</span>
-                                             <span class="s2">&quot;&#39;ngeso&#39;, &#39;cpo&#39;, &#39;ons&#39; or &#39;all&#39;.&quot;</span><span class="p">)</span>
+                                             <span class="s2">&quot;&#39;ngeso&#39;, &#39;cpo&#39;, &#39;ons&#39;, &#39;eurostat&#39; or &#39;all&#39;.&quot;</span><span class="p">)</span>
     <span class="n">parser</span><span class="o">.</span><span class="n">add_argument</span><span class="p">(</span><span class="s2">&quot;--load-cpo-zip&quot;</span><span class="p">,</span> <span class="n">dest</span><span class="o">=</span><span class="s2">&quot;cpo_zip&quot;</span><span class="p">,</span> <span class="n">action</span><span class="o">=</span><span class="s2">&quot;store&quot;</span><span class="p">,</span> <span class="nb">type</span><span class="o">=</span><span class="nb">str</span><span class="p">,</span>
                         <span class="n">required</span><span class="o">=</span><span class="kc">False</span><span class="p">,</span> <span class="n">default</span><span class="o">=</span><span class="kc">None</span><span class="p">,</span> <span class="n">metavar</span><span class="o">=</span><span class="s2">&quot;&lt;/path/to/zip-file&gt;&quot;</span><span class="p">,</span>
                         <span class="n">help</span><span class="o">=</span><span class="s2">&quot;Load the Code Point Open data from a local zip file.&quot;</span><span class="p">)</span>
@@ -392,7 +408,7 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
     <span class="n">options</span> <span class="o">=</span> <span class="n">parser</span><span class="o">.</span><span class="n">parse_args</span><span class="p">()</span>
     <span class="k">def</span> <span class="nf">handle_options</span><span class="p">(</span><span class="n">options</span><span class="p">):</span>
         <span class="k">if</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span> <span class="ow">is</span> <span class="ow">not</span> <span class="kc">None</span><span class="p">:</span>
-            <span class="n">valid_options</span> <span class="o">=</span> <span class="p">[</span><span class="s2">&quot;ngeso&quot;</span><span class="p">,</span> <span class="s2">&quot;cpo&quot;</span><span class="p">,</span> <span class="s2">&quot;ons&quot;</span><span class="p">,</span> <span class="s2">&quot;all&quot;</span><span class="p">]</span>
+            <span class="n">valid_options</span> <span class="o">=</span> <span class="p">[</span><span class="s2">&quot;ngeso&quot;</span><span class="p">,</span> <span class="s2">&quot;cpo&quot;</span><span class="p">,</span> <span class="s2">&quot;ons&quot;</span><span class="p">,</span> <span class="s2">&quot;eurostat&quot;</span><span class="p">,</span> <span class="s2">&quot;all&quot;</span><span class="p">]</span>
             <span class="n">options</span><span class="o">.</span><span class="n">setup</span> <span class="o">=</span> <span class="nb">list</span><span class="p">(</span><span class="nb">map</span><span class="p">(</span><span class="nb">str</span><span class="o">.</span><span class="n">lower</span><span class="p">,</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span><span class="p">))</span>
             <span class="k">if</span> <span class="nb">any</span><span class="p">(</span><span class="n">s</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">valid_options</span> <span class="k">for</span> <span class="n">s</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span><span class="p">):</span>
                 <span class="k">raise</span> <span class="ne">ValueError</span><span class="p">(</span><span class="sa">f</span><span class="s2">&quot;Invalid value for `--setup` - valid values are </span><span class="si">{</span><span class="n">valid_options</span><span class="si">}</span><span class="s2">&quot;</span><span class="p">)</span>
@@ -422,7 +438,7 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
     <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="s2">&quot;Time taken: </span><span class="si">%s</span><span class="s2"> seconds&quot;</span><span class="p">,</span> <span class="nb">round</span><span class="p">(</span><span class="n">TIME</span><span class="o">.</span><span class="n">time</span><span class="p">()</span> <span class="o">-</span> <span class="n">timerstart</span><span class="p">,</span> <span class="mi">1</span><span class="p">))</span>
     <span class="k">for</span> <span class="p">(</span><span class="n">lat</span><span class="p">,</span> <span class="n">lon</span><span class="p">),</span> <span class="n">llsoa</span> <span class="ow">in</span> <span class="nb">zip</span><span class="p">(</span><span class="n">sample_latlons</span><span class="p">,</span> <span class="n">results</span><span class="p">):</span>
         <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="s2">&quot;</span><span class="si">%s</span><span class="s2">, </span><span class="si">%s</span><span class="s2"> :    </span><span class="si">%s</span><span class="s2">&quot;</span><span class="p">,</span> <span class="n">lat</span><span class="p">,</span> <span class="n">lon</span><span class="p">,</span> <span class="n">llsoa</span><span class="p">)</span>
-    <span class="n">sample_file</span> <span class="o">=</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">join</span><span class="p">(</span><span class="n">SCRIPT_DIR</span><span class="p">,</span> <span class="s2">&quot;sample_latlons.txt&quot;</span><span class="p">)</span>
+    <span class="n">sample_file</span> <span class="o">=</span> <span class="n">SCRIPT_DIR</span><span class="o">.</span><span class="n">joinpath</span><span class="p">(</span><span class="s2">&quot;sample_latlons.txt&quot;</span><span class="p">)</span>
     <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">sample_file</span><span class="p">)</span> <span class="k">as</span> <span class="n">fid</span><span class="p">:</span>
         <span class="n">sample_latlons</span> <span class="o">=</span> <span class="p">[</span><span class="nb">tuple</span><span class="p">(</span><span class="nb">map</span><span class="p">(</span><span class="nb">float</span><span class="p">,</span> <span class="n">line</span><span class="o">.</span><span class="n">strip</span><span class="p">()</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s2">&quot;,&quot;</span><span class="p">)))</span>
                           <span class="k">for</span> <span class="n">line</span> <span class="ow">in</span> <span class="n">fid</span> <span class="k">if</span> <span class="n">line</span><span class="o">.</span><span class="n">strip</span><span class="p">()][:</span><span class="mi">10</span><span class="p">]</span>
@@ -443,7 +459,7 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
     <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="s2">&quot;Time taken: </span><span class="si">%s</span><span class="s2"> seconds&quot;</span><span class="p">,</span> <span class="nb">round</span><span class="p">(</span><span class="n">TIME</span><span class="o">.</span><span class="n">time</span><span class="p">()</span> <span class="o">-</span> <span class="n">timerstart</span><span class="p">,</span> <span class="mi">1</span><span class="p">))</span>
     <span class="k">for</span> <span class="n">constituency</span><span class="p">,</span> <span class="p">(</span><span class="n">lat</span><span class="p">,</span> <span class="n">lon</span><span class="p">)</span> <span class="ow">in</span> <span class="nb">zip</span><span class="p">(</span><span class="n">sample_constituencies</span><span class="p">,</span> <span class="n">results</span><span class="p">):</span>
         <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="s2">&quot;</span><span class="si">%s</span><span class="s2"> :    </span><span class="si">%s</span><span class="s2">, </span><span class="si">%s</span><span class="s2">&quot;</span><span class="p">,</span> <span class="n">constituency</span><span class="p">,</span> <span class="n">lat</span><span class="p">,</span> <span class="n">lon</span><span class="p">)</span>
-    <span class="n">sample_file</span> <span class="o">=</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">join</span><span class="p">(</span><span class="n">SCRIPT_DIR</span><span class="p">,</span> <span class="s2">&quot;sample_postcodes.txt&quot;</span><span class="p">)</span>
+    <span class="n">sample_file</span> <span class="o">=</span> <span class="n">SCRIPT_DIR</span><span class="o">.</span><span class="n">joinpath</span><span class="p">(</span><span class="s2">&quot;sample_postcodes.txt&quot;</span><span class="p">)</span>
     <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">sample_file</span><span class="p">)</span> <span class="k">as</span> <span class="n">fid</span><span class="p">:</span>
         <span class="n">postcodes</span> <span class="o">=</span> <span class="p">[</span><span class="n">line</span><span class="o">.</span><span class="n">strip</span><span class="p">()</span> <span class="k">for</span> <span class="n">line</span> <span class="ow">in</span> <span class="n">fid</span> <span class="k">if</span> <span class="n">line</span><span class="o">.</span><span class="n">strip</span><span class="p">()][:</span><span class="mi">10</span><span class="p">]</span>
     <span class="n">timerstart</span> <span class="o">=</span> <span class="n">TIME</span><span class="o">.</span><span class="n">time</span><span class="p">()</span>
@@ -478,9 +494,11 @@ <h1>Source code for geocode.geocode</h1><div class="highlight"><pre>
         <span class="n">ngeso_setup</span> <span class="o">=</span>  <span class="s2">&quot;ngeso&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span> <span class="ow">or</span> <span class="s2">&quot;all&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span>
         <span class="n">cpo_setup</span> <span class="o">=</span> <span class="s2">&quot;cpo&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span> <span class="ow">or</span> <span class="s2">&quot;all&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span>
         <span class="n">ons_setup</span> <span class="o">=</span> <span class="s2">&quot;ons&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span> <span class="ow">or</span> <span class="s2">&quot;all&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span>
+        <span class="n">eurostat_setup</span> <span class="o">=</span> <span class="s2">&quot;eurostat&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span> <span class="ow">or</span> <span class="s2">&quot;all&quot;</span> <span class="ow">in</span> <span class="n">options</span><span class="o">.</span><span class="n">setup</span>
         <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="s2">&quot;Running forced setup&quot;</span><span class="p">)</span>
         <span class="k">with</span> <span class="n">Geocoder</span><span class="p">()</span> <span class="k">as</span> <span class="n">geocoder</span><span class="p">:</span>
-            <span class="n">geocoder</span><span class="o">.</span><span class="n">force_setup</span><span class="p">(</span><span class="n">ngeso_setup</span><span class="o">=</span><span class="n">ngeso_setup</span><span class="p">,</span> <span class="n">cpo_setup</span><span class="o">=</span><span class="n">cpo_setup</span><span class="p">,</span> <span class="n">ons_setup</span><span class="o">=</span><span class="n">ons_setup</span><span class="p">)</span>
+            <span class="n">geocoder</span><span class="o">.</span><span class="n">force_setup</span><span class="p">(</span><span class="n">ngeso_setup</span><span class="o">=</span><span class="n">ngeso_setup</span><span class="p">,</span> <span class="n">cpo_setup</span><span class="o">=</span><span class="n">cpo_setup</span><span class="p">,</span> <span class="n">ons_setup</span><span class="o">=</span><span class="n">ons_setup</span><span class="p">,</span>
+                                 <span class="n">eurostat_setup</span><span class="o">=</span><span class="n">eurostat_setup</span><span class="p">)</span>
     <span class="k">if</span> <span class="n">options</span><span class="o">.</span><span class="n">debug</span><span class="p">:</span>
         <span class="n">debug</span><span class="p">()</span>
 
diff --git a/README.md b/README.md
index efd655a..66e6741 100644
--- a/README.md
+++ b/README.md
@@ -2,7 +2,7 @@
 
 Geocode various geographical entities including postcodes and LLSOAs. Reverse-geocode to LLSOA or GSP/GNode.
 
-*Latest Version: 0.11.1*
+*Latest Version: 0.12.0*
 
 ## What is this repository for?
 
@@ -10,7 +10,8 @@ Geocode various geographical entities including postcodes and LLSOAs. Reverse-ge
 * Use ONS & NRS LLSOA Population Weighted Centroids to geocode Lower Layer Super Output Areas.
 * Use GIS data from data.gov.uk to geocode GB constituencies based on geospatial centroid.
 * Use GIS boundaries data from ONS and NRS to reverse-geocode lat/lon to LLSOA.
-* Use GIS data from National Grid ESO's data portal to reverse-geocode to GSP / GNode. 
+* Use GIS data from National Grid ESO's data portal to reverse-geocode to GSP / GNode.
+* Use GIS boudnaries from the Europa/Eurostats API to reverse-geocode to NUTS regions.
 
 ## Benefits
 * Prioritises Code Point Open for postcode lookup to save expensive GMaps API bills.
@@ -37,16 +38,16 @@ This will print the helper for the limited command line interface which provides
 usage: geocode.py [-h] [--clear-cache] [--debug] [--setup SETUP [SETUP ...]]
                   [--load-cpo-zip </path/to/zip-file>] [--load-gmaps-key <gmaps-api-key>]
 
-This is a command line interface (CLI) for the Geocode module version 0.11.1.
+This is a command line interface (CLI) for the Geocode module version 0.12.0.
 
 optional arguments:
   -h, --help            show this help message and exit
   --clear-cache         Specify to delete the cache files.
   --debug               Geocode some sample postcodes/addresses/LLSOAs.
   --setup SETUP [SETUP ...]
-                        Force download all datasets to local cache (useful if running inside a
-                        Docker container i.e. run this as part of image build). Possible values
-                        are 'ngeso', 'cpo', 'ons' or 'all'.
+                        Force download all datasets to local cache (useful if running
+                        inside a Docker container i.e. run this as part of image build).
+                        Possible values are 'ngeso', 'cpo', 'ons', 'eurostat' or 'all'.
   --load-cpo-zip </path/to/zip-file>
                         Load the Code Point Open data from a local zip file.
   --load-gmaps-key <gmaps-api-key>
@@ -84,82 +85,77 @@ from geocode import Geocoder
 def main():
     with Geocoder() as geocoder:
         # Geocode some postcodes / addresses...
+        print("GEOCODE POSTCODES / ADDRESSES:")
         postcodes = ["S3 7RH", "S3 7", "S3", None, None, "S3 7RH"]
         addresses = [None, None, None, "Hicks Building, Sheffield", "Hicks", "Hicks Building"]
         results = geocoder.geocode(postcodes, "postcode", address=addresses)
         for postcode, address, (lat, lon, status) in zip(postcodes, addresses, results):
-            print(f"Postcode: `{postcode}`    Address: `{address}`")
-            if status == 0:
-                print("    Failed to geocode!")
-            else:
-                print(f"    {lat:.3f}, {lon:.3f}    ->  {geocoder.status_codes[status]}")
+            print(f"    Postcode + Address: `{postcode}` + `{address}`  ->  {lat:.3f}, {lon:.3f} "
+                  f"({geocoder.status_codes[status]})")
         # Geocode some LLSOAs...
+        print("GEOCODE LLSOAs:")
         llsoas = ["E01033264", "E01033262"]
         results = geocoder.geocode_llsoa(llsoas)
         for llsoa, (lat, lon) in zip(llsoas, results):
-            print(f"LLSOA: `{llsoa}`")
-            print(f"    {lat:.3f}, {lon:.3f}")
+            print(f"    LLSOA: `{llsoa}`  ->  {lat:.3f}, {lon:.3f}")
+        # Geocode some Constituencies...
+        print("GEOCODE CONSTITUENCIES:")
+        constituencies = ["Sheffield Central", "Sheffield Hallam"]
+        results = geocoder.geocode_constituency(constituencies)
+        for constituency, (lat, lon) in zip(constituencies, results):
+            print(f"    Constituency: `{constituency}`  ->  {lat:.3f}, {lon:.3f}")
         # Reverse-geocode some lat/lons to LLSOAs...
+        print("REVERSE-GEOCODE TO LLSOA:")
         latlons = [(53.384, -1.467), (53.388, -1.470)]
         results = geocoder.reverse_geocode_llsoa(latlons)
         for llsoa, (lat, lon) in zip(results, latlons):
-            print(f"LATLON: {lat:.3f}, {lon:.3f}:")
-            print(f"    `{llsoa}`")
-        # Reverse-geocode some lat/lons to GSP/GNode...
+            print(f"    LATLON: {lat:.3f}, {lon:.3f}  ->  `{llsoa}`")
+        # Reverse-geocode some lat/lons to GSP...
+        print("REVERSE-GEOCODE TO GSP:")
         latlons = [(53.384, -1.467), (53.388, -1.470)]
         results = geocoder.reverse_geocode_gsp(latlons)
         for (lat, lon), region_id in zip(latlons, results):
-            print(f"LATLON: {lat:.3f}, {lon:.3f}:")
-            print(f"    {region_id}")
-        # Geocode some Constituencies...
-        constituencies = ["Sheffield Central", "Sheffield Hallam"]
-        results = geocoder.geocode_constituency(constituencies)
-        for constituency, (lat, lon) in zip(constituencies, results):
-            print(f"Constituency: `{constituency}`")
-            print(f"    {lat:.3f}, {lon:.3f}")
+            print(f"    LATLON: {lat:.3f}, {lon:.3f}  ->  {region_id}")
+        # Reverse-geocode some lat/lons to 2021 NUTS2...
+        print("REVERSE-GEOCODE TO NUTS2:")
+        latlons = [(51.3259, -1.9613), (47.9995, 0.2335), (50.8356, 8.7343)]
+        results = geocoder.reverse_geocode_nuts(latlons, year=2021, level=2)
+        for (lat, lon), nuts2 in zip(latlons, results):
+            print(f"    LATLON: {lat:.3f}, {lon:.3f}  ->  {nuts2}")
 
 if __name__ == "__main__":
     log_fmt = "%(asctime)s [%(levelname)s] [%(filename)s:%(funcName)s] - %(message)s"
     fmt = os.environ.get("GEOCODE_LOGGING_FMT", log_fmt)
     datefmt = os.environ.get("GEOCODE_LOGGING_DATEFMT", "%Y-%m-%dT%H:%M:%SZ")
-    logging.basicConfig(format=fmt, datefmt=datefmt, level=os.environ.get("LOGLEVEL", "DEBUG"))
+    logging.basicConfig(format=fmt, datefmt=datefmt, level=os.environ.get("LOGLEVEL", "WARNING"))
     main()
 ```
 
 ```
->> python geocode_example.py
-Postcode: `S3 7RH`    Address: `None`
-    53.381, -1.486    ->  Full match with Code Point Open
-Postcode: `S3 7`    Address: `None`
-    53.383, -1.481    ->  Partial match with Code Point Open
-Postcode: `S3`    Address: `None`
-    53.387, -1.474    ->  Partial match with Code Point Open
-Postcode: `None`    Address: `Hicks Building, Sheffield`
-    Failed to geocode!
-Postcode: `None`    Address: `Hicks`
-    Failed to geocode!
-Postcode: `S3 7RH`    Address: `Hicks Building`
-    Failed to geocode!
-LLSOA: `E01033264`
-    53.384, -1.467
-LLSOA: `E01033262`
-    53.388, -1.470
-LATLON: 53.384, -1.467:
-    `E01033264`
-LATLON: 53.388, -1.470:
-    `E01033262`
-LATLON: 53.384, -1.467:
-    179
-        {'ng_id': 310, 'ggd_id': 310, 'gnode_id': 310, 'gnode_name': 'PIT1', 'gnode_lat': 53.400459999999995, 'gnode_lon': -1.44215, 'gsp_id': 293, 'gsp_name': 'PITS_3', 'gsp_lat': 53.400459999999995, 'gsp_lon': -1.44215, 'dc_id': <NA>, 'dc_name': <NA>, 'dc_lat': <NA>, 'dc_lon': <NA>, 'region_id': 179, 'region_name': 'Pitsmoor', 'has_pv': 1, 'pes_id': 23, 'pes_name': '_M'}
-        {'ng_id': 311, 'ggd_id': 311, 'gnode_id': 311, 'gnode_name': 'PIT2', 'gnode_lat': 53.400459999999995, 'gnode_lon': -1.44215, 'gsp_id': 293, 'gsp_name': 'PITS_3', 'gsp_lat': 53.400459999999995, 'gsp_lon': -1.44215, 'dc_id': <NA>, 'dc_name': <NA>, 'dc_lat': <NA>, 'dc_lon': <NA>, 'region_id': 179, 'region_name': 'Pitsmoor', 'has_pv': 1, 'pes_id': 23, 'pes_name': '_M'}
-LATLON: 53.388, -1.470:
-    179
-        {'ng_id': 310, 'ggd_id': 310, 'gnode_id': 310, 'gnode_name': 'PIT1', 'gnode_lat': 53.400459999999995, 'gnode_lon': -1.44215, 'gsp_id': 293, 'gsp_name': 'PITS_3', 'gsp_lat': 53.400459999999995, 'gsp_lon': -1.44215, 'dc_id': <NA>, 'dc_name': <NA>, 'dc_lat': <NA>, 'dc_lon': <NA>, 'region_id': 179, 'region_name': 'Pitsmoor', 'has_pv': 1, 'pes_id': 23, 'pes_name': '_M'}
-        {'ng_id': 311, 'ggd_id': 311, 'gnode_id': 311, 'gnode_name': 'PIT2', 'gnode_lat': 53.400459999999995, 'gnode_lon': -1.44215, 'gsp_id': 293, 'gsp_name': 'PITS_3', 'gsp_lat': 53.400459999999995, 'gsp_lon': -1.44215, 'dc_id': <NA>, 'dc_name': <NA>, 'dc_lat': <NA>, 'dc_lon': <NA>, 'region_id': 179, 'region_name': 'Pitsmoor', 'has_pv': 1, 'pes_id': 23, 'pes_name': '_M'}
-Constituency: `Sheffield Central`
-    53.376, -1.474
-Constituency: `Sheffield Hallam`
-    53.382, -1.590
+>> python example.py
+GEOCODE POSTCODES / ADDRESSES:
+    Postcode + Address: `S3 7RH` + `None`  ->  53.381, -1.486 (Full match with GMaps)
+    Postcode + Address: `S3 7` + `None`  ->  nan, nan (Failed)
+    Postcode + Address: `S3` + `None`  ->  nan, nan (Failed)
+    Postcode + Address: `None` + `Hicks Building, Sheffield`  ->  53.381, -1.486 (Full match with GMaps)
+    Postcode + Address: `None` + `Hicks`  ->  nan, nan (Failed)
+    Postcode + Address: `S3 7RH` + `Hicks Building`  ->  53.381, -1.486 (Full match with GMaps)
+GEOCODE LLSOAs:
+    LLSOA: `E01033264`  ->  53.384, -1.467
+    LLSOA: `E01033262`  ->  53.388, -1.470
+GEOCODE CONSTITUENCIES:
+    Constituency: `Sheffield Central`  ->  53.376, -1.464
+    Constituency: `Sheffield Hallam`  ->  53.396, -1.604
+REVERSE-GEOCODE TO LLSOA:
+    LATLON: 53.384, -1.467  ->  E01033264
+    LATLON: 53.388, -1.470  ->  E01033262
+REVERSE-GEOCODE TO GSP:
+    LATLON: 53.384, -1.467  ->  ('PITS_3', '_M')
+    LATLON: 53.388, -1.470  ->  ('NEEP_3', '_M')
+REVERSE-GEOCODE TO NUTS2:
+    LATLON: 51.326, -1.961  ->  UKK1
+    LATLON: 47.999, 0.234  ->  FRG0
+    LATLON: 50.836, 8.734  ->  DE72
 ```
 
 In the above example, `postcodes` and `addresses` are lists of strings, but it should be fine to use any iterator such as Numpy arrays or Pandas DataFrame columns, although the `geocode()` method will still return a list of tuples.
@@ -360,3 +356,7 @@ To update your locally cached boundary definitions, clear your local cache:
 ```
 geocode --clear-cache
 ```
+
+### Eurostat
+
+The Geocode library makes use of the eurostat API to download boundaries for NUTS (Nomenclature of territorial units for statistics) regions. For more information, see [here](https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts). For license conditions, see the FAQ section of the Europa website [here](https://ec.europa.eu/eurostat/web/gisco/faq).
\ No newline at end of file
diff --git a/Tests/test_geocode.py b/Tests/test_geocode.py
index efd7a55..66d06c9 100644
--- a/Tests/test_geocode.py
+++ b/Tests/test_geocode.py
@@ -10,7 +10,10 @@
 """
 
 import sys
+import os
 import unittest
+from pathlib import Path
+
 from numpy.testing import assert_almost_equal, assert_equal
 
 # sys.path.append("../geocode/")
@@ -18,10 +21,21 @@
 
 class geocodeTestCase(unittest.TestCase):
     """Tests for `geocode.py`."""
-    #def setUp(self):
-    #    with Geocoder() as geo:
-    #        geo.cache_manager.clear(delete_gmaps_cache=True)
-    #        geo.force_setup()
+    def test_clear_cache(self):
+        """Test the `cache_manager.clear()` method."""
+        with Geocoder() as geo:
+            geo.cache_manager.clear(delete_gmaps_cache=False)
+        cache_dir = geo.cache_manager.cache_dir
+        assert cache_dir.is_dir()
+        assert len([c for c in cache_dir.glob("*.p") if "gmaps" not in c.name]) == 0
+
+    def test_force_setup(self):
+        """Test the `force_setup()` method."""
+        with Geocoder() as geo:
+            geo.force_setup()
+        cache_dir = geo.cache_manager.cache_dir
+        assert cache_dir.is_dir()
+        assert len([c for c in cache_dir.glob("*.p") if "gmaps" not in c.name]) == 13
 
     def test_geocode_llsoa(self):
         """
@@ -56,6 +70,25 @@ def test_reverse_geocode_llsoa(self):
             self.assertEqual(geo.reverse_geocode_llsoa(latlons), llsoas)
             self.assertEqual(geo.reverse_geocode_llsoa(datazone_latlons, dz=True), datazones)
 
+    def test_reverse_geocode_nuts(self):
+        """
+        Test the `reverse_geocode_nuts()` function with several test cases.
+        """
+        nuts0 = ["UK", "FR", "DE"]
+        nuts1 = ["UKK", "FRG", "DE7"]
+        nuts2 = ["UKK1", "FRG0", "DE72"]
+        nuts3 = ["UKK15", "FRG04", "DE724"]
+        latlons = [
+            (51.3259, -1.9613),
+            (47.9995, 0.2335),
+            (50.8356, 8.7343)
+        ]
+        with Geocoder() as geo:
+            self.assertEqual(geo.reverse_geocode_nuts(latlons, level=3), nuts3)
+            self.assertEqual(geo.reverse_geocode_nuts(latlons, level=2), nuts2)
+            self.assertEqual(geo.reverse_geocode_nuts(latlons, level=1), nuts1)
+            self.assertEqual(geo.reverse_geocode_nuts(latlons, level=0), nuts0)
+
     def test_reverse_geocode_gsp(self):
         """
         Test the `reverse_geocode_gsp` function with several test cases.
diff --git a/docs/build/doctrees/environment.pickle b/docs/build/doctrees/environment.pickle
index c293395430e5a0dba6e92ecc6f2a15232a6bb134..06d8e5635ffd01eeb88e7ec99a6ec90b7546e74b 100644
GIT binary patch
delta 9299
zcmd5?eQ;FO70+(I34tIC*<DC@gb+5cSvDUa1{Q-D5SN5NB1uS=&1UnmyU)#T-1i<B
ztpZx(WL%4K7HS8A!B(q=h`Nqr?bk4NM#u4^;`BpCTWkO5bSloYI9972d+vR2zu&Fc
zaT=M)z4x8_JLjJByZ79E?_N&-QGek>eg5J6<J+HqCV#m;x4~+Q)#_{Jkp}E3YodaQ
zGH)?$TsA(Mmxn!3Xj!%>_H^Y$$^6J%&b%tj2Y9#0ud0p4Hm-g&cd5{wo0GGB%R2>c
z+E%-Ygr3Ax$EQK~vhEta&;za3DxC;E>popKyku?BO+wz<O}a;5%i2cW6v%6PhlSIq
zCl69};bCmo9*N!&H4Bem(>E^s5G!X?cnqsZVOwZO4)9|1o@f%2lnQ4NB)o`)Jm?l?
zuqJlSVeO2<f~KWKvhWjvI1Qt=@|DDSVTu&QaYZeXCxtvB9iDbi!daWimgb3Bm_&G_
zIN%p0mtW#Xq7cpS7(TJBTc4(+>V)}WKFkL_e6|kSBvFF9zHuc{6{bl+1NsBDv;%?z
z8&=i&#?`RiO=!Hs)EH8+ze`}>rdZ*1tTMtt3A*%KtRjU-ixlRWerI!C5-AvUW2{nC
z=!Op(nstRx(Ri1sG!J7@X5pX2aDf<-!^lUZ>A~IA@Im8}B{bxZu}?eXRjT%Dh*0QJ
z0^l{~J09#?3G5rxu5bY>&0v4e>UfXU@fgfqXDCjrM|?6WA*>;hYA!>ZH8Xs7UHMA7
z7?ZPRW@^NPPR`oOuf%2stwqqG6f)x($Xp04(#*pDz#7ttdE2Cs&cb=@Q)ae@86;A;
zxS6DFWR5%N@rA-c_Xr==$b_4iCgQOtIo(6C>`eG(3D`^+@bV!~Xo%0)&*D;rg@eox
z2n)I%xPL=M=2DWFI@w|0@_!6_GPjz99l(7v><+MZ!S7wOxJunX>3%TDj3=HNClR6W
z^$-)Q^0leCCX=t@Y@HJ=#zQnwlMCe-MKc$Q`J3EdPT^RKHTSB=VE4xAG?rdwMv=|Z
zub5bsrP56enJk@W>mW-CI?2*6nCSQ#Om0y>rPyzpH*MN-bD-B=nMT~lOjoms`#Te>
z5;uRdHIulH*g6`-{f&v9ZQ^osBFms-e?^_SHf_to*EE`@%kx>kdD%QxGO;R;rk0JF
zJT77DXz*xYqGy{&2yIEv+;*m`*~B@RSe3Xn*EMAlcLQ5TgSb{EdbWwf1y}|1_i=D)
zYe`xL-_CS7o4{L{Se3wQuiu<Wpv>0MATZ2Ce+vj)0=H~4vNHH_rsvr_KFY+ZJnr3g
zbtaFGuyr(ee3*&;7Vv0-xAqxWJpP91c{Yz1m{^s^AGBA)#eHRQ&Eyxun`})DHeX|+
zzZGnraT(y(JJQoyUvT9=C|@$eOQG>emj;brur)Pk{EUg7Z5rVv$C|WVaqZWj&%(s2
z^qE>};F4=$_A0DpYiiI}!$i+EeR(+%({_0Jq`{*7cGJ)FJkCk-Bczv!O)SRsnS0}2
zwvGmmolNv>^N2svm9VRSDQj<>V0xa-<2Vzm@;KemoXO+eY#j|A4>QrT%_Drd+me>f
z&of=kCN9Rrs>HptXHzC|PqTG2h?`-eXPY?snX;?jSRa%nk@2bWXQubrT)xl5Cbq_E
zc(bF-LdJ`#XKKk`>G#;`8nj+yqT`Nbel}~J?8ULkV^ZbN;=f^&c8@Z@K(p;C$46IY
zH~2N|DyMN*W_K6UZhhJbryQA&2idNmJmSKW>98Mf637ifHy$K9@UntbWiWM<-9jXI
z_Faf}-D0F{*pHVEE-B;+^MkTKFoasP?<i4*;=P^aedH{%pmS_wFqCqvbr!_EH|PxT
zv%TANy^z~oRo+X+$dMyO9>3(4{GlLqpcgiE+p1d&@d#@_0)%bYO57@tFY9*UMijLR
z4tSuu$=ycV?X)zYi|xWLY)=RsSfbZ^XeAFu^2w<v2~mg-o9_<}@wX~tQcog<nM9z0
z9$49PgDwogp3Z_W8jGxe^L;Df!=B|RoY%`M!df`ezXpEPy9$M`_1=u-hJ8*fPweZ%
z@*lKp@8hxkTpy8?^?R_C`uU}1fm*eN`2&1d<X!4{C(er`!ks;p(As+y$)`;;pM2x(
zQk0kFh+a^ps|#JaBN{0jQBRLk;>#OxOFkN27nC@Sm1G_#=)_LQ2@>?9V=Cl2=hgx;
z>w*OLIV}rj(6Cv#4U(`aJO!^f&2v+v2jL56U8>+Vc-d70EiDU6g%sEwaQptoRM=s7
zalbuPa5v<-)~5=_p>%-0<fj>$pwcH@3sRvG9J4u9a1>GNQU%ALec+l@!Eu0rjj4i%
z&<5SIS}Db;&YvU`UNFK-zR*x~Lc)tSS(u`XCkV_9i;}{iDu+@xr|SZ%yFzyvcDn0y
zAH!|#tIIwmBz{JUf05#IIPb1G^a(c2LV8zZu5z)*?V!t(44Q!ua*J|u?v0K|2k(vI
zr4jEAdwfx$7=?KDq-nCGlY+Pox0jD5EzV8)Gp7i`gS(qA6W>VT8y}AfIy5Qhsa<>`
z1E{0mQyF<$p!VT<f`1Ph4=JlLlcK!fXq|_0M6-H3npaOF)6yBm?ZQr4USU!T>Ue5K
z%W&Ld*O}pM&q4jvjj^3WD|Ck8aA<@Z#Y71B2RZ-9XecZp3z83l{AzfPF9+b4z+d=k
zgN56v9K_F=qryhVM#CWi-t|6c;N@^g#62o*p;dhb-siolbL*ZR9bK;LckJQp)V3MG
zyGB<Fulioiqho#F-$Wa+H~cT?7bc9^xrCkiKvir!a6xAf`M|K%JK`P{t-}FX^tn84
zkB>LkLVs}Q919oMg~{M$FojBr5(0O~n0SaC34LAw_Z*mni(_SVWL2AqT?-d*OMy^O
zY{pncsjo{O4e))#fsk9WaI|tdVd}UES{}$Z1#mwK5Fovs=(Ab4dJETJ;TjW=4P&^d
zl86SItq~D6+j>&gqXL#xQcoof^-koc!PZa@XO0_jjJ$cZT-|2H6;TR1VS!WymABR9
zwcpg$3El&haJys%;{ng10qV>E*Ub5Y(i|@3U*{<A@ehx2l8@(z8b^^9ZbTL(Zjk3f
zLHrx$(EWN#gQd}FO=u4g@4w0<!EhnhuwpyRdu3!0f2gttpwZLdY3wv~EY@2{z)nkI
zSmOH-KkSwRlE_IRg^ww*3jTV)R6%B{3TKH>Z&xEBF4#!{e-anQRv(-{AI1*P>!r}<
zG_zCzDZcrjJM0y$US2X+QCU@O8jc2)4kIEQm`=`H_n2DatIDeNt7KMJkzK8s4nK9*
zK<xf;cWH@*TNSq$y97t|$%T6BiT;8`$po?rA%Bo$t2uV*)WtmQp7+mODlo<ExkO%)
zM=|QExCG`sz7Af0{H(f?i_7C33e+BdjEPUYpv$dsYys=Ep~#LqnPaP?722oEm51$)
z9WXL&QXU+RE;u?}IhVrrI)-rLHNY8{Nf8V?1Qb-jnVprT0)+$|x1eM(%=DJ464CJ;
zl#p*hBfQaHpi1v_Oi)dfUZwo{Q1$y9r>MAN|1xF#(~cj*k+U_*!xN5Y;oR9OW$fo1
zKS!t$&oNEK%GkejymD-2u>yV7@gY>sRC=k^kJ<`!h+VW@xo8brUCInL3<kYXs@qni
zYg?Xy!OjK+r(#y<+Ux{V6t@!<p3x(F%iHX5W@gb2GFIEQy0$(7${)hF#JF{*uI(5R
z;+u&-s9KBa+MY$3jPIVQ7Rj3kgr*Hfda+#=9|~8sa!HKuWow_RrY|fC&4u(`CBor}
z_;R&PNy(>jcyM~D0;bcZH0gwq|Chs^TZUK1E$3D!^28*mb}GoEEm5Y2xFaX;s->tX
H4_f~NK-#+#

delta 9837
zcmd5?Yj7LY6=wamV-wr4V+V+{wnHK(vFtnxA#QBQ%Tbhvonpt1tE}X;rHiZ;t5p-X
zq#_w`l2$<4qcUZHCS*t(I*?#1?vzp}eLxuqWJ2kK5DKNV({?D7(qU*Bn4WvD9xGWg
zUL`-`A9r`pJ@>o!obR4<*L(EJyIJ3RGb{J@-0#4t(cIds?8WBA;fAbvSJE5kC|$+`
zhe`(v+0BQuj_2fHPY7D)RD>U|IW#THITGHm;CS|&V6My4FYgJhZ~Mci?^Jzv<F83q
z_7zvKr}n<r=0l<R`Vv!A>7Eejf;Hwk!4HGx9$^H2VQv+UKz7qA;bCZPS}Gg^t!alX
zn5(%ow;b}@85ns85gYb}_J&O4VZ3#YkmGpi3y~9e3Fcdse!pAxh3*W+b>dQT5|zkP
z*!b-dIfYj^(6lN*o}}WVOzaO9xfRbKo7Bgq8KIrWXhW^xqk}n`;?Q*;!JgnWyWj1W
zHJjVz(?SqRwMr9VnuS!w=F;RLst0GA^9rKIz#GlWr$-+U6)iggwTr9E6J^v|nBLS}
zL&fAV+A!Pxig!`l&_{4?#Wh*v2{_eq^-5|-ut@SaoQmp^hUAcQ#Lp}b7OL`qtjZpT
z9O5W{nuro0PoS_b6d6+HP=>scL$i4$RT@+!uQM|AbvBgz7B8tBK>mmw!F=?drpj^}
zFEM5)!=r3U62=R>v7)p*vW8$D&MJGf5E+5+l9fUkytib3MO6-(&JBv=2bONEV2=3=
zF_Ik9oM=)4ecTV+7&*Eg>dhoS=$Pa~e#hn^=kb!#{<qR~{3Trn+t}spo_?1@_K9$M
z*$gu!Js1;7o_<l2>~!1nHC#tM5!v90n*hVCX}PI}?e~umgzGP#s6Ri!H{{cp?vl^d
zU5pxfUrgdYvVfuSz)HJi()rRvj61U<eVQP@iWw;<8cCkR9t_0i`Q-(EIm0h6)8GR0
z^0GiQp2$1c!#!X(0>tq64kmBLAZ*-4pYNSOQbx?gJ={Bkvf@zs<)kDUm2%7Om|)q}
zpYNk`=((;qz;<8UNBbrqO#+nKM_-?+kEn|Rx>F+)AF&Lc@5bz{O73(Sj5)GkJ$Iy<
z;7Qyk+=V~P)bJBHkV70d3=u4VAGK8Ael3T>m6fZKxHx764<6cB!^w>A)*}ocZTdG|
z%Eb5x&NVY-l05+>H_Q#tEE)Cnt40jdZ?d?r&l!Y(v6gD)=kqtzCj0qiBUq{Y{Gt&6
zy*;nt1K&H(Gl(e?q7NF9_Nbzaq$UJs4VFj4!5m2?riO&Th;k+5CLCWZ`|eDee)Y!1
z$@<KQe<B&Y94@rXPG;&8LqI7a;$tHKdSfzVYPBxCWx^Dct!_*R9dimg46?D%AVe&5
z9@#zv=Bz1z`*)RmF>*4-4-S(N)Ko!qtr37rK+T$qq=p7jHSFu12^Y3aOW>$CQ%m9z
z(`^tU#?i-Z^~pi=B_mj=nCUPAaG5Yu2h$}nkS3=c%+#PbL3d6Cq9Ez5v>b}pG$zxt
ze=0~ZpcFYRU<BY2&~wo_ZQI(0gm@Xx)POi0V+J8&@lwBTUNRx4j9{gTmq(2NTn>at
z>rDxSyq&25aYEiS2oWQsbz>Fm950@Dz+5nbn2MIyi~w8?v?%R~v`o*N_$p@b>@PM5
zp`+#6<Txobf|ZJpJR<;?10j1Yi7#WTGBqF`C)XQ<h!H|x#%}4WoERr7j3B0>WvLN>
zOF+w|coDk_yql}j9rvzG4T{q=FcpZ2$9?O@8E~erbRtQPso+g^=<hNDaG8))1Fhj}
z(vfs;rUu1Hf~i1EBx%pqhUC<Ccq&Mf(Q}s(fXjp)5e|>lO!#N!mzi1;=j9oL5V0J0
zxMN{5Gfx@8N|oJyVgz6^W@2B`C;uQF&||!c5Bz%%|1>v63_V2rL&O79`j)}xCok(U
z;>+OQYA<RTUz>axk@=&+l6o2Zv-(9XV`rLW)Ug>!QkfuO8DQzUri2#cELp7hVVMqE
zS4#nv;OkWxJkfP^VHv{TgGI0rc%!RR=z+!EmBL9_-@QSohr`|Lgl_nKcU@IC1-Zf>
zEO59qNpmS4R<rDe(jConugyp9`%r<*!@F2TlGejElLaW+Ocp}_j%89KyEn789;e+*
zn(+P*X~q`UvxHsz_Fyip8`82TiB21-_|jq>7X|aNFsHhh7OtHe)CK?Du|ep78+$sj
zy|>4b*MX41X86}w1$@%8356SayRbdldkeP3JNvLDJA1Kxj<t1tGTtBT>*WKRV023*
zeAd^_9~Z+1_DbmKt%2LW+~|`1{$Lj2#o1*z&E8Nje<)HG)~!5eNYb3lSrMfYr{uHA
zxB8hE<F(yGjG#mQn8{h*VLl9Qv)z=(YQo#$S=*xOZ{TDmUd!5!x8$4H9MnSlz|$j^
z5M6H<tk~69tP}KrV4nx?_2%J+t}FE&eK5Xjsa~)P{<Et^FK|GMv_dZ!fM+BV{y$Tq
z7rIcmTrY6rj4Sm51vTc<H;06ks&LSLy<W`^FWGNcty}sYG%}bsBx_EkKXeGK#R~CW
zl<7R3tztc#?hx*1(hYYxL}3}6bj%jcz!^uQ@Dlvb@rBZJ)M>BK<~Ou?725je1<vBF
zi6nLfi+Pz}#G)wv+<;D+q(Ls=umhoyko`a?nCFuv)qw@-w^4{6nNvbe>tsseR%uYy
zAr{Bi^XGi{cmKAPXXx5!;~W_Yk@NKF_2?6wATkT~{hdDjBhrWa2lmRd0udi1xhTiu
zZFDdnv6!8YgCgJ26M_}Ro5>>9PBCO-KIdt|50wKgf&iTZyR)8N7d|skFF=8_VWu|h
zb$JFw*N|6HHF1lVCKAaF3!SrJy|Xn?6zT8|dsT%L6=~}5$|4c}22^E8^f@K3>>j3H
z^>H7l!%-L}UnFi-Tpn$rw&L|g75yGb9k%&BI4RTX+|agl?WP^Jwwv49H?`fkDX>;-
z5x1*;xk%pQkiDAN&epkBRTY)5!>tSsq5<aJlIp>$scxO*a?4m2SE6%t@<obnmtDo#
z#71#g@r%1%Znx-BG|?`Lw1Q3tsPb22)h9PbiWIeC9hBc!8K~=&WvYrbc$Y`ioU+)l
zwr!1dEsj(C8p?c{<aSdpMyFCdVi*1(M%nEXha`v7<&k~nx`vb%QUljXZl7FK)Ok}!
z`_}EYHSHZJQ+#GGrYUp33tn{9))du?zTt>S@_J=YzbUTX07snPf~_tayzhEG+o!4U
zBw0qU;Sb4=v%+r;y&~jzC7%yoP&UCr@87~V-}-JIto~{t+_>k#aOG`P)8VlL+ru3}
zsU-Zf`|i#P?|7g)Z&q{#1X<Kgwnz*=d+dD9B;Da3{4)>UKly0*^;2U)c9W$Y4v$tB
z%q+m;%3`s?<D*slVT)xOoE^Q2KXh7p;DgbcVy3p!LNI!(LFk)RNSMTJ@j}a3J?~a6
zyRmyFeh;sTbPrqhLvU<v5!1cHa`-?9etx=^BY)d+KU^7}#~&WFd>>YY>$u_x%P1-?
z8lSy@jSE|Tim;{0wW<(b!PAzv*qXeXXAN-4+m??|QW4X)xg`yMLA1o~8H|ki*jgqe
zs$ky|uEtfYGX!fR)rnV($!HdOYJs&8non1>)0J791gnKA<@U>Ve+1lW6|DQI5RWsL
zCxUi|U_FX5KjvILB0Z|$+-UVIdXvh+E_nL%Ty~;yx8l1A2;=jGE@&7RvE4X6hxam0
pa@kP_tH)}w>(uxhY|oC@=dfAs9;>UQH~wILB#bbp`uukD{{iE|vk(9P

diff --git a/docs/build/doctrees/modules.doctree b/docs/build/doctrees/modules.doctree
index b64da015b8f3f052e3137e27d25b3f2aee4a66cf..fa971041cf99d761507ed2b9ed621c13e8c6be95 100644
GIT binary patch
literal 67043
zcmeHw3zQsJd8S9w%xFe$ziqu*F~%)<G}6dmyuvcrmS0E~#x}+viFQwS%~Y3qx|@DT
z8pi};Vn9%M59=mCvYSmHgh$A3_AJhBc$_869xz}?fK3Q&LV#p>E}SHLl7)pO`~CM-
z_uj6$RW&m$IEUyQ_0;1&{`=qmegC@m)uV6Oc)`d8^e^1hY5BFr+39k%+Gtk2PCr~*
zZB}{<UZdN8qQCo>`>*U5!-;nJOsCmvSG<0B0ZLT-TD{u#8vT>~a3dA(*6N*UQGdMS
zRl2oiLsj=j{V{+1WWVT7grnVBx9(}L!tuMC)m|MPQxTL{bH8_%x{h8v)@M+lnTx`)
zRvG<v`~H@2qND)zSBGN)zGLB<dadF0kM(BD;bghnZP#XdT{KVQSY4}j2<Guhz1->G
zA+ARo<pqAVw%Th&H?!r+skS%A%iy9`?Nj}kelHyFmglKp-{0hK_1F8`{MG*U@cLU1
zo_KM`Yj;ir<%ODe;;`2_)or#;+~GAV!0d^`%}VD4hw{V%AM<p(gQ4s=vZ^=BaMNtn
zDmrc^Z*4ZZ7;a~BlA)*QZ{{=NkA`dIJC3+ZwFeQxo&)%v&EY%d?;xD`JAnv)7f^i_
z{$Gv%*Wmx|ll`8*mknpe-|kNXrf~+-y4ha6?sZG`8X%SQASV<sPtc1}t?Mleyrt=$
zz;u=8+vS!&@S5HNtsbTmWHr$ArLpsm;oP4voT;)sqHt79+qghC8{@rb<ArI6xEqCD
z00bPC=;jpc-(}q&v+nP+?lFpaKK*?)DbXWz(vBYuQQCpNV??r5nl&0xGs}qTX@%}Y
zOSQB@w{qNgPD-MdZp+frjE;b)O~@EQ`Wy?!J*d{W*C=?Oq+ayT7!#AY&YTpDiTp(-
zbmC?SueSQXrvsPFI{vl(Ba(gI$Ez?FnndRpgrl`a)jK;#;Yo}bj)UXdH71jakZZNZ
zJS1kgRmhQcp~$3Hq1Gq}D%jl*H;c<@S*q`k9qabX1l9*J5&;$$0sv|aM7=I*q=dcz
z@*HnUof;%ux7XG=GwjI5i<NT4_e#}T`(WXhx{ED<IQ_yF3H@^Y_%F4}UB6zNoxTTu
zPZkapjx`(JRAGLh-0GB0d5fjFTHVW9Z?}+pP*e=Aa--4gmM!%X9yosERb8l!gN17t
zZmzj$s&JzmN2${)SJ0s<Qn==(n;`R;M!GZ{Wi=x91SYonf38nF(4~;V8`CB;94SZv
zj4Q?Dtd+0pt(fE!j+_BuTjKYt)N^k(;(`H*HgR^Levpj^9MEceP&cPNCf5E8bwAJh
z3ty6=i3F{jkTla;JP1|V0ns()3(Yx6R@3smoOV)GBaEik^T~kCKx!}&-r-tVM<OM9
zVB+S-$+#(o>p*^cF+A3KxPdxjv&8N_9c}*q@qa)x{(qz4`#)4*PB1VBV2s8?CI>X9
zdSJj9Y=n*E0&^{a$zfbeJtcw{tf!XL1Zfan$RH#wx!`rZwxwpn+IeF%#b=z<6-4^v
z7BpxSyX#!88evmLLZoCPe?qk(DMWi>qJPSX{&1u*z(K(thUeg!_`~S*VuoIVqt$Cv
zy1fC(F>+9Xfh)87Ln+<22D^W*(f!w?-f_r?;Ha|~|9f#lLX0oC7Ijr08_gL}Q?b6~
zHL76j(pgwT+H4Ao2GlPrTH}An4Vfmpm?2kfE|hDHez>Z&2m>2^K`%gUt!~Z+eO91b
zZqIw&ez?g}82xY^N_Sdbr8Woi_Z*DI#e>{p>mP%ZTRg}m^fB;x@N$mfm*G!uHn>qo
zz{qhmmB+22oOTW+xD&5!JO^SfMNcKICxy8+MEh1(zvApl7hO8mf(EK-SQARgz}tFV
zudW!<xFUoQ964!?<n7LpWI&p!$Aagey#Xwn44a8#e$NS)lVhOCw<w9wY<cZ+x7pTs
zA{;rSXy}lIj*mFI4M&=T;^0G6Wvi4%nFjZ?vlsW^Cc?2=gE_yNLS94GmmCzoF~!_p
zIC~`JJ6G=2yQS0Rde1Y*7+%!rwOY+~x71lEx4W<CHObLFM|}Eaeds13=sY0Qckm2e
zynsJo>#Olc$+fF(!%>Q><Aa*Vl89zaiG+>TV-p?!NRTjxKF(ow67k=2=uM%NtY6lw
za`xmBO1vjO6GdI<>?VPt#8A=|Ob{y=lT^wx+^{656iYY;E#TRZj4Pp03IHvKM!CT`
z5~pw!J;fp}PNUrJ>`Lb*heipWu#vWGbrOrQWOec~=P)v0%#<P-I58W4G@Nj``QyVt
zsuENV-SLRCkKidZGoTrT^2o#LZ*X?w9?JjFhmxu^-sbEn_@+IebWOp>Y0j}u+p(G=
zz4#75GU#{#93J&X%wYO-%{xN@T3Dv?vQ%zVe8VyjHeX`QhC-Rrz`eXA;*I-+Gc5li
z*o<NbBnO)TRyc+m_}#FUwfJQ40Ya^E6}b}9#`#GW6htT+PE9|YBo7o8HX`nQ!!HiU
zgasW%8<`03$HKL92<bAnjDs;*>NYDSIDFksKlluKXOCYCqn*8oN|K-bCy1{x{~92E
zhV-Ns4iFJL&7R^|EBY!-8eVZr8FBZ6g@e+pKll&_{~@NHi!emTMUlw<iXOq5`l3a8
zh6V~ROlTO<=J7hh)a`OzgD?0xHO7MgV*-t)Iq(guH#3S;4A(HPiaR}URrYdM#h3~*
zhyI3|*&2Kme_TOPMH0csP|^_Re@71@w@Nw)<T`n0TlHF{)`jat1<sV)4bzKav#NPg
zeBO1Em4jC<*7o5Bk;FW@VqHXScVk&)sKG42@!x4`I3Y)5S`iVU0+z)PWdq4h6<W<s
zm*fO4j~>16g*SH!&2|BcFdZyu^{|M6h%G&7d#AAm<n0x_U#KJI(CwjGq1>oW)|-ua
zemgDg(eu)mqUvDtYDS71=})4at4rSVF%hfGi2$=!nhVAnBv!n(?%GvowVM^1k7{9d
zvC!x(v=$*j3%hnnWID?zBtR~04R!(pE@qS(XFErjq!W;2Wk~VW>A5DBSW6wR+iS^K
z>xZO$k*0z4RRjMM;Wb)Y#odLXAV;-6P(ZxVZ&tB%u+4sDf}t-XC}#1ui<N|@P4J5T
z=j4O%xZt6--d?dW?{%7dh==ZP_pmTfX_@!nDC51?YlHL3a54!mtddw22sY?O%Q|hu
zk{pa@=Jem{WPoCD1=<d-#E&7&B~kEzZty(S$oV9SV5dIr+g1dnT>YA^_h!3gQOIx$
z*J*_m+Vy89+FzzHmfOF609E0lWi`5ecTX;(Tl8ew#6+!ffNR)|IkMC4{+3`9b1*dH
z&swsV#3-UC_j3Fg(UXlPcg0QMgeoFOr#&%;_%Gc<bVc~y$bri`dH=@UlWX!sPx&VA
zcisIhVe&W_%bPq?6vj@TDxxMYdSXuAe{&DfHF+P+H+i3T_vD&9(Nn(3`>MOYB}^U%
zV_B1TQH<V-o577{IXR8?V{<Nl=pM0aF8?9lT#m1C*@@23h3F~YTwdbtZwYhB!N@t6
z>BV_uL9+mY7ty4!74qOoI2F<mN;BMaKP)|T5MCT-R`pHKk=!^QN6_H2LkYn_z{{fv
z!6yK3@DzRwV<&hTPpkj}^BsfFQu!M(Rly(A)nQy|tLUL&cD=42e2&_=G!F3MZkD}(
zu!2G#O!FHBc|+TO0aZ;*>fwZ|4mZy<>dkVsP_EYt)pEB?URHrwihQ-Yqx`H)WfhwR
z$A}QO21lLr80ROQFYkbSrn>}mK-3Q4Od@KfU8UpbCWosep~9rTi~!tm<l+=hKlnVR
zS(|u`9PFA!9&8!o{)-*;606Wb8~ujk2<fm@=%7QlEJs7TA8#~yYSy`My||>zt3?D)
zHZ_p5$x|=IBl3bj%^YSCGLRn7QXRCbgxlqI&Mp_Fl;1Wj9+>%Exg4{Bwu%^It-Dx4
zcI;GvUr!bG>^X%<(tIb@*HT%H$@z2SV*e|e{lp+cMjz{O)ks(xmgB9E%Eh<gUu|<?
zM`(1z=AMiyu-E3Qr2Q5VPyLY+KrBY_UN!dI-qO2JH4D!Q=i#?lant`B;{%Bqg#HEt
zkhdoAf3$SqKeBw_*Gd8wQ)GA^tM=?+j{!Sfns9%~4Y$jD{Zo!b!$BV7d?fpmw`QN;
zboc2p=y;#$28%J#>BWV86LDYgn*`&~?VHVa@Pz5Ml5~l-O$6J3ZF3H=F2}anFe+`E
zOK_#zHcIqKyXHb0V9lnHD1$6~qnP1To90!by+RpwOr3vAL0Or1+;<Ox%kxOGg_N$>
zJI!)ow%Kb`%akjOFdCjmk!ZxIZ$6BtDrk++^rQ-A7rIR|E?tO;OaIj2P#hBM#2f&H
z+v!iDI*GHz&WFMvJM2mU8CJt4!pXUsSFa|1GBL_tYpYP)L7>kiL`^$qa`GOTA1c0*
zfy&SxDLgTiaYG&?1=n-TUPphnOAd$Lh=7@GRu>0=wLgp~CVv?b-WI9DLXj=O&DK34
zXr6G~*<IMXw{XAj6=Vt^H$iAQ1!O9ifm`uKO2rn-cqLTFVu}Wup+|*blB`-^k_%5S
z=6HG${s1}98ik^zHRzR2oW{DO@BuFfH%sWfg_d9L)H*u+YbAE1u81C(bKJ1{GFK%n
zMhNO0_8Nmdxv0L}$ZhPH#YJvD+sPzB3vHWeoUad7!66AqdIqcMik7^|FBKV_9FxHU
z-Yc(V4N-+0JHIi<{d$zvahyCD>Pd`S^fZ`nV|G3+$WivM?oANhz7!CC6_s4Jn2!43
zH7Lw*X|3xfp#IlUJa`>`vSA#wJF+{4bWvus@by<eVSUdK*d&f}yV4m~;V(;J$WptP
zWK)HFGJL5O97{njmBR-~!4X!@A*x=)MOC55X@N<m`&eS~VxbEOC0^4^XtJp=*KRHp
zYSZ2{^fE#bSOk`8SBS0P0rafMTw2RA@jzJF3QtUoe-1UWle!R3pWy?H6HoL?Z!bnX
znfx%uaUj`2-kMN;!)i4XL;f{`GP&q-xZzu9I}<l3R1I*0jj1GVh$(l|GKobYH(YK@
zpmrO+bUyl2%EV;mqeIg7BGu0Gb(u16;j3F>9O5NduK3<%kv=(^=N7+Ml8)sW)0@MT
zL=P0w7~kv|GwD2;%t%EDzPS>Wm|Bk&d-PUk9rQkve^3g{rBudr+EAoV=VfeEP6TLT
z(cK<C`V8(&Cy(HEltEFlC+)zUZI;aNqSF~#o#I8Ml`XRmh8pvHX$tn_;;TpxiC_z^
zC1`gXxiSp_BIYR_!^r-|H+JT16uD6l8}kg#z_DcSgR^wi0^Gr(x_TH_=Yn5kzmL!_
zG4YB9GWQGiEomO;qC2X2wWu->i8R_+Mu>B-a2r$*!DzUsK?*k`Ho|#v+^D27jDM*(
zoC%_Hmeer$Fk|fVG1W*|A9C!|z>vs3G80WUnrLIlFaw6MC5&KLTnWfi`>~12-&9!0
z9R>Z@sOqx3B^Qg}*xND|7s0V;dV;qSj6)BOZFkaj@+x2i=s=|G)c|}sk+P4YK3@g=
zEUxrOSv;4>u-Q)MP_!7D#4Oaxql9FZtC_Mc`25@Hb5DfI$l4H`jCfk5iOnqdb9kC2
zIpceJ9u-3<ELt>NenX0sgFh$ew+4SE1|hYWDHtp*$x>j?JL^+O>kARC_d=mb?}bM)
zf&&AOiotWVXGVZZ#$pm<R^QEo5*f4l7FF6V84p|@+@S%$;nXtEUu;cxVknB+@U@Cf
z9~vb^gu;Yi(e&?ZVJ3!eKURdn)WrqkoU}7hg|u6L3=|}Uqv(}xw5IxdQLIV88$B=~
zy(AZ;A_#LKD!Jsa?(^jhK3#^kb)S#qP;AVJqsw%2Q)5okOBt@{X@z1=&sz!zH|D^2
zum^>6#GGzM@f;Y3f*<g*ZTdu+Sm!II+}?Z~s;D>2SrJ0e?%g!Y2W9ov;D?MQoibfI
zNhH8*OW$fT5Khwn%+RWElBC?0f@G*p(jzHECe171Bt@_Ve}fUhUb-|50m4a=nm-cw
zF-{VrVuzDN@41uo1YI#FiCqofN!l&l3AS+O_KtRlvV-0ctGo=is5oGAad6mgDOyd8
z^@uhYLy)eJ^;9t%Q)=%Ansp^)vPDK8=sl{D&>OZ7BoG{v1+v-3QH=YOP{we8;>$q#
z4ZbjO`7wot+)mJcK~<OXm0T*Epbup%E}Wofeu57Yj6-*VqC`TIA`|ljA_xb(pLYV{
z<#<2eMt$!6{2i`z@5f>c>GoWnz=Ba8q(mar&$|epS^S+_viLiUg0-B1c#mf#no^V-
z{D?*XttbW}m69n4ETzYMp7Y1*<m~)ogy^9=I~&)!grbED>FlhdO4}uaft%gV&P#1+
zmAvw<jHv<g$n5Y8iZ0<o48Bm5!Susbxlkd3^*d21!Ret_y1rRCr{@N%WhVj4mFH(O
zNRx{$$LTqUwsjMn#80f#!_k$~>3PvohAVnnA*W|{DIk<`U_7`Rg>yJPRTR&G@i?c)
zrq3j&C#KxKynL!SEzDgN0?~#b-SOF$f_hqR6RSfc%+l<B;rOi0(5i5JqzsojW~h$O
zzfFNlnpDE^iC_!HQ8S0*BTaxv_{TUtjFKIW554D(&mYhgb9~s<@ExD)l#SM`l<TE(
z4__U(k!<dEfd)QTmdT){y4lP!oS|n9T4H9~MI>{?FAiOnFevF1CB9co$<*H;6DcI-
zGiG#;KCK!Fn?rYxOa#TQ60ISu#W+XWbHgEuo<kb)8xu?4Py@{E_<S8zb2vU<$yi)C
zJ`v`Fza$uk?)b=L$Y@p5=7vNUD1vps1^PH(Uycj3egYzeeFk<hu5=eDkqx9X^b4sh
z7`a6f#ZXR95w^3qMfYZLizrdLP&`VcUAS3_r@~?EYu>I^x=2rrT^Jfh(~9nc9l#_8
zx{dxM3JA+pU~bb-Lj;fmb%lTwj>*s+sF^%ikq*=}RoX5YHPwN-nWR_-r->0yW=APi
z*4>_KXu=`71(gyUB6_70(8@VP$E}u=y&N+9#46B(4ASJH%W;TKqHWzgCt(!p5OH+n
zbcni38LsGQg&d;CmIA_~sN~XsI>v(*3g>W$UWej2FdpX++4Py@5XF?cN;yOh9#j%F
z%v2nL%<d`O@p)bfiZ#cFgiuDu=h_Ud3dcvv-t*%4{MQt?EXOB;E!cyaIUFCU<RhUS
z<M=R&bvQoso;yBYqATY3u&d!aK2y?1kh$P)k7t)CKDg*3>o3FAIWS;n$wrIv0^f@9
zFRBk+gZ*HxiCqDbhfI`YCF$=BzEF(E9G$<XsZX?HGCDfnSB-@JupOOffH6x!m&vyl
z&KNd|W=7nED$89W!}YOVLilM<Obkv;x-jV640RM$T`E9wtXTaSVa3f+tk+3H4AU5i
zWI+-4mAB>?b}fAj+c{9iH@uGhR`<Dka(%%odNMXG8D%UXiCx^8V>{&TaS7X&a6ppP
zPwcd%7Z<*3#7w~j1p3f@*DDf5P;%B{1mu9z`fUh<<v6WX44ONwbzJFA>!2)?-l{r*
zM+KuKb8Fd(5}Pn3zAFS)>KDJy&ElrYzB=*#=(5^Sv>co!V%QpV=})4-j&nTEt4Jn~
z^A!Tb^Wkxx%YzQ-alV!+ZI|o?5$s;FI}%NBdTL@+HnZ=ksTLunwI_-In7aOrT*DP!
z=&zwtf)`4!bZS~TFZ4fJEhoE{niu-}4ASJHx=)lh7x787t@}hIU1GgZj;@?u=+jFX
zuIOonywJZ~3J70DC6@`JV?1~Yg>!hJe}&>XFdpZH+Vq*^g~pUCDlasRB^8T>+3P^e
z*?lH!cu+X|8uqOz$kzN(5>Xlb(br^XRrsS)GM^WJ^vZQ^(voJ9Scr>Y3m(OYa`>ZC
zn@1u&#vf(0E3IF1FE#d=7QN^G=nUZI{wTW|zCU^`4qqzajMtLbL&WWl%8w>=chso9
z3~%&G%@j1+gc(Jusc7!&yU7Av+TjDw6^Ah;_5vVSS30J=&**dBt{MsFhv#z|NQhlj
zTP0AFc;==gb79^`a0d_45Qc8v>uGJ~a1xiKA{Ym(`<nsta;*C%>T~P<46bzRet>AC
zaql8ZlB~2tp+qIb@~wo`EK5*_?(n+od*{Q)URXF<t9Xr$SGc+Dl?%71&3pI;0`yVz
zZTG#faF}*6ya%T$mFsHn!ezFiI3xIFf_rQ5DgnL}(I#`m+4tv@{K>q3Btq~~7L>;C
z#b)8{hT1Pj-78cuA<@r+2Wi;9nJOu2NHgpoq~nSn(gg2Ag<=&8i)?d=!?-A%@tdFt
zoQ#Si0@~PE0(8J-bm1<nOwxXqyaYRei<c=EuoXIX(ZM%OU7WgykLV9h7I`I^#=STt
zTHy~myQkSEJECQg>%Ca0p~1p^e!1n<7bov&))(iSjeU3Hz^|IPVE_<JtV&Y%h2q(T
zT4SoP$bQc*l+V(2`RwlL$w@YN+S8~rw4onE>)7%N`vK9_ZoxO!>CM_EMGuZ(QugdI
zre)8bLU9kx%$}c@so|PboqF-lNFjRVCJ6Swq8(4H4W_mXkc>1&HW;L2?XU9(xunbo
zsPxw0_wi@YpJp3RiMP6S19;cQZ7DUzDO!TLhUnXZ?$f!T`=kZkw|BFhf_@GJS%kLm
zTR%S%f)k=nYH1$m_TWuKBRh_)0iC{|QL)>+psktS%gRMtGlj1gtj6`mFlDo?J0{J4
zmXHh=e`<hn$rL;ZRLQTg!%Y@6lF`>`N;e{0z8FWdlQp&*2x;_vK0nsZQc@bxTfNEP
zQ#y<pWYFk(vjHc;1VRQWm?}e%K{i!HNDOCJp8E%S154OQ&CGY>sQeu^Ursc#y?;30
zB{v<8&Em8;I<|Ra58J8zR#?6of4@u2yeqmV15@uvhe>9LlZl-%F!gwbRz+Y+`V`XZ
z8ES_3VUV;7E7D^ZfvE^KiQX&H5FkE{m9A*yk_-l>9DAno{ezAJS{|4>L02qMf?b`b
zz|?%FwVc3|R(-_+Q(Bvjz?4=>%3c;@(4U(=XU{3?G6tqzp&A)(U`ii`h$8LO51+~#
zn2O*IUPD7T&w(j@IEjI&2*!cHR2?ucCouJEsLuma@4(eSU`husJ}{+~&=#E~Dk1nA
zgw-s8soP{=>cIX32lih#^DMT(^uEF!$8g4>{*@@(-7+RPRk-<{yHtQ`g<4ab61<<r
zu{C&~7>*R#ZjDCgMGiGkrJ<kVv=pd%cZBd^9>MX4=q86ckx)S?{AU6%`w<+HfKXGC
zm92<eM{sa;@+at%tCP>s6;mhdYWV5|+tGDP_UEte#g4e`=my(ZhGNn8x*BXaUb|rT
ztihsE#l*i<+{5hr@6fy_+W8sv%r{gc!`y%_ZlE!VcX^df1bpxV8pqI;jlMlw%$O22
zP6X(H#`zimUyjDvv_Wc|%W$P@oOl|NYA2Os<)VnhEJWd73t~x4czVN8I<kp$3Oh+x
z(>6GNriSNjz9hmvZ{p`T6>4*Z^67G|UdHb6C{;(ue9r4ud|t9$K12FwvWU|#>l8Al
z|1T-mvAZL?Y?Ujg%JVojmYoq*ysy&iQZX4RzBK$<0!i_4unThr#O|a&iL4yATr<Qi
z<6;O~sKRXmOw9i<RpB(c$)O4*R8VMp3BVmkF0mDwg+X~5m856XlvHllOjm`sp--*~
zZ>KA!3fa~0RpGVbEb5ZTEq80d!6kpI`Z5&Yb!y4qY7^hsc4@=KjOR{mvZ_mcIPhGt
z5L1OEn)gKh$*2mCt44-rV-5^2Mhjk%h%$x^D*=66lxFYDlc<bW$azxtCh|_JadYKl
zn&8gV^ERqx;Z%`&-eASe)N_5ad=q^A4WO5|Ch(6h9r!maANXjKa-P>><gYu=p!yv*
z+^!7vcW^A4`T1GSp!$HjPuCe#ywBt?q6~H!6P;dMWV1w^7c3BrL$~9iY!;g#6HT}X
z#sL%V0l>T*6Yk5X&w~bE$CYlv4TwEyy+zq9)&R7ifkY+jiw6m-S%L;9?)6T4ZJ1s{
zqab$xXKm84Ef6!cZm)_%d&u&Nj+3Gdc!f!)9biRgWZ3%+saV2rJ5ll>_+K<?gv!O>
ze+EJf)?~>L!#AVp**mbs?Fw(j`%ihWAv0orOqI4v=7h+POLb4q#(<7S+Ml*p?=hjY
z3dK4;!zc1$=`Q@JsGAGsLoEpg<pmkyBo?gQGm4Z;`sOCRO4*5Icy+JngINw&ZseSi
zi8$o8yqBWV;L<&7k6!7PNNOtl8XBxLbNG!3;XbQHQ$5?3t@sdoib0rMR8Ppxn*u+B
zw)KQ;5?Zk-@ElD!Q{eAh%5X(bX({kyoXRA&<*kYDmo5c>2Xf#$IE=zM65x-ccn*BW
zCBWNMnUnw@Q|`bL*ij)gTuN<-Y#WYbOG2(aUX(HcRF3ixHnVb39>Sn#TJ{l9h#@QW
zC$f!9(IhQdA(`Z@!Wt7KZw8>cAW7zB3{48AH}TO~tN5%vdZo9QDwyPfP!dXUT$Bmr
z?^vz6p%lXDoeaw4qPlR(%Pj9f+nJa}p=w~~Y8z8Y%o0=Xx%_L9tcQu2l%g;jE~^<3
zxeZx5`@A7#N;0$0ZTU`y5(d#;X>RiGkYl;V^tCW4(F27!M(JP4MU&#|wl6aTBuzEF
z)tLw7{zdu;QAvYS;4Y^TOtTGD`ZRwr*EGA^qZ6Yxq_-zCKFEZ$gKM_&BK8FQAuF0J
zA4uSdj1Or!OY?20amgz-Ws7!*jE@Mm;0I(s?>Le*<3qY@k;y%n@xcY7<5MKEGI(5a
zKj7wZNp>~-xa2;|&nZFBxnq;h5CUL2DUL`@EF(U7wX0njqm2DJ>7tmiNCPtVo#JjL
zjP3&Nb;)5;b;ii#ZK{#*)Udu%iXCQ5bW-_IWEYsBJoA+ipE6!Sda8R9TMwy$<_<-k
zK-C<f$OkhP7ZIvx27;Fmj6;u5ZF3Q9@{S4-oC8s+=L7ELM5$hj`aDYYIIi?4RV*9G
zNYxJ4K(uI;L@Tt)Erj7L(W=|@eH8=&YmqDSG$<KzV7vX1<t4t8RAH$VYE4me@HQI9
z*5FsgaHNhhMS!IP*(QJIo5?B7<%tOCLytlIZXS%t7}R^H(ss#y2{9=3txI@R&UgA|
z3$moS)xrRYSf@ZK5#|mseKgltg*EzNR7wai(JNh{r228uF;4=}=z$_trmFtbYS~Gr
za>u{VGf0z*>b_jwK+_k|w(iT3+=&e|adhPjG<|<5!xcTH1)9{^;*vJF^W2G&%o)c2
zX(=H57?oU#PRDrg9Td(HY8u&`3*%7a1KzVupGhuNOu1{-UJDsGQqgXhs~RMo-Lu=x
zq4UT8Fa_o0kch3%NlImObKaGqRpI7H;Vf0iP~DsxK?E)WlWTp#&52+Oej6i#X_WQM
zC#mNnsUG9zFba0KIrN^pIrkCqvh)jfHGDT`#<rQ+ikI8bdFDJn3p?4-VT~-q)!CnD
zaB=ViZ%q}d4<EX}-Sc9-o=akfCjFj-w~8g1aI0fFbVX(Ah>V_(uNny#&X?ySCLqT5
z*=*yEA!Nl9GG3zwo!bj~6{=?8U9nF8i>$ci#R<9)V>8^iRn><j4bY`I;lI9gu;0qz
zFw!R++yA}S-IFVMLiA*8|F?$NK#q&wm?Qg7?*5jr`5y-(S@TF4Ys^!6apA2-ToXJ@
zpby<!O?L7)GbLtKM34^ns}(@L9DnsIsL%b?Z{SMzR};jS^jM=^_+y2j=CexV!bI>0
z&slud7cALlWvJ9rX_e!;l1(bY5B`=$gs&9Ez@+J6it#w__52HJ@?XCjO~BCo*VS8G
zs?EZQ^j|NaO4}tvLP{?)hvG%jNM5Rg8&6M5^Ky6O!iey4uSBH;FPC2F#>2{axz|}O
zJLy?^x%(NU$whUaE3cRPY_zR=xg@P(y<Cp2oL=tHr3_c}v_f9)$)$ksvK$x>?nL1n
zUTz7+b6`Bq%eCn<$;*u?w?D<p9b`)-lf(RVAsKBLlFbNO<t2ajvJ~XY)!>U+IY%<Z
z4b_JW#ao)qx>s(nkLvGgV1p8*a5Wp)^T=Z#f*|TfC%4zC_$Y)v-gvm#@c6b)mMC%>
zpKk#uy1+>0RdUva0D3crKQ3K@Ug_<n3ZQW@Ole1|do_E^D|(;Rt{YB~F!62%XL3<p
z{N!bZ_oMAh%%G4p5c08cmBb7&<%*dJLZLCtUI(OX_nA)okEbBpBz|{$v_`%qy*(K`
zA`>BF@aX?ROSl|8qH+;Dl9rLQ;D#DJ`rj#|m5z-F9!0PPKPCfi$C1m^5FmC%kOp^T
zQpE(1m|}1QkLW!Q9!*||D;7LrSHlk;?YEtaMz7m($BqU+V~uNL8KI-66GKM~IAL%e
zEc{Y;v4t(^rjH_Z^|-vbpEj&JK38v+yHf@B=j2Ib=^~%JuGi>#?Q;G2{;9&uRN=r>
z;X0IA^vdnQA_w;GzYazA@4udYXHWyb?8yv!a$x3U;Sio3n1N=w$%F~vcH|omZZ&#w
zJi;hB`O+}{rII^Ld|w4J(ba+}lQRajcB)2(_swnmKx1NBZ~IxDRWW~8>$_3>f&!M2
zUL47ZYDU;9<5Jy0H7~j2Zr7tKJh~foU`Gy0TL4!k8Qtmb=}t6Lw%y3*;vv)t@5%p~
z+97JBpt>no-8|sz#{Zt(hyPX8f%z%c6=~=mbQQo*Yp9k5lJ5rTYH;E;_`myP-!N(#
z7<@S1_U0gYynaAiZnF2?avPt~_rrCXM8X{kc3NJgHdm|k&*9LO#e)mYYOju(pIJq*
zAUvql%biX?*o>A)SH7HMw;DH(#}$H)TcdfzIhw#VkD{ko#Emk-w)%R5voD=ma$#=0
zz#Q-L{GskQoWof<>X_!SQ1^akUs+N2giUm#yswqU*_0x;DQg@bb&ewg&a#P(&pNwF
z5F1W<j+>OfaQBi?5}dC|;pPN{+Q><3QvTLCjtrAxX^rHbezM%2OlQUUWiMAVW(A#5
zq#Z97C=LKL=s!bQUny+AqqPB>k5PEh)G?#1xTJk?XR6I@7rCr2=Xfqh)f@_XkGm(A
zVCOw$<Fp%`-6U|D7>Z%xT$GIDSPe5S%3HMj0yjJvd2&EIZ&ccuQNnT*^zCu;n5>{*
z?i`SlIeAb11>t4*=(SSn)aFZdPRXT4Z(-JJ_n+uTM%aEyxFh>5oHm>wgd+=m+YAfp
zRAp-rJJw|F@JeS#S;m^lsC=Wdn*^E_Ls>xG2xl`2XEPSg-tO!pSj6iULfX6CJ-HZN
z5J)DZJ?ZQw57G`Oq#dx3_AzH4D~YtvxqDhEr2Un%n><LnP9g0&3u)hT_K_cH$=drz
z?w-<-W-8+h+Rm(yQ+waQqcT#jWRJ!$R9m`z4K*g`fZv#!{fdj1tl2L?)f}4r8h1~w
zY2ZC&)9lZ7c9TH6Vkm}Yzc$AV2-N4%_;<MB{29>rx5mwHvc^B|9FS8Bik=K5wT1{^
zQ7*eNDYxeC&!wtjN0zML&$xSXQ7-Q(8|6Oc>?VP7lTdfPj=JkD+53dEC#USqiMro%
z_q0-|`}fXn@}X`<S(G3_+?Z7N5ob>;iMmg_ds->feaYEPPSoiPJ7a77Z#sKgN!0y=
zyQh^x-RetRx-%#0bcUU=8Fq`ar<FwAPIpf$g}S}YZt|h-fO2QRu(&Z9_PNfU@}n-<
z&b-UrQ#$HQv(ag1=CO!QI^jydwXle?CI+8!_7OaV2=0L2_#qzo&oh)oXE*Mln3>0n
zQp86>ab)v4XHUTrGLQK<%{fkW4bqE?%;RWjGnk=hC0lSr3}rowp|r)BZCc=Di(Ect
zD@@wDZLAR+x1#hOCyMl;<yse6sTHN)26zMRNwWF+VdrR^=Bwx_7ICYASxX{)%Gpy^
zybZ+PSx?JY5{boEvb6h(a}*iyl}_64pf0dSuZq5Fq#Zno^+XxbLWTfk66}0szv%D6
z#Z~xkJ5X}L^70a2k?iaf9|Xz2{csekm$hm?j5eE&UCd9+sgD31$f>^)V^~g3{ZD`j
zo>M<|DXpW*e35=a@`OyOk2af*8;X`uFNFs1y-na7%=3i-Bp0>K<nCH>M*XoD-+y0R
zLj4N1sH}isGcbl>ZKOYmTh;3HLWyQ{;K9(gdXAe5xh3H`0Vc3N^mO{G@}NYf(-)}H
zcFBq)V(%E(I?>OuCE?FhI{i(#Fe0*9528{+Iz7G8i7qvpHTp(e7$?yKMbRuF>OQMw
zGa<?%XPHia4}&zhsII8<rqjO|ZR`20B&cH3={dS`rqkD#GF;Kqpkg*V&-l^!E0+So
z!>Hs+7u7Ky_$bUYxOcWyuT^SY<j_R9H@3zPUWMX0FivaeyDeM9rcaa^WyemYm&>0L
zuLyW)>GU?N$p!=$;9G;|q#&Ql(cXm5V&xp^^v`EY<b@)(Fk#vB`fhVfm$CO#g@qnA
zRh{(;$fT#eRgrBuGj(9<x|6Z{s%2Pre-eP{;v<=JDV0iyqyNB16PH|1uk`j(#nEbv
zSVj`rwVFdFgg#}p=Y~%t)qjG)mt1r?Eb%nj&cqT5Ndqil<0y$GV#>|rrwfJ8a0!(`
zx;7N)O!4nh(4Uzpy0o?w&lBto7=ZBxw)^QJeH5Rb3^E8?-NMiHQ!LY%B^u+0Zjoj@
zHCgrM$~ZE+Lz`XGqgX~tUF2GInMwA^w5PGrPZ*vs`>P5?OB&Lvm1LvKa=|Cq=wehF
zmS(vEZD(R5g`@#CvT>BeMlt2)bT%3yDMh?@zIp`1+T%@U*dOI$Sa*B$UCc~+d-A&&
zvh~wG3lV!qZqLxF_%24eIOpZNn1@o>MCQ|p?_wg@f^8TP;%Hmb5Foybk#2M38^wGV
z!z?_<cQN#ye;0FxSde`e!>%%a7vr3NasH0~87{s2AM!Y%Sx20_Gs+XSJ>8YyXVCo$
zfp&HcejZamvMU#m{`=%~YrQ?nVj|Q-I4;65NXim<=%k1y{jL6b0qit@oirw6R;w<e
zI)l|$Yt>gd1?E_%*=tv@A;mOxJrZuL^xAFgN?j7qkA>^ylX|V;^^b*XaYjQIALMlU
z;g)K%((BghooO8F*ywr<1d_%cDA#-Vrc#AnHysNn=Dh~K#OXHMs6T=4F1=2-ywJkU
zI^GSE-5RHGJ1K5^?K<u@beomZEY8}Zx9DyZaJ3iv;e~U}X1CGodfLk)BjJXQSLq^K
zu5@;;*<R@5$wt}<t5vEr7Zy;pgC@7~CtlrSPfBNMUBAR89U!AeQGwF+*@N|Muhneh
zkShAJsDETtIO&(`>UK2T;+?GlSKV^A*O7%r!_lfY+nYyw8_$&64d4nT+Fp49PbRwU
za>b)T^<N)utTkG_Zpmv@npKp-R|viC+}<1ZBK@x)ZZlt&ymq_UF8MeRqmFdIknNaV
z!Pf|zn!T>s$|!3EF+J;z)kN??P;$7Dw;RDAYi~1OSe3)gG=>fk<pV`^jHw@Pz#(Y!
zJsgY%j90z>TmA?!68aae@*4dR8#h*(?P{s&wLmE7uf{MuhVOCbFkYahRxPzC9C;RB
zm4zG2y>7D<JqlM*m~un4x6mpn3IvAP16b@{yHT#w9ByJyfaXSr=pcI3>Uh0svvhWW
zXlpav&1T0d=`|U-RJx_pwN8!G>ZVe))}b@vFd_lvxa16fI$V#zbifPD&_kmeM6uFb
z4QIvm!Qq$E6*}p7Z88Qi)zn%QbH5IU%&A+@>K0-zkWgbD6j`q=z!AQC^EKhdmPcn9
z&v{_V3MWZ=GFt}h#3kY7ez)7|9Nf1Lk2}+?#jf9MOgG!}`}*M}RFVOJlBKS<fO+a-
zp4Q^|-mCXN3;5jFo~u~R)03$AXUgX3YiXk0Q(k+z=5^;-AK?Fees`gc&L+C$S+G;J
z>-UeW4Yz~w7D|dw)c{HzoSo;v#oydnY;?<KOMY$Muj4-mc?})1(GS;m7RunbSM<Q(
z9U$r=cF%81><mAS+U*s;40tf{obiBVhC<Xn)oY=Pjeriv$N>LbAe;mBG%K~n9B4pF
zhzm((bb7N^x!M9zF}^P3Y>AUHAf7-o$i%IJmlrX+zK2hN&V`p$Q0i34Tfkw0rE;~}
z_L$^&eYk~MrYY}u^*Juz&V}oF3#|DLYOnShr~Aixv*B7UXTaDX%^JR22eWibJbupz
zX$yb^sdXW4(I#AkPXCpVgvQfGc>$=wW3)D>ZpMhEs_gkVJ>=KLpbzf%XZ?5k?+o5V
zmyh<s33^W-JyB_1;+X&LW9<F0p1+IM8fKxJ-s#Wyk3zY9(*JyLa5Z$|^XTUsX~)Ou
z=hb8Qd5nJEIgX$A(9g%$;^))!^X&=ze4l=9#%b%p?es%OAqR9+azIB^26TjFKqnmq
zbP`fPyU7N$+iXDl^aQkzPe7{+e3@j*zu$kvZ}-9tI=!(0Z$LbWF(=VDe?nqKmgAV!
zpYU>ET+rvjY_qONm=(R)dZ<cI7%fi7URk9JjrO4`+FjC3m)LEIHy4DH5PFqnk3<@h
z*Q0%9p#jHl70;i*1O-oF+K#O1Jl4M#jNB^YB&8Avpaat79V;1VJYJ_BGiqB`YgFnu
XY6;XnPZA66C&aXnXk<q8W~cu@Qgq`G

delta 7280
zcmb7JX>?RY7N#nlUef*2oq*C=NeCoJ(gAEeBqXvgvIWALut@?9$qT^_A)=zP$f6Py
zxCszJNI(QsBC!Pxg21pIaTIh$o#QB*9055qFpdiY>fBqeUN4Y8{OPK?Rd=iU?ydW(
zFRm~wZU{3j46mJcx>b|$K_hz%vW*?XD(B9DGGjWmzO6=67#j<3S>l7qV2j9bwwqQo
zt;F-#^wK#!l16yvmOehaXx2pH$&XmV6eWao(duCkeYG`AQYOGNSCf2al}|bA0F90q
zq2`pMj9L6F@WJI|h_pF<PemtdY>lnQ6C0dWoW6l=k|NkZSe4KQ<|lUZ9Zpmj>k6N>
z>*{Nh)Gr+F>lq7<k!I*U*#yJVZdpx+R;|>}IR4mSRHrX2F*q)XQ`;Kg@YHnJJSrTf
zJTQy5F~EqG9-fNS8kf?B*2I^Dl_aPg4RC#>3tsI!IO!3jR+}C#R_8wD_0*_6X=tAA
zshlRqzOeKej7^55E(!2TNg6AJcQS7AxhS$Zn|lqAJTetdC!6>z1H8IEpI2)WFN92t
zuN%f*M+rT6y#dbEb?5JDt8f(6GZ7uL6T8+GJ}QL)TgwdKm^B?<>zd9BQR<2qKE?nQ
z8+zKl#x>d=oHIa|ZaGRBF(nRVxDn$}L1UV4Q@7U4a+Sz&>SZL!y}9~e3BP23G42<6
zSU6HM>?~ZY`Gb}wsAUlQ(2u;mn6e~kQ}N+&rYMGEbI=(FQ;!+p$%hB?Y_Z1NU0d@I
zoXKhjhpy+w<shNjI>$4k((9R4I-_)2%?m0ILwnKG^N+SQ$K!lev6{7Ju#1n@^hg<~
zVD$|r%YtVIyZQ7G@kK)X^<xUp3x|dVSM&^&qK?9!5aDEF(Re3|hp4PH{LRkt5Gi?u
zb4{+3#h{Qhc0csVO@f&0LY|DG*0+Vc%ox5g#NqXt!`NsJo6x419^H6Th)s=m+T?55
zCKN+9trs?}he0`;Ig+GJwSJpCWY3hk7`XUG9KRwoS>*v-Lc!oY63Z`zL$^nkP$Std
zNOnihG}apX?IYAET#RS_2carCLcOg);3h&~M|;b!+6WPF##6>S7)5kkBBB-@zZjsX
zmrEHJpdbcsoOBxU$b|gbw&{IwN~Mw7k9s*q{ORxS4Zz-u>Ok>32pluV`to|WWo!&I
z^%XJ1M;f8wR53*M-NMR{*n`hB`r~LCbzRCx6e-OqHi`h)*Ds40DFB8VDFCdiNv%i%
zbx_cMoKj3Y8vci)Q;lhGp?|G;D^cUr2BCe`s5Y>BK$Uxw6ptAE4efkjQO&b8ol3$>
zj3uVYHR|tokvgdy3I=Agy+Vq2s3So}jE@8*s8TgyIvTfSO?C|YJa94{rnhfzc1q{c
zgXx4(1QRv^Y_WGAX6O71Wi)(i1XeL2>H=yB26`Qwd(>$@fm9hp?B<=EDBB-LECQMa
zyYwhE!Hp5^fWqi!JckBXl+DzFbdCZLvD8!?Ph}@y*N|svpka<5#SKO%Nu2qNAvoKu
zO-h_?LT5c?Y!IqVg9STHioW@{Z|@jkL7__-E{{X(q=4GbNC`lVMjF(nkgWo0qIVb6
zpmO*xY#zF+W$Q}cJ9$RNUKWIN9TV6~XhkMlCTzjWM5#o15@$<I0%r$Dc4L)rV`MC!
zZ4xkhk5men?T2ysk19)vMT42wBw+S=ex3O^QRCF3s3`+4pi1{)DfwO@_76HM9bR9O
zkOWn#42IW*(M1HWiV^+UzlcMJ)>V^$7B(T!;$db=^cIP$R%XIgbhtV6<oBCIcD{;Q
zC9rsW2p@{|ql&viraGHpQIB-)G>g>SrA5gc`3Y3mAgHL>I({L$in|>{x%mfy%=!X@
zj;k+Cge!}lzrOI#OK2TNXEBATXb4{L`i|p<eQpxC^2#$1TPSelM7sjGqLGHHMbs0x
z61}@{<(oAjp7B@BklZC+)B$VZ(~0MJ6KcjfU>DM59k7$SE@dsslepSp7P$Igat>dM
zlbF9&o5cnG73mbfI%|drQzqJW5R-;1fWuSHTdoixPW^#4S}SY5dF)B|_fqyfB8%S^
z=cPmICP_(9rOsezr4$Wj$xvS8Vo`8ts>>`dAO%hr#d5_G3M`gNYjj+tSOl&*SweAz
z{?V<kAX!gXy@J_Ls2^5rWbn*zi$GP5MWAX_i*lI)71LwLGRjm6_@<ZAKxV243uP*^
z%v3jFQ9+7D-iN>7h6L`z@67^CD{*~|%cgHKZqd*1Tw$gS<gU}6LDP&FLk=1ffD(-~
zD7{WS0VUB3fs)TLvxxBnXq*Nkb2Ba`oIw6;FJzW3<$pu<h^W)JL?Y@F3$#-U_%&_f
zvIYL^=;Hq)5!oC^iDeHngq1e|rr!S$5%iK0`GQMIT@%`W%7Tbo&zotFV)%QwZ+M$_
ziO|v59U{x~lv8Ax2FAUp8gJ88<txm6a1|o57b?tv5wo9g=Mn)OxJNKzHIVG}r@4F*
zLUsw2!C0()d^DxGsim>7t~8#N5|^IpJQ0-Ym<KAnDOjA*1^OJ#>lfYpH5C)17m`~+
zHDdP&0YwjzbttYALS<$|p#doF2%gy%A)p8m0*Y_8D3>EU)p;I>s!YS*?8-_aNhlUa
z1Yer@5h0gmo;E2nTQ2&Bne7r1F|+BLjC(XZLS(kR5;xjuQgH@oPymiJ(%^WTdIFB3
z7YfJvssW7M2Nesn>oN+`=e@5)S+yfY*(p|$yRN+y2`itl!mKA_d6ZS4#%2|$>49st
z9M%`EJvmHCAzlqNkye44Sx>!U9!AtSbqFfTfG$s0xd%xV_kx<Qo_6E6Y{ddUY9^5o
zvQelEpvLJdpC4vmi)gK7fu?%C6|vIBQq7Y()u1Kw;pl=PY%Q&bs?k>vB1-PBz5-{K
zjP}(pe3^0N!lQ#jOPxxArp18F^Q^G?r&KoASF^Y?V+-Nn^S#)CU}eohS*)(CnOH?>
zomY~Whc4H`v?X2GX<{);?agp<$#)zX@6cyLeS3sHKG!PZIepoEwjg7?Fl<>mFUP$O
zlw*@<WuWXY)Fvge2BC3?EW1slIsG7T9}A#<d8QKSH<k|n#>uwuQRlb}RaE?`HXT%H
z)Dxf*y%11Aeboej7k=cuZF-vHxi-j{Hv`gN8mZ(WtGpTWY!Y2=mJ3NMVtIF)0F{TN
zi;J-wmaZ74<PnPos%$j7u@e}pH(BNpHBPNU3SN#~o-Ri5_Hd4*Igqz<24A9iI47X*
zdHDCrOeH%&Km$TO+<<dh*leVp@KE$ZJ*01%Z(p`Bz9-;flTCcoyt1lLX$ojKqPe)Q
zpA9^#Ggy;w^KEjomz6XKFE7~Ox7A~nBgCe;*+_1t`4+8tfLX2)IZnN*`AKgZy9Z2!
zo;9&pFOi=&g`dg~0j0m$;1*_$#sCG)2zAt<=#EBFPdF-ip^ny9jq_cpH8b8>@t<X?
zf@fI_nay_Cxo#=%qlgRZ5nP~MSS%ar){o+a8hvaqeO%peUO|$wSvj1o%i;^PMdgaT
zs_Z;l?`jVrsc4-<OO_(h5&`KO9^hyO-d<R!cyTMYkrQGscfr05!;~yy)3!1NRp8~-
zW8+^f?-Lc{2`K_GQ^<*&8eB-Ic~`2smz;#(g#2R?mk}#J2v3ld1XUU!k`nAV_y{y?
z8b=q|`E^d%{5F2X3nu+0JFkeM9h*SWUcVV{d2N-vntoM8adjIt1m9@86mg@SN3y)p
zc<V?qNiWJ!tahc@8oKcgAydhbbf+;p)Hh^Cr$)9Um?JyY?ATDvj!n*%Wz@(^?K_nI
zQ!4*Pflq&F2f15?oBET=B%w{OAPb;qb1W~24Dym*R4PJDnjLIXWy>a^2r>y9s`v_E
z=C(TgWp8xjV<UsDng7Oc+I0OvWv5W2Kd5X&+4P{ow!x_FwfIZ?Tm1Y{JyY4&aA-%0
z{xNcbkQw+GIV$F4QhgT>5R>Znr~m^`M-t-Thyo4ET}q0?mOzI15V4bX|62-;w7heR
zdLkS}FEkt}?<{=N%2*p{*sWdJ4mk7HWtJ|+^nsfO`Fq;YlOjuzJuQ$Nc7GtXw~IhJ
zva36H*eQs()h_aGKCaQCs1T+$j<K~T9xaL#IMjH;JdFr(Y7uIqkKAxL^v+!O6e;;0
z<AyzfQv7#$B?(DTrN+Qr9s`D^F+qXDRuGpSH_PoJZm<c&4V8iQv`8-lUlaR)uaePQ
zP)Ts)6xv1PEJw189M&wv%YZ|<fm-n6V3{X`OnO7vV;AM%CoPKRD1m^A4FPpIY<ur1
z8fc-zmY@_ryejP>@bZc|39l*WC&J6*5b)AB9QW&SyQl^g+CA;@4iA^A!97R`z>7v2
zyb>KcyhQJAc+Gr2k+JS@exC*}H=O$5zbp^OiSSw9cs|78PxAwz^S)thhJUKqL0=|p
zE)sfh9Gq+KW{V;J!+aX(A0-3zu7*lCnj;M#;1H;+*q_bas02}&=@6*gNhXWj-UN;V
z1<C;8(NNhHHw1J0#sfRe9}+cAZ5Ec(Gv9L`J>lLf#osG5KNV8&nNQE^r$|tODpl)Q
zeQhu>ujo8FI4@yn;)^*pq2T<(zfi*LmP6$AR#BnoY;p*6UP9G6I@?0gLI1Zl3J##N
z6V2#|?+9m(W#CirjlX0@cWzNCSw+CahLSoq)khXnBa`|Cl#NOKtV1OA7mxNPPxJP5
zVpD|$@<m`$l<*LnfQKc=+^j0tm^nel%xl>g6oAICN#H%6hQI5NS8^oDASr}CN7Kj>
zOvwD7c;JCL_Y@HWb?!-RQpU_)+<ts^{0P}8ADZduv9%5eJj)Ywq=`eFz-rTMT_sTH
zjN$sp7(;8cG!Q>D(&Fdu)D!7h^zMwGc=&aD0@UB`qdk#d@bg{t^NC&vpU?JtBj?U&
z*#GC`J?+`rFxP)@bxCo__L8mYkJMdURc)DFRPLewq$=sDTBQ|vCEN3Ge%+-o)}d?+
zzFDjEw~;!O9!u%{T%C>2Un;$Ds;lwwLZvq;^*#J>ROu~1wPB^N(zl-49ScR3>T>mM
o`l0ewO{OY?sxQ?_$%C#BH4}G9bvxEQXJEH^YV^8i!r9FK1L>V*^Z)<=

diff --git a/docs/build/html/.buildinfo b/docs/build/html/.buildinfo
index 7803520..3a339db 100644
--- a/docs/build/html/.buildinfo
+++ b/docs/build/html/.buildinfo
@@ -1,4 +1,4 @@
 # Sphinx build info version 1
 # This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done.
-config: ece5e90ad69b60bc3c014a5b8ba25dcf
+config: ea54a8d4e80d8f2678387db897b0d1fd
 tags: 645f666f9bcd5a90fca523b33c5a78b7
diff --git a/docs/build/html/_modules/index.html b/docs/build/html/_modules/index.html
index 1bb5074..4ac44dd 100644
--- a/docs/build/html/_modules/index.html
+++ b/docs/build/html/_modules/index.html
@@ -3,7 +3,7 @@
 <head>
   <meta charset="utf-8" />
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
-  <title>Overview: module code &mdash; Geocode 0.11.1 documentation</title>
+  <title>Overview: module code &mdash; Geocode 0.12.0 documentation</title>
       <link rel="stylesheet" href="../_static/pygments.css" type="text/css" />
       <link rel="stylesheet" href="../_static/css/theme.css" type="text/css" />
   <!--[if lt IE 9]>
diff --git a/docs/build/html/_static/documentation_options.js b/docs/build/html/_static/documentation_options.js
index 5c6f190..9390130 100644
--- a/docs/build/html/_static/documentation_options.js
+++ b/docs/build/html/_static/documentation_options.js
@@ -1,6 +1,6 @@
 var DOCUMENTATION_OPTIONS = {
     URL_ROOT: document.getElementById("documentation_options").getAttribute('data-url_root'),
-    VERSION: '0.11.1',
+    VERSION: '0.12.0',
     LANGUAGE: 'en',
     COLLAPSE_INDEX: false,
     BUILDER: 'html',
diff --git a/docs/build/html/genindex.html b/docs/build/html/genindex.html
index 322cd91..25c7ebb 100644
--- a/docs/build/html/genindex.html
+++ b/docs/build/html/genindex.html
@@ -3,7 +3,7 @@
 <head>
   <meta charset="utf-8" />
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
-  <title>Index &mdash; Geocode 0.11.1 documentation</title>
+  <title>Index &mdash; Geocode 0.12.0 documentation</title>
       <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
       <link rel="stylesheet" href="_static/css/theme.css" type="text/css" />
   <!--[if lt IE 9]>
@@ -113,11 +113,13 @@ <h2 id="R">R</h2>
   <td style="width: 33%; vertical-align: top;"><ul>
       <li><a href="modules.html#geocode.Geocoder.reverse_geocode">reverse_geocode() (geocode.Geocoder method)</a>
 </li>
-  </ul></td>
-  <td style="width: 33%; vertical-align: top;"><ul>
       <li><a href="modules.html#geocode.Geocoder.reverse_geocode_gsp">reverse_geocode_gsp() (geocode.Geocoder method)</a>
 </li>
+  </ul></td>
+  <td style="width: 33%; vertical-align: top;"><ul>
       <li><a href="modules.html#geocode.Geocoder.reverse_geocode_llsoa">reverse_geocode_llsoa() (geocode.Geocoder method)</a>
+</li>
+      <li><a href="modules.html#geocode.Geocoder.reverse_geocode_nuts">reverse_geocode_nuts() (geocode.Geocoder method)</a>
 </li>
   </ul></td>
 </tr></table>
diff --git a/docs/build/html/index.html b/docs/build/html/index.html
index f3f2a28..f33482d 100644
--- a/docs/build/html/index.html
+++ b/docs/build/html/index.html
@@ -4,7 +4,7 @@
   <meta charset="utf-8" /><meta name="generator" content="Docutils 0.17.1: http://docutils.sourceforge.net/" />
 
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
-  <title>Welcome to Geocode’s documentation! &mdash; Geocode 0.11.1 documentation</title>
+  <title>Welcome to Geocode’s documentation! &mdash; Geocode 0.12.0 documentation</title>
       <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
       <link rel="stylesheet" href="_static/css/theme.css" type="text/css" />
   <!--[if lt IE 9]>
diff --git a/docs/build/html/modules.html b/docs/build/html/modules.html
index 3105484..52d39f5 100644
--- a/docs/build/html/modules.html
+++ b/docs/build/html/modules.html
@@ -4,7 +4,7 @@
   <meta charset="utf-8" /><meta name="generator" content="Docutils 0.17.1: http://docutils.sourceforge.net/" />
 
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
-  <title>Modules &mdash; Geocode 0.11.1 documentation</title>
+  <title>Modules &mdash; Geocode 0.12.0 documentation</title>
       <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
       <link rel="stylesheet" href="_static/css/theme.css" type="text/css" />
   <!--[if lt IE 9]>
@@ -56,6 +56,7 @@
 <li class="toctree-l4"><a class="reference internal" href="#geocode.Geocoder.reverse_geocode"><code class="docutils literal notranslate"><span class="pre">Geocoder.reverse_geocode()</span></code></a></li>
 <li class="toctree-l4"><a class="reference internal" href="#geocode.Geocoder.reverse_geocode_gsp"><code class="docutils literal notranslate"><span class="pre">Geocoder.reverse_geocode_gsp()</span></code></a></li>
 <li class="toctree-l4"><a class="reference internal" href="#geocode.Geocoder.reverse_geocode_llsoa"><code class="docutils literal notranslate"><span class="pre">Geocoder.reverse_geocode_llsoa()</span></code></a></li>
+<li class="toctree-l4"><a class="reference internal" href="#geocode.Geocoder.reverse_geocode_nuts"><code class="docutils literal notranslate"><span class="pre">Geocoder.reverse_geocode_nuts()</span></code></a></li>
 </ul>
 </li>
 </ul>
@@ -104,13 +105,13 @@ <h1><a class="toc-backref" href="#id1">Modules</a><a class="headerlink" href="#m
 <h2><a class="toc-backref" href="#id2">geocode</a><a class="headerlink" href="#geocode" title="Permalink to this heading"></a></h2>
 <dl class="py class">
 <dt class="sig sig-object py" id="geocode.Geocoder">
-<em class="property"><span class="pre">class</span><span class="w"> </span></em><span class="sig-prename descclassname"><span class="pre">geocode.</span></span><span class="sig-name descname"><span class="pre">Geocoder</span></span><span class="sig-paren">(</span><em class="sig-param"><span class="n"><span class="pre">cache_dir</span></span><span class="p"><span class="pre">:</span></span><span class="w"> </span><span class="n"><span class="pre">Optional</span><span class="p"><span class="pre">[</span></span><span class="pre">str</span><span class="p"><span class="pre">]</span></span></span><span class="w"> </span><span class="o"><span class="pre">=</span></span><span class="w"> </span><span class="default_value"><span class="pre">None</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">gmaps_key_file</span></span><span class="p"><span class="pre">:</span></span><span class="w"> </span><span class="n"><span class="pre">Optional</span><span class="p"><span class="pre">[</span></span><span class="pre">str</span><span class="p"><span class="pre">]</span></span></span><span class="w"> </span><span class="o"><span class="pre">=</span></span><span class="w"> </span><span class="default_value"><span class="pre">None</span></span></em><span class="sig-paren">)</span><a class="reference internal" href="_modules/geocode/geocode.html#Geocoder"><span class="viewcode-link"><span class="pre">[source]</span></span></a><a class="headerlink" href="#geocode.Geocoder" title="Permalink to this definition"></a></dt>
+<em class="property"><span class="pre">class</span><span class="w"> </span></em><span class="sig-prename descclassname"><span class="pre">geocode.</span></span><span class="sig-name descname"><span class="pre">Geocoder</span></span><span class="sig-paren">(</span><em class="sig-param"><span class="n"><span class="pre">cache_dir</span></span><span class="p"><span class="pre">:</span></span><span class="w"> </span><span class="n"><span class="pre">Optional</span><span class="p"><span class="pre">[</span></span><span class="pre">Path</span><span class="p"><span class="pre">]</span></span></span><span class="w"> </span><span class="o"><span class="pre">=</span></span><span class="w"> </span><span class="default_value"><span class="pre">None</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">gmaps_key_file</span></span><span class="p"><span class="pre">:</span></span><span class="w"> </span><span class="n"><span class="pre">Optional</span><span class="p"><span class="pre">[</span></span><span class="pre">Path</span><span class="p"><span class="pre">]</span></span></span><span class="w"> </span><span class="o"><span class="pre">=</span></span><span class="w"> </span><span class="default_value"><span class="pre">None</span></span></em><span class="sig-paren">)</span><a class="reference internal" href="_modules/geocode/geocode.html#Geocoder"><span class="viewcode-link"><span class="pre">[source]</span></span></a><a class="headerlink" href="#geocode.Geocoder" title="Permalink to this definition"></a></dt>
 <dd><p>Bases: <code class="xref py py-class docutils literal notranslate"><span class="pre">object</span></code></p>
 <p>Geocode addresses, postcodes, LLSOAs or Constituencies or reverse-geocode latitudes and
 longitudes.</p>
 <dl class="py method">
 <dt class="sig sig-object py" id="geocode.Geocoder.force_setup">
-<span class="sig-name descname"><span class="pre">force_setup</span></span><span class="sig-paren">(</span><em class="sig-param"><span class="n"><span class="pre">ngeso_setup</span></span><span class="o"><span class="pre">=</span></span><span class="default_value"><span class="pre">True</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">cpo_setup</span></span><span class="o"><span class="pre">=</span></span><span class="default_value"><span class="pre">True</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">ons_setup</span></span><span class="o"><span class="pre">=</span></span><span class="default_value"><span class="pre">True</span></span></em><span class="sig-paren">)</span><a class="reference internal" href="_modules/geocode/geocode.html#Geocoder.force_setup"><span class="viewcode-link"><span class="pre">[source]</span></span></a><a class="headerlink" href="#geocode.Geocoder.force_setup" title="Permalink to this definition"></a></dt>
+<span class="sig-name descname"><span class="pre">force_setup</span></span><span class="sig-paren">(</span><em class="sig-param"><span class="n"><span class="pre">ngeso_setup</span></span><span class="o"><span class="pre">=</span></span><span class="default_value"><span class="pre">True</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">cpo_setup</span></span><span class="o"><span class="pre">=</span></span><span class="default_value"><span class="pre">True</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">ons_setup</span></span><span class="o"><span class="pre">=</span></span><span class="default_value"><span class="pre">True</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">eurostat_setup</span></span><span class="o"><span class="pre">=</span></span><span class="default_value"><span class="pre">True</span></span></em><span class="sig-paren">)</span><a class="reference internal" href="_modules/geocode/geocode.html#Geocoder.force_setup"><span class="viewcode-link"><span class="pre">[source]</span></span></a><a class="headerlink" href="#geocode.Geocoder.force_setup" title="Permalink to this definition"></a></dt>
 <dd><p>Download all data and setup caches.</p>
 </dd></dl>
 
@@ -245,6 +246,22 @@ <h2><a class="toc-backref" href="#id2">geocode</a><a class="headerlink" href="#g
 </dl>
 </dd></dl>
 
+<dl class="py method">
+<dt class="sig sig-object py" id="geocode.Geocoder.reverse_geocode_nuts">
+<span class="sig-name descname"><span class="pre">reverse_geocode_nuts</span></span><span class="sig-paren">(</span><em class="sig-param"><span class="n"><span class="pre">latlons</span></span><span class="p"><span class="pre">:</span></span><span class="w"> </span><span class="n"><span class="pre">List</span><span class="p"><span class="pre">[</span></span><span class="pre">Tuple</span><span class="p"><span class="pre">[</span></span><span class="pre">float</span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="pre">float</span><span class="p"><span class="pre">]</span></span><span class="p"><span class="pre">]</span></span></span></em>, <em class="sig-param"><span class="n"><span class="pre">level</span></span><span class="p"><span class="pre">:</span></span><span class="w"> </span><span class="n"><span class="pre">Literal</span><span class="p"><span class="pre">[</span></span><span class="m"><span class="pre">0</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">1</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">2</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">3</span></span><span class="p"><span class="pre">]</span></span></span></em>, <em class="sig-param"><span class="n"><span class="pre">year</span></span><span class="p"><span class="pre">:</span></span><span class="w"> </span><span class="n"><span class="pre">Literal</span><span class="p"><span class="pre">[</span></span><span class="m"><span class="pre">2003</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">2006</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">2010</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">2013</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">2016</span></span><span class="p"><span class="pre">,</span></span><span class="w"> </span><span class="m"><span class="pre">2021</span></span><span class="p"><span class="pre">]</span></span></span><span class="w"> </span><span class="o"><span class="pre">=</span></span><span class="w"> </span><span class="default_value"><span class="pre">2021</span></span></em><span class="sig-paren">)</span> <span class="sig-return"><span class="sig-return-icon">&#x2192;</span> <span class="sig-return-typehint"><span class="pre">List</span><span class="p"><span class="pre">[</span></span><span class="pre">str</span><span class="p"><span class="pre">]</span></span></span></span><a class="reference internal" href="_modules/geocode/geocode.html#Geocoder.reverse_geocode_nuts"><span class="viewcode-link"><span class="pre">[source]</span></span></a><a class="headerlink" href="#geocode.Geocoder.reverse_geocode_nuts" title="Permalink to this definition"></a></dt>
+<dd><p>Function to reverse geocode a collection of latlons into NUTS boundaries.</p>
+<dl class="field-list simple">
+<dt class="field-odd">Parameters</dt>
+<dd class="field-odd"><ul class="simple">
+<li><p><strong>latlons</strong> (<em>iterable</em><em> of </em><em>strings</em>) – Specific latlons to geocode to llsoa boundaries.</p></li>
+<li><p><strong>level</strong> (<em>int</em>) – Specify the NUTS level, must be one of [0,1,2,3].</p></li>
+<li><p><strong>year</strong> (<em>int</em>) – Specify the year of NUTS regulation, must be one of [2003,2006,2010,2013,2016,2021],
+defaults to 2021.</p></li>
+</ul>
+</dd>
+</dl>
+</dd></dl>
+
 </dd></dl>
 
 </section>
diff --git a/docs/build/html/objects.inv b/docs/build/html/objects.inv
index a0dbfe1fdeb98307ddb7f505bf4c7fe79464ab4b..0a405ad88d05d17dc46d2671a5614a8a9ddc9612 100644
GIT binary patch
delta 274
zcmV+t0qy>f1C|4jcYi@n#4rp+_dP`rY+H%#TmZ2^8i@tVD3h49QW{sWgXoTHa6L{!
z(so3rfEY5XM$UiV&X;QeH5zE{=5G+gxYJG%iCfX5i4Le~xp&t$qOR;l+&p8@APG1Q
z8SH^uG-r`Zp^TCmJ#ohoydDcJoFi1FENbPXiU|xImrgJsaZFJLER9DQz;-Ey0@~7W
z0X+pL4;Z~sp@aWg1bBf!Aj=kDL+7Vl=uol;TMr-k%oyyizgxYbzA3Yb@4SImL2PzT
z9l&MkSe-C&pO2F~0UA2A7gtr9<%eFwSRXHDvxjNqrRy^Py6^f)i<qjF%OCdpF$L;R
YyAY_L+f#6^&CKI`zgHHjPi6hTc(6}>ApigX

delta 290
zcmV+-0p0$V1CRrdcYnc7#2^rc?|ll4Ufaa$zJT#yvWW+egoTkdp-jUJ#@>7lU(Y8I
z+HUMtjV^nF%>3U!!vGdgqk-mbUV|9Mopy>y+=?DebU;nZy}P~<b!9i=<{5(qNx*T)
zU=Q4)`4_na$|$MP6L%cJ>#@+nIYL#+qE=3-n846+X$J!m7k|xwrST{O*e;EsfVMPT
zKu?2{2aMjR(80MD0bU>w$g&5pq4QBLbU5Q3*kAujy`jD-vx)D#fmcCnc1|6@W$GB6
zFma!^)tAhm&Z3p#37p0r1fy8IKi=L+Bw_5qv*L8AS5@lehg!o}9}~0L!_@Lpb?JXy
ocYUNqOx4Qe4;%aam<H-ks}QK5+f#6k&CKJR-zy8%C$;gX6~4EV^8f$<

diff --git a/docs/build/html/search.html b/docs/build/html/search.html
index 24f504f..904f543 100644
--- a/docs/build/html/search.html
+++ b/docs/build/html/search.html
@@ -3,7 +3,7 @@
 <head>
   <meta charset="utf-8" />
   <meta name="viewport" content="width=device-width, initial-scale=1.0" />
-  <title>Search &mdash; Geocode 0.11.1 documentation</title>
+  <title>Search &mdash; Geocode 0.12.0 documentation</title>
       <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
       <link rel="stylesheet" href="_static/css/theme.css" type="text/css" />
     
diff --git a/docs/build/html/searchindex.js b/docs/build/html/searchindex.js
index 8484e7f..966e884 100644
--- a/docs/build/html/searchindex.js
+++ b/docs/build/html/searchindex.js
@@ -1 +1 @@
-Search.setIndex({"docnames": ["index", "modules"], "filenames": ["index.rst", "modules.rst"], "titles": ["Welcome to Geocode\u2019s documentation!", "Modules"], "terms": {"modul": 0, "index": 0, "search": 0, "page": 0, "class": 1, "cache_dir": 1, "option": 1, "str": 1, "none": 1, "gmaps_key_fil": 1, "sourc": 1, "base": 1, "object": 1, "address": 1, "postcod": 1, "llsoa": 1, "constitu": 1, "revers": 1, "latitud": 1, "longitud": 1, "force_setup": 1, "ngeso_setup": 1, "true": 1, "cpo_setup": 1, "ons_setup": 1, "download": 1, "all": 1, "data": 1, "setup": 1, "cach": 1, "entity_id": 1, "entiti": 1, "kwarg": 1, "select": 1, "gsp": 1, "region": 1, "boundari": 1, "lad": 1, "paramet": 1, "string": 1, "specifi": 1, "type": 1, "from": 1, "i": 1, "e": 1, "iter": 1, "The": 1, "specif": 1, "geocode_constitu": 1, "function": 1, "collect": 1, "latlon": 1, "geocode_llsoa": 1, "llsoa_boundari": 1, "geocode_local_author": 1, "local": 1, "author": 1, "district": 1, "geocode_postcod": 1, "method": 1, "cpo": 1, "get_dno_region": 1, "get": 1, "dno": 1, "licens": 1, "area": 1, "eso": 1, "portal": 1, "return": 1, "dno_region": 1, "dict": 1, "whose": 1, "kei": 1, "ar": 1, "id": 1, "valu": 1, "tupl": 1, "contain": 1, "region_boundari": 1, "region_bound": 1, "shape": 1, "polygon": 1, "multipolygon": 1, "bound": 1, "xmin": 1, "ymin": 1, "xmax": 1, "ymax": 1, "dno_nam": 1, "name": 1, "longnam": 1, "get_gsp_region": 1, "20220314": 1, "gnode": 1, "api": 1, "get_llsoa_boundari": 1, "load": 1, "either": 1, "avail": 1, "els": 1, "fetch": 1, "raw": 1, "england": 1, "wale": 1, "packag": 1, "scotland": 1, "reverse_geocod": 1, "set": 1, "list": 1, "A": 1, "reverse_geocode_gsp": 1, "reverse_geocode_llsoa": 1, "dz": 1, "boolean": 1, "indic": 1, "whether": 1, "consid": 1, "datazon": 1, "pass": 1, "underli": 1}, "objects": {"geocode": [[1, 0, 1, "", "Geocoder"]], "geocode.Geocoder": [[1, 1, 1, "", "force_setup"], [1, 1, 1, "", "geocode"], [1, 1, 1, "", "geocode_constituency"], [1, 1, 1, "", "geocode_llsoa"], [1, 1, 1, "", "geocode_local_authority"], [1, 1, 1, "", "geocode_postcode"], [1, 1, 1, "", "get_dno_regions"], [1, 1, 1, "", "get_gsp_regions"], [1, 1, 1, "", "get_llsoa_boundaries"], [1, 1, 1, "", "reverse_geocode"], [1, 1, 1, "", "reverse_geocode_gsp"], [1, 1, 1, "", "reverse_geocode_llsoa"]]}, "objtypes": {"0": "py:class", "1": "py:method"}, "objnames": {"0": ["py", "class", "Python class"], "1": ["py", "method", "Python method"]}, "titleterms": {"welcom": 0, "geocod": [0, 1], "": 0, "document": 0, "content": [0, 1], "indic": 0, "tabl": 0, "modul": 1}, "envversion": {"sphinx.domains.c": 2, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 8, "sphinx.domains.index": 1, "sphinx.domains.javascript": 2, "sphinx.domains.math": 2, "sphinx.domains.python": 3, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "sphinx.ext.viewcode": 1, "sphinx": 57}, "alltitles": {"Welcome to Geocode\u2019s documentation!": [[0, "welcome-to-geocode-s-documentation"]], "Contents:": [[0, null]], "Indices and tables": [[0, "indices-and-tables"]], "Modules": [[1, "modules"]], "Contents": [[1, "contents"]], "geocode": [[1, "geocode"]]}, "indexentries": {"geocoder (class in geocode)": [[1, "geocode.Geocoder"]], "force_setup() (geocode.geocoder method)": [[1, "geocode.Geocoder.force_setup"]], "geocode() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode"]], "geocode_constituency() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_constituency"]], "geocode_llsoa() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_llsoa"]], "geocode_local_authority() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_local_authority"]], "geocode_postcode() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_postcode"]], "get_dno_regions() (geocode.geocoder method)": [[1, "geocode.Geocoder.get_dno_regions"]], "get_gsp_regions() (geocode.geocoder method)": [[1, "geocode.Geocoder.get_gsp_regions"]], "get_llsoa_boundaries() (geocode.geocoder method)": [[1, "geocode.Geocoder.get_llsoa_boundaries"]], "reverse_geocode() (geocode.geocoder method)": [[1, "geocode.Geocoder.reverse_geocode"]], "reverse_geocode_gsp() (geocode.geocoder method)": [[1, "geocode.Geocoder.reverse_geocode_gsp"]], "reverse_geocode_llsoa() (geocode.geocoder method)": [[1, "geocode.Geocoder.reverse_geocode_llsoa"]]}})
\ No newline at end of file
+Search.setIndex({"docnames": ["index", "modules"], "filenames": ["index.rst", "modules.rst"], "titles": ["Welcome to Geocode\u2019s documentation!", "Modules"], "terms": {"modul": 0, "index": 0, "search": 0, "page": 0, "class": 1, "cache_dir": 1, "option": 1, "str": 1, "none": 1, "gmaps_key_fil": 1, "sourc": 1, "base": 1, "object": 1, "address": 1, "postcod": 1, "llsoa": 1, "constitu": 1, "revers": 1, "latitud": 1, "longitud": 1, "force_setup": 1, "ngeso_setup": 1, "true": 1, "cpo_setup": 1, "ons_setup": 1, "download": 1, "all": 1, "data": 1, "setup": 1, "cach": 1, "entity_id": 1, "entiti": 1, "kwarg": 1, "select": 1, "gsp": 1, "region": 1, "boundari": 1, "lad": 1, "paramet": 1, "string": 1, "specifi": 1, "type": 1, "from": 1, "i": 1, "e": 1, "iter": 1, "The": 1, "specif": 1, "geocode_constitu": 1, "function": 1, "collect": 1, "latlon": 1, "geocode_llsoa": 1, "llsoa_boundari": 1, "geocode_local_author": 1, "local": 1, "author": 1, "district": 1, "geocode_postcod": 1, "method": 1, "cpo": 1, "get_dno_region": 1, "get": 1, "dno": 1, "licens": 1, "area": 1, "eso": 1, "portal": 1, "return": 1, "dno_region": 1, "dict": 1, "whose": 1, "kei": 1, "ar": 1, "id": 1, "valu": 1, "tupl": 1, "contain": 1, "region_boundari": 1, "region_bound": 1, "shape": 1, "polygon": 1, "multipolygon": 1, "bound": 1, "xmin": 1, "ymin": 1, "xmax": 1, "ymax": 1, "dno_nam": 1, "name": 1, "longnam": 1, "get_gsp_region": 1, "20220314": 1, "gnode": 1, "api": 1, "get_llsoa_boundari": 1, "load": 1, "either": 1, "avail": 1, "els": 1, "fetch": 1, "raw": 1, "england": 1, "wale": 1, "packag": 1, "scotland": 1, "reverse_geocod": 1, "set": 1, "list": 1, "A": 1, "reverse_geocode_gsp": 1, "reverse_geocode_llsoa": 1, "dz": 1, "boolean": 1, "indic": 1, "whether": 1, "consid": 1, "datazon": 1, "pass": 1, "underli": 1, "path": 1, "eurostat_setup": 1, "reverse_geocode_nut": 1, "float": 1, "level": 1, "liter": 1, "0": 1, "1": 1, "2": 1, "3": 1, "year": 1, "2003": 1, "2006": 1, "2010": 1, "2013": 1, "2016": 1, "2021": 1, "nut": 1, "int": 1, "must": 1, "one": 1, "regul": 1, "default": 1}, "objects": {"geocode": [[1, 0, 1, "", "Geocoder"]], "geocode.Geocoder": [[1, 1, 1, "", "force_setup"], [1, 1, 1, "", "geocode"], [1, 1, 1, "", "geocode_constituency"], [1, 1, 1, "", "geocode_llsoa"], [1, 1, 1, "", "geocode_local_authority"], [1, 1, 1, "", "geocode_postcode"], [1, 1, 1, "", "get_dno_regions"], [1, 1, 1, "", "get_gsp_regions"], [1, 1, 1, "", "get_llsoa_boundaries"], [1, 1, 1, "", "reverse_geocode"], [1, 1, 1, "", "reverse_geocode_gsp"], [1, 1, 1, "", "reverse_geocode_llsoa"], [1, 1, 1, "", "reverse_geocode_nuts"]]}, "objtypes": {"0": "py:class", "1": "py:method"}, "objnames": {"0": ["py", "class", "Python class"], "1": ["py", "method", "Python method"]}, "titleterms": {"welcom": 0, "geocod": [0, 1], "": 0, "document": 0, "content": [0, 1], "indic": 0, "tabl": 0, "modul": 1}, "envversion": {"sphinx.domains.c": 2, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 8, "sphinx.domains.index": 1, "sphinx.domains.javascript": 2, "sphinx.domains.math": 2, "sphinx.domains.python": 3, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "sphinx.ext.viewcode": 1, "sphinx": 57}, "alltitles": {"Welcome to Geocode\u2019s documentation!": [[0, "welcome-to-geocode-s-documentation"]], "Contents:": [[0, null]], "Indices and tables": [[0, "indices-and-tables"]], "Modules": [[1, "modules"]], "Contents": [[1, "contents"]], "geocode": [[1, "geocode"]]}, "indexentries": {"geocoder (class in geocode)": [[1, "geocode.Geocoder"]], "force_setup() (geocode.geocoder method)": [[1, "geocode.Geocoder.force_setup"]], "geocode() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode"]], "geocode_constituency() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_constituency"]], "geocode_llsoa() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_llsoa"]], "geocode_local_authority() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_local_authority"]], "geocode_postcode() (geocode.geocoder method)": [[1, "geocode.Geocoder.geocode_postcode"]], "get_dno_regions() (geocode.geocoder method)": [[1, "geocode.Geocoder.get_dno_regions"]], "get_gsp_regions() (geocode.geocoder method)": [[1, "geocode.Geocoder.get_gsp_regions"]], "get_llsoa_boundaries() (geocode.geocoder method)": [[1, "geocode.Geocoder.get_llsoa_boundaries"]], "reverse_geocode() (geocode.geocoder method)": [[1, "geocode.Geocoder.reverse_geocode"]], "reverse_geocode_gsp() (geocode.geocoder method)": [[1, "geocode.Geocoder.reverse_geocode_gsp"]], "reverse_geocode_llsoa() (geocode.geocoder method)": [[1, "geocode.Geocoder.reverse_geocode_llsoa"]], "reverse_geocode_nuts() (geocode.geocoder method)": [[1, "geocode.Geocoder.reverse_geocode_nuts"]]}})
\ No newline at end of file
diff --git a/docs/source/conf.py b/docs/source/conf.py
index b9493ea..7320a73 100644
--- a/docs/source/conf.py
+++ b/docs/source/conf.py
@@ -22,7 +22,7 @@
 author = "Jamie Taylor"
 
 # The full version, including alpha/beta/rc tags
-release = "0.11.1"
+release = "0.12.0"
 
 
 # -- General configuration ---------------------------------------------------
diff --git a/example.py b/example.py
index da052bf..914453c 100644
--- a/example.py
+++ b/example.py
@@ -6,43 +6,47 @@
 def main():
     with Geocoder() as geocoder:
         # Geocode some postcodes / addresses...
+        print("GEOCODE POSTCODES / ADDRESSES:")
         postcodes = ["S3 7RH", "S3 7", "S3", None, None, "S3 7RH"]
         addresses = [None, None, None, "Hicks Building, Sheffield", "Hicks", "Hicks Building"]
         results = geocoder.geocode(postcodes, "postcode", address=addresses)
         for postcode, address, (lat, lon, status) in zip(postcodes, addresses, results):
-            print(f"Postcode: `{postcode}`    Address: `{address}`")
-            if status == 0:
-                print("    Failed to geocode!")
-            else:
-                print(f"    {lat:.3f}, {lon:.3f}    ->  {geocoder.status_codes[status]}")
+            print(f"    Postcode + Address: `{postcode}` + `{address}`  ->  {lat:.3f}, {lon:.3f} "
+                  f"({geocoder.status_codes[status]})")
         # Geocode some LLSOAs...
+        print("GEOCODE LLSOAs:")
         llsoas = ["E01033264", "E01033262"]
         results = geocoder.geocode_llsoa(llsoas)
         for llsoa, (lat, lon) in zip(llsoas, results):
-            print(f"LLSOA: `{llsoa}`")
-            print(f"    {lat:.3f}, {lon:.3f}")
+            print(f"    LLSOA: `{llsoa}`  ->  {lat:.3f}, {lon:.3f}")
+        # Geocode some Constituencies...
+        print("GEOCODE CONSTITUENCIES:")
+        constituencies = ["Sheffield Central", "Sheffield Hallam"]
+        results = geocoder.geocode_constituency(constituencies)
+        for constituency, (lat, lon) in zip(constituencies, results):
+            print(f"    Constituency: `{constituency}`  ->  {lat:.3f}, {lon:.3f}")
         # Reverse-geocode some lat/lons to LLSOAs...
+        print("REVERSE-GEOCODE TO LLSOA:")
         latlons = [(53.384, -1.467), (53.388, -1.470)]
         results = geocoder.reverse_geocode_llsoa(latlons)
         for llsoa, (lat, lon) in zip(results, latlons):
-            print(f"LATLON: {lat:.3f}, {lon:.3f}:")
-            print(f"    `{llsoa}`")
-        # Reverse-geocode some lat/lons to GSP/GNode...
+            print(f"    LATLON: {lat:.3f}, {lon:.3f}  ->  {llsoa}")
+        # Reverse-geocode some lat/lons to GSP...
+        print("REVERSE-GEOCODE TO GSP:")
         latlons = [(53.384, -1.467), (53.388, -1.470)]
         results = geocoder.reverse_geocode_gsp(latlons)
         for (lat, lon), region_id in zip(latlons, results):
-            print(f"LATLON: {lat:.3f}, {lon:.3f}:")
-            print(f"    {region_id}")
-        # Geocode some Constituencies...
-        constituencies = ["Sheffield Central", "Sheffield Hallam"]
-        results = geocoder.geocode_constituency(constituencies)
-        for constituency, (lat, lon) in zip(constituencies, results):
-            print(f"Constituency: `{constituency}`")
-            print(f"    {lat:.3f}, {lon:.3f}")
+            print(f"    LATLON: {lat:.3f}, {lon:.3f}  ->  {region_id}")
+        # Reverse-geocode some lat/lons to 2021 NUTS2...
+        print("REVERSE-GEOCODE TO NUTS2:")
+        latlons = [(51.3259, -1.9613), (47.9995, 0.2335), (50.8356, 8.7343)]
+        results = geocoder.reverse_geocode_nuts(latlons, year=2021, level=2)
+        for (lat, lon), nuts2 in zip(latlons, results):
+            print(f"    LATLON: {lat:.3f}, {lon:.3f}  ->  {nuts2}")
 
 if __name__ == "__main__":
     log_fmt = "%(asctime)s [%(levelname)s] [%(filename)s:%(funcName)s] - %(message)s"
     fmt = os.environ.get("GEOCODE_LOGGING_FMT", log_fmt)
     datefmt = os.environ.get("GEOCODE_LOGGING_DATEFMT", "%Y-%m-%dT%H:%M:%SZ")
-    logging.basicConfig(format=fmt, datefmt=datefmt, level=os.environ.get("LOGLEVEL", "DEBUG"))
+    logging.basicConfig(format=fmt, datefmt=datefmt, level=os.environ.get("LOGLEVEL", "WARNING"))
     main()
\ No newline at end of file
diff --git a/geocode/bng2latlon.py b/geocode/bng2latlon.py
index 3e8329f..78cd3c8 100644
--- a/geocode/bng2latlon.py
+++ b/geocode/bng2latlon.py
@@ -8,9 +8,10 @@
 """
 
 import sys
-import os
 import argparse
 import time as TIME
+from pathlib import Path
+
 import pandas as pd
 
 from geocode import Geocoder
@@ -21,17 +22,17 @@ def parse_options():
     parser = argparse.ArgumentParser(description=("This is a command line interface (CLI) for "
                                                   "the bng2latlon.py module"),
                                      epilog="Jamie Taylor & Ethan Jones, 2020-06-12")
-    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=str,
+    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=Path,
                         required=True, metavar="</path/to/file>",
                         help="Specify a CSV file containing a list of eastings and northings to "
                              "be converted. The file must contain the columns 'eastings' "
                              "and 'northings' (it can contain others, all of which will be kept).")
-    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=str,
+    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=Path,
                         required=True, metavar="</path/to/file>", help="Specify an output file.")
     options = parser.parse_args()
-    if not os.path.isfile(options.infile):
+    if not options.infile.is_file():
         raise GenericException(f"The input file '{options.infile}' does not exist.")
-    if os.path.isfile(options.outfile):
+    if options.outfile.is_file():
         check = query_yes_no(f"The outfile '{options.outfile}' already exists, are you sure you "
                              "wish to overwrite?", "no")
         if not check:
diff --git a/geocode/cache_manager.py b/geocode/cache_manager.py
index 95700fa..c70d805 100644
--- a/geocode/cache_manager.py
+++ b/geocode/cache_manager.py
@@ -10,16 +10,17 @@
 import pickle
 import glob
 import logging
+from pathlib import Path
 import errno
 from typing import Any, Optional
 
-SCRIPT_DIR = os.path.dirname(os.path.realpath(__file__))
+SCRIPT_DIR = Path(os.path.dirname(os.path.realpath(__file__)))
 
 from . version import __version__
 
 class CacheManager:
     """Cache Python variables to files using Pickle."""
-    def __init__(self, cache_dir: Optional[str] = None):
+    def __init__(self, cache_dir: Optional[Path] = None):
         """
         Cache Python variables to files using Pickle.
 
@@ -29,15 +30,15 @@ def __init__(self, cache_dir: Optional[str] = None):
             Path to a directory to use for writing cache files. If not set, or set to None, default
             location will be `os.path.join(os.path.dirname(os.path.realpath(__file__)), "cache")`.
         """
-        self.cache_dir = os.path.join(SCRIPT_DIR, "cache") if cache_dir is None else cache_dir
-        if not os.path.isdir(self.cache_dir):
+        self.cache_dir = SCRIPT_DIR.joinpath("cache") if cache_dir is None else cache_dir
+        if not self.cache_dir.is_dir():
             raise FileNotFoundError(errno.ENOENT, os.strerror(errno.ENOENT), self.cache_dir)
         self.version_string = __version__.replace(".", "-")
 
     def _get_filename(self, label):
         """Generate a filename using `label` and the current version string."""
         file_name = f"{label}_{self.version_string}.p"
-        return os.path.join(self.cache_dir, file_name)
+        return self.cache_dir.joinpath(file_name)
 
     def retrieve(self, label: str) -> Any:
         """
@@ -55,7 +56,7 @@ def retrieve(self, label: str) -> Any:
             cache was not found.
         """
         cache_file = self._get_filename(label)
-        if os.path.isfile(cache_file):
+        if cache_file.is_file():
             with open(cache_file, "rb") as pickle_fid:
                 return pickle.load(pickle_fid)
         return None
@@ -93,15 +94,15 @@ def clear(self,
         """
         logging.debug("Deleting cache files (delete_gmaps_cache=%s, old_versions_only=%s)",
                       delete_gmaps_cache, old_versions_only)
-        cache_files = glob.glob(os.path.join(self.cache_dir, "*.p"))
+        cache_files = self.cache_dir.glob("*.p")
         for cache_file in cache_files:
-            if not delete_gmaps_cache and "gmaps" in cache_file:
+            if not delete_gmaps_cache and "gmaps" in cache_file.name:
                 continue
             if old_versions_only:
-                if self.version_string in cache_file:
+                if self.version_string in cache_file.name:
                     continue
             try:
-                os.remove(cache_file)
+                cache_file.unlink()
             except FileNotFoundError:
                 pass
             except:
diff --git a/geocode/cpo.py b/geocode/cpo.py
index b2e3b8d..fd6f116 100644
--- a/geocode/cpo.py
+++ b/geocode/cpo.py
@@ -11,21 +11,22 @@
 import pickle
 import zipfile
 import logging
+from pathlib import Path
+from typing import Optional, Iterable, Tuple, Union, List, Dict
+
 import pandas as pd
 import numpy as np
-from typing import Optional, Iterable, Tuple, Union, List, Dict
 
 from . utilities import GenericException, bng2latlon
 
-SCRIPT_DIR = os.path.dirname(os.path.realpath(__file__))
+SCRIPT_DIR = Path(os.path.dirname(os.path.realpath(__file__)))
 
 class CodePointOpen:
     """The Code Point Open data manager for the Geocode class."""
     def __init__(self, cache_manager):
         """The Code Point Open data manager for the Geocode class."""
         self.cache_manager = cache_manager
-        self.cpo_zipfile = os.path.join(SCRIPT_DIR, "code_point_open", "codepo_gb.zip")
-        self.cpo_cache_file = "code_point_open"
+        self.cpo_zipfile = SCRIPT_DIR.joinpath("code_point_open", "codepo_gb.zip")
         self.cpo = None
         self.cache = None
 
@@ -39,9 +40,10 @@ def _load(self, force_reload : bool = False):
         """Load the OS Code Point Open Database, either from raw zip file or local cache."""
         if self.cpo is not None and not force_reload:
             return
-        cpo_cache_contents = self.cache_manager.retrieve(self.cpo_cache_file)
+        cpo_cache_name = "code_point_open"
+        cpo_cache_contents = self.cache_manager.retrieve(cpo_cache_name)
         if cpo_cache_contents is not None and not force_reload:
-            logging.debug("Loading Code Point Open data from cache ('%s')", self.cpo_cache_file)
+            logging.debug("Loading Code Point Open data from cache ('%s')", cpo_cache_name)
             self.cpo = cpo_cache_contents
             return
         logging.info("Extracting the Code Point Open data (this only needs to be done once)")
@@ -72,8 +74,8 @@ def _load(self, force_reload : bool = False):
         cpo.loc[nn_indices, "latitude"] = lats
         cpo["outward_postcode"] = cpo["Postcode"].str.slice(0, -3).str.strip()
         cpo["inward_postcode"] = cpo["Postcode"].str.slice(-3).str.strip()
-        self.cache_manager.write(self.cpo_cache_file, cpo)
-        logging.info("Code Point Open extracted and pickled to '%s'", self.cpo_cache_file)
+        self.cache_manager.write(cpo_cache_name, cpo)
+        logging.info("Code Point Open extracted and pickled to cache")
         self.cpo = cpo
         return
     
diff --git a/geocode/eurostat.py b/geocode/eurostat.py
new file mode 100644
index 0000000..3c90eb4
--- /dev/null
+++ b/geocode/eurostat.py
@@ -0,0 +1,136 @@
+"""
+Manage data from Eurostat, primarily Nomenclature of territorial units for statistics (NUTS)
+regions.
+
+- Jamie Taylor <jamie.taylor@sheffield.ac.uk>
+- First Authored: 2023-03-14
+"""
+
+import os
+import sys
+import zipfile
+import json
+import csv
+import logging
+from pathlib import Path
+from typing import Optional, Iterable, Tuple, Union, List, Dict, Literal
+
+import pandas as pd
+import geopandas as gpd
+import shapefile
+
+from . import utilities as utils
+
+SCRIPT_DIR = Path(os.path.dirname(os.path.realpath(__file__)))
+
+class Eurostat:
+    """
+    Manage data from Eurostat.
+    """
+    def __init__(self, cache_manager):
+        """
+        Manage data from Eurostat.
+        """
+        self.cache_manager = cache_manager
+        self.data_dir = SCRIPT_DIR.joinpath("eurostat")
+        self.nuts_regions = {
+            (l, y): None for y in [2003, 2006, 2010, 2013, 2016, 2021] for l in range(4)
+        }
+
+    def force_setup(self):
+        """
+        Function to setup all lookup files.
+        """
+        for l in range(0, 4):
+            self._load_nuts_boundaries(l)
+
+    def load_nuts_boundaries(self,
+                             level: Literal[0, 1, 2, 3],
+                             year: Literal[2003, 2006, 2010, 2013, 2016, 2021] = 2021
+                             ) -> gpd.GeoDataFrame:
+        """
+        Load the NUTS boundaries, either from local cache if available, else fetch from Eurostat
+        API.
+
+        Parameters
+        ----------
+        `level` : int
+            Specify the NUTS level, must be one of [0,1,2,3].
+        `year` : int
+            Specify the year of NUTS regulation, must be one of [2003,2006,2010,2013,2016,2021],
+            defaults to 2021.
+
+        Returns
+        -------
+        Geopandas GeoDataFrame
+            Contains columns NUTS_ID, LEVL_CODE, CNTR_CODE, NAME_LATN, NUTS_NAME, MOUNT_TYPE,
+            URBN_TYPE, COAST_TYPE, FID, id, geometry
+        """
+        cache_label = f"nuts_{year}_{level}"
+        nuts_boundaries_cache_contents = self.cache_manager.retrieve(cache_label)
+        if nuts_boundaries_cache_contents is not None:
+            logging.debug("Loading %s NUTS%s boundaries from cache ('%s')",
+                          year, level, cache_label)
+            return nuts_boundaries_cache_contents
+        logging.info("Extracting the %s NUTS%s boundary data from Eurostat (this only needs to be "
+                     "done once)", year, level)
+        eurostat_url = ("https://gisco-services.ec.europa.eu/distribution/v2/nuts/geojson/"
+                        f"NUTS_RG_01M_{year}_4326_LEVL_{level}.geojson")
+        success, api_response = utils.fetch_from_api(eurostat_url)
+        if success:
+            raw = json.loads(api_response.text)
+            nuts_regions = gpd.GeoDataFrame.from_features(raw["features"],
+                                                          crs=raw["crs"]["properties"]["name"])
+            nuts_regions.geometry = nuts_regions.buffer(0)
+        else:
+            raise utils.GenericException("Encountered an error while extracting %s NUTS%s region "
+                                         "data from Eurostat API.")
+        self.cache_manager.write(cache_label, nuts_regions)
+        logging.info("%s NUTS%s boundaries extracted and pickled to file ('%s')",
+                     year, level, cache_label)
+        return nuts_regions
+
+    def _load_nuts_boundaries(self, level, year=2021):
+        """
+        For backwards compatibility pending https://github.com/SheffieldSolar/Geocode/issues/6
+
+        Load the NUTS boundaries, either from local cache if available, else fetch from Eurostat
+        API.
+        """
+        nuts_gdf = self.load_nuts_boundaries(level, year)
+        nuts_gdf["bounds"] = nuts_gdf.bounds.apply(tuple, axis=1)
+        nuts_dict = nuts_gdf[["NUTS_ID", "geometry", "bounds"]].set_index("NUTS_ID")\
+                                                               .to_dict("index")
+        for r in nuts_dict:
+            nuts_dict[r] = tuple(nuts_dict[r].values())
+        return nuts_dict
+
+    def reverse_geocode_nuts(self,
+                             latlons: List[Tuple[float, float]],
+                             level: Literal[0, 1, 2, 3],
+                             year: Literal[2003, 2006, 2010, 2013, 2016, 2021] = 2021
+                            ) -> List[str]:
+        """
+        Reverse-geocode latitudes and longitudes to NUTS regions.
+
+        Parameters
+        ----------
+        `latlons` : list of tuples
+            A list of tuples containing (latitude, longitude).
+        `level` : int
+            Specify the NUTS level, must be one of [0,1,2,3].
+        `year` : int
+            Specify the year of NUTS regulation, must be one of [2003,2006,2010,2013,2016,2021],
+            defaults to 2021.
+
+        Returns
+        -------
+        list of strings
+            The NUTS_ID codes that the input latitudes and longitudes fall within. Any lat/lons which
+            do not fall inside a NUTS boundary will return None.
+        """
+        if self.nuts_regions[(level, year)] is None:
+            self.nuts_regions[(level, year)] = self._load_nuts_boundaries(level=level, year=year)
+        results = utils.reverse_geocode(latlons, self.nuts_regions[(level, year)])
+        return results
+
diff --git a/geocode/geocode.py b/geocode/geocode.py
index d23b942..a4a288a 100644
--- a/geocode/geocode.py
+++ b/geocode/geocode.py
@@ -11,42 +11,33 @@
 import os
 import sys
 import logging
+from pathlib import Path
 import time as TIME
 import argparse
-from typing import Optional
 from shutil import copyfile
+from typing import Optional, Iterable, Tuple, Union, List, Dict, Literal
 
-from typing import Optional, Iterable, Tuple, Union, List, Dict
 import pyproj
 
-try:
-    from shapely.geometry import shape, Point
-    from shapely.ops import unary_union
-    SHAPELY_AVAILABLE = True
-except ImportError:
-    logging.warning("Failed to import Shapely library - you will not be able to reverse-geocode! "
-                    "See notes in the README about installing Shapely on Windows machines.")
-    SHAPELY_AVAILABLE = False
-
-SCRIPT_DIR = os.path.dirname(os.path.realpath(__file__))
-# sys.path.append(SCRIPT_DIR)
-
 from . utilities import GenericException
 from . cpo import CodePointOpen
 from . ngeso import NationalGrid
 from . ons_nrs import ONS_NRS
+from . eurostat import Eurostat
 from . gmaps import GMaps
 from . cache_manager import CacheManager
 from . version import __version__
 
+SCRIPT_DIR = Path(os.path.dirname(os.path.realpath(__file__)))
+
 class Geocoder:
     """
     Geocode addresses, postcodes, LLSOAs or Constituencies or reverse-geocode latitudes and
     longitudes.
     """
     def __init__(self,
-                 cache_dir: Optional[str] = None,
-                 gmaps_key_file: Optional[str] = None) -> None:
+                 cache_dir: Optional[Path] = None,
+                 gmaps_key_file: Optional[Path] = None) -> None:
         """
         Geocode addresses, postcodes, LLSOAs or Constituencies or reverse-geocode latitudes and
         longitudes.
@@ -63,6 +54,7 @@ def __init__(self,
         self.cpo = CodePointOpen(self.cache_manager)
         self.ngeso = NationalGrid(self.cache_manager)
         self.ons_nrs = ONS_NRS(self.cache_manager)
+        self.eurostat = Eurostat(self.cache_manager)
         self.gmaps = GMaps(self.cache_manager, gmaps_key_file)
         self.status_codes = {
             0: "Failed",
@@ -78,9 +70,9 @@ def __enter__(self):
 
     def __exit__(self, type, value, traceback):
         """Context manager."""
-        pass
+        self.gmaps.flush_cache()
 
-    def force_setup(self, ngeso_setup=True, cpo_setup=True, ons_setup=True):
+    def force_setup(self, ngeso_setup=True, cpo_setup=True, ons_setup=True, eurostat_setup=True):
         """Download all data and setup caches."""
         if ngeso_setup:
             self.ngeso.force_setup()
@@ -88,6 +80,8 @@ def force_setup(self, ngeso_setup=True, cpo_setup=True, ons_setup=True):
             self.cpo.force_setup()
         if ons_setup:
             self.ons_nrs.force_setup()
+        if eurostat_setup:
+            self.eurostat.force_setup()
 
     def get_dno_regions(self):
         """
@@ -137,12 +131,30 @@ def reverse_geocode_llsoa(self, latlons, dz=True):
         ----------
         `latlons` : iterable of strings
             Specific latlons to geocode to llsoa boundaries.
-        
         `dz` : Boolean
             Indication whether to consider datazones
         """
         return self.reverse_geocode(latlons, "llsoa", datazones=dz)
 
+    def reverse_geocode_nuts(self, latlons: List[Tuple[float, float]],
+                             level: Literal[0, 1, 2, 3],
+                             year: Literal[2003, 2006, 2010, 2013, 2016, 2021] = 2021
+                            ) -> List[str]:
+        """
+        Function to reverse geocode a collection of latlons into NUTS boundaries.
+        
+        Parameters
+        ----------
+        `latlons` : iterable of strings
+            Specific latlons to geocode to llsoa boundaries.
+        `level` : int
+            Specify the NUTS level, must be one of [0,1,2,3].
+        `year` : int
+            Specify the year of NUTS regulation, must be one of [2003,2006,2010,2013,2016,2021],
+            defaults to 2021.
+        """
+        return self.reverse_geocode(latlons, "nuts", level=level, year=year)
+
     def geocode_constituency(self, constituencies):
         """
         Function to geocode a collection of constituencies into latlons.
@@ -242,6 +254,10 @@ def reverse_geocode(self, latlons, entity, **kwargs):
         elif entity == "llsoa":
             datazones = kwargs.get("datazones", False)
             return self.ons_nrs.reverse_geocode_llsoa(latlons=latlons, datazones=datazones)
+        elif entity == "nuts":
+            level = kwargs.get("level")
+            year = kwargs.get("year", 2021)
+            return self.eurostat.reverse_geocode_nuts(latlons=latlons, level=level, year=year)
         else:
             raise GenericException(f"Entity '{entity}' is not supported.")
 
@@ -313,7 +329,7 @@ def parse_options():
                         required=False, help="Force download all datasets to local cache (useful "
                                              "if running inside a Docker container i.e. run this "
                                              "as part of image build). Possible values are "
-                                             "'ngeso', 'cpo', 'ons' or 'all'.")
+                                             "'ngeso', 'cpo', 'ons', 'eurostat' or 'all'.")
     parser.add_argument("--load-cpo-zip", dest="cpo_zip", action="store", type=str,
                         required=False, default=None, metavar="</path/to/zip-file>",
                         help="Load the Code Point Open data from a local zip file.")
@@ -323,7 +339,7 @@ def parse_options():
     options = parser.parse_args()
     def handle_options(options):
         if options.setup is not None:
-            valid_options = ["ngeso", "cpo", "ons", "all"]
+            valid_options = ["ngeso", "cpo", "ons", "eurostat", "all"]
             options.setup = list(map(str.lower, options.setup))
             if any(s not in valid_options for s in options.setup):
                 raise ValueError(f"Invalid value for `--setup` - valid values are {valid_options}")
@@ -353,7 +369,7 @@ def debug():
     logging.info("Time taken: %s seconds", round(TIME.time() - timerstart, 1))
     for (lat, lon), llsoa in zip(sample_latlons, results):
         logging.info("%s, %s :    %s", lat, lon, llsoa)
-    sample_file = os.path.join(SCRIPT_DIR, "sample_latlons.txt")
+    sample_file = SCRIPT_DIR.joinpath("sample_latlons.txt")
     with open(sample_file) as fid:
         sample_latlons = [tuple(map(float, line.strip().split(",")))
                           for line in fid if line.strip()][:10]
@@ -374,7 +390,7 @@ def debug():
     logging.info("Time taken: %s seconds", round(TIME.time() - timerstart, 1))
     for constituency, (lat, lon) in zip(sample_constituencies, results):
         logging.info("%s :    %s, %s", constituency, lat, lon)
-    sample_file = os.path.join(SCRIPT_DIR, "sample_postcodes.txt")
+    sample_file = SCRIPT_DIR.joinpath("sample_postcodes.txt")
     with open(sample_file) as fid:
         postcodes = [line.strip() for line in fid if line.strip()][:10]
     timerstart = TIME.time()
@@ -409,9 +425,11 @@ def main():
         ngeso_setup =  "ngeso" in options.setup or "all" in options.setup
         cpo_setup = "cpo" in options.setup or "all" in options.setup
         ons_setup = "ons" in options.setup or "all" in options.setup
+        eurostat_setup = "eurostat" in options.setup or "all" in options.setup
         logging.info("Running forced setup")
         with Geocoder() as geocoder:
-            geocoder.force_setup(ngeso_setup=ngeso_setup, cpo_setup=cpo_setup, ons_setup=ons_setup)
+            geocoder.force_setup(ngeso_setup=ngeso_setup, cpo_setup=cpo_setup, ons_setup=ons_setup,
+                                 eurostat_setup=eurostat_setup)
     if options.debug:
         debug()
 
diff --git a/geocode/gmaps.py b/geocode/gmaps.py
index 3962395..a708b22 100644
--- a/geocode/gmaps.py
+++ b/geocode/gmaps.py
@@ -10,13 +10,16 @@
 import sys
 import logging
 import pickle
+from pathlib import Path
+from typing import Optional, Iterable, Tuple, Union, List, Dict
+
 import numpy as np
 import pandas as pd
-from typing import Optional, Iterable, Tuple, Union, List, Dict
+import googlemaps
 
 from . utilities import GenericException
 
-SCRIPT_DIR = os.path.dirname(os.path.realpath(__file__))
+SCRIPT_DIR = Path(os.path.dirname(os.path.realpath(__file__)))
 
 class GMaps:
     """The Gmaps data manager for the Geocode class."""
@@ -27,16 +30,17 @@ def __init__(self, cache_manager, gmaps_key_file=None):
         self.gmaps_client = None
         self.cache_file = "gmaps_cache"
         self.gmaps_key_file = gmaps_key_file if gmaps_key_file is not None \
-                                  else os.path.join(self.cache_manager.cache_dir, "key.txt")
+                                  else self.cache_manager.cache_dir.joinpath("key.txt")
         self._load_cache()
+        self.cache_modified = False
 
     def __enter__(self):
         """Context manager."""
         return self
 
-    def __exit__(self, *args):
+    def __exit__(self, type, value, traceback):
         """Context manager - flush GMaps cache on exit."""
-        self.flush_cache(self.cache, self.cache_file)
+        self.flush_cache()
 
     def _load_key(self):
         """Load the user's GMaps API key from installation directory."""
@@ -50,9 +54,10 @@ def _load_key(self):
             return None
         return key
 
-    def flush_cache(self, gmaps_cache_contents):
+    def flush_cache(self):
         """Flush any new GMaps API queries to the GMaps cache."""
-        self.cache_manager.write(self.cache_file, self.cache)
+        if self.cache_modified:
+            self.cache_manager.write(self.cache_file, self.cache)
 
     def geocode_postcode(self, postcode: [str],
                          address: Optional[str] = None) -> Union[Tuple[float, float], List[Tuple[float, float]]]:
@@ -112,14 +117,16 @@ class attribute *self.status_codes* (a dict) to get a string representation.
         postcode = postcode if postcode is not None else ""
         address = address if address is not None else ""
         search_term = f"{address}{sep}{postcode}"
-        logging.debug("Querying Google Maps Geocoder API for '%s'", search_term)
         if search_term in self.cache:
+            logging.debug("Loading GMaps Geocoder API result from cache: '%s'", search_term)
             geocode_result = self.cache[search_term]
         else:
+            logging.debug("Querying Google Maps Geocoder API for '%s'", search_term)
             if self.gmaps_key is None:
                 return pd.Series({"latitude": np.nan, "longitude": np.nan, "match_status": 0})
             geocode_result = self.gmaps_client.geocode(search_term, region="uk")
-            self.gmaps_cache[search_term] = geocode_result
+            self.cache[search_term] = geocode_result
+            self.cache_modified = True
         if not geocode_result or len(geocode_result) > 1:
             return pd.Series({"latitude": np.nan, "longitude": np.nan, "match_status": 0})
         geometry = geocode_result[0]["geometry"]
@@ -135,7 +142,6 @@ def _load_cache(self):
         """Load the cache of prior addresses/postcodes for better performance."""
         self.cache = self.cache_manager.retrieve(self.cache_file)
         if self.cache is None:
-            self.cache = pd.DataFrame({"input_postcode": [], "input_address": [], "latitude": [],
-                                       "longitude": [], "match_status": []})
+            self.cache = {}
         return
 
diff --git a/geocode/latlons2gsp.py b/geocode/latlons2gsp.py
index ba54eed..d6fb789 100644
--- a/geocode/latlons2gsp.py
+++ b/geocode/latlons2gsp.py
@@ -8,9 +8,10 @@
 """
 
 import sys
-import os
 import argparse
 import time as TIME
+from pathlib import Path
+
 import pandas as pd
 
 from geocode import Geocoder
@@ -21,17 +22,17 @@ def parse_options():
     parser = argparse.ArgumentParser(description=("This is a command line interface (CLI) for "
                                                   "the latlons2gsp.py module"),
                                      epilog="Jamie Taylor & Ethan Jones, 2020-05-29")
-    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=str,
+    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=Path,
                         required=True, metavar="</path/to/file>",
                         help="Specify a CSV file containing a list of latitudes and longitudes to "
                              "be reverse-geocoded. The file must contain the columns 'latitude' "
                              "and 'longitude' (it can contain others, all of which will be kept).")
-    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=str,
+    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=Path,
                         required=True, metavar="</path/to/file>", help="Specify an output file.")
     options = parser.parse_args()
-    if not os.path.isfile(options.infile):
+    if not options.infile.is_file():
         raise GenericException(f"The input file '{options.infile}' does not exist.")
-    if os.path.isfile(options.outfile):
+    if options.outfile.is_file():
         check = query_yes_no(f"The outfile '{options.outfile}' already exists, are you sure you "
                              "wish to overwrite?", "no")
         if not check:
diff --git a/geocode/latlons2llsoa.py b/geocode/latlons2llsoa.py
index b70b50a..54be08f 100644
--- a/geocode/latlons2llsoa.py
+++ b/geocode/latlons2llsoa.py
@@ -8,9 +8,10 @@
 """
 
 import sys
-import os
 import argparse
 import time as TIME
+from pathlib import Path
+
 import pandas as pd
 
 from geocode import Geocoder
@@ -21,19 +22,19 @@ def parse_options():
     parser = argparse.ArgumentParser(description=("This is a command line interface (CLI) for "
                                                   "the latlons2llsoa.py module"),
                                      epilog="Jamie Taylor & Ethan Jones, 2020-04-16")
-    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=str,
+    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=Path,
                         required=True, metavar="</path/to/file>",
                         help="Specify a CSV file containing a list of latitudes and longitudes to "
                              "be reverse-geocoded. The file must contain the columns 'latitude' "
                              "and 'longitude' (it can contain others, all of which will be kept).")
-    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=str,
+    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=Path,
                         required=True, metavar="</path/to/file>", help="Specify an output file.")
     parser.add_argument("--datazones", dest="datazones", action="store_true",
                         required=False, help="Specify to use Data Zones in Scotland.")
     options = parser.parse_args()
-    if not os.path.isfile(options.infile):
+    if not options.infile.is_file():
         raise GenericException(f"The input file '{options.infile}' does not exist.")
-    if os.path.isfile(options.outfile):
+    if options.outfile.is_file():
         check = query_yes_no(f"The outfile '{options.outfile}' already exists, are you sure you "
                              "wish to overwrite?", "no")
         if not check:
diff --git a/geocode/ngeso.py b/geocode/ngeso.py
index ce38491..bf55ed5 100644
--- a/geocode/ngeso.py
+++ b/geocode/ngeso.py
@@ -11,22 +11,21 @@
 import pickle
 import logging
 import json
+from pathlib import Path
+from typing import Optional, Iterable, Tuple, Union, List, Dict
+
 import pandas as pd
 import geopandas as gpd
-
-from typing import Optional, Iterable, Tuple, Union, List, Dict
 try:
     from shapely.geometry import shape, Point
     from shapely.ops import unary_union
-    SHAPELY_AVAILABLE = True
 except ImportError:
     logging.warning("Failed to import Shapely library - you will not be able to reverse-geocode! "
                     "See notes in the README about installing Shapely on Windows machines.")
-    SHAPELY_AVAILABLE = False
 
 from . import utilities as utils
 
-SCRIPT_DIR = os.path.dirname(os.path.realpath(__file__))
+SCRIPT_DIR = Path(os.path.dirname(os.path.realpath(__file__)))
 
 class NationalGrid:
     """The NGESO data manager for the Geocode class."""
@@ -156,10 +155,6 @@ def reverse_geocode_gsp(self,
         -----
         Return format needs some work, maybe switch to DataFrames in future release.
         """
-        #if not SHAPELY_AVAILABLE:
-        #    raise utilities.GenericException("Geocode was unable to import the Shapely library, follow the "
-        #                           "installation instructions at "
-        #                           "https://github.com/SheffieldSolar/Geocode")
         logging.debug("Reverse geocoding %s latlons to 20220314 GSP", len(latlons))
         if self.gsp_regions is None:
             self.gsp_regions = self._load_gsp_boundaries_20220314()
@@ -198,10 +193,6 @@ def reverse_geocode_gsp_20181031(self,
         -----
         Return format needs some work, maybe switch to DataFrames in future release.
         """
-        if not SHAPELY_AVAILABLE:
-            raise utils.GenericException("Geocode was unable to import the Shapely library, follow the "
-                                   "installation instructions at "
-                                   "https://github.com/SheffieldSolar/Geocode")
         if self.gsp_regions_20181031 is None:
             self.gsp_regions_20181031 = self.load_gsp_boundaries_20181031()
         lats = [l[0] for l in latlons]
diff --git a/geocode/ons_nrs.py b/geocode/ons_nrs.py
index 52e71bf..af32f39 100644
--- a/geocode/ons_nrs.py
+++ b/geocode/ons_nrs.py
@@ -12,22 +12,21 @@
 import json
 import csv
 import logging
+from pathlib import Path
+from typing import Optional, Iterable, Tuple, Union, List, Dict
+
 import pandas as pd
 import shapefile
-
-from typing import Optional, Iterable, Tuple, Union, List, Dict
 try:
     from shapely.geometry import shape, Point
     from shapely.ops import unary_union
-    SHAPELY_AVAILABLE = True
 except ImportError:
     logging.warning("Failed to import Shapely library - you will not be able to reverse-geocode! "
                     "See notes in the README about installing Shapely on Windows machines.")
-    SHAPELY_AVAILABLE = False
 
 from . import utilities as utils
 
-SCRIPT_DIR = os.path.dirname(os.path.realpath(__file__))
+SCRIPT_DIR = Path(os.path.dirname(os.path.realpath(__file__)))
 
 class ONS_NRS:
     """
@@ -39,15 +38,11 @@ def __init__(self, cache_manager):
         (NRS).
         """
         self.cache_manager = cache_manager
-        data_dir = os.path.join(SCRIPT_DIR, "ons")
-        self.nrs_zipfile = os.path.join(data_dir, "nrs.zip")
-        self.dz_lookup_cache_file = "datazone_lookup"
-        self.constituency_lookup_file = os.path.join(data_dir, "constituency_centroids_Dec2020.psv")
-        self.constituency_cache_file = os.path.join(data_dir, "constituency_centroids.psv")
-        self.lad_lookup_file = os.path.join(data_dir, "lad_centroids_May2021.psv")
-        self.lad_cache_file = "lad_centroids"
-        self.pc_llsoa_lookup_cache_file = "pc_llsoa_lookup"
-        self.pc_llsoa_zipfile = os.path.join(data_dir, "PCD_OA_LSOA_MSOA_LAD_MAY22_UK_LU.zip")
+        data_dir = SCRIPT_DIR.joinpath("ons")
+        self.nrs_zipfile = data_dir.joinpath("nrs.zip")
+        self.constituency_lookup_file = data_dir.joinpath("constituency_centroids_Dec2020.psv")
+        self.lad_lookup_file = data_dir.joinpath("lad_centroids_May2021.psv")
+        self.pc_llsoa_zipfile = data_dir.joinpath("PCD_OA_LSOA_MSOA_LAD_MAY22_UK_LU.zip")
         self.llsoa_lookup = None
         self.llsoa_regions = None
         self.llsoa_reverse_lookup = None
@@ -142,10 +137,9 @@ def _load_llsoa_boundaries(self):
 
     def _load_datazone_lookup(self):
         """Load a lookup of Scottish LLSOA <-> Datazone."""
-        datazone_lookup_cache_contents = self.cache_manager.retrieve(self.dz_lookup_cache_file)
+        datazone_lookup_cache_contents = self.cache_manager.retrieve("datazone_lookup")
         if datazone_lookup_cache_contents is not None:
-            logging.debug("Loading LLSOA<->Datazone lookup from cache ('%s')",
-                          self.dz_lookup_cache_file)
+            logging.debug("Loading LLSOA<->Datazone lookup from cache ('%s')", "datazone_lookup")
             return datazone_lookup_cache_contents
         with zipfile.ZipFile(self.nrs_zipfile, "r") as nrs_zip:
             with nrs_zip.open("OA_DZ_IZ_2011.csv", "r") as fid:
@@ -153,15 +147,14 @@ def _load_datazone_lookup(self):
         dz_lookup.set_index("OutputArea2011Code", inplace=True)
         dz_lookup.drop(columns=["IntermediateZone2011Code"], inplace=True)
         dz_lookup = dz_lookup.to_dict()["DataZone2011Code"]
-        self.cache_manager.write(self.dz_lookup_cache_file, dz_lookup)
+        self.cache_manager.write("datazone_lookup", dz_lookup)
         return dz_lookup
 
     def _load_constituency_lookup(self):
         """Load a lookup of UK constituency -> Geospatial Centroid."""
-        constituency_lookup_cache_contents = self.cache_manager.retrieve(self.constituency_cache_file)
+        constituency_lookup_cache_contents = self.cache_manager.retrieve("constituency_centroids")
         if constituency_lookup_cache_contents is not None:
-            logging.debug("Loading Constituency lookup from cache ('%s')",
-                      self.constituency_cache_file)
+            logging.debug("Loading Constituency lookup from cache ('%s')", "constituency_centroids")
             return constituency_lookup_cache_contents
         logging.info("Extracting the Constituency Centroids data (this only needs to be done "
                      "once)")
@@ -171,9 +164,8 @@ def _load_constituency_lookup(self):
                 _, name, longitude, latitude = line.strip().split("|")
                 match_str = name.strip().replace(" ", "").replace(",", "").lower()
                 constituency_lookup[match_str] = (float(latitude), float(longitude))
-        self.cache_manager.write(self.constituency_cache_file, constituency_lookup)
-        logging.info("Constituency lookup extracted and pickled to '%s'",
-                     self.constituency_cache_file)
+        self.cache_manager.write("constituency_centroids", constituency_lookup)
+        logging.info("Constituency lookup extracted and pickled to '%s'", "constituency_centroids")
         return constituency_lookup
 
     def geocode_llsoa(self,
@@ -231,10 +223,6 @@ def reverse_geocode_llsoa(self,
             The LLSOA codes that the input latitudes and longitudes fall within. Any lat/lons which
             do not fall inside an LLSOA boundary will return None.
         """
-        if not SHAPELY_AVAILABLE:
-            raise utils.GenericException("Geocode was unable to import the Shapely library, follow the "
-                                   "installation instructions at "
-                                   "https://github.com/SheffieldSolar/Geocode")
         if self.llsoa_regions is None:
             self.llsoa_regions = self._load_llsoa_boundaries()
         results = utils.reverse_geocode(latlons, self.llsoa_regions)
@@ -335,10 +323,10 @@ def postcode2llsoa(self, pcs: pd.DataFrame) -> pd.DataFrame:
 
     def _load_lad_lookup(self):
         """Load a lookup of UK Local Authority District -> Geospatial Centroid."""
-        lad_lookup_cache_contents = self.cache_manager.retrieve(self.lad_cache_file)
+        lad_lookup_cache_contents = self.cache_manager.retrieve("lad_centroids")
         if lad_lookup_cache_contents is not None:
             logging.debug("Loading Local Authority District lookup from cache ('%s')",
-                          self.lad_cache_file)
+                          "lad_centroids")
             return lad_lookup_cache_contents
         logging.info("Extracting the Local Authority District Centroids data (this only needs to "
                      "be done once)")
@@ -348,22 +336,21 @@ def _load_lad_lookup(self):
                 _, name, longitude, latitude = line.strip().split("|")
                 match_str = name.strip().replace(" ", "").replace(",", "").lower()
                 lad_lookup[match_str] = (float(latitude), float(longitude))
-        self.cache_manager.write(self.lad_cache_file, lad_lookup)
+        self.cache_manager.write("lad_centroids", lad_lookup)
         logging.info("Local Authority District lookup extracted and pickled to '%s'",
-                     self.lad_cache_file)
+                     "lad_centroids")
         return lad_lookup
 
     def _load_postcode_llsoa_lookup(self):
         """Load a lookup of postcode <-> LLSOA."""
-        postcode_llsoa_lookup_cache_contents = self.cache_manager.retrieve(self.pc_llsoa_lookup_cache_file)
+        postcode_llsoa_lookup_cache_contents = self.cache_manager.retrieve("pc_llsoa_lookup")
         if postcode_llsoa_lookup_cache_contents is not None:
-            logging.debug("Loading postcode<->LLSOA lookup from cache ('%s')",
-                          self.pc_llsoa_lookup_cache_file)
+            logging.debug("Loading postcode<->LLSOA lookup from cache ('%s')", "pc_llsoa_lookup")
             return postcode_llsoa_lookup_cache_contents
         pc_llsoa_lookup = pd.read_csv(self.pc_llsoa_zipfile, dtype=str)
         pc_llsoa_lookup["postcode"] = pc_llsoa_lookup.pcds.str.strip().str.upper()\
                                                                       .str.replace(" ", "")
-        self.cache_manager.write(self.pc_llsoa_lookup_cache_file, pc_llsoa_lookup)
+        self.cache_manager.write("pc_llsoa_lookup", pc_llsoa_lookup)
         return pc_llsoa_lookup
 
     def _lad_centroid(self, local_authority):
diff --git a/geocode/postcodes2latlon.py b/geocode/postcodes2latlon.py
index b8f8f69..c311794 100644
--- a/geocode/postcodes2latlon.py
+++ b/geocode/postcodes2latlon.py
@@ -8,9 +8,10 @@
 """
 
 import sys
-import os
 import argparse
 import time as TIME
+from pathlib import Path
+
 import pandas as pd
 
 from geocode import Geocoder
@@ -21,17 +22,17 @@ def parse_options():
     parser = argparse.ArgumentParser(description=("This is a command line interface (CLI) for "
                                                   "the postcodes2latlon.py module"),
                                      epilog="Jamie Taylor & Ethan Jones, 2020-06-12")
-    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=str,
+    parser.add_argument("-f", "--input-file", dest="infile", action="store", type=Path,
                         required=True, metavar="</path/to/file>",
                         help="Specify a CSV file containing a list of postcodes to be geocoded. "
                              "The file must contain the column 'postcode' (it can contain others, "
                              "all of which will be kept).")
-    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=str,
+    parser.add_argument("-o", "--output-file", dest="outfile", action="store", type=Path,
                         required=True, metavar="</path/to/file>", help="Specify an output file.")
     options = parser.parse_args()
-    if not os.path.isfile(options.infile):
+    if not options.infile.is_file():
         raise GenericException(f"The input file '{options.infile}' does not exist.")
-    if os.path.isfile(options.outfile):
+    if options.outfile.is_file():
         check = query_yes_no(f"The outfile '{options.outfile}' already exists, are you sure you "
                              "wish to overwrite?", "no")
         if not check:
diff --git a/geocode/utilities.py b/geocode/utilities.py
index 56dfa81..47f2bd5 100644
--- a/geocode/utilities.py
+++ b/geocode/utilities.py
@@ -10,10 +10,9 @@
 import logging
 import requests
 import json
-
 from typing import Optional, Iterable, Tuple, Union, List, Dict
-import pyproj
 
+import pyproj
 try:
     from shapely.geometry import shape, Point
     from shapely.ops import unary_union
@@ -238,8 +237,8 @@ def reverse_geocode(coords: List[Tuple[float, float]],
     `Shapely.MultiPolygon.contains()` is not great.
     """
     if not SHAPELY_AVAILABLE:
-        raise utilities.GenericException("Geocode was unable to import the Shapely library, follow the "
-                                         "installation instructions at "
+        raise utilities.GenericException("Geocode was unable to import the Shapely library, follow "
+                                         "the installation instructions at "
                                          "https://github.com/SheffieldSolar/Geocode")
     results = []
     tot = len(coords)
diff --git a/geocode/version.py b/geocode/version.py
index cb73775..7ff9620 100644
--- a/geocode/version.py
+++ b/geocode/version.py
@@ -1 +1 @@
-__version__ = "0.11.1"
\ No newline at end of file
+__version__ = "0.12.0"
\ No newline at end of file
diff --git a/setup.py b/setup.py
index 05e2d4c..1c8abb1 100644
--- a/setup.py
+++ b/setup.py
@@ -24,7 +24,7 @@
     # Versions should comply with PEP440.  For a discussion on single-sourcing
     # the version across setup.py and the project code, see
     # https://packaging.python.org/en/latest/single_source_version.html
-    version="0.11.1",
+    version="0.12.0",
 
     description="Geocode postcodes, addresses or LLSOA using the Code Point Open database and GMaps API.",
     long_description=long_description,
@@ -32,7 +32,7 @@
 
     # The project's main homepage.
     url="https://github.com/SheffieldSolar/Geocode",
-    download_url="https://github.com/SheffieldSolar/Geocode/archive/refs/tags/0.11.1.tar.gz",
+    download_url="https://github.com/SheffieldSolar/Geocode/archive/refs/tags/0.12.0.tar.gz",
 
     # Author details
     author="Jamie Taylor",
