{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Using xesmf to efficiently regrid data to another resolution\n", "**Annette Hirsch, CLEX CMS**\n", "\n", "Often we want to compare to datasets but find that they are on different grids and first require to be interpolated to a common grid. \n", "\n", "There exist a number of different tools for doing this however in this post we'll talk about the regridding tools from xesmf." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "First load the relevant packages" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import xesmf as xe\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Open the data. In this example, we will open a monthly surface temperature dataset from the ACCESS-1.3 CMIP5 model." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "url = 'http://dapds00.nci.org.au/thredds/dodsC/rr3/CMIP5/output1/CSIRO-BOM/ACCESS1-3/historical/mon/atmos/Amon/r1i1p1/latest/tas/tas_Amon_ACCESS1-3_historical_r1i1p1_185001-200512.nc'\n", "ds = xr.open_dataset(url)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n", "Dimensions: (bnds: 2, lat: 145, lon: 192, time: 1872)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1850-01-16T12:00:00 ... 2005-12-16T12:00:00\n", " * lat (lat) float64 -90.0 -88.75 -87.5 -86.25 ... 86.25 87.5 88.75 90.0\n", " * lon (lon) float64 0.0 1.875 3.75 5.625 ... 352.5 354.4 356.2 358.1\n", " height float64 1.5\n", "Dimensions without coordinates: bnds\n", "Data variables:\n", " time_bnds (time, bnds) datetime64[ns] 1850-01-01 1850-02-01 ... 2006-01-01\n", " lat_bnds (lat, bnds) float64 -90.0 -89.38 -89.38 ... 89.38 89.38 90.0\n", " lon_bnds (lon, bnds) float64 -0.9375 0.9375 0.9375 ... 357.2 357.2 359.1\n", " tas (time, lat, lon) float32 ...\n", "Attributes: (12/29)\n", " institution: CSIRO (Commonwealth Scientific and Indus...\n", " institute_id: CSIRO-BOM\n", " experiment_id: historical\n", " source: ACCESS1-3 2011. Atmosphere: AGCM v1.0 (N...\n", " model_id: ACCESS1.3\n", " forcing: GHG, Oz, SA, Sl, Vl, BC, OC, (GHG = CO2,...\n", " ... ...\n", " title: ACCESS1-3 model output prepared for CMIP...\n", " parent_experiment: pre-industrial control\n", " modeling_realm: atmos\n", " realization: 1\n", " cmor_version: 2.8.0\n", " DODS_EXTRA.Unlimited_Dimension: time
array(['1850-01-16T12:00:00.000000000', '1850-02-15T00:00:00.000000000',\n", " '1850-03-16T12:00:00.000000000', ..., '2005-10-16T12:00:00.000000000',\n", " '2005-11-16T00:00:00.000000000', '2005-12-16T12:00:00.000000000'],\n", " dtype='datetime64[ns]')
array([-90. , -88.75, -87.5 , -86.25, -85. , -83.75, -82.5 , -81.25, -80. ,\n", " -78.75, -77.5 , -76.25, -75. , -73.75, -72.5 , -71.25, -70. , -68.75,\n", " -67.5 , -66.25, -65. , -63.75, -62.5 , -61.25, -60. , -58.75, -57.5 ,\n", " -56.25, -55. , -53.75, -52.5 , -51.25, -50. , -48.75, -47.5 , -46.25,\n", " -45. , -43.75, -42.5 , -41.25, -40. , -38.75, -37.5 , -36.25, -35. ,\n", " -33.75, -32.5 , -31.25, -30. , -28.75, -27.5 , -26.25, -25. , -23.75,\n", " -22.5 , -21.25, -20. , -18.75, -17.5 , -16.25, -15. , -13.75, -12.5 ,\n", " -11.25, -10. , -8.75, -7.5 , -6.25, -5. , -3.75, -2.5 , -1.25,\n", " 0. , 1.25, 2.5 , 3.75, 5. , 6.25, 7.5 , 8.75, 10. ,\n", " 11.25, 12.5 , 13.75, 15. , 16.25, 17.5 , 18.75, 20. , 21.25,\n", " 22.5 , 23.75, 25. , 26.25, 27.5 , 28.75, 30. , 31.25, 32.5 ,\n", " 33.75, 35. , 36.25, 37.5 , 38.75, 40. , 41.25, 42.5 , 43.75,\n", " 45. , 46.25, 47.5 , 48.75, 50. , 51.25, 52.5 , 53.75, 55. ,\n", " 56.25, 57.5 , 58.75, 60. , 61.25, 62.5 , 63.75, 65. , 66.25,\n", " 67.5 , 68.75, 70. , 71.25, 72.5 , 73.75, 75. , 76.25, 77.5 ,\n", " 78.75, 80. , 81.25, 82.5 , 83.75, 85. , 86.25, 87.5 , 88.75,\n", " 90. ])
array([ 0. , 1.875, 3.75 , 5.625, 7.5 , 9.375, 11.25 , 13.125,\n", " 15. , 16.875, 18.75 , 20.625, 22.5 , 24.375, 26.25 , 28.125,\n", " 30. , 31.875, 33.75 , 35.625, 37.5 , 39.375, 41.25 , 43.125,\n", " 45. , 46.875, 48.75 , 50.625, 52.5 , 54.375, 56.25 , 58.125,\n", " 60. , 61.875, 63.75 , 65.625, 67.5 , 69.375, 71.25 , 73.125,\n", " 75. , 76.875, 78.75 , 80.625, 82.5 , 84.375, 86.25 , 88.125,\n", " 90. , 91.875, 93.75 , 95.625, 97.5 , 99.375, 101.25 , 103.125,\n", " 105. , 106.875, 108.75 , 110.625, 112.5 , 114.375, 116.25 , 118.125,\n", " 120. , 121.875, 123.75 , 125.625, 127.5 , 129.375, 131.25 , 133.125,\n", " 135. , 136.875, 138.75 , 140.625, 142.5 , 144.375, 146.25 , 148.125,\n", " 150. , 151.875, 153.75 , 155.625, 157.5 , 159.375, 161.25 , 163.125,\n", " 165. , 166.875, 168.75 , 170.625, 172.5 , 174.375, 176.25 , 178.125,\n", " 180. , 181.875, 183.75 , 185.625, 187.5 , 189.375, 191.25 , 193.125,\n", " 195. , 196.875, 198.75 , 200.625, 202.5 , 204.375, 206.25 , 208.125,\n", " 210. , 211.875, 213.75 , 215.625, 217.5 , 219.375, 221.25 , 223.125,\n", " 225. , 226.875, 228.75 , 230.625, 232.5 , 234.375, 236.25 , 238.125,\n", " 240. , 241.875, 243.75 , 245.625, 247.5 , 249.375, 251.25 , 253.125,\n", " 255. , 256.875, 258.75 , 260.625, 262.5 , 264.375, 266.25 , 268.125,\n", " 270. , 271.875, 273.75 , 275.625, 277.5 , 279.375, 281.25 , 283.125,\n", " 285. , 286.875, 288.75 , 290.625, 292.5 , 294.375, 296.25 , 298.125,\n", " 300. , 301.875, 303.75 , 305.625, 307.5 , 309.375, 311.25 , 313.125,\n", " 315. , 316.875, 318.75 , 320.625, 322.5 , 324.375, 326.25 , 328.125,\n", " 330. , 331.875, 333.75 , 335.625, 337.5 , 339.375, 341.25 , 343.125,\n", " 345. , 346.875, 348.75 , 350.625, 352.5 , 354.375, 356.25 , 358.125])
array(1.5)
array([['1850-01-01T00:00:00.000000000', '1850-02-01T00:00:00.000000000'],\n", " ['1850-02-01T00:00:00.000000000', '1850-03-01T00:00:00.000000000'],\n", " ['1850-03-01T00:00:00.000000000', '1850-04-01T00:00:00.000000000'],\n", " ...,\n", " ['2005-10-01T00:00:00.000000000', '2005-11-01T00:00:00.000000000'],\n", " ['2005-11-01T00:00:00.000000000', '2005-12-01T00:00:00.000000000'],\n", " ['2005-12-01T00:00:00.000000000', '2006-01-01T00:00:00.000000000']],\n", " dtype='datetime64[ns]')
array([[-90. , -89.375],\n", " [-89.375, -88.125],\n", " [-88.125, -86.875],\n", " ...,\n", " [ 86.875, 88.125],\n", " [ 88.125, 89.375],\n", " [ 89.375, 90. ]])
array([[ -0.9375, 0.9375],\n", " [ 0.9375, 2.8125],\n", " [ 2.8125, 4.6875],\n", " ...,\n", " [353.4375, 355.3125],\n", " [355.3125, 357.1875],\n", " [357.1875, 359.0625]])
[52116480 values with dtype=float32]
<xarray.DataArray 'lat' ()>\n", "array(1.25)\n", "Coordinates:\n", " height float64 1.5
array(1.25)
array(1.5)
<xarray.DataArray 'lon' ()>\n", "array(1.875)\n", "Coordinates:\n", " height float64 1.5
array(1.875)
array(1.5)
<xarray.Dataset>\n", "Dimensions: (x: 144, x_b: 145, y: 72, y_b: 73)\n", "Coordinates:\n", " lon (y, x) float64 -178.8 -176.2 -173.8 -171.2 ... 173.8 176.2 178.8\n", " lat (y, x) float64 -88.75 -88.75 -88.75 -88.75 ... 88.75 88.75 88.75\n", " lon_b (y_b, x_b) float64 -180.0 -177.5 -175.0 ... 175.0 177.5 180.0\n", " lat_b (y_b, x_b) float64 -90.0 -90.0 -90.0 -90.0 ... 90.0 90.0 90.0 90.0\n", "Dimensions without coordinates: x, x_b, y, y_b\n", "Data variables:\n", " *empty*
array([[-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " ...,\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75]])
array([[-88.75, -88.75, -88.75, ..., -88.75, -88.75, -88.75],\n", " [-86.25, -86.25, -86.25, ..., -86.25, -86.25, -86.25],\n", " [-83.75, -83.75, -83.75, ..., -83.75, -83.75, -83.75],\n", " ...,\n", " [ 83.75, 83.75, 83.75, ..., 83.75, 83.75, 83.75],\n", " [ 86.25, 86.25, 86.25, ..., 86.25, 86.25, 86.25],\n", " [ 88.75, 88.75, 88.75, ..., 88.75, 88.75, 88.75]])
array([[-180. , -177.5, -175. , ..., 175. , 177.5, 180. ],\n", " [-180. , -177.5, -175. , ..., 175. , 177.5, 180. ],\n", " [-180. , -177.5, -175. , ..., 175. , 177.5, 180. ],\n", " ...,\n", " [-180. , -177.5, -175. , ..., 175. , 177.5, 180. ],\n", " [-180. , -177.5, -175. , ..., 175. , 177.5, 180. ],\n", " [-180. , -177.5, -175. , ..., 175. , 177.5, 180. ]])
array([[-90. , -90. , -90. , ..., -90. , -90. , -90. ],\n", " [-87.5, -87.5, -87.5, ..., -87.5, -87.5, -87.5],\n", " [-85. , -85. , -85. , ..., -85. , -85. , -85. ],\n", " ...,\n", " [ 85. , 85. , 85. , ..., 85. , 85. , 85. ],\n", " [ 87.5, 87.5, 87.5, ..., 87.5, 87.5, 87.5],\n", " [ 90. , 90. , 90. , ..., 90. , 90. , 90. ]])
<xarray.DataArray 'tas' (time: 1872, lat: 145, lon: 192)>\n", "[52116480 values with dtype=float32]\n", "Coordinates:\n", " * time (time) datetime64[ns] 1850-01-16T12:00:00 ... 2005-12-16T12:00:00\n", " * lat (lat) float64 -90.0 -88.75 -87.5 -86.25 ... 86.25 87.5 88.75 90.0\n", " * lon (lon) float64 0.0 1.875 3.75 5.625 7.5 ... 352.5 354.4 356.2 358.1\n", " height float64 1.5\n", "Attributes:\n", " standard_name: air_temperature\n", " long_name: Near-Surface Air Temperature\n", " units: K\n", " cell_methods: time: mean\n", " cell_measures: area: areacella\n", " history: 2012-02-05T23:49:51Z altered by CMOR: Treated scalar d...\n", " associated_files: baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation...
[52116480 values with dtype=float32]
array(['1850-01-16T12:00:00.000000000', '1850-02-15T00:00:00.000000000',\n", " '1850-03-16T12:00:00.000000000', ..., '2005-10-16T12:00:00.000000000',\n", " '2005-11-16T00:00:00.000000000', '2005-12-16T12:00:00.000000000'],\n", " dtype='datetime64[ns]')
array([-90. , -88.75, -87.5 , -86.25, -85. , -83.75, -82.5 , -81.25, -80. ,\n", " -78.75, -77.5 , -76.25, -75. , -73.75, -72.5 , -71.25, -70. , -68.75,\n", " -67.5 , -66.25, -65. , -63.75, -62.5 , -61.25, -60. , -58.75, -57.5 ,\n", " -56.25, -55. , -53.75, -52.5 , -51.25, -50. , -48.75, -47.5 , -46.25,\n", " -45. , -43.75, -42.5 , -41.25, -40. , -38.75, -37.5 , -36.25, -35. ,\n", " -33.75, -32.5 , -31.25, -30. , -28.75, -27.5 , -26.25, -25. , -23.75,\n", " -22.5 , -21.25, -20. , -18.75, -17.5 , -16.25, -15. , -13.75, -12.5 ,\n", " -11.25, -10. , -8.75, -7.5 , -6.25, -5. , -3.75, -2.5 , -1.25,\n", " 0. , 1.25, 2.5 , 3.75, 5. , 6.25, 7.5 , 8.75, 10. ,\n", " 11.25, 12.5 , 13.75, 15. , 16.25, 17.5 , 18.75, 20. , 21.25,\n", " 22.5 , 23.75, 25. , 26.25, 27.5 , 28.75, 30. , 31.25, 32.5 ,\n", " 33.75, 35. , 36.25, 37.5 , 38.75, 40. , 41.25, 42.5 , 43.75,\n", " 45. , 46.25, 47.5 , 48.75, 50. , 51.25, 52.5 , 53.75, 55. ,\n", " 56.25, 57.5 , 58.75, 60. , 61.25, 62.5 , 63.75, 65. , 66.25,\n", " 67.5 , 68.75, 70. , 71.25, 72.5 , 73.75, 75. , 76.25, 77.5 ,\n", " 78.75, 80. , 81.25, 82.5 , 83.75, 85. , 86.25, 87.5 , 88.75,\n", " 90. ])
array([ 0. , 1.875, 3.75 , 5.625, 7.5 , 9.375, 11.25 , 13.125,\n", " 15. , 16.875, 18.75 , 20.625, 22.5 , 24.375, 26.25 , 28.125,\n", " 30. , 31.875, 33.75 , 35.625, 37.5 , 39.375, 41.25 , 43.125,\n", " 45. , 46.875, 48.75 , 50.625, 52.5 , 54.375, 56.25 , 58.125,\n", " 60. , 61.875, 63.75 , 65.625, 67.5 , 69.375, 71.25 , 73.125,\n", " 75. , 76.875, 78.75 , 80.625, 82.5 , 84.375, 86.25 , 88.125,\n", " 90. , 91.875, 93.75 , 95.625, 97.5 , 99.375, 101.25 , 103.125,\n", " 105. , 106.875, 108.75 , 110.625, 112.5 , 114.375, 116.25 , 118.125,\n", " 120. , 121.875, 123.75 , 125.625, 127.5 , 129.375, 131.25 , 133.125,\n", " 135. , 136.875, 138.75 , 140.625, 142.5 , 144.375, 146.25 , 148.125,\n", " 150. , 151.875, 153.75 , 155.625, 157.5 , 159.375, 161.25 , 163.125,\n", " 165. , 166.875, 168.75 , 170.625, 172.5 , 174.375, 176.25 , 178.125,\n", " 180. , 181.875, 183.75 , 185.625, 187.5 , 189.375, 191.25 , 193.125,\n", " 195. , 196.875, 198.75 , 200.625, 202.5 , 204.375, 206.25 , 208.125,\n", " 210. , 211.875, 213.75 , 215.625, 217.5 , 219.375, 221.25 , 223.125,\n", " 225. , 226.875, 228.75 , 230.625, 232.5 , 234.375, 236.25 , 238.125,\n", " 240. , 241.875, 243.75 , 245.625, 247.5 , 249.375, 251.25 , 253.125,\n", " 255. , 256.875, 258.75 , 260.625, 262.5 , 264.375, 266.25 , 268.125,\n", " 270. , 271.875, 273.75 , 275.625, 277.5 , 279.375, 281.25 , 283.125,\n", " 285. , 286.875, 288.75 , 290.625, 292.5 , 294.375, 296.25 , 298.125,\n", " 300. , 301.875, 303.75 , 305.625, 307.5 , 309.375, 311.25 , 313.125,\n", " 315. , 316.875, 318.75 , 320.625, 322.5 , 324.375, 326.25 , 328.125,\n", " 330. , 331.875, 333.75 , 335.625, 337.5 , 339.375, 341.25 , 343.125,\n", " 345. , 346.875, 348.75 , 350.625, 352.5 , 354.375, 356.25 , 358.125])
array(1.5)
<xarray.DataArray 'tas' (time: 1872, y: 72, x: 144)>\n", "array([[[241.42548074, 241.46530151, 241.51575264, ..., 241.29264361,\n", " 241.3331604 , 241.37937428],\n", " [249.29878775, 249.49865723, 249.70174442, ..., 248.53185656,\n", " 248.81828308, 249.07307976],\n", " [258.90998768, 259.34362793, 259.62011489, ..., 257.00554583,\n", " 257.77041626, 258.39012844],\n", " ...,\n", " [240.0256449 , 240.15310669, 240.28757747, ..., 239.80319719,\n", " 239.92158508, 239.9827772 ],\n", " [239.95589676, 239.94976807, 239.9454768 , ..., 239.91784213,\n", " 239.93392944, 239.96497476],\n", " [239.91920767, 239.95617676, 239.98541397, ..., 239.76462219,\n", " 239.8212738 , 239.87550099]],\n", "\n", " [[235.35084386, 235.43763733, 235.53132127, ..., 235.12648029,\n", " 235.19128418, 235.26917861],\n", " [247.80141914, 248.36372375, 248.88822159, ..., 246.27575141,\n", " 246.73565674, 247.26041022],\n", " [257.63035803, 257.96582031, 258.17488678, ..., 256.06124283,\n", " 256.6065979 , 257.11943182],\n", "...\n", " [250.55788953, 250.56117249, 250.62096828, ..., 251.29435661,\n", " 250.8780365 , 250.63096473],\n", " [252.6061043 , 252.39373779, 252.23463555, ..., 253.47047825,\n", " 253.12278748, 252.8264588 ],\n", " [254.88175595, 254.8631897 , 254.86619789, ..., 255.01445768,\n", " 254.97032166, 254.92613684]],\n", "\n", " [[246.75074031, 246.7674408 , 246.81818953, ..., 246.66889692,\n", " 246.68983459, 246.7283168 ],\n", " [253.50811955, 253.73731995, 253.98962857, ..., 252.68888363,\n", " 252.96813965, 253.25622345],\n", " [262.74184175, 263.19360352, 263.53787576, ..., 261.00175357,\n", " 261.69522095, 262.25131706],\n", " ...,\n", " [245.47192229, 245.0090332 , 244.5782707 , ..., 246.80896796,\n", " 246.32383728, 245.88781399],\n", " [246.79688785, 246.40936279, 246.01888102, ..., 247.81054793,\n", " 247.51338196, 247.15659795],\n", " [245.98595642, 245.89811707, 245.81729825, ..., 246.27256695,\n", " 246.18641663, 246.08195717]]])\n", "Coordinates:\n", " * time (time) datetime64[ns] 1850-01-16T12:00:00 ... 2005-12-16T12:00:00\n", " height float64 1.5\n", " lon (y, x) float64 -178.8 -176.2 -173.8 -171.2 ... 173.8 176.2 178.8\n", " lat (y, x) float64 -88.75 -88.75 -88.75 -88.75 ... 88.75 88.75 88.75\n", "Dimensions without coordinates: y, x\n", "Attributes:\n", " regrid_method: bilinear
array([[[241.42548074, 241.46530151, 241.51575264, ..., 241.29264361,\n", " 241.3331604 , 241.37937428],\n", " [249.29878775, 249.49865723, 249.70174442, ..., 248.53185656,\n", " 248.81828308, 249.07307976],\n", " [258.90998768, 259.34362793, 259.62011489, ..., 257.00554583,\n", " 257.77041626, 258.39012844],\n", " ...,\n", " [240.0256449 , 240.15310669, 240.28757747, ..., 239.80319719,\n", " 239.92158508, 239.9827772 ],\n", " [239.95589676, 239.94976807, 239.9454768 , ..., 239.91784213,\n", " 239.93392944, 239.96497476],\n", " [239.91920767, 239.95617676, 239.98541397, ..., 239.76462219,\n", " 239.8212738 , 239.87550099]],\n", "\n", " [[235.35084386, 235.43763733, 235.53132127, ..., 235.12648029,\n", " 235.19128418, 235.26917861],\n", " [247.80141914, 248.36372375, 248.88822159, ..., 246.27575141,\n", " 246.73565674, 247.26041022],\n", " [257.63035803, 257.96582031, 258.17488678, ..., 256.06124283,\n", " 256.6065979 , 257.11943182],\n", "...\n", " [250.55788953, 250.56117249, 250.62096828, ..., 251.29435661,\n", " 250.8780365 , 250.63096473],\n", " [252.6061043 , 252.39373779, 252.23463555, ..., 253.47047825,\n", " 253.12278748, 252.8264588 ],\n", " [254.88175595, 254.8631897 , 254.86619789, ..., 255.01445768,\n", " 254.97032166, 254.92613684]],\n", "\n", " [[246.75074031, 246.7674408 , 246.81818953, ..., 246.66889692,\n", " 246.68983459, 246.7283168 ],\n", " [253.50811955, 253.73731995, 253.98962857, ..., 252.68888363,\n", " 252.96813965, 253.25622345],\n", " [262.74184175, 263.19360352, 263.53787576, ..., 261.00175357,\n", " 261.69522095, 262.25131706],\n", " ...,\n", " [245.47192229, 245.0090332 , 244.5782707 , ..., 246.80896796,\n", " 246.32383728, 245.88781399],\n", " [246.79688785, 246.40936279, 246.01888102, ..., 247.81054793,\n", " 247.51338196, 247.15659795],\n", " [245.98595642, 245.89811707, 245.81729825, ..., 246.27256695,\n", " 246.18641663, 246.08195717]]])
array(['1850-01-16T12:00:00.000000000', '1850-02-15T00:00:00.000000000',\n", " '1850-03-16T12:00:00.000000000', ..., '2005-10-16T12:00:00.000000000',\n", " '2005-11-16T00:00:00.000000000', '2005-12-16T12:00:00.000000000'],\n", " dtype='datetime64[ns]')
array(1.5)
array([[-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " ...,\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75],\n", " [-178.75, -176.25, -173.75, ..., 173.75, 176.25, 178.75]])
array([[-88.75, -88.75, -88.75, ..., -88.75, -88.75, -88.75],\n", " [-86.25, -86.25, -86.25, ..., -86.25, -86.25, -86.25],\n", " [-83.75, -83.75, -83.75, ..., -83.75, -83.75, -83.75],\n", " ...,\n", " [ 83.75, 83.75, 83.75, ..., 83.75, 83.75, 83.75],\n", " [ 86.25, 86.25, 86.25, ..., 86.25, 86.25, 86.25],\n", " [ 88.75, 88.75, 88.75, ..., 88.75, 88.75, 88.75]])