Skip to content
Snippets Groups Projects
Commit 25b237cf authored by Anna Lanteri's avatar Anna Lanteri
Browse files

hope this fixes

parent f57a232b
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
scenario comparison/scenarios_vs_historical.ipynb
+ 2
− 95
View file @ 25b237cf
......@@ -939,53 +939,23 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 48,
=======
"execution_count": 19,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "84377e1d-1620-42a7-aa7f-a67ad35011f1",
"metadata": {
"tags": []
},
"outputs": [
{
"ename": "ValueError",
"evalue": "Length mismatch: Expected axis has 1 elements, new values have 2 elements",
"output_type": "error",
"traceback": [
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
"\u001b[0;31mValueError\u001b[0m Traceback (most recent call last)",
"Cell \u001b[0;32mIn[19], line 3\u001b[0m\n\u001b[1;32m 1\u001b[0m cams_nox_file \u001b[38;5;241m=\u001b[39m \u001b[38;5;124m'\u001b[39m\u001b[38;5;124m/home/b/b309253/data-exploration/emissions/cams-glob-ant-anthro-nox.csv\u001b[39m\u001b[38;5;124m'\u001b[39m\n\u001b[1;32m 2\u001b[0m cams_nox_data \u001b[38;5;241m=\u001b[39m pd\u001b[38;5;241m.\u001b[39mread_csv(cams_nox_file, sep\u001b[38;5;241m=\u001b[39m\u001b[38;5;124m'\u001b[39m\u001b[38;5;130;01m\\t\u001b[39;00m\u001b[38;5;124m'\u001b[39m)\n\u001b[0;32m----> 3\u001b[0m \u001b[43mcams_nox_data\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mcolumns\u001b[49m \u001b[38;5;241m=\u001b[39m [\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mDate\u001b[39m\u001b[38;5;124m'\u001b[39m, \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mCAMS Road\u001b[39m\u001b[38;5;124m'\u001b[39m]\n\u001b[1;32m 4\u001b[0m cams_nox_data \u001b[38;5;241m=\u001b[39m cams_nox_data\u001b[38;5;241m.\u001b[39mset_index(pd\u001b[38;5;241m.\u001b[39mdate_range(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m1/1/2000\u001b[39m\u001b[38;5;124m\"\u001b[39m, periods\u001b[38;5;241m=\u001b[39m\u001b[38;5;241m24\u001b[39m, freq\u001b[38;5;241m=\u001b[39m\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mY\u001b[39m\u001b[38;5;124m'\u001b[39m))\u001b[38;5;241m.\u001b[39mdrop(columns \u001b[38;5;241m=\u001b[39m [\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mDate\u001b[39m\u001b[38;5;124m'\u001b[39m])\n",
"File \u001b[0;32m/sw/spack-levante/mambaforge-22.9.0-2-Linux-x86_64-kptncg/lib/python3.10/site-packages/pandas/core/generic.py:6002\u001b[0m, in \u001b[0;36mNDFrame.__setattr__\u001b[0;34m(self, name, value)\u001b[0m\n\u001b[1;32m 6000\u001b[0m \u001b[38;5;28;01mtry\u001b[39;00m:\n\u001b[1;32m 6001\u001b[0m \u001b[38;5;28mobject\u001b[39m\u001b[38;5;241m.\u001b[39m\u001b[38;5;21m__getattribute__\u001b[39m(\u001b[38;5;28mself\u001b[39m, name)\n\u001b[0;32m-> 6002\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;43mobject\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[38;5;21;43m__setattr__\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mname\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mvalue\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 6003\u001b[0m \u001b[38;5;28;01mexcept\u001b[39;00m \u001b[38;5;167;01mAttributeError\u001b[39;00m:\n\u001b[1;32m 6004\u001b[0m \u001b[38;5;28;01mpass\u001b[39;00m\n",
"File \u001b[0;32m/sw/spack-levante/mambaforge-22.9.0-2-Linux-x86_64-kptncg/lib/python3.10/site-packages/pandas/_libs/properties.pyx:69\u001b[0m, in \u001b[0;36mpandas._libs.properties.AxisProperty.__set__\u001b[0;34m()\u001b[0m\n",
"File \u001b[0;32m/sw/spack-levante/mambaforge-22.9.0-2-Linux-x86_64-kptncg/lib/python3.10/site-packages/pandas/core/generic.py:730\u001b[0m, in \u001b[0;36mNDFrame._set_axis\u001b[0;34m(self, axis, labels)\u001b[0m\n\u001b[1;32m 725\u001b[0m \u001b[38;5;250m\u001b[39m\u001b[38;5;124;03m\"\"\"\u001b[39;00m\n\u001b[1;32m 726\u001b[0m \u001b[38;5;124;03mThis is called from the cython code when we set the `index` attribute\u001b[39;00m\n\u001b[1;32m 727\u001b[0m \u001b[38;5;124;03mdirectly, e.g. `series.index = [1, 2, 3]`.\u001b[39;00m\n\u001b[1;32m 728\u001b[0m \u001b[38;5;124;03m\"\"\"\u001b[39;00m\n\u001b[1;32m 729\u001b[0m labels \u001b[38;5;241m=\u001b[39m ensure_index(labels)\n\u001b[0;32m--> 730\u001b[0m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43m_mgr\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mset_axis\u001b[49m\u001b[43m(\u001b[49m\u001b[43maxis\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mlabels\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 731\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39m_clear_item_cache()\n",
"File \u001b[0;32m/sw/spack-levante/mambaforge-22.9.0-2-Linux-x86_64-kptncg/lib/python3.10/site-packages/pandas/core/internals/managers.py:225\u001b[0m, in \u001b[0;36mBaseBlockManager.set_axis\u001b[0;34m(self, axis, new_labels)\u001b[0m\n\u001b[1;32m 223\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21mset_axis\u001b[39m(\u001b[38;5;28mself\u001b[39m, axis: AxisInt, new_labels: Index) \u001b[38;5;241m-\u001b[39m\u001b[38;5;241m>\u001b[39m \u001b[38;5;28;01mNone\u001b[39;00m:\n\u001b[1;32m 224\u001b[0m \u001b[38;5;66;03m# Caller is responsible for ensuring we have an Index object.\u001b[39;00m\n\u001b[0;32m--> 225\u001b[0m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43m_validate_set_axis\u001b[49m\u001b[43m(\u001b[49m\u001b[43maxis\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mnew_labels\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 226\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39maxes[axis] \u001b[38;5;241m=\u001b[39m new_labels\n",
"File \u001b[0;32m/sw/spack-levante/mambaforge-22.9.0-2-Linux-x86_64-kptncg/lib/python3.10/site-packages/pandas/core/internals/base.py:70\u001b[0m, in \u001b[0;36mDataManager._validate_set_axis\u001b[0;34m(self, axis, new_labels)\u001b[0m\n\u001b[1;32m 67\u001b[0m \u001b[38;5;28;01mpass\u001b[39;00m\n\u001b[1;32m 69\u001b[0m \u001b[38;5;28;01melif\u001b[39;00m new_len \u001b[38;5;241m!=\u001b[39m old_len:\n\u001b[0;32m---> 70\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mValueError\u001b[39;00m(\n\u001b[1;32m 71\u001b[0m \u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mLength mismatch: Expected axis has \u001b[39m\u001b[38;5;132;01m{\u001b[39;00mold_len\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m elements, new \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 72\u001b[0m \u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mvalues have \u001b[39m\u001b[38;5;132;01m{\u001b[39;00mnew_len\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m elements\u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 73\u001b[0m )\n",
"\u001b[0;31mValueError\u001b[0m: Length mismatch: Expected axis has 1 elements, new values have 2 elements"
]
}
],
"outputs": [],
"source": [
"cams_nox_file = '/home/b/b309253/data-exploration/emissions/cams-glob-ant-anthro-nox.csv'\n",
<<<<<<< HEAD
"cams_nox_data = pd.read_csv(cams_nox_file, sep=',')\n",
"cams_nox_data.columns = ['Date', 'CAMS Road']\n",
"cams_nox_data['CAMS Road'] = cams_nox_data['CAMS Road']\n",
=======
"cams_nox_data = pd.read_csv(cams_nox_file, sep='\\t')\n",
"cams_nox_data.columns = ['Date', 'CAMS Road']\n",
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"cams_nox_data = cams_nox_data.set_index(pd.date_range(\"1/1/2000\", periods=24, freq='Y')).drop(columns = ['Date'])"
]
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 49,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "34c63c20-e320-4323-8469-7e672d3cde19",
"metadata": {
"tags": []
......@@ -995,26 +965,17 @@
"ceds_nox_file = '/home/b/b309253/data-exploration/emissions/ceds-anthro-nox-yearly.csv'\n",
"ceds_nox_data = pd.read_csv(ceds_nox_file, sep=',')\n",
"ceds_nox_data.columns = ['Date', 'CEDS tranport']\n",
<<<<<<< HEAD
"ceds_nox_data['CEDS tranport'] = ceds_nox_data['CEDS tranport']\n",
=======
"ceds_nox_data\n",
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"ceds_nox_data = ceds_nox_data.set_index(pd.date_range(\"1/1/1950\", periods=65, freq='Y')).drop(columns = ['Date'])"
]
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 50,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "5da24dc3-098e-4ed7-9cec-e6cbb37697bb",
"metadata": {
"tags": []
},
<<<<<<< HEAD
"outputs": [
{
"data": {
......@@ -1037,9 +998,6 @@
"output_type": "display_data"
}
],
=======
"outputs": [],
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"source": [
"ax = ceds_nox_data.plot()\n",
"cams_nox_data.plot(ax=ax)\n",
......@@ -1059,11 +1017,7 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 26,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "febfd5ef-20f2-46ad-a352-0206af7e59ee",
"metadata": {
"tags": []
......@@ -1078,11 +1032,7 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 27,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "31f6de82-eca3-4ed3-aabc-b6124d6033c6",
"metadata": {
"tags": []
......@@ -1097,16 +1047,11 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 28,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "f0c7b549-97d3-4f6e-bdcb-cb03896d2799",
"metadata": {
"tags": []
},
<<<<<<< HEAD
"outputs": [
{
"data": {
......@@ -1129,9 +1074,6 @@
"output_type": "display_data"
}
],
=======
"outputs": [],
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"source": [
"ax = ceds_nh3_data.plot()\n",
"cams_nh3_data.plot(ax=ax)\n",
......@@ -1151,11 +1093,7 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 29,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "ea3c0a21-9aa3-48b4-a018-de01399de37a",
"metadata": {
"tags": []
......@@ -1170,11 +1108,7 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 30,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "b6dc224e-6598-421e-a5c4-088b4d5a6c48",
"metadata": {
"tags": []
......@@ -1189,16 +1123,11 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 31,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"id": "bc450025-f98f-4e40-ae97-f7ac5dd030c9",
"metadata": {
"tags": []
},
<<<<<<< HEAD
"outputs": [
{
"data": {
......@@ -1221,9 +1150,6 @@
"output_type": "display_data"
}
],
=======
"outputs": [],
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"source": [
"ax = ceds_so2_data.plot()\n",
"cams_so2_data.plot(ax=ax)\n",
......@@ -1278,7 +1204,6 @@
"metadata": {
"tags": []
},
<<<<<<< HEAD
"outputs": [
{
"data": {
......@@ -1301,9 +1226,6 @@
"output_type": "display_data"
}
],
=======
"outputs": [],
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"source": [
"ax = ceds_co_data.plot()\n",
"cams_co_data.plot(ax=ax)\n",
......@@ -1323,16 +1245,11 @@
},
{
"cell_type": "code",
<<<<<<< HEAD
"execution_count": 36,
=======
"execution_count": null,
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"execution_count": 51,
"id": "4375409e-3a29-4e3b-82d0-2eb6449c4b08",
"metadata": {
"tags": []
},
<<<<<<< HEAD
"outputs": [
{
"data": {
......@@ -1345,22 +1262,13 @@
"output_type": "display_data"
}
],
=======
"outputs": [],
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"source": [
"fig4, axs4 = plt.subplots(2,2)\n",
"fig4.set_figwidth(12)\n",
"\n",
"\n",
"NOx_hist_emissions[220:].plot(ax= axs4[0,0])\n",
<<<<<<< HEAD
"ceds_nox_data.plot(ax= axs4[0,0])\n",
"cams_nox_data.plot(ax= axs4[0,0])\n",
=======
"ceds_data.plot(ax= axs4[0,0])\n",
"cams_data.plot(ax= axs4[0,0])\n",
>>>>>>> d1a035698f335f4838e24def47fa9bbc9e5fed3d
"emissions_NOx_road.plot(title = 'Yearly total road emissions of NOx across scenarios', ax= axs4[0,0])\n",
"NH3_hist_emissions[220:].plot(ax= axs4[0,1])\n",
"ceds_nh3_data.plot(ax= axs4[0,1])\n",
......@@ -1378,7 +1286,6 @@
"axs4[0,0].set_ylabel('Tg(NO)')\n",
"axs4[1,0].set_ylabel('Tg(NO)')\n",
"\n",
"\n",
"fig4.tight_layout()"
]
},
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Imprint | Privacy Policy