From a471b8d92ee8e48e6873d995ea4d21a02f80b35d Mon Sep 17 00:00:00 2001 From: Alba Vilanova Cortezon Date: Wed, 31 Aug 2022 18:25:17 +0200 Subject: [PATCH] Fix Providentia interpolation --- CHANGELOG.md | 9 +- nes/nc_projections/points_nes.py | 3 +- .../4.3.Providentia_Interpolation.ipynb | 1491 +++++++---------- 3 files changed, 615 insertions(+), 888 deletions(-) diff --git a/CHANGELOG.md b/CHANGELOG.md index 1677780..57f4dc2 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -1,7 +1,7 @@ # NES CHANGELOG -### 0.9.0 -* Release date: 2022/08/12 +### 1.0.0 +* Release date: 2022/08/31 * Changes and new features: * First beta release * Open: @@ -11,8 +11,8 @@ * Lambert Conformal Conic * Mercator * Points - * Points * Points in GHOST format + * Points in Providentia format * Parallelization: * Balanced / Unbalanced * By time axis @@ -25,13 +25,16 @@ * Lambert Conformal Conic * Mercator * Points + * Shapefile * Write: * NetCDF * Grib2 + * Shapefile * Interpolation: * Vertical interpolation * Horizontal interpolation * Nearest Neighbours + * Providentia interpolation * Statistics: * Daily_mean * Daily_max diff --git a/nes/nc_projections/points_nes.py b/nes/nc_projections/points_nes.py index 1847fde..1f4f5d7 100644 --- a/nes/nc_projections/points_nes.py +++ b/nes/nc_projections/points_nes.py @@ -612,8 +612,7 @@ class PointsNes(Nes): print(var_name) variables[var_name] = {} if var_info['dimensions'] == ('time', 'lat', 'lon') and len(var_info['data'].shape) == 2: - variables[var_name]['data'] = var_info['data'].reshape((var_info['data'].shape[1], - var_info['data'].shape[0])) + variables[var_name]['data'] = var_info['data'].T variables[var_name]['dimensions'] = ('station', 'time') else: variables[var_name]['data'] = var_info['data'] diff --git a/tutorials/4.Interpolation/4.3.Providentia_Interpolation.ipynb b/tutorials/4.Interpolation/4.3.Providentia_Interpolation.ipynb index 9c9e928..43456b9 100644 --- a/tutorials/4.Interpolation/4.3.Providentia_Interpolation.ipynb +++ b/tutorials/4.Interpolation/4.3.Providentia_Interpolation.ipynb @@ -41,7 +41,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 3, @@ -216,6 +216,135 @@ "source_grid.load()" ] }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
geometry
FID
0POLYGON ((-25.10000 29.90000, -24.90000 29.900...
1POLYGON ((-24.90000 29.90000, -24.70000 29.900...
2POLYGON ((-24.70000 29.90000, -24.50000 29.900...
3POLYGON ((-24.50000 29.90000, -24.30000 29.900...
4POLYGON ((-24.30000 29.90000, -24.10000 29.900...
......
74056POLYGON ((44.10000 71.90000, 44.30000 71.90000...
74057POLYGON ((44.30000 71.90000, 44.50000 71.90000...
74058POLYGON ((44.50000 71.90000, 44.70000 71.90000...
74059POLYGON ((44.70000 71.90000, 44.90000 71.90000...
74060POLYGON ((44.90000 71.90000, 45.10000 71.90000...
\n", + "

74061 rows × 1 columns

\n", + "
" + ], + "text/plain": [ + " geometry\n", + "FID \n", + "0 POLYGON ((-25.10000 29.90000, -24.90000 29.900...\n", + "1 POLYGON ((-24.90000 29.90000, -24.70000 29.900...\n", + "2 POLYGON ((-24.70000 29.90000, -24.50000 29.900...\n", + "3 POLYGON ((-24.50000 29.90000, -24.30000 29.900...\n", + "4 POLYGON ((-24.30000 29.90000, -24.10000 29.900...\n", + "... ...\n", + "74056 POLYGON ((44.10000 71.90000, 44.30000 71.90000...\n", + "74057 POLYGON ((44.30000 71.90000, 44.50000 71.90000...\n", + "74058 POLYGON ((44.50000 71.90000, 44.70000 71.90000...\n", + "74059 POLYGON ((44.70000 71.90000, 44.90000 71.90000...\n", + "74060 POLYGON ((44.90000 71.90000, 45.10000 71.90000...\n", + "\n", + "[74061 rows x 1 columns]" + ] + }, + "execution_count": 8, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "source_grid.create_shapefile()\n", + "source_grid.shapefile" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [], + "source": [ + "source_grid.shapefile['sconco3'] = source_grid.variables['sconco3']['data'][0, 0, :, :].ravel()" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [], + "source": [ + "source_grid.write_shapefile('model')" + ] + }, { "cell_type": "markdown", "metadata": {}, @@ -225,7 +354,7 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": 11, "metadata": {}, "outputs": [ { @@ -239,84 +368,81 @@ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 8, + "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "#dst_grid_path = '/gpfs/projects/bsc32/AC_cache/obs/ghost/EBAS/1.4/hourly/sconco3/sconco3_201804.nc'\n", - "dst_grid_path = '/gpfs/projects/bsc32/AC_cache/recon/exp_interp/1.3.3/cams61_chimere_ph2-eu-000/hourly/sconco3/EBAS/sconco3_201804.nc'\n", + "dst_grid_path = '/gpfs/projects/bsc32/AC_cache/obs/ghost/EBAS/1.4/hourly/sconco3/sconco3_201804.nc'\n", "dst_grid = open_netcdf(path=dst_grid_path, info=True)\n", "dst_grid" ] }, { "cell_type": "code", - "execution_count": 9, + "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "{'data': masked_array(data=[-64.24006 , -54.84846497, -22.10333333, -31.66861111,\n", - " 47.76666641, 46.677778 , 48.721111 , 47.529167 ,\n", - " 47.05407 , 46.693611 , 47.348056 , 47.973056 ,\n", - " 48.878611 , 48.106111 , 48.371111 , 48.334722 ,\n", - " 48.050833 , 47.838611 , 47.040277 , 47.06694444,\n", - " 49.877778 , 50.629421 , 50.503333 , 41.695833 ,\n", - " 32.27000046, 80.05000305, 46.5475 , 46.813056 ,\n", - " 47.479722 , 47.049722 , 47.0675 , 47.18961391,\n", - " -30.17254 , 16.86403 , 35.0381 , 49.73508444,\n", - " 49.573394 , 49.066667 , 54.925556 , 52.802222 ,\n", - " 47.914722 , 53.166667 , 50.65 , 54.4368 ,\n", - " 47.80149841, 47.4165 , -70.666 , 54.746495 ,\n", - " 81.6 , 55.693588 , 72.58000183, 56.290424 ,\n", - " 59.5 , 58.383333 , 39.54694 , 42.72056 ,\n", - " 39.87528 , 37.23722 , 43.43917 , 41.27417 ,\n", - " 42.31917 , 38.47278 , 39.08278 , 41.23889 ,\n", - " 41.39389 , 42.63472 , 37.05194 , 28.309 ,\n", - " 59.779167 , 60.53002 , 66.320278 , 67.97333333,\n", - " 48.5 , 49.9 , 47.266667 , 43.616667 ,\n", - " 47.3 , 46.65 , 45. , 45.8 ,\n", - " 48.633333 , 42.936667 , 48.70861111, 44.56944444,\n", + "{'data': masked_array(data=[-64.24006 , -54.84846497, 47.76666641, 46.677778 ,\n", + " 48.721111 , 47.529167 , 47.05407 , 47.348056 ,\n", + " 47.973056 , 48.878611 , 48.106111 , 48.371111 ,\n", + " 48.334722 , 48.050833 , 47.838611 , 47.040277 ,\n", + " 47.06694444, 49.877778 , 50.629421 , 50.503333 ,\n", + " 41.695833 , 32.27000046, 80.05000305, 46.5475 ,\n", + " 46.813056 , 47.479722 , 47.049722 , 47.0675 ,\n", + " 47.18961391, -30.17254 , 16.86403 , 35.0381 ,\n", + " 49.73508444, 49.573394 , 49.066667 , 54.925556 ,\n", + " 52.802222 , 47.914722 , 53.166667 , 50.65 ,\n", + " 54.4368 , 47.80149841, 47.4165 , -70.666 ,\n", + " 54.746495 , 81.6 , 55.693588 , 72.58000183,\n", + " 56.290424 , 59.5 , 58.383333 , 39.54694 ,\n", + " 42.72056 , 39.87528 , 37.23722 , 43.43917 ,\n", + " 41.27417 , 42.31917 , 38.47278 , 39.08278 ,\n", + " 41.23889 , 41.39389 , 42.63472 , 37.05194 ,\n", + " 28.309 , 59.779167 , 60.53002 , 66.320278 ,\n", + " 67.97333333, 48.5 , 49.9 , 47.266667 ,\n", + " 43.616667 , 47.3 , 46.65 , 45. ,\n", + " 45.8 , 48.633333 , 42.936667 , 44.56944444,\n", " 46.81472778, 45.772223 , 55.313056 , 54.443056 ,\n", " 50.596389 , 54.334444 , 57.734444 , 52.503889 ,\n", " 55.858611 , 53.398889 , 50.792778 , 52.293889 ,\n", " 51.781784 , 52.298333 , 55.79216 , 52.950556 ,\n", - " 51.778056 , 60.13922 , -75.62 , 51.149617 ,\n", - " 38.366667 , 35.316667 , 46.966667 , 46.91 ,\n", - " -0.20194 , 51.939722 , 53.32583 , 45.8 ,\n", - " 44.183333 , 37.571111 , 35.5182 , 42.805462 ,\n", - " -69.005 , 39.0319 , 24.2883 , 24.466941 ,\n", - " 36.53833389, 33.293917 , 55.37611111, 56.161944 ,\n", - " 57.135278 , 41.536111 , 36.0722 , 52.083333 ,\n", - " 53.333889 , 51.541111 , 52.3 , 51.974444 ,\n", - " 58.38853 , 65.833333 , 62.783333 , 78.90715 ,\n", - " 59. , 69.45 , 59.2 , 60.372386 ,\n", - " -72.0117 , 59.2 , -41.40819168, -77.83200073,\n", - " -45.0379982 , 51.814408 , 50.736444 , 54.753894 ,\n", - " 54.15 , 43.4 , 71.58616638, 63.85 ,\n", - " 67.883333 , 57.394 , 57.1645 , 57.9525 ,\n", - " 56.0429 , 60.0858 , 57.816667 , 64.25 ,\n", - " 59.728 , 45.566667 , 46.428611 , 46.299444 ,\n", - " 48.933333 , 49.15 , 49.05 , 47.96 ,\n", - " 71.32301331, 40.12498 , 19.53623009, -89.99694824,\n", - " 41.05410004, 21.5731 , -34.35348 ],\n", + " 51.778056 , 60.13922 , 51.149617 , 38.366667 ,\n", + " 35.316667 , 46.966667 , 46.91 , -0.20194 ,\n", + " 51.939722 , 53.32583 , 45.8 , 44.183333 ,\n", + " 37.571111 , 42.805462 , -69.005 , 39.0319 ,\n", + " 24.2883 , 24.466941 , 36.53833389, 33.293917 ,\n", + " 55.37611111, 56.161944 , 57.135278 , 36.0722 ,\n", + " 52.083333 , 53.333889 , 51.541111 , 52.3 ,\n", + " 51.974444 , 58.38853 , 65.833333 , 62.783333 ,\n", + " 78.90715 , 59. , 69.45 , 59.2 ,\n", + " 60.372386 , -72.0117 , 59.2 , -41.40819168,\n", + " -77.83200073, -45.0379982 , 51.814408 , 50.736444 ,\n", + " 54.753894 , 54.15 , 43.4 , 71.58616638,\n", + " 63.85 , 67.883333 , 57.394 , 57.1645 ,\n", + " 57.9525 , 56.0429 , 60.0858 , 57.816667 ,\n", + " 64.25 , 59.728 , 45.566667 , 46.428611 ,\n", + " 46.299444 , 48.933333 , 49.15 , 49.05 ,\n", + " 47.96 , 71.32301331, 40.12498 , 19.53623009,\n", + " -89.99694824, 41.05410004, 21.5731 , -34.35348 ],\n", " mask=False,\n", " fill_value=1e+20),\n", " 'dimensions': ('station',),\n", " 'standard_name': 'latitude',\n", - " 'units': 'decimal degrees North',\n", " 'long_name': 'latitude',\n", + " 'units': 'decimal degrees North',\n", " 'description': 'Geodetic latitude of measuring instrument, in decimal degrees North, following the stated horizontal datum.',\n", " 'axis': 'Y'}" ] }, - "execution_count": 9, + "execution_count": 12, "metadata": {}, "output_type": "execute_result" } @@ -327,16 +453,15 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "{'data': masked_array(data=[-5.66247800e+01, -6.83106918e+01, -6.56008333e+01,\n", - " -6.38819444e+01, 1.67666664e+01, 1.29722220e+01,\n", - " 1.59422220e+01, 9.92666700e+00, 1.29579400e+01,\n", - " 1.39150000e+01, 1.58822220e+01, 1.30161110e+01,\n", + "{'data': masked_array(data=[-5.66247800e+01, -6.83106918e+01, 1.67666664e+01,\n", + " 1.29722220e+01, 1.59422220e+01, 9.92666700e+00,\n", + " 1.29579400e+01, 1.58822220e+01, 1.30161110e+01,\n", " 1.50466670e+01, 1.59194440e+01, 1.55466670e+01,\n", " 1.67305560e+01, 1.66766670e+01, 1.44413890e+01,\n", " 1.43300000e+01, 1.54936111e+01, 5.20361100e+00,\n", @@ -360,49 +485,47 @@ " 7.13333300e+00, 4.63333300e+00, 4.08333300e+00,\n", " 1.83333000e-01, 6.83333300e+00, -7.50000000e-01,\n", " 6.46666700e+00, 2.06666700e+00, -4.50000000e-01,\n", - " 1.41944000e-01, 2.15888889e+00, 5.27897222e+00,\n", - " 2.61000833e+00, 2.96488600e+00, -3.20416700e+00,\n", - " -7.87000000e+00, -3.71305600e+00, -8.07500000e-01,\n", - " -4.77444400e+00, -3.03305600e+00, -3.20500000e+00,\n", - " -1.75333300e+00, 1.79444000e-01, 1.46305600e+00,\n", - " -4.69146200e+00, 2.92778000e-01, -3.24290000e+00,\n", - " 1.12194400e+00, 1.08223000e+00, -1.18531900e+00,\n", - " -2.61800000e+01, -1.43822800e+00, 2.30833330e+01,\n", - " 2.56666670e+01, 1.95833330e+01, 1.63200000e+01,\n", - " 1.00318100e+02, -1.02444440e+01, -9.89944000e+00,\n", - " 8.63333300e+00, 1.07000000e+01, 1.26597220e+01,\n", - " 1.26305000e+01, 1.25656450e+01, 3.95905556e+01,\n", + " 1.41944000e-01, 5.27897222e+00, 2.61000833e+00,\n", + " 2.96488600e+00, -3.20416700e+00, -7.87000000e+00,\n", + " -3.71305600e+00, -8.07500000e-01, -4.77444400e+00,\n", + " -3.03305600e+00, -3.20500000e+00, -1.75333300e+00,\n", + " 1.79444000e-01, 1.46305600e+00, -4.69146200e+00,\n", + " 2.92778000e-01, -3.24290000e+00, 1.12194400e+00,\n", + " 1.08223000e+00, -1.18531900e+00, -1.43822800e+00,\n", + " 2.30833330e+01, 2.56666670e+01, 1.95833330e+01,\n", + " 1.63200000e+01, 1.00318100e+02, -1.02444440e+01,\n", + " -9.89944000e+00, 8.63333300e+00, 1.07000000e+01,\n", + " 1.26597220e+01, 1.25656450e+01, 3.95905556e+01,\n", " 1.41822200e+02, 1.53983300e+02, 1.23010872e+02,\n", " 1.26330002e+02, 1.26163111e+02, 2.10305556e+01,\n", - " 2.11730560e+01, 2.59055560e+01, 2.06938900e+01,\n", - " 1.42184000e+01, 6.56666700e+00, 6.27722200e+00,\n", - " 5.85361100e+00, 4.50000000e+00, 4.92361100e+00,\n", - " 8.25200000e+00, 1.39166670e+01, 8.88333300e+00,\n", - " 1.18866800e+01, 1.15333330e+01, 3.00333330e+01,\n", - " 5.20000000e+00, 1.10781420e+01, 2.53510000e+00,\n", - " 9.51666700e+00, 1.74870804e+02, 1.66660004e+02,\n", - " 1.69684006e+02, 2.19724190e+01, 1.57395000e+01,\n", - " 1.75342640e+01, 2.20666670e+01, 2.19500000e+01,\n", - " 1.28918823e+02, 1.53333330e+01, 2.10666670e+01,\n", - " 1.19140000e+01, 1.47825000e+01, 1.24030000e+01,\n", - " 1.31480000e+01, 1.75052800e+01, 1.55666670e+01,\n", - " 1.97666670e+01, 1.54720000e+01, 1.48666670e+01,\n", - " 1.50033330e+01, 1.45386110e+01, 1.95833330e+01,\n", - " 2.02833330e+01, 2.22666670e+01, 1.78605560e+01,\n", - " -1.56611465e+02, -1.05236800e+02, -1.55576157e+02,\n", - " -2.47999992e+01, -1.24151001e+02, 1.03515700e+02,\n", - " 1.84896800e+01],\n", + " 2.11730560e+01, 2.59055560e+01, 1.42184000e+01,\n", + " 6.56666700e+00, 6.27722200e+00, 5.85361100e+00,\n", + " 4.50000000e+00, 4.92361100e+00, 8.25200000e+00,\n", + " 1.39166670e+01, 8.88333300e+00, 1.18866800e+01,\n", + " 1.15333330e+01, 3.00333330e+01, 5.20000000e+00,\n", + " 1.10781420e+01, 2.53510000e+00, 9.51666700e+00,\n", + " 1.74870804e+02, 1.66660004e+02, 1.69684006e+02,\n", + " 2.19724190e+01, 1.57395000e+01, 1.75342640e+01,\n", + " 2.20666670e+01, 2.19500000e+01, 1.28918823e+02,\n", + " 1.53333330e+01, 2.10666670e+01, 1.19140000e+01,\n", + " 1.47825000e+01, 1.24030000e+01, 1.31480000e+01,\n", + " 1.75052800e+01, 1.55666670e+01, 1.97666670e+01,\n", + " 1.54720000e+01, 1.48666670e+01, 1.50033330e+01,\n", + " 1.45386110e+01, 1.95833330e+01, 2.02833330e+01,\n", + " 2.22666670e+01, 1.78605560e+01, -1.56611465e+02,\n", + " -1.05236800e+02, -1.55576157e+02, -2.47999992e+01,\n", + " -1.24151001e+02, 1.03515700e+02, 1.84896800e+01],\n", " mask=False,\n", " fill_value=1e+20),\n", " 'dimensions': ('station',),\n", " 'standard_name': 'longitude',\n", - " 'units': 'decimal degrees East',\n", " 'long_name': 'longitude',\n", + " 'units': 'decimal degrees East',\n", " 'description': 'Geodetic longitude of measuring instrument, in decimal degrees East, following the stated horizontal datum.',\n", " 'axis': 'X'}" ] }, - "execution_count": 10, + "execution_count": 13, "metadata": {}, "output_type": "execute_result" } @@ -413,7 +536,7 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": 14, "metadata": {}, "outputs": [ { @@ -433,70 +556,59 @@ }, { "cell_type": "code", - "execution_count": 12, - "metadata": {}, - "outputs": [], - "source": [ - "interpolated_source_grid.netcdf" - ] - }, - { - "cell_type": "code", - "execution_count": 13, + "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "{'data': masked_array(data=[-64.24006 , -54.84846497, -22.10333333, -31.66861111,\n", - " 47.76666641, 46.677778 , 48.721111 , 47.529167 ,\n", - " 47.05407 , 46.693611 , 47.348056 , 47.973056 ,\n", - " 48.878611 , 48.106111 , 48.371111 , 48.334722 ,\n", - " 48.050833 , 47.838611 , 47.040277 , 47.06694444,\n", - " 49.877778 , 50.629421 , 50.503333 , 41.695833 ,\n", - " 32.27000046, 80.05000305, 46.5475 , 46.813056 ,\n", - " 47.479722 , 47.049722 , 47.0675 , 47.18961391,\n", - " -30.17254 , 16.86403 , 35.0381 , 49.73508444,\n", - " 49.573394 , 49.066667 , 54.925556 , 52.802222 ,\n", - " 47.914722 , 53.166667 , 50.65 , 54.4368 ,\n", - " 47.80149841, 47.4165 , -70.666 , 54.746495 ,\n", - " 81.6 , 55.693588 , 72.58000183, 56.290424 ,\n", - " 59.5 , 58.383333 , 39.54694 , 42.72056 ,\n", - " 39.87528 , 37.23722 , 43.43917 , 41.27417 ,\n", - " 42.31917 , 38.47278 , 39.08278 , 41.23889 ,\n", - " 41.39389 , 42.63472 , 37.05194 , 28.309 ,\n", - " 59.779167 , 60.53002 , 66.320278 , 67.97333333,\n", - " 48.5 , 49.9 , 47.266667 , 43.616667 ,\n", - " 47.3 , 46.65 , 45. , 45.8 ,\n", - " 48.633333 , 42.936667 , 48.70861111, 44.56944444,\n", + "{'data': masked_array(data=[-64.24006 , -54.84846497, 47.76666641, 46.677778 ,\n", + " 48.721111 , 47.529167 , 47.05407 , 47.348056 ,\n", + " 47.973056 , 48.878611 , 48.106111 , 48.371111 ,\n", + " 48.334722 , 48.050833 , 47.838611 , 47.040277 ,\n", + " 47.06694444, 49.877778 , 50.629421 , 50.503333 ,\n", + " 41.695833 , 32.27000046, 80.05000305, 46.5475 ,\n", + " 46.813056 , 47.479722 , 47.049722 , 47.0675 ,\n", + " 47.18961391, -30.17254 , 16.86403 , 35.0381 ,\n", + " 49.73508444, 49.573394 , 49.066667 , 54.925556 ,\n", + " 52.802222 , 47.914722 , 53.166667 , 50.65 ,\n", + " 54.4368 , 47.80149841, 47.4165 , -70.666 ,\n", + " 54.746495 , 81.6 , 55.693588 , 72.58000183,\n", + " 56.290424 , 59.5 , 58.383333 , 39.54694 ,\n", + " 42.72056 , 39.87528 , 37.23722 , 43.43917 ,\n", + " 41.27417 , 42.31917 , 38.47278 , 39.08278 ,\n", + " 41.23889 , 41.39389 , 42.63472 , 37.05194 ,\n", + " 28.309 , 59.779167 , 60.53002 , 66.320278 ,\n", + " 67.97333333, 48.5 , 49.9 , 47.266667 ,\n", + " 43.616667 , 47.3 , 46.65 , 45. ,\n", + " 45.8 , 48.633333 , 42.936667 , 44.56944444,\n", " 46.81472778, 45.772223 , 55.313056 , 54.443056 ,\n", " 50.596389 , 54.334444 , 57.734444 , 52.503889 ,\n", " 55.858611 , 53.398889 , 50.792778 , 52.293889 ,\n", " 51.781784 , 52.298333 , 55.79216 , 52.950556 ,\n", - " 51.778056 , 60.13922 , -75.62 , 51.149617 ,\n", - " 38.366667 , 35.316667 , 46.966667 , 46.91 ,\n", - " -0.20194 , 51.939722 , 53.32583 , 45.8 ,\n", - " 44.183333 , 37.571111 , 35.5182 , 42.805462 ,\n", - " -69.005 , 39.0319 , 24.2883 , 24.466941 ,\n", - " 36.53833389, 33.293917 , 55.37611111, 56.161944 ,\n", - " 57.135278 , 41.536111 , 36.0722 , 52.083333 ,\n", - " 53.333889 , 51.541111 , 52.3 , 51.974444 ,\n", - " 58.38853 , 65.833333 , 62.783333 , 78.90715 ,\n", - " 59. , 69.45 , 59.2 , 60.372386 ,\n", - " -72.0117 , 59.2 , -41.40819168, -77.83200073,\n", - " -45.0379982 , 51.814408 , 50.736444 , 54.753894 ,\n", - " 54.15 , 43.4 , 71.58616638, 63.85 ,\n", - " 67.883333 , 57.394 , 57.1645 , 57.9525 ,\n", - " 56.0429 , 60.0858 , 57.816667 , 64.25 ,\n", - " 59.728 , 45.566667 , 46.428611 , 46.299444 ,\n", - " 48.933333 , 49.15 , 49.05 , 47.96 ,\n", - " 71.32301331, 40.12498 , 19.53623009, -89.99694824,\n", - " 41.05410004, 21.5731 , -34.35348 ],\n", + " 51.778056 , 60.13922 , 51.149617 , 38.366667 ,\n", + " 35.316667 , 46.966667 , 46.91 , -0.20194 ,\n", + " 51.939722 , 53.32583 , 45.8 , 44.183333 ,\n", + " 37.571111 , 42.805462 , -69.005 , 39.0319 ,\n", + " 24.2883 , 24.466941 , 36.53833389, 33.293917 ,\n", + " 55.37611111, 56.161944 , 57.135278 , 36.0722 ,\n", + " 52.083333 , 53.333889 , 51.541111 , 52.3 ,\n", + " 51.974444 , 58.38853 , 65.833333 , 62.783333 ,\n", + " 78.90715 , 59. , 69.45 , 59.2 ,\n", + " 60.372386 , -72.0117 , 59.2 , -41.40819168,\n", + " -77.83200073, -45.0379982 , 51.814408 , 50.736444 ,\n", + " 54.753894 , 54.15 , 43.4 , 71.58616638,\n", + " 63.85 , 67.883333 , 57.394 , 57.1645 ,\n", + " 57.9525 , 56.0429 , 60.0858 , 57.816667 ,\n", + " 64.25 , 59.728 , 45.566667 , 46.428611 ,\n", + " 46.299444 , 48.933333 , 49.15 , 49.05 ,\n", + " 47.96 , 71.32301331, 40.12498 , 19.53623009,\n", + " -89.99694824, 41.05410004, 21.5731 , -34.35348 ],\n", " mask=False,\n", " fill_value=1e+20)}" ] }, - "execution_count": 13, + "execution_count": 15, "metadata": {}, "output_type": "execute_result" } @@ -507,16 +619,15 @@ }, { "cell_type": "code", - "execution_count": 14, + "execution_count": 16, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "{'data': masked_array(data=[-5.66247800e+01, -6.83106918e+01, -6.56008333e+01,\n", - " -6.38819444e+01, 1.67666664e+01, 1.29722220e+01,\n", - " 1.59422220e+01, 9.92666700e+00, 1.29579400e+01,\n", - " 1.39150000e+01, 1.58822220e+01, 1.30161110e+01,\n", + "{'data': masked_array(data=[-5.66247800e+01, -6.83106918e+01, 1.67666664e+01,\n", + " 1.29722220e+01, 1.59422220e+01, 9.92666700e+00,\n", + " 1.29579400e+01, 1.58822220e+01, 1.30161110e+01,\n", " 1.50466670e+01, 1.59194440e+01, 1.55466670e+01,\n", " 1.67305560e+01, 1.66766670e+01, 1.44413890e+01,\n", " 1.43300000e+01, 1.54936111e+01, 5.20361100e+00,\n", @@ -540,43 +651,41 @@ " 7.13333300e+00, 4.63333300e+00, 4.08333300e+00,\n", " 1.83333000e-01, 6.83333300e+00, -7.50000000e-01,\n", " 6.46666700e+00, 2.06666700e+00, -4.50000000e-01,\n", - " 1.41944000e-01, 2.15888889e+00, 5.27897222e+00,\n", - " 2.61000833e+00, 2.96488600e+00, -3.20416700e+00,\n", - " -7.87000000e+00, -3.71305600e+00, -8.07500000e-01,\n", - " -4.77444400e+00, -3.03305600e+00, -3.20500000e+00,\n", - " -1.75333300e+00, 1.79444000e-01, 1.46305600e+00,\n", - " -4.69146200e+00, 2.92778000e-01, -3.24290000e+00,\n", - " 1.12194400e+00, 1.08223000e+00, -1.18531900e+00,\n", - " -2.61800000e+01, -1.43822800e+00, 2.30833330e+01,\n", - " 2.56666670e+01, 1.95833330e+01, 1.63200000e+01,\n", - " 1.00318100e+02, -1.02444440e+01, -9.89944000e+00,\n", - " 8.63333300e+00, 1.07000000e+01, 1.26597220e+01,\n", - " 1.26305000e+01, 1.25656450e+01, 3.95905556e+01,\n", + " 1.41944000e-01, 5.27897222e+00, 2.61000833e+00,\n", + " 2.96488600e+00, -3.20416700e+00, -7.87000000e+00,\n", + " -3.71305600e+00, -8.07500000e-01, -4.77444400e+00,\n", + " -3.03305600e+00, -3.20500000e+00, -1.75333300e+00,\n", + " 1.79444000e-01, 1.46305600e+00, -4.69146200e+00,\n", + " 2.92778000e-01, -3.24290000e+00, 1.12194400e+00,\n", + " 1.08223000e+00, -1.18531900e+00, -1.43822800e+00,\n", + " 2.30833330e+01, 2.56666670e+01, 1.95833330e+01,\n", + " 1.63200000e+01, 1.00318100e+02, -1.02444440e+01,\n", + " -9.89944000e+00, 8.63333300e+00, 1.07000000e+01,\n", + " 1.26597220e+01, 1.25656450e+01, 3.95905556e+01,\n", " 1.41822200e+02, 1.53983300e+02, 1.23010872e+02,\n", " 1.26330002e+02, 1.26163111e+02, 2.10305556e+01,\n", - " 2.11730560e+01, 2.59055560e+01, 2.06938900e+01,\n", - " 1.42184000e+01, 6.56666700e+00, 6.27722200e+00,\n", - " 5.85361100e+00, 4.50000000e+00, 4.92361100e+00,\n", - " 8.25200000e+00, 1.39166670e+01, 8.88333300e+00,\n", - " 1.18866800e+01, 1.15333330e+01, 3.00333330e+01,\n", - " 5.20000000e+00, 1.10781420e+01, 2.53510000e+00,\n", - " 9.51666700e+00, 1.74870804e+02, 1.66660004e+02,\n", - " 1.69684006e+02, 2.19724190e+01, 1.57395000e+01,\n", - " 1.75342640e+01, 2.20666670e+01, 2.19500000e+01,\n", - " 1.28918823e+02, 1.53333330e+01, 2.10666670e+01,\n", - " 1.19140000e+01, 1.47825000e+01, 1.24030000e+01,\n", - " 1.31480000e+01, 1.75052800e+01, 1.55666670e+01,\n", - " 1.97666670e+01, 1.54720000e+01, 1.48666670e+01,\n", - " 1.50033330e+01, 1.45386110e+01, 1.95833330e+01,\n", - " 2.02833330e+01, 2.22666670e+01, 1.78605560e+01,\n", - " -1.56611465e+02, -1.05236800e+02, -1.55576157e+02,\n", - " -2.47999992e+01, -1.24151001e+02, 1.03515700e+02,\n", - " 1.84896800e+01],\n", + " 2.11730560e+01, 2.59055560e+01, 1.42184000e+01,\n", + " 6.56666700e+00, 6.27722200e+00, 5.85361100e+00,\n", + " 4.50000000e+00, 4.92361100e+00, 8.25200000e+00,\n", + " 1.39166670e+01, 8.88333300e+00, 1.18866800e+01,\n", + " 1.15333330e+01, 3.00333330e+01, 5.20000000e+00,\n", + " 1.10781420e+01, 2.53510000e+00, 9.51666700e+00,\n", + " 1.74870804e+02, 1.66660004e+02, 1.69684006e+02,\n", + " 2.19724190e+01, 1.57395000e+01, 1.75342640e+01,\n", + " 2.20666670e+01, 2.19500000e+01, 1.28918823e+02,\n", + " 1.53333330e+01, 2.10666670e+01, 1.19140000e+01,\n", + " 1.47825000e+01, 1.24030000e+01, 1.31480000e+01,\n", + " 1.75052800e+01, 1.55666670e+01, 1.97666670e+01,\n", + " 1.54720000e+01, 1.48666670e+01, 1.50033330e+01,\n", + " 1.45386110e+01, 1.95833330e+01, 2.02833330e+01,\n", + " 2.22666670e+01, 1.78605560e+01, -1.56611465e+02,\n", + " -1.05236800e+02, -1.55576157e+02, -2.47999992e+01,\n", + " -1.24151001e+02, 1.03515700e+02, 1.84896800e+01],\n", " mask=False,\n", " fill_value=1e+20)}" ] }, - "execution_count": 14, + "execution_count": 16, "metadata": {}, "output_type": "execute_result" } @@ -585,115 +694,204 @@ "interpolated_source_grid.lon" ] }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "## 3. Compare outputs" - ] - }, - { - "cell_type": "markdown", - "metadata": {}, - "source": [ - "### From NES" - ] - }, { "cell_type": "code", - "execution_count": 15, + "execution_count": 17, "metadata": {}, "outputs": [ { "data": { + "text/html": [ + "
\n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
geometry
FID
0POINT (-56.62478 -64.24006)
1POINT (-68.31069 -54.84846)
2POINT (16.76667 47.76667)
3POINT (12.97222 46.67778)
4POINT (15.94222 48.72111)
......
163POINT (-155.57616 19.53623)
164POINT (-24.80000 -89.99695)
165POINT (-124.15100 41.05410)
166POINT (103.51570 21.57310)
167POINT (18.48968 -34.35348)
\n", + "

168 rows × 1 columns

\n", + "
" + ], "text/plain": [ - "" + " geometry\n", + "FID \n", + "0 POINT (-56.62478 -64.24006)\n", + "1 POINT (-68.31069 -54.84846)\n", + "2 POINT (16.76667 47.76667)\n", + "3 POINT (12.97222 46.67778)\n", + "4 POINT (15.94222 48.72111)\n", + ".. ...\n", + "163 POINT (-155.57616 19.53623)\n", + "164 POINT (-24.80000 -89.99695)\n", + "165 POINT (-124.15100 41.05410)\n", + "166 POINT (103.51570 21.57310)\n", + "167 POINT (18.48968 -34.35348)\n", + "\n", + "[168 rows x 1 columns]" ] }, - "execution_count": 15, + "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "interpolated_source_grid" + "interpolated_source_grid.create_shapefile()\n", + "interpolated_source_grid.shapefile" ] }, { "cell_type": "code", - "execution_count": 16, + "execution_count": 18, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "array([[53.66018016, 53.72655496, 53.72658317, ..., 43.57299226,\n", - " 47.24854741, 32.18078968],\n", - " [33.18998267, 31.79456739, 27.61967959, ..., 39.97857755,\n", - " 40.50389612, 37.44205106],\n", - " [39.16614765, 38.60188681, 33.58534075, ..., 46.94782439,\n", - " 42.2582876 , 44.41170698],\n", - " ...,\n", - " [47.78859072, 53.85937227, 40.06637886, ..., 46.7489719 ,\n", - " 57.41357903, 39.15433177],\n", - " [52.09374241, 39.10817074, 44.6270592 , ..., 34.05485287,\n", - " 45.79934222, 46.92515764],\n", - " [47.54277469, 47.60884789, 38.51539661, ..., 50.60927108,\n", - " 45.68748656, 45.30077914]])" - ] - }, - "execution_count": 16, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ - "interpolated_source_grid.variables['sconco3']['data']" + "interpolated_source_grid.shapefile['sconco3'] = interpolated_source_grid.variables['sconco3']['data'][:, 0].ravel()" ] }, { "cell_type": "code", - "execution_count": 17, + "execution_count": 19, + "metadata": {}, + "outputs": [], + "source": [ + "interpolated_source_grid.write_shapefile('interpolated_points')" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## 3. Compare outputs" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### From NES" + ] + }, + { + "cell_type": "code", + "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "43.88671143497364" + "" ] }, - "execution_count": 17, + "execution_count": 20, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "interpolated_source_grid.variables['sconco3']['data'][100, 100]" + "interpolated_source_grid" ] }, { "cell_type": "code", - "execution_count": 18, + "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "51.778056" + "array([[53.66018016, 53.83203477, 53.63669104, ..., 53.2461802 ,\n", + " 53.21874343, 53.10930358],\n", + " [53.72655496, 53.84765007, 53.58983273, ..., 53.39846363,\n", + " 53.19921385, 52.96872593],\n", + " [39.71929631, 40.01980784, 39.58707525, ..., 54.41383702,\n", + " 53.30565688, 51.58627296],\n", + " ...,\n", + " [50.31659288, 50.36346616, 50.3400049 , ..., 50.5349958 ,\n", + " 50.66003373, 50.60927108],\n", + " [38.61718492, 39.64453277, 40.28125328, ..., 47.66795862,\n", + " 46.42186248, 45.68748656],\n", + " [44.8632812 , 44.78906232, 44.64843684, ..., 45.16406056,\n", + " 45.16796676, 45.30077914]])" ] }, - "execution_count": 18, + "execution_count": 21, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "interpolated_source_grid.lat['data'][100]" + "interpolated_source_grid.variables['sconco3']['data']" ] }, { "cell_type": "code", - "execution_count": 19, + "execution_count": 24, "metadata": {}, "outputs": [ { @@ -711,7 +909,7 @@ }, { "cell_type": "code", - "execution_count": 20, + "execution_count": 25, "metadata": {}, "outputs": [ { @@ -1069,10 +1267,10 @@ " fill: currentColor;\n", "}\n", "
<xarray.Dataset>\n",
-       "Dimensions:                 (time: 720, station: 175, model_latitude: 211, model_longitude: 351, grid_edge: 1125)\n",
+       "Dimensions:                 (time: 720, station: 168, model_latitude: 211, model_longitude: 351, grid_edge: 1125)\n",
        "Coordinates:\n",
        "  * time                    (time) datetime64[ns] 2018-04-01 ... 2018-04-30T2...\n",
-       "  * station                 (station) float64 0.0 1.0 2.0 ... 172.0 173.0 174.0\n",
+       "  * station                 (station) float64 0.0 1.0 2.0 ... 165.0 166.0 167.0\n",
        "Dimensions without coordinates: model_latitude, model_longitude, grid_edge\n",
        "Data variables:\n",
        "    lat                     (station) float64 -64.24 -54.85 ... 21.57 -34.35\n",
@@ -1083,10 +1281,10 @@
        "    grid_edge_latitude      (grid_edge) float64 29.9 30.1 30.3 ... 29.9 29.9\n",
        "    sconco3                 (station, time) float64 ...\n",
        "Attributes:\n",
-       "    Conventions:  CF-1.7
  • Conventions :
    CF-1.7
    • " ], "text/plain": [ "\n", - "Dimensions: (time: 720, station: 175, model_latitude: 211, model_longitude: 351, grid_edge: 1125)\n", + "Dimensions: (time: 720, station: 168, model_latitude: 211, model_longitude: 351, grid_edge: 1125)\n", "Coordinates:\n", " * time (time) datetime64[ns] 2018-04-01 ... 2018-04-30T2...\n", - " * station (station) float64 0.0 1.0 2.0 ... 172.0 173.0 174.0\n", + " * station (station) float64 0.0 1.0 2.0 ... 165.0 166.0 167.0\n", "Dimensions without coordinates: model_latitude, model_longitude, grid_edge\n", "Data variables:\n", " lat (station) float64 ...\n", @@ -1205,7 +1398,7 @@ " Conventions: CF-1.7" ] }, - "execution_count": 20, + "execution_count": 25, "metadata": {}, "output_type": "execute_result" } @@ -1223,7 +1416,7 @@ }, { "cell_type": "code", - "execution_count": 21, + "execution_count": 26, "metadata": {}, "outputs": [], "source": [ @@ -1232,7 +1425,7 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 27, "metadata": {}, "outputs": [ { @@ -1246,10 +1439,10 @@ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 22, + "execution_count": 27, "metadata": {}, "output_type": "execute_result" } @@ -1261,7 +1454,7 @@ }, { "cell_type": "code", - "execution_count": 23, + "execution_count": 28, "metadata": {}, "outputs": [ { @@ -1281,7 +1474,7 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 29, "metadata": {}, "outputs": [ { @@ -1300,7 +1493,7 @@ " dtype=float32)" ] }, - "execution_count": 24, + "execution_count": 29, "metadata": {}, "output_type": "execute_result" } @@ -1311,599 +1504,131 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 32, "metadata": {}, "outputs": [ { "data": { + "text/html": [ + "
    \n", + "\n", + "\n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + "
    geometry
    FID
    0POINT (-56.62478 -64.24006)
    1POINT (-68.31069 -54.84846)
    2POINT (-65.60083 -22.10333)
    3POINT (-63.88194 -31.66861)
    4POINT (16.76667 47.76667)
    ......
    170POINT (-155.57616 19.53623)
    171POINT (-24.80000 -89.99695)
    172POINT (-124.15100 41.05410)
    173POINT (103.51570 21.57310)
    174POINT (18.48968 -34.35348)
    \n", + "

    175 rows × 1 columns

    \n", + "
    " + ], "text/plain": [ - "38.09361" + " geometry\n", + "FID \n", + "0 POINT (-56.62478 -64.24006)\n", + "1 POINT (-68.31069 -54.84846)\n", + "2 POINT (-65.60083 -22.10333)\n", + "3 POINT (-63.88194 -31.66861)\n", + "4 POINT (16.76667 47.76667)\n", + ".. ...\n", + "170 POINT (-155.57616 19.53623)\n", + "171 POINT (-24.80000 -89.99695)\n", + "172 POINT (-124.15100 41.05410)\n", + "173 POINT (103.51570 21.57310)\n", + "174 POINT (18.48968 -34.35348)\n", + "\n", + "[175 rows x 1 columns]" ] }, - "execution_count": 25, + "execution_count": 32, "metadata": {}, "output_type": "execute_result" } ], "source": [ - "tool_interpolated_grid.variables['sconco3']['data'][100, 100]" + "tool_interpolated_grid.create_shapefile()\n", + "tool_interpolated_grid.shapefile" ] }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 33, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "51.778056" - ] - }, - "execution_count": 26, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ - "tool_interpolated_grid.lat['data'][100]" + "tool_interpolated_grid.shapefile['sconco3'] = tool_interpolated_grid.variables['sconco3']['data'][:, 0].ravel()" ] }, { "cell_type": "code", - "execution_count": 27, + "execution_count": 34, "metadata": {}, - "outputs": [ - { - "data": { - "text/html": [ - "
    \n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "\n", - "
    <xarray.Dataset>\n",
-       "Dimensions:                 (grid_edge: 1125, station: 175, model_latitude: 211, model_longitude: 351, time: 720)\n",
-       "Coordinates:\n",
-       "  * time                    (time) datetime64[ns] 2018-04-01 ... 2018-04-30T2...\n",
-       "Dimensions without coordinates: grid_edge, station, model_latitude, model_longitude\n",
-       "Data variables:\n",
-       "    grid_edge_latitude      (grid_edge) float64 29.9 30.1 30.3 ... 29.9 29.9\n",
-       "    grid_edge_longitude     (grid_edge) float64 -25.1 -25.1 ... -24.9 -25.1\n",
-       "    latitude                (station) float64 -64.24 -54.85 ... 21.57 -34.35\n",
-       "    longitude               (station) float64 -56.62 -68.31 ... 103.5 18.49\n",
-       "    model_centre_latitude   (model_latitude, model_longitude) float64 30.0 .....\n",
-       "    model_centre_longitude  (model_latitude, model_longitude) float64 -25.0 ....\n",
-       "    sconco3                 (station, time) float32 ...\n",
-       "    station_reference       (station) object 'AR0001R_UVP' ... 'ZA0001G_UVP'\n",
-       "Attributes:\n",
-       "    title:          Inverse distance weighting (4 neighbours) interpolated ca...\n",
-       "    institution:    Barcelona Supercomputing Center\n",
-       "    source:         Experiment cams61_chimere_ph2\n",
-       "    creator_name:   Dene R. Bowdalo\n",
-       "    creator_email:  dene.bowdalo@bsc.es\n",
-       "    conventions:    CF-1.7\n",
-       "    data_version:   1.0\n",
-       "    history:        Thu Feb 11 10:19:01 2021: ncks -O --dfl_lvl 1 /gpfs/proje...\n",
-       "    NCO:            4.7.2
    " - ], - "text/plain": [ - "\n", - "Dimensions: (grid_edge: 1125, station: 175, model_latitude: 211, model_longitude: 351, time: 720)\n", - "Coordinates:\n", - " * time (time) datetime64[ns] 2018-04-01 ... 2018-04-30T2...\n", - "Dimensions without coordinates: grid_edge, station, model_latitude, model_longitude\n", - "Data variables:\n", - " grid_edge_latitude (grid_edge) float64 ...\n", - " grid_edge_longitude (grid_edge) float64 ...\n", - " latitude (station) float64 ...\n", - " longitude (station) float64 ...\n", - " model_centre_latitude (model_latitude, model_longitude) float64 ...\n", - " model_centre_longitude (model_latitude, model_longitude) float64 ...\n", - " sconco3 (station, time) float32 ...\n", - " station_reference (station) object ...\n", - "Attributes:\n", - " title: Inverse distance weighting (4 neighbours) interpolated ca...\n", - " institution: Barcelona Supercomputing Center\n", - " source: Experiment cams61_chimere_ph2\n", - " creator_name: Dene R. Bowdalo\n", - " creator_email: dene.bowdalo@bsc.es\n", - " conventions: CF-1.7\n", - " data_version: 1.0\n", - " history: Thu Feb 11 10:19:01 2021: ncks -O --dfl_lvl 1 /gpfs/proje...\n", - " NCO: 4.7.2" - ] - }, - "execution_count": 27, - "metadata": {}, - "output_type": "execute_result" - } - ], + "outputs": [], "source": [ - "xr.open_dataset('/gpfs/projects/bsc32/AC_cache/recon/exp_interp/1.3.3/cams61_chimere_ph2-eu-000/hourly/sconco3/EBAS/sconco3_201804.nc')" + "tool_interpolated_grid.write_shapefile('providentia_it_points')" ] } ], -- GitLab