diff --git a/documentation/images/examples_ozon.png b/documentation/images/examples_ozon.png
index bf6bd1bfd6373f34cd1b6ddc0e790fab61dadcd6..c6fccdfe502a7e0e13776c807889ce01f4295187 100644
Binary files a/documentation/images/examples_ozon.png and b/documentation/images/examples_ozon.png differ
diff --git a/examples/00_download_and_visualization.ipynb b/examples/00_download_and_visualization.ipynb
index d1b32aae0d0c5eba94b76a6f9c68965d65b3c5c1..d868cff2e39e402bcad158300a393e7a7582ca60 100644
--- a/examples/00_download_and_visualization.ipynb
+++ b/examples/00_download_and_visualization.ipynb
@@ -13,7 +13,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -49,7 +49,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -73,13 +73,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
     "\n",
     "time = TimeSample(dt(2016,3,1), dt(2016,3,8), \"daily\")\n",
-    "metadata = Metadata.construct(\"mole_fraction_of_ozone_in_air\", time, \"mean\")\n",
+    "metadata = Metadata.construct(\"mole_fraction_of_ozone_in_air\", time, \"mean\", True)\n",
     "my_grid = RegularGrid(2.0, 2.5)"
    ]
   },
@@ -98,7 +98,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -164,7 +164,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -188,7 +188,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -282,7 +282,7 @@
     "time = ds.time[timestep]\n",
     "data = ds.sel(time=time)\n",
     "\n",
-    "plot_cells(data[\"mean\"], not_na_coords, all_na_coords, discrete=False, plot_stations=True)\n",
+    "plot_cells(data[\"mean\"], not_na_coords, all_na_coords, discrete=False, plot_stations=False)\n",
     "plt.show()\n",
     "\n",
     "plot_cells(data[\"n\"], not_na_coords, all_na_coords, discrete=True)\n",
@@ -290,9 +290,9 @@
     "\n",
     "n_observations = ds[\"n\"].sum([\"latitude\", \"longitude\"])\n",
     "plt.plot(ds.time, n_observations)\n",
-    "plt.title(\"Number of Stations per Time Point\")\n",
+    "plt.title(\"Number of Contributions per Time Point\")\n",
     "plt.xlabel(\"Time\")\n",
-    "plt.ylabel(\"Number of Stations\")\n",
+    "plt.ylabel(\"Number of contributions\")\n",
     "print(np.unique(ds[\"n\"]))"
    ]
   },
@@ -315,7 +315,7 @@
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "interactive",
+   "display_name": ".venv",
    "language": "python",
    "name": "python3"
   },
@@ -329,7 +329,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.7"
+   "version": "3.11.9"
   }
  },
  "nbformat": 4,
diff --git a/examples/01_produce_data_one_week.ipynb b/examples/01_produce_data_one_week.ipynb
index db0c98679594361096750df567f9f822997a559a..ff1e9bbccf5254e8b3e02a72ac2c3ddb26685750 100644
--- a/examples/01_produce_data_one_week.ipynb
+++ b/examples/01_produce_data_one_week.ipynb
@@ -121,7 +121,10 @@
     "\n",
     "The obtained data are stored in the cache directory. Before submitting a request, toargridding checks his cache, if the data have already been downloaded.\n",
     "\n",
-    "This function also creates a file with the extension \".contributors\"."
+    "This function also creates a file with the extension \".contributors\".\n",
+    "\n",
+    "\n",
+    "We also calculate first the average of all timeseries at each station before we average all stations per grid cell (option `averaage_TS_at_station`). Otherwise the results of all timeseries in each grid cell are averaged."
    ]
   },
   {
@@ -136,7 +139,8 @@
     "    time=time_sampling,\n",
     "    variables=variable,\n",
     "    stats=statistics,\n",
-    "    contributors_path=result_basepath\n",
+    "    contributors_path=result_basepath,\n",
+    "    average_TS_at_station=True\n",
     ")"
    ]
   },
@@ -180,7 +184,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.11.9"
+   "version": "3.9.21"
   }
  },
  "nbformat": 4,
diff --git a/examples/02_produce_data_manyStations.ipynb b/examples/02_produce_data_manyStations.ipynb
index 50638aa3fb5bc83a73cc38946dc093af9018eac0..74c42fee464ab9e3d0bfb689a54cf3a6109c399d 100644
--- a/examples/02_produce_data_manyStations.ipynb
+++ b/examples/02_produce_data_manyStations.ipynb
@@ -182,7 +182,8 @@
     "        time=config.time,\n",
     "        variables=config.variables,\n",
     "        stats=config.stats,\n",
-    "        contributors_path=result_basepath\n",
+    "        contributors_path=result_basepath,\n",
+    "        average_TS_at_station=True\n",
     "    )\n",
     "\n",
     "    for dataset, metadata in zip(datasets, metadatas):\n",