add bias

Author: malmans2

notebooks/wp5/steric_sea_level_ora5.ipynb

@@ -527,15 +527,28 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "for label, ds in zip([\"ORAS5\", \"ARGO\"], [ds_reanalysis_trend, ds_argo_trend]):\n",
-    "    da = ds.to_dataarray()\n",
-    "    facet = plot.projected_map(da, col=\"variable\", robust=True)\n",
-    "    for ax in facet.axs.flatten():\n",
-    "        ax.set_extent(\n",
-    "            [lon_slice.start, lon_slice.stop, lat_slice.start, lat_slice.stop]\n",
-    "        )\n",
-    "    facet.fig.suptitle(label)\n",
-    "    plt.show()"
+    "ds_reanalysis_trend_interp = diagnostics.regrid(\n",
+    "    ds_reanalysis_trend, ds_argo_trend, method=\"nearest_s2d\"\n",
+    ")\n",
+    "ds_trend = xr.concat(\n",
+    "    [\n",
+    "        ds_reanalysis_trend_interp.expand_dims(product=[\"ORAS5\"]),\n",
+    "        ds_argo_trend.expand_dims(product=[\"ARGO\"]),\n",
+    "    ],\n",
+    "    \"product\",\n",
+    ")\n",
+    "da = ds_trend.to_dataarray()\n",
+    "seconds_per_year = 365 * 24 * 360\n",
+    "with xr.set_options(keep_attrs=True):\n",
+    "    meters_to_millimeters = 1000\n",
+    "    da = da * meters_to_millimeters * seconds_per_year\n",
+    "    da.attrs[\"units\"] = \"mm/yr\"\n",
+    "    bias = da.diff(\"product\")\n",
+    "plot.projected_map(da, col=\"variable\", row=\"product\", robust=True)\n",
+    "\n",
+    "\n",
+    "plot.projected_map(bias, col=\"variable\", row=\"product\", robust=True)\n",
+    "_ = plt.suptitle(f\"Bias {'-'.join(da['product'].values.tolist())}\")"
    ]
   },
   {