{ "cells": [ { "attachments": {}, "cell_type": "markdown", "id": "13c9f12f", "metadata": {}, "source": [ "# 7. 進階運算與統計方法" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## `flox`: 進階`xarray.groupby` \n", "\n", "`xarray.groupby` 有一些限制:\n", "1. 只能針對單一變數分群。\n", "2. 若在平行運算的陣列 (dask array) 下,為了尋找各個值適合的組別,必須耗費較多計算資源。 \n", "\n", "`flox`套件則可以解決這個問題。以下是針對兩個變數分群的範例。`flox.xarray_reduce`的用法請參考[官方說明](https://flox.readthedocs.io/en/latest/generated/flox.xarray.xarray_reduce.html#flox.xarray.xarray_reduce)。\n", "\n", "**Example 1:** 將每日OLR資料轉化成(year, pentad, lat, lon)格式 (同第五章的範例)。" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/waynetsai/.local/lib/python3.10/site-packages/ecmwflibs/__init__.py:83: UserWarning: dlopen(/Users/waynetsai/.local/lib/python3.10/site-packages/ecmwflibs/_ecmwflibs.cpython-310-darwin.so, 0x0002): tried: '/Users/waynetsai/.local/lib/python3.10/site-packages/ecmwflibs/_ecmwflibs.cpython-310-darwin.so' (mach-o file, but is an incompatible architecture (have 'x86_64', need 'arm64')), '/System/Volumes/Preboot/Cryptexes/OS/Users/waynetsai/.local/lib/python3.10/site-packages/ecmwflibs/_ecmwflibs.cpython-310-darwin.so' (no such file), '/Users/waynetsai/.local/lib/python3.10/site-packages/ecmwflibs/_ecmwflibs.cpython-310-darwin.so' (mach-o file, but is an incompatible architecture (have 'x86_64', need 'arm64'))\n", " warnings.warn(str(e))\n" ] }, { "data": { "text/html": [ "
<xarray.DataArray 'olr' (year: 24, dayofyear_bins: 73, lat: 50, lon: 82)>\n", "array([[[[284.27383, 285.73615, 288.25266, ..., 231.88802, 241.6569 ,\n", " 249.84024],\n", " [286.65228, 288.05078, 288.27435, ..., 215.95383, 230.82944,\n", " 246.1481 ],\n", " [290.7838 , 290.25864, 287.82904, ..., 198.27682, 214.98921,\n", " 225.9622 ],\n", " ...,\n", " [247.29614, 246.06245, 247.94983, ..., 279.0298 , 277.99088,\n", " 276.57684],\n", " [244.55948, 242.81021, 239.5767 , ..., 273.2099 , 271.63373,\n", " 270.99365],\n", " [236.1903 , 229.96805, 221.43008, ..., 266.5147 , 266.81476,\n", " 266.7093 ]],\n", "\n", " [[257.097 , 249.06888, 260.7727 , ..., 251.24788, 258.4394 ,\n", " 267.57504],\n", " [265.49664, 263.90622, 270.4431 , ..., 237.29634, 248.2174 ,\n", " 250.10226],\n", " [261.7508 , 266.22888, 276.74783, ..., 222.44157, 235.64157,\n", " 227.53247],\n", "...\n", " [246.52588, 246.06047, 249.06099, ..., 256.54034, 259.21564,\n", " 260.4043 ],\n", " [244.2377 , 247.02632, 246.57187, ..., 258.28467, 255.90445,\n", " 257.71054],\n", " [237.547 , 235.2325 , 224.96394, ..., 255.28891, 253.66435,\n", " 255.53586]],\n", "\n", " [[286.86404, 287.83688, 290.30945, ..., 194.66504, 212.79057,\n", " 210.36751],\n", " [288.31 , 288.92453, 290.12537, ..., 204.64694, 225.9674 ,\n", " 223.37724],\n", " [284.80383, 288.20563, 289.7766 , ..., 215.14029, 236.13554,\n", " 232.43852],\n", " ...,\n", " [236.44717, 236.95541, 234.64787, ..., 268.50082, 268.91937,\n", " 267.2533 ],\n", " [235.04623, 231.87955, 229.08992, ..., 267.69647, 267.49542,\n", " 266.0851 ],\n", " [233.4418 , 225.93411, 211.10928, ..., 268.23108, 265.9177 ,\n", " 264.8466 ]]]], dtype=float32)\n", "Coordinates:\n", " * lon (lon) float32 79.5 80.5 81.5 82.5 ... 158.5 159.5 160.5\n", " * lat (lat) float32 -19.5 -18.5 -17.5 -16.5 ... 27.5 28.5 29.5\n", " * year (year) int64 1998 1999 2000 2001 ... 2018 2019 2020 2021\n", " * dayofyear_bins (dayofyear_bins) object [1, 6) [6, 11) ... [361, 366)\n", "Attributes:\n", " standard_name: toa_outgoing_longwave_flux\n", " long_name: NOAA Climate Data Record of Daily Mean Upward Longwave Fl...\n", " units: W m-2\n", " cell_methods: time: mean area: mean
<xarray.DataArray 'olr' (time: 8760, lat: 20, lon: 32)>\n", "array([[[210.37125, 199.29874, 197.44077, ..., 268.05853, 268.25916,\n", " 266.9823 ],\n", " [204.61313, 213.93433, 230.92719, ..., 284.24072, 282.60632,\n", " 278.35385],\n", " [236.69272, 247.62692, 259.35587, ..., 294.86554, 292.45862,\n", " 292.01804],\n", " ...,\n", " [297.27472, 302.5027 , 305.18634, ..., 193.11395, 206.42505,\n", " 197.59212],\n", " [291.31464, 291.11176, 298.43726, ..., 155.4393 , 146.58469,\n", " 162.1478 ],\n", " [286.89136, 290.6887 , 291.20618, ..., 171.58977, 183.45123,\n", " 195.39035]],\n", "\n", " [[259.8478 , 225.8486 , 211.98647, ..., 269.50348, 269.39746,\n", " 272.64 ],\n", " [276.39136, 272.1738 , 268.95856, ..., 291.94778, 289.79297,\n", " 284.37097],\n", " [260.77368, 265.89886, 275.41898, ..., 301.7583 , 299.66028,\n", " 292.9348 ],\n", "...\n", " [289.9028 , 281.8377 , 285.37466, ..., 169.81645, 173.29109,\n", " 213.98285],\n", " [286.0285 , 282.36426, 284.31454, ..., 111.96492, 154.03249,\n", " 224.03793],\n", " [290.05353, 295.02466, 291.87067, ..., 125.69296, 153.77583,\n", " 191.87067]],\n", "\n", " [[240.00652, 211.87251, 198.56546, ..., 264.23984, 258.3333 ,\n", " 255.25803],\n", " [224.23218, 243.02698, 246.38681, ..., 279.6334 , 274.43097,\n", " 277.21472],\n", " [234.04726, 253.98969, 263.0449 , ..., 289.64203, 291.71185,\n", " 293.91864],\n", " ...,\n", " [286.8454 , 287.29822, 284.08643, ..., 165.07278, 165.68568,\n", " 180.77968],\n", " [291.43097, 289.64923, 295.36426, ..., 178.16243, 158.17258,\n", " 182.52911],\n", " [294.91678, 297.96875, 296.73392, ..., 174.35318, 157.13422,\n", " 152.69089]]], dtype=float32)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1998-01-01 1998-01-02 ... 2021-12-31\n", " * lon (lon) float32 81.25 83.75 86.25 88.75 ... 151.2 153.8 156.2 158.8\n", " * lat (lat) float32 28.75 26.25 23.75 21.25 ... -13.75 -16.25 -18.75\n", "Attributes:\n", " standard_name: toa_outgoing_longwave_flux\n", " long_name: NOAA Climate Data Record of Daily Mean Upward Longwave Fl...\n", " units: W m-2\n", " cell_methods: time: mean area: mean
<xarray.DataArray 'olr' (time: 1752, lat: 50, lon: 82)>\n", "array([[[284.27383, 285.73615, 288.25266, ..., 231.88799, 241.6569 ,\n", " 249.84024],\n", " [286.65228, 288.0508 , 288.27435, ..., 215.95383, 230.82944,\n", " 246.1481 ],\n", " [290.7838 , 290.2586 , 287.82904, ..., 198.27681, 214.98921,\n", " 225.9622 ],\n", " ...,\n", " [247.29614, 246.06245, 247.9498 , ..., 279.0298 , 277.99088,\n", " 276.57684],\n", " [244.55948, 242.81021, 239.57669, ..., 273.2099 , 271.63373,\n", " 270.99365],\n", " [236.1903 , 229.96805, 221.43008, ..., 266.5147 , 266.81476,\n", " 266.7093 ]],\n", "\n", " [[257.097 , 249.0689 , 260.7727 , ..., 251.24785, 258.4394 ,\n", " 267.57504],\n", " [265.49664, 263.9062 , 270.4431 , ..., 237.29636, 248.2174 ,\n", " 250.10226],\n", " [261.7508 , 266.22888, 276.74783, ..., 222.4416 , 235.64157,\n", " 227.53247],\n", "...\n", " [246.52588, 246.06046, 249.06099, ..., 256.54034, 259.21564,\n", " 260.40427],\n", " [244.2377 , 247.02632, 246.57185, ..., 258.28464, 255.90445,\n", " 257.71054],\n", " [237.54703, 235.23251, 224.96394, ..., 255.28891, 253.66435,\n", " 255.53586]],\n", "\n", " [[286.864 , 287.83688, 290.30945, ..., 194.66502, 212.79057,\n", " 210.36748],\n", " [288.31 , 288.92453, 290.12537, ..., 204.64694, 225.9674 ,\n", " 223.37724],\n", " [284.80383, 288.20563, 289.7766 , ..., 215.14029, 236.13554,\n", " 232.43852],\n", " ...,\n", " [236.44717, 236.95544, 234.64787, ..., 268.50082, 268.91937,\n", " 267.25333],\n", " [235.04623, 231.87952, 229.08992, ..., 267.69647, 267.49542,\n", " 266.0851 ],\n", " [233.44182, 225.93411, 211.10928, ..., 268.23108, 265.9177 ,\n", " 264.84662]]], dtype=float32)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1998-01-01 1998-01-06 ... 2021-12-27\n", " * lon (lon) float32 79.5 80.5 81.5 82.5 83.5 ... 157.5 158.5 159.5 160.5\n", " * lat (lat) float32 -19.5 -18.5 -17.5 -16.5 -15.5 ... 26.5 27.5 28.5 29.5\n", "Attributes:\n", " standard_name: toa_outgoing_longwave_flux\n", " long_name: NOAA Climate Data Record of Daily Mean Upward Longwave Fl...\n", " units: W m-2\n", " cell_methods: time: mean area: mean
<xarray.DataArray 'olr' (time: 1752, lat: 50, lon: 82)>\n", "array([[[ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan]],\n", "\n", " [[ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", "...\n", " [258.01865, 258.54132, 260.59457, ..., 265.62262, 266.46997,\n", " 267.72852],\n", " [254.52092, 255.43341, 253.74515, ..., 262.07703, 262.12787,\n", " 264.36252],\n", " [249.05486, 243.50095, 230.53357, ..., 257.8312 , 259.1329 ,\n", " 259.55786]],\n", "\n", " [[285.45612, 286.57214, 287.6652 , ..., 247.677 , 257.28473,\n", " 259.3609 ],\n", " [288.13184, 288.3128 , 289.0145 , ..., 255.27068, 264.8253 ,\n", " 264.77902],\n", " [288.0848 , 289.45697, 289.93057, ..., 261.4938 , 269.81818,\n", " 268.2691 ],\n", " ...,\n", " [248.90976, 249.10083, 250.33572, ..., 262.9018 , 263.4146 ,\n", " 263.58185],\n", " [247.00436, 247.02068, 245.05011, ..., 261.76514, 261.24374,\n", " 261.8905 ],\n", " [243.24814, 237.90851, 224.16135, ..., 258.3669 , 257.77066,\n", " 257.97882]]], dtype=float32)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1998-01-01 1998-01-06 ... 2021-12-27\n", " * lon (lon) float32 79.5 80.5 81.5 82.5 83.5 ... 157.5 158.5 159.5 160.5\n", " * lat (lat) float32 -19.5 -18.5 -17.5 -16.5 -15.5 ... 26.5 27.5 28.5 29.5\n", "Attributes:\n", " standard_name: toa_outgoing_longwave_flux\n", " long_name: NOAA Climate Data Record of Daily Mean Upward Longwave Fl...\n", " units: W m-2\n", " cell_methods: time: mean area: mean
<xarray.DataArray 'data' (year: 41, lat: 20, lon: 32)>\n", "array([[[1.0133333e+00, 9.1833329e-01, 7.9833335e-01, ...,\n", " 3.6333334e+00, 2.7116668e+00, 2.5383332e+00],\n", " [2.0333333e-01, 4.2833331e-01, 2.5000000e-01, ...,\n", " 2.4483333e+00, 2.3216665e+00, 1.7766665e+00],\n", " [5.4999996e-02, 1.3666667e-01, 1.5833335e-01, ...,\n", " 2.2033336e+00, 1.8683332e+00, 1.2850000e+00],\n", " ...,\n", " [1.1466666e+00, 1.4949999e+00, 1.8333334e+00, ...,\n", " 3.6483333e+00, 4.2433333e+00, 2.9916668e+00],\n", " [7.1833330e-01, 1.3866667e+00, 1.4466667e+00, ...,\n", " 4.0750003e+00, 3.2466669e+00, 2.8599999e+00],\n", " [3.9666665e-01, 8.2333332e-01, 1.1883334e+00, ...,\n", " 2.9016666e+00, 2.7083333e+00, 1.7416667e+00]],\n", "\n", " [[4.6500000e-01, 2.8333333e-01, 3.4333333e-01, ...,\n", " 3.8716667e+00, 3.8799999e+00, 3.4633334e+00],\n", " [6.0833335e-01, 1.7000000e-01, 1.4666666e-01, ...,\n", " 3.9750001e+00, 3.6300001e+00, 3.1483333e+00],\n", " [1.2800001e+00, 3.2166669e-01, 1.2000000e-01, ...,\n", " 2.4466667e+00, 2.1550000e+00, 1.1616668e+00],\n", "...\n", " [7.9235015e+00, 7.0794005e+00, 5.7172127e+00, ...,\n", " 9.7668800e+00, 1.0939357e+01, 1.0846936e+01],\n", " [9.5492640e+00, 7.4848499e+00, 3.3095343e+00, ...,\n", " 1.2228072e+01, 1.2528176e+01, 1.1302447e+01],\n", " [7.2531962e+00, 5.3915715e+00, 2.0546637e+00, ...,\n", " 1.0486613e+01, 1.3482674e+01, 1.1387763e+01]],\n", "\n", " [[2.0096309e+00, 1.1509016e+00, 4.8980752e-01, ...,\n", " 2.0202532e+00, 2.7056713e+00, 2.6676538e+00],\n", " [1.2896461e+00, 1.1470729e+00, 2.6983759e-01, ...,\n", " 1.3239046e+00, 1.6525987e+00, 1.6284056e+00],\n", " [8.9463830e-01, 3.9747801e-01, 1.5476857e-01, ...,\n", " 1.1533502e+00, 1.1403440e+00, 9.9328095e-01],\n", " ...,\n", " [1.2956619e+00, 6.6536534e-01, 5.8920252e-01, ...,\n", " 4.1320643e+00, 4.3791137e+00, 4.1250496e+00],\n", " [1.7519032e+00, 1.2151116e+00, 9.6764684e-01, ...,\n", " 5.1612849e+00, 5.1797566e+00, 2.8651981e+00],\n", " [3.1523111e+00, 2.0467207e+00, 1.0773114e+00, ...,\n", " 3.6295269e+00, 4.3244128e+00, 2.4604805e+00]]], dtype=float32)\n", "Coordinates:\n", " * lon (lon) float32 81.25 83.75 86.25 88.75 ... 151.2 153.8 156.2 158.8\n", " * lat (lat) float32 28.75 26.25 23.75 21.25 ... -13.75 -16.25 -18.75\n", " * year (year) int64 1979 1980 1981 1982 1983 ... 2015 2016 2017 2018 2019\n", "Attributes:\n", " long_name: GPCP pentad precipitation (mm/day)\n", " units: mm/day
<xarray.DataArray (year: 41)>\n", "array([ 0.6, 0. , -0.1, 2.2, -0.9, -1.1, -0.4, 1.2, 1.1, -1.8, -0.1,\n", " 0.4, 1.5, -0.1, 0.1, 1.1, -1. , -0.5, 2.4, -1.6, -1.7, -0.7,\n", " -0.3, 1.1, 0.4, 0.7, -0.8, 0.9, -1.6, -0.7, 1.6, -1.6, -1. ,\n", " -0.2, -0.3, 0.7, 2.6, -0.6, -1. , 0.8, 0.5])\n", "Coordinates:\n", " * year (year) int64 1979 1980 1981 1982 1983 ... 2015 2016 2017 2018 2019
<xarray.DataArray (lat: 20, lon: 32)>\n", "array([[ 3.62929452e-01, 3.42366935e-01, 1.36141385e-01,\n", " 2.86181885e-01, 1.81635007e-01, 1.49191933e-01,\n", " 2.91984747e-01, 4.77790062e-01, 4.72934980e-01,\n", " 2.59950214e-01, 3.45925729e-01, 3.45836394e-01,\n", " 2.41839289e-01, 3.26920710e-01, 3.12002374e-01,\n", " 2.78955901e-01, 2.83943380e-01, 3.95803450e-01,\n", " 4.20742051e-01, 4.71375054e-01, 5.46175408e-01,\n", " 4.34086465e-01, 3.36992606e-01, 2.34458518e-01,\n", " 3.30061209e-01, 2.83089068e-01, 6.62670164e-02,\n", " -6.41038383e-02, -8.02357126e-02, -1.32435603e-01,\n", " -6.14116577e-02, -3.57627814e-01],\n", " [ 3.83909987e-01, 3.03894582e-01, 2.37785616e-01,\n", " 2.14941992e-01, 1.29919746e-01, 9.02064326e-02,\n", " 1.65142443e-01, 2.02511225e-01, 2.68816770e-01,\n", " 1.56273306e-01, 2.46169180e-01, 3.13843048e-01,\n", " 3.71282366e-01, 3.76425080e-01, 3.56976452e-01,\n", " 3.31840158e-01, 1.84695348e-01, 4.19843103e-01,\n", " 4.71277802e-01, 4.27318315e-01, 3.51365941e-01,\n", " 2.37002913e-01, 2.34139483e-01, 2.04142657e-01,\n", " 1.27191143e-01, -7.16084696e-02, -1.17030947e-01,\n", "...\n", " -1.63141768e-01, 2.02073967e-02, 8.64961743e-02,\n", " -6.46165929e-02, -1.30967877e-01, -1.11179142e-01,\n", " -2.19488532e-01, -3.89394397e-01, -5.41138832e-01,\n", " -5.55411086e-01, -4.62892783e-01, -3.87886292e-01,\n", " -2.86447708e-01, -1.87858258e-01, -4.88526454e-02,\n", " -6.03390847e-02, -5.92937090e-02, -1.09022023e-02,\n", " -1.62940949e-01, -1.06521328e-01, -1.63654736e-01,\n", " -2.66138362e-01, -3.63495740e-01, -3.56191141e-01,\n", " -3.43018577e-01, -4.02329768e-01],\n", " [ 2.18229674e-01, 1.50764053e-01, 4.34983952e-02,\n", " -2.31485281e-02, -1.24680834e-01, -2.82829671e-01,\n", " -2.27022137e-01, -1.80625263e-01, -6.68313177e-02,\n", " -1.22718872e-01, -7.70160205e-02, -1.92838402e-01,\n", " -2.63437000e-01, -4.77678670e-01, -6.13706400e-01,\n", " -5.25689479e-01, -3.34162312e-01, -3.46265743e-01,\n", " -2.73299118e-01, -3.36175549e-01, -1.75612517e-01,\n", " -7.18667328e-02, -2.45879382e-02, 2.50092725e-03,\n", " -1.46619847e-01, -7.41479082e-02, -3.23609599e-02,\n", " -1.90437806e-01, -2.59385838e-01, -3.08636754e-01,\n", " -2.93624430e-01, -3.09437350e-01]])\n", "Coordinates:\n", " * lon (lon) float32 81.25 83.75 86.25 88.75 ... 151.2 153.8 156.2 158.8\n", " * lat (lat) float32 28.75 26.25 23.75 21.25 ... -13.75 -16.25 -18.75
<xarray.DataArray 'mslp' (time: 31, lat: 25, lon: 33)>\n", "array([[[1013.425, 1015.375, 1017.25 , ..., 1015.8 , 1015.5 ,\n", " 1015.45 ],\n", " [1019.35 , 1020.625, 1021.65 , ..., 1009.975, 1012.825,\n", " 1015.25 ],\n", " [1023.575, 1024.4 , 1025.175, ..., 1008.025, 1013.475,\n", " 1017.35 ],\n", " ...,\n", " [1010.3 , 1009.675, 1008.85 , ..., 1006.8 , 1007.1 ,\n", " 1007.05 ],\n", " [1010.15 , 1009.875, 1009.05 , ..., 1006.775, 1007.25 ,\n", " 1007.125],\n", " [1010.075, 1009.9 , 1009.175, ..., 1007.15 , 1007.35 ,\n", " 1007.275]],\n", "\n", " [[1003.85 , 1006.025, 1008.05 , ..., 983.225, 983.4 ,\n", " 984.375],\n", " [1010.625, 1012.975, 1015.025, ..., 974.25 , 978.325,\n", " 981.95 ],\n", " [1015.5 , 1017.675, 1019.55 , ..., 976.95 , 981.475,\n", " 985.325],\n", "...\n", " [1011.025, 1010.875, 1010.8 , ..., 1008.775, 1008.8 ,\n", " 1008.825],\n", " [1010.625, 1010.525, 1010.65 , ..., 1009.075, 1008.9 ,\n", " 1008.65 ],\n", " [1011.025, 1010.85 , 1010.825, ..., 1008.475, 1008.325,\n", " 1008.75 ]],\n", "\n", " [[1024.175, 1025.275, 1026.675, ..., 1019.775, 1019.725,\n", " 1019.45 ],\n", " [1029.125, 1029.8 , 1030.6 , ..., 1006.275, 1007.75 ,\n", " 1010. ],\n", " [1034.4 , 1034.875, 1035.45 , ..., 995.475, 997.3 ,\n", " 1000.925],\n", " ...,\n", " [1011.825, 1011.275, 1010.975, ..., 1009.225, 1009.4 ,\n", " 1009.475],\n", " [1011.55 , 1011. , 1010.475, ..., 1008.975, 1009.125,\n", " 1008.975],\n", " [1011.575, 1011.275, 1010.925, ..., 1008.9 , 1009.2 ,\n", " 1009.025]]], dtype=float32)\n", "Coordinates:\n", " * lat (lat) float32 60.0 57.5 55.0 52.5 50.0 ... 10.0 7.5 5.0 2.5 0.0\n", " * lon (lon) float32 80.0 82.5 85.0 87.5 90.0 ... 152.5 155.0 157.5 160.0\n", " * time (time) datetime64[ns] 2021-12-01 2021-12-02 ... 2021-12-31\n", " X (lon) float32 80.0 82.5 85.0 87.5 90.0 ... 152.5 155.0 157.5 160.0\n", " Y (lat) float32 60.0 57.5 55.0 52.5 50.0 ... 10.0 7.5 5.0 2.5 0.0
<xarray.DataArray (lat: 21, lon: 33)>\n", "<Quantity([[ 2.27227035e-06 -6.35822645e-06 -1.46749070e-05 -2.08053115e-05\n", " -2.30204190e-05 -2.14179636e-05 -1.73882539e-05 -1.29046760e-05\n", " -9.84929500e-06 -9.54206099e-06 -1.12507049e-05 -1.23807539e-05\n", " -1.01248504e-05 -4.97617877e-06 -9.28601335e-07 -2.13257232e-06\n", " -8.04519296e-06 -1.36745794e-05 -1.37081362e-05 -8.49913312e-06\n", " -2.77764758e-06 -1.18851542e-06 -3.59719131e-06 -6.00788240e-06\n", " -5.68558831e-06 -3.57765352e-06 -2.98596786e-06 -5.61280440e-06\n", " -9.78352706e-06 -1.28709753e-05 -1.36409031e-05 -1.30405734e-05\n", " -1.21284570e-05]\n", " [ 1.59641488e-06 1.27951490e-07 -2.96147587e-06 -6.09627077e-06\n", " -7.14075885e-06 -6.30349598e-06 -5.16875042e-06 -4.56615673e-06\n", " -4.01070010e-06 -3.29133037e-06 -3.14996178e-06 -4.03313727e-06\n", " -5.06772766e-06 -5.30085544e-06 -5.33611754e-06 -7.06711890e-06\n", " -1.08910896e-05 -1.43131536e-05 -1.42338578e-05 -1.03047389e-05\n", " -5.71612173e-06 -3.96175474e-06 -5.53658936e-06 -7.88921794e-06\n", " -8.64479278e-06 -8.03751306e-06 -8.04710263e-06 -9.82336098e-06\n", " -1.24495122e-05 -1.41660453e-05 -1.44190610e-05 -1.42504704e-05\n", " -1.44654921e-05]\n", " [-2.25853323e-06 1.13172266e-06 2.50715062e-06 3.24109081e-06\n", " 4.54796451e-06 5.22032636e-06 3.42384104e-06 8.61968466e-08\n", "...\n", " -8.52059086e-06 -1.10716107e-05 -1.19881288e-05 -1.13950856e-05\n", " -9.98141214e-06]\n", " [ 3.40814844e-07 -1.87193747e-06 -2.55199434e-06 -1.83217317e-06\n", " -7.07930811e-07 -2.45485863e-07 -7.17219789e-07 -1.21135839e-06\n", " -6.77305687e-07 1.98063110e-07 -8.42058040e-07 -4.36596747e-06\n", " -7.21982209e-06 -6.11138032e-06 -2.24744337e-06 6.50249030e-08\n", " -5.85517420e-07 -8.93727851e-07 1.92617077e-06 5.34523983e-06\n", " 4.91827661e-06 6.98802749e-07 -2.36916696e-06 -8.49899076e-07\n", " 3.18171194e-06 5.20927594e-06 3.40469449e-06 -5.15668196e-07\n", " -4.15454945e-06 -6.51466556e-06 -7.57468620e-06 -7.70031620e-06\n", " -7.32678885e-06]\n", " [ 3.86059996e-06 4.83366406e-07 -1.67842807e-06 -2.14714011e-06\n", " -1.31075583e-06 -6.41593796e-07 -9.82999754e-07 -1.46979384e-06\n", " -7.01365099e-07 5.12796672e-07 -4.05116393e-07 -4.63963087e-06\n", " -8.87224706e-06 -8.92474419e-06 -4.83957896e-06 -1.08197146e-06\n", " -7.43141058e-07 -2.15767543e-06 -1.95926642e-06 -1.83126193e-07\n", " -9.75771821e-09 -3.08505634e-06 -6.80966442e-06 -7.53604087e-06\n", " -4.54671714e-06 -6.88500940e-08 3.41230334e-06 4.97175514e-06\n", " 4.77222735e-06 3.07068630e-06 7.25767405e-07 -1.38476001e-06\n", " -2.89335635e-06]], '1 / second')>\n", "Coordinates:\n", " * lon (lon) float32 80.0 82.5 85.0 87.5 90.0 ... 152.5 155.0 157.5 160.0\n", " * lat (lat) float32 30.0 27.5 25.0 22.5 20.0 ... -12.5 -15.0 -17.5 -20.0\n", " level float32 850.0
<xarray.DataArray (lat: 21, lon: 33)>\n", "array([[ 2.27227035e-06, -6.35822645e-06, -1.46749070e-05,\n", " -2.08053115e-05, -2.30204190e-05, -2.14179636e-05,\n", " -1.73882539e-05, -1.29046760e-05, -9.84929500e-06,\n", " -9.54206099e-06, -1.12507049e-05, -1.23807539e-05,\n", " -1.01248504e-05, -4.97617877e-06, -9.28601335e-07,\n", " -2.13257232e-06, -8.04519296e-06, -1.36745794e-05,\n", " -1.37081362e-05, -8.49913312e-06, -2.77764758e-06,\n", " -1.18851542e-06, -3.59719131e-06, -6.00788240e-06,\n", " -5.68558831e-06, -3.57765352e-06, -2.98596786e-06,\n", " -5.61280440e-06, -9.78352706e-06, -1.28709753e-05,\n", " -1.36409031e-05, -1.30405734e-05, -1.21284570e-05],\n", " [ 1.59641488e-06, 1.27951490e-07, -2.96147587e-06,\n", " -6.09627077e-06, -7.14075885e-06, -6.30349598e-06,\n", " -5.16875042e-06, -4.56615673e-06, -4.01070010e-06,\n", " -3.29133037e-06, -3.14996178e-06, -4.03313727e-06,\n", " -5.06772766e-06, -5.30085544e-06, -5.33611754e-06,\n", " -7.06711890e-06, -1.08910896e-05, -1.43131536e-05,\n", " -1.42338578e-05, -1.03047389e-05, -5.71612173e-06,\n", " -3.96175474e-06, -5.53658936e-06, -7.88921794e-06,\n", " -8.64479278e-06, -8.03751306e-06, -8.04710263e-06,\n", "...\n", " -7.17219789e-07, -1.21135839e-06, -6.77305687e-07,\n", " 1.98063110e-07, -8.42058040e-07, -4.36596747e-06,\n", " -7.21982209e-06, -6.11138032e-06, -2.24744337e-06,\n", " 6.50249030e-08, -5.85517420e-07, -8.93727851e-07,\n", " 1.92617077e-06, 5.34523983e-06, 4.91827661e-06,\n", " 6.98802749e-07, -2.36916696e-06, -8.49899076e-07,\n", " 3.18171194e-06, 5.20927594e-06, 3.40469449e-06,\n", " -5.15668196e-07, -4.15454945e-06, -6.51466556e-06,\n", " -7.57468620e-06, -7.70031620e-06, -7.32678885e-06],\n", " [ 3.86059996e-06, 4.83366406e-07, -1.67842807e-06,\n", " -2.14714011e-06, -1.31075583e-06, -6.41593796e-07,\n", " -9.82999754e-07, -1.46979384e-06, -7.01365099e-07,\n", " 5.12796672e-07, -4.05116393e-07, -4.63963087e-06,\n", " -8.87224706e-06, -8.92474419e-06, -4.83957896e-06,\n", " -1.08197146e-06, -7.43141058e-07, -2.15767543e-06,\n", " -1.95926642e-06, -1.83126193e-07, -9.75771821e-09,\n", " -3.08505634e-06, -6.80966442e-06, -7.53604087e-06,\n", " -4.54671714e-06, -6.88500940e-08, 3.41230334e-06,\n", " 4.97175514e-06, 4.77222735e-06, 3.07068630e-06,\n", " 7.25767405e-07, -1.38476001e-06, -2.89335635e-06]])\n", "Coordinates:\n", " * lon (lon) float32 80.0 82.5 85.0 87.5 90.0 ... 152.5 155.0 157.5 160.0\n", " * lat (lat) float32 30.0 27.5 25.0 22.5 20.0 ... -12.5 -15.0 -17.5 -20.0\n", " level float32 850.0\n", "Attributes:\n", " units: 1 / second