diff --git a/Jupyter_notebooks/1.2-rotated_grids.ipynb b/Jupyter_notebooks/1.2-rotated_grids.ipynb index 9ace309cd1f519f0977b568452fa71f07a45de01..31b19f7ff6d02f55466f9b20c3f8c67356701e8e 100644 --- a/Jupyter_notebooks/1.2-rotated_grids.ipynb +++ b/Jupyter_notebooks/1.2-rotated_grids.ipynb @@ -415,14 +415,14 @@ " rotated_pole |S1 b''\n", "Attributes:\n", " Conventions: CF-1.7\n", - " comment: Generated on marenostrum4
array(['2021-08-03T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13.,\n", - " 14., 15., 16., 17., 18., 19., 20., 21., 22., 23.])
array([[16.350338, 16.43293 , 16.515146, ..., 16.515146, 16.43293 , 16.350338],\n", + " comment: Generated on marenostrum4
array(['2021-08-03T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13.,\n", + " 14., 15., 16., 17., 18., 19., 20., 21., 22., 23.])
array([[16.350338, 16.43293 , 16.515146, ..., 16.515146, 16.43293 , 16.350338],\n", " [16.527426, 16.610239, 16.692677, ..., 16.692677, 16.610243, 16.527426],\n", " [16.704472, 16.787508, 16.870167, ..., 16.870167, 16.78751 , 16.704472],\n", " ...,\n", " [58.32095 , 58.472683, 58.62431 , ..., 58.62431 , 58.472683, 58.32095 ],\n", " [58.426285, 58.578201, 58.730026, ..., 58.730026, 58.578201, 58.426285],\n", - " [58.530792, 58.682899, 58.834919, ..., 58.834919, 58.682903, 58.530792]])
array([[-22.181265, -22.016672, -21.851799, ..., 41.851795, 42.016666,\n", + " [58.530792, 58.682899, 58.834919, ..., 58.834919, 58.682903, 58.530792]])
array([[-22.181265, -22.016672, -21.851799, ..., 41.851795, 42.016666,\n", " 42.181259],\n", " [-22.278179, -22.113186, -21.947905, ..., 41.947899, 42.113174,\n", " 42.278172],\n", @@ -434,10 +434,10 @@ " [-67.901878, -67.722473, -67.541939, ..., 87.541939, 87.722458,\n", " 87.901871],\n", " [-68.228035, -68.04982 , -67.870514, ..., 87.870506, 88.04982 ,\n", - " 88.228035]])
array([-27. , -26.799999, -26.6 , ..., 26.6 , 26.800001,\n", - " 27. ])
array([-35. , -34.799999, -34.599998, ..., 34.600002, 34.799999,\n", - " 35. ])
array([['2021-08-03T00:00:00.000000000', '2021-08-07T00:00:00.000000000']],\n", - " dtype='datetime64[ns]')
[2282904 values with dtype=float32]
array(b'', dtype='|S1')
array([-27. , -26.799999, -26.6 , ..., 26.6 , 26.800001,\n", + " 27. ])
array([-35. , -34.799999, -34.599998, ..., 34.600002, 34.799999,\n", + " 35. ])
array([['2021-08-03T00:00:00.000000000', '2021-08-07T00:00:00.000000000']],\n", + " dtype='datetime64[ns]')
[2282904 values with dtype=float32]
array(b'', dtype='|S1')
array(['2021-08-03T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13.,\n", - " 14., 15., 16., 17., 18., 19., 20., 21., 22., 23.])
array([[16.350338, 16.43293 , 16.515146, ..., 16.515146, 16.43293 , 16.350338],\n", + " comment: Generated on marenostrum4
array(['2021-08-03T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13.,\n", + " 14., 15., 16., 17., 18., 19., 20., 21., 22., 23.])
array([[16.350338, 16.43293 , 16.515146, ..., 16.515146, 16.43293 , 16.350338],\n", " [16.527426, 16.610239, 16.692677, ..., 16.692677, 16.610243, 16.527426],\n", " [16.704472, 16.787508, 16.870167, ..., 16.870167, 16.78751 , 16.704472],\n", " ...,\n", " [58.32095 , 58.472683, 58.62431 , ..., 58.62431 , 58.472683, 58.32095 ],\n", " [58.426285, 58.578201, 58.730026, ..., 58.730026, 58.578201, 58.426285],\n", - " [58.530792, 58.682899, 58.834919, ..., 58.834919, 58.682903, 58.530792]])
array([[-22.181265, -22.016672, -21.851799, ..., 41.851795, 42.016666,\n", + " [58.530792, 58.682899, 58.834919, ..., 58.834919, 58.682903, 58.530792]])
array([[-22.181265, -22.016672, -21.851799, ..., 41.851795, 42.016666,\n", " 42.181259],\n", " [-22.278179, -22.113186, -21.947905, ..., 41.947899, 42.113174,\n", " 42.278172],\n", @@ -1217,10 +1217,10 @@ " [-67.901878, -67.722473, -67.541939, ..., 87.541939, 87.722458,\n", " 87.901871],\n", " [-68.228035, -68.04982 , -67.870514, ..., 87.870506, 88.04982 ,\n", - " 88.228035]])
array([-27. , -26.799999, -26.6 , ..., 26.6 , 26.800001,\n", - " 27. ])
array([-35. , -34.799999, -34.599998, ..., 34.600002, 34.799999,\n", - " 35. ])
array([['2021-08-03T00:00:00.000000000', '2021-08-07T00:00:00.000000000']],\n", - " dtype='datetime64[ns]')
[2282904 values with dtype=float32]
array(b'', dtype='|S1')
array([-27. , -26.799999, -26.6 , ..., 26.6 , 26.800001,\n", + " 27. ])
array([-35. , -34.799999, -34.599998, ..., 34.600002, 34.799999,\n", + " 35. ])
array([['2021-08-03T00:00:00.000000000', '2021-08-07T00:00:00.000000000']],\n", + " dtype='datetime64[ns]')
[2282904 values with dtype=float32]
array(b'', dtype='|S1')
array(['2022-06-26T00:00:00.000000000', '2022-06-26T01:00:00.000000000',\n", + " Lambert_conformal int32 -2147483647
array(['2022-06-26T00:00:00.000000000', '2022-06-26T01:00:00.000000000',\n", " '2022-06-26T02:00:00.000000000', '2022-06-26T03:00:00.000000000',\n", " '2022-06-26T04:00:00.000000000', '2022-06-26T05:00:00.000000000',\n", " '2022-06-26T06:00:00.000000000', '2022-06-26T07:00:00.000000000',\n", @@ -434,9 +434,9 @@ " '2022-06-27T18:00:00.000000000', '2022-06-27T19:00:00.000000000',\n", " '2022-06-27T20:00:00.000000000', '2022-06-27T21:00:00.000000000',\n", " '2022-06-27T22:00:00.000000000', '2022-06-27T23:00:00.000000000'],\n", - " dtype='datetime64[ns]')
[190244 values with dtype=float32]
[190244 values with dtype=float32]
array([-2125847.534765, -2113847.610112, -2101847.010684, ..., 3574152.262286,\n", - " 3586150.20523 , 3598150.638118])
array([-2067138.479486, -2055138.337429, -2043137.546656, ..., 2672862.222453,\n", - " 2684862.70005 , 2696862.03252 ])
array([0.], dtype=float32)
[9131712 values with dtype=float32]
array(-2147483647, dtype=int32)
[190244 values with dtype=float32]
[190244 values with dtype=float32]
array([-2125847.534765, -2113847.610112, -2101847.010684, ..., 3574152.262286,\n", + " 3586150.20523 , 3598150.638118])
array([-2067138.479486, -2055138.337429, -2043137.546656, ..., 2672862.222453,\n", + " 2684862.70005 , 2696862.03252 ])
array([0.], dtype=float32)
[9131712 values with dtype=float32]
array(-2147483647, dtype=int32)
<xarray.Dataset>\n", + "Dimensions: (time: 1, y: 236, x: 210)\n", + "Coordinates:\n", + " * time (time) float64 0.0\n", + " lat (y, x) float32 -43.52 -43.52 -43.52 -43.52 ... 49.6 49.6 49.6\n", + " lon (y, x) float32 -18.91 -18.46 -18.01 -17.56 ... 74.22 74.67 75.12\n", + " * x (x) float64 -1.01e+05 -5.101e+04 ... 1.03e+07 1.035e+07\n", + " * y (y) float64 -5.382e+06 -5.332e+06 ... 6.318e+06 6.368e+06\n", + "Data variables:\n", + " var_aux (time, y, x) float32 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0 0.0\n", + " mercator int32 -2147483647\n", + " cell_area (y, x) float32 1.316e+09 1.316e+09 ... 1.051e+09 1.051e+09
array([0.])
array([[-43.5182 , -43.5182 , -43.5182 , ..., -43.5182 , -43.5182 ,\n", + " -43.5182 ],\n", + " [-43.191082, -43.191082, -43.191082, ..., -43.191082, -43.191082,\n", + " -43.191082],\n", + " [-42.8622 , -42.8622 , -42.8622 , ..., -42.8622 , -42.8622 ,\n", + " -42.8622 ],\n", + " ...,\n", + " [ 49.016712, 49.016712, 49.016712, ..., 49.016712, 49.016712,\n", + " 49.016712],\n", + " [ 49.31089 , 49.31089 , 49.31089 , ..., 49.31089 , 49.31089 ,\n", + " 49.31089 ],\n", + " [ 49.60332 , 49.60332 , 49.60332 , ..., 49.60332 , 49.60332 ,\n", + " 49.60332 ]], dtype=float32)
array([[-18.9089 , -18.459013, -18.009129, ..., 74.217415, 74.6673 ,\n", + " 75.11718 ],\n", + " [-18.9089 , -18.459013, -18.009129, ..., 74.217415, 74.6673 ,\n", + " 75.11718 ],\n", + " [-18.9089 , -18.459013, -18.009129, ..., 74.217415, 74.6673 ,\n", + " 75.11718 ],\n", + " ...,\n", + " [-18.9089 , -18.459013, -18.009129, ..., 74.217415, 74.6673 ,\n", + " 75.11718 ],\n", + " [-18.9089 , -18.459013, -18.009129, ..., 74.217415, 74.6673 ,\n", + " 75.11718 ],\n", + " [-18.9089 , -18.459013, -18.009129, ..., 74.217415, 74.6673 ,\n", + " 75.11718 ]], dtype=float32)
array([-1.010145e+05, -5.101450e+04, -1.014500e+03, ..., 1.024899e+07,\n", + " 1.029899e+07, 1.034899e+07])
array([-5382460., -5332460., -5282460., ..., 6267540., 6317540., 6367540.])
array([[[0., 0., ..., 0., 0.],\n", + " [0., 0., ..., 0., 0.],\n", + " ...,\n", + " [0., 0., ..., 0., 0.],\n", + " [0., 0., ..., 0., 0.]]], dtype=float32)
array(-2147483647, dtype=int32)
array([[1.315942e+09, 1.315937e+09, 1.315942e+09, ..., 1.315931e+09,\n", + " 1.315954e+09, 1.315931e+09],\n", + " [1.330202e+09, 1.330197e+09, 1.330202e+09, ..., 1.330191e+09,\n", + " 1.330214e+09, 1.330191e+09],\n", + " [1.344550e+09, 1.344544e+09, 1.344550e+09, ..., 1.344538e+09,\n", + " 1.344561e+09, 1.344538e+09],\n", + " ...,\n", + " [1.076385e+09, 1.076381e+09, 1.076385e+09, ..., 1.076376e+09,\n", + " 1.076394e+09, 1.076376e+09],\n", + " [1.063687e+09, 1.063683e+09, 1.063687e+09, ..., 1.063678e+09,\n", + " 1.063696e+09, 1.063678e+09],\n", + " [1.051047e+09, 1.051043e+09, 1.051047e+09, ..., 1.051038e+09,\n", + " 1.051056e+09, 1.051038e+09]], dtype=float32)
<xarray.Dataset>\n", + "Dimensions: (time: 1, lev: 1, y: 236, x: 210)\n", + "Coordinates:\n", + " * time (time) datetime64[ns] 2000-01-01\n", + " * lev (lev) float64 0.0\n", + " lat (y, x) float64 -43.52 -43.52 -43.52 -43.52 ... 49.6 49.6 49.6\n", + " lon (y, x) float64 -18.91 -18.46 -18.01 -17.56 ... 74.22 74.67 75.12\n", + " * y (y) float64 -5.382e+06 -5.332e+06 ... 6.318e+06 6.368e+06\n", + " * x (x) float64 -1.01e+05 -5.101e+04 ... 1.03e+07 1.035e+07\n", + "Data variables:\n", + " var_aux (time, lev, y, x) float32 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0\n", + " cell_area (time, lev, y, x) float32 1.316e+09 1.316e+09 ... 1.051e+09\n", + " mercator |S1 b''\n", + "Attributes:\n", + " Conventions: CF-1.7
array(['2000-01-01T00:00:00.000000000'], dtype='datetime64[ns]')
array([0.])
array([[-43.5182 , -43.5182 , -43.5182 , ..., -43.5182 , -43.5182 ,\n", + " -43.5182 ],\n", + " [-43.191082, -43.191082, -43.191082, ..., -43.191082, -43.191082,\n", + " -43.191082],\n", + " [-42.862202, -42.862202, -42.862202, ..., -42.862202, -42.862202,\n", + " -42.862202],\n", + " ...,\n", + " [ 49.016712, 49.016712, 49.016712, ..., 49.016712, 49.016712,\n", + " 49.016712],\n", + " [ 49.31089 , 49.31089 , 49.31089 , ..., 49.31089 , 49.31089 ,\n", + " 49.31089 ],\n", + " [ 49.603321, 49.603321, 49.603321, ..., 49.603321, 49.603321,\n", + " 49.603321]])
array([[-18.908899, -18.459013, -18.009129, ..., 74.217415, 74.667297,\n", + " 75.11718 ],\n", + " [-18.908899, -18.459013, -18.009129, ..., 74.217415, 74.667297,\n", + " 75.11718 ],\n", + " [-18.908899, -18.459013, -18.009129, ..., 74.217415, 74.667297,\n", + " 75.11718 ],\n", + " ...,\n", + " [-18.908899, -18.459013, -18.009129, ..., 74.217415, 74.667297,\n", + " 75.11718 ],\n", + " [-18.908899, -18.459013, -18.009129, ..., 74.217415, 74.667297,\n", + " 75.11718 ],\n", + " [-18.908899, -18.459013, -18.009129, ..., 74.217415, 74.667297,\n", + " 75.11718 ]])
array([-5382460., -5332460., -5282460., ..., 6267540., 6317540., 6367540.])
array([-1.010145e+05, -5.101450e+04, -1.014500e+03, ..., 1.024899e+07,\n", + " 1.029899e+07, 1.034899e+07])
array([[[[0., ..., 0.],\n", + " ...,\n", + " [0., ..., 0.]]]], dtype=float32)
array([[[[1.315942e+09, ..., 1.315931e+09],\n", + " ...,\n", + " [1.051047e+09, ..., 1.051038e+09]]]], dtype=float32)
array(b'', dtype='|S1')
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([0.])
array([41.15, 41.25, 41.35, 41.45, 41.55, 41.65, 41.75, 41.85, 41.95, 42.05])
array([1.85, 1.95, 2.05, 2.15, 2.25, 2.35, 2.45, 2.55, 2.65, 2.75])
array(b'', dtype='|S1')
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([0.])
array([41.15, 41.25, 41.35, 41.45, 41.55, 41.65, 41.75, 41.85, 41.95, 42.05])
array([1.85, 1.95, 2.05, 2.15, 2.25, 2.35, 2.45, 2.55, 2.65, 2.75])
array(b'', dtype='|S1')
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([0.])
array([-27. , -26.8, -26.6, ..., 26.6, 26.8, 27. ])
array([-35. , -34.8, -34.6, ..., 34.6, 34.8, 35. ])
array([[16.350341, 16.432934, 16.515147, ..., 16.515147, 16.432934, 16.350341],\n", + " Conventions: CF-1.7
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([0.])
array([-27. , -26.8, -26.6, ..., 26.6, 26.8, 27. ])
array([-35. , -34.8, -34.6, ..., 34.6, 34.8, 35. ])
array([[16.350341, 16.432934, 16.515147, ..., 16.515147, 16.432934, 16.350341],\n", " [16.527429, 16.610244, 16.692679, ..., 16.692679, 16.610244, 16.527429],\n", " [16.704473, 16.78751 , 16.870167, ..., 16.870167, 16.78751 , 16.704473],\n", " ...,\n", " [58.320945, 58.472682, 58.624312, ..., 58.624312, 58.472682, 58.320945],\n", " [58.426281, 58.578207, 58.730029, ..., 58.730029, 58.578207, 58.426281],\n", - " [58.530794, 58.682904, 58.834914, ..., 58.834914, 58.682904, 58.530794]])
array([[-22.181266, -22.01667 , -21.851799, ..., 41.851799, 42.01667 ,\n", + " [58.530794, 58.682904, 58.834914, ..., 58.834914, 58.682904, 58.530794]])
array([[-22.181266, -22.01667 , -21.851799, ..., 41.851799, 42.01667 ,\n", " 42.181266],\n", " [-22.278178, -22.113181, -21.947906, ..., 41.947906, 42.113181,\n", " 42.278178],\n", @@ -889,7 +889,7 @@ " [-67.901883, -67.722467, -67.541948, ..., 87.541948, 87.722467,\n", " 87.901883],\n", " [-68.22804 , -68.049824, -67.870513, ..., 87.870513, 88.049824,\n", - " 88.22804 ]])
array(b'', dtype='|S1')
array(b'', dtype='|S1')
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11.,\n", + " Conventions: CF-1.7
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11.,\n", " 12., 13., 14., 15., 16., 17., 18., 19., 20., 21., 22., 23.,\n", " 24., 25., 26., 27., 28., 29., 30., 31., 32., 33., 34., 35.,\n", " 36., 37., 38., 39., 40., 41., 42., 43., 44., 45., 46., 47.,\n", @@ -1539,7 +1539,7 @@ " 96., 97., 98., 99., 100., 101., 102., 103., 104., 105., 106., 107.,\n", " 108., 109., 110., 111., 112., 113., 114., 115., 116., 117., 118., 119.,\n", " 120., 121., 122., 123., 124., 125., 126., 127., 128., 129., 130., 131.,\n", - " 132., 133.])
array(['ES0266A', 'ES0392A', 'ES0395A', 'ES0559A', 'ES0567A', 'ES0584A',\n", + " 132., 133.])
array(['ES0266A', 'ES0392A', 'ES0395A', 'ES0559A', 'ES0567A', 'ES0584A',\n", " 'ES0586A', 'ES0691A', 'ES0692A', 'ES0694A', 'ES0700A', 'ES0704A',\n", " 'ES0963A', 'ES0971A', 'ES0991A', 'ES1018A', 'ES1117A', 'ES1120A',\n", " 'ES1122A', 'ES1123A', 'ES1124A', 'ES1125A', 'ES1126A', 'ES1135A',\n", @@ -1561,7 +1561,7 @@ " 'ES2011A', 'ES2012A', 'ES2017A', 'ES2027A', 'ES2033A', 'ES2034A',\n", " 'ES2035A', 'ES2071A', 'ES2079A', 'ES2043A', 'ES2090A', 'ES0554A',\n", " 'ES0977A', 'ES1398A', 'ES1200A', 'ES2087A', 'ES2091A', 'ES2088A',\n", - " 'ES1908A', 'ES9994A'], dtype=object)
array(['urban-centre', 'urban-suburban', 'urban-centre', 'urban-centre',\n", + " 'ES1908A', 'ES9994A'], dtype=object)
array(['urban-centre', 'urban-suburban', 'urban-centre', 'urban-centre',\n", " 'urban-centre', 'urban-suburban', 'urban-centre', 'urban-centre',\n", " 'urban-centre', 'urban-suburban', 'urban-centre', 'urban-suburban',\n", " 'urban-suburban', 'urban-suburban', 'urban-suburban', 'urban-centre',\n", @@ -1593,7 +1593,7 @@ " 'urban-suburban', 'urban-suburban', 'rural', 'urban-centre', 'rural',\n", " 'rural', 'urban-suburban', 'urban-centre', 'urban-suburban', 'nan',\n", " 'nan', 'nan', 'nan', 'nan', 'nan', 'nan', 'nan', 'nan', 'nan', 'nan'],\n", - " dtype=object)
array([41.379322, 41.727704, 41.567824, 41.387424, 41.384906, 41.482016,\n", + " dtype=object)
array([41.379322, 41.727704, 41.567824, 41.387424, 41.384906, 41.482016,\n", " 41.413621, 41.403716, 41.37076 , 41.392157, 41.515609, 41.552421,\n", " 41.473378, 41.450745, 41.884897, 41.556115, 41.111995, 41.115881,\n", " 41.193743, 41.155004, 41.159532, 41.73028 , 41.475573, 41.577627,\n", @@ -1615,7 +1615,7 @@ " 41.400771, 41.415319, 40.552819, 41.451311, 41.242375, 41.544105,\n", " 41.567703, 41.120064, 41.559607, 41.230073, 41.418413, 40.289167,\n", " 41.604722, 41.53667 , 40.281389, 41.929283, 40.5799 , 41.77106 ,\n", - " 41.239069, 42.358363])
array([2.08614 , 1.838531, 2.014598, 2.164918, 2.119574, 2.188297, 2.015986,\n", + " 41.239069, 42.358363])
array([2.08614 , 1.838531, 2.014598, 2.164918, 2.119574, 2.188297, 2.015986,\n", " 2.204736, 2.114771, 2.009802, 2.124997, 2.26525 , 1.982017, 1.975021,\n", " 2.874243, 2.007401, 1.151879, 1.191975, 1.236904, 1.217734, 1.239709,\n", " 1.825306, 1.92319 , 1.625893, 2.222245, 2.842138, 1.120172, 0.581651,\n", @@ -1634,7 +1634,7 @@ " 2.816547, 0.546797, 1.999635, 1.990521, 0.529983, 2.248236, 1.859334,\n", " 0.829933, 1.637614, 1.254472, 1.995963, 0.547183, 2.123899, 0.289167,\n", " 1.604722, 2.183611, 0.281389, 2.257302, 0.5535 , 2.250647, 1.856564,\n", - " 1.459455])
array(['2017-01-01T00:00:00.000000000', '2017-01-02T00:00:00.000000000',\n", + " Conventions: CF-1.7
array(['2017-01-01T00:00:00.000000000', '2017-01-02T00:00:00.000000000',\n", " '2017-01-03T00:00:00.000000000', ..., '2017-12-29T00:00:00.000000000',\n", " '2017-12-30T00:00:00.000000000', '2017-12-31T00:00:00.000000000'],\n", - " dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13.,\n", + " dtype='datetime64[ns]')
array([ 0., 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13.,\n", " 14., 15., 16., 17., 18., 19., 20., 21., 22., 23., 24., 25., 26., 27.,\n", " 28., 29., 30., 31., 32., 33., 34., 35., 36., 37., 38., 39., 40., 41.,\n", " 42., 43., 44., 45., 46., 47., 48., 49., 50., 51., 52., 53., 54., 55.,\n", " 56., 57., 58., 59., 60., 61., 62., 63., 64., 65., 66., 67., 68., 69.,\n", - " 70., 71., 72., 73., 74., 75., 76., 77., 78., 79., 80., 81., 82.])
array(['Barcelona (Eixample)', 'Badalona (guàrdia urbana)',\n", + " 70., 71., 72., 73., 74., 75., 76., 77., 78., 79., 80., 81., 82.])
array(['Barcelona (Eixample)', 'Badalona (guàrdia urbana)',\n", " 'Badalona (Assamblea de Catalunya)',\n", " 'Barcelona (Pl. de la Universitat)', 'Barcelona (Poblenou)',\n", " 'Barcelona (Zona Universitària)', 'Barcelona (el Port Vell)',\n", @@ -2652,7 +2652,7 @@ " 'Montseny (la Castanya)', 'Caldes de Montbui (Ajuntament)',\n", " "Sant Feliu de Llobregat (Eugeni d'Ors)",\n", " "La Seu d'Urgell (CC Les Monges)", 'Vic (Centre Cívic Santa Anna)'],\n", - " dtype=object)
array(['ES1438A', 'ES1928A', 'ES2027A', 'ES0559A', 'ES0691A', 'ES0567A',\n", + " dtype=object)
array(['ES1438A', 'ES1928A', 'ES2027A', 'ES0559A', 'ES0691A', 'ES0567A',\n", " 'ES1870A', 'ES1852A', 'ES1900A', 'ES1843A', 'ES1839A', 'ES0963A',\n", " 'ES0266A', 'ES1891A', 'ES1135A', 'ES1397A', 'ES0392A', 'ES1861A',\n", " 'ES1126A', 'ES1390A', 'ES1413A', 'ES0586A', 'ES1931A', 'ES0704A',\n", @@ -2665,23 +2665,23 @@ " 'ES1480A', 'ES1856A', 'ES1851A', 'ES1964A', 'ES1225A', 'ES1312A',\n", " 'ES1262A', 'ES1910A', 'ES1983A', 'ES1899A', 'ES1903A', 'ES1396A',\n", " 'ES1348A', 'ES1123A', 'ES1117A', 'ES2033A', 'ES2011A', 'ES1982A',\n", - " 'ES1778A', 'ES1680A', 'ES1362A', 'ES9994A', 'ES1874A'], dtype=object)
array([[19.6 , 27.2 , 35.7 , ..., 24.6 , 27.4 , 17.3 ],\n", + " 'ES1778A', 'ES1680A', 'ES1362A', 'ES9994A', 'ES1874A'], dtype=object)
array([[19.6 , 27.2 , 35.7 , ..., 24.6 , 27.4 , 17.3 ],\n", " [ nan, 20.86, nan, ..., 21. , nan, 12.5 ],\n", " [20. , 23. , 32. , ..., 24. , 15. , 13. ],\n", " ...,\n", " [ nan, 22.06, 35.84, ..., 14.09, nan, nan],\n", " [ nan, nan, nan, ..., nan, nan, nan],\n", - " [ nan, nan, nan, ..., nan, nan, nan]])
array([nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", + " [ nan, nan, nan, ..., nan, nan, nan]])
array([nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", - " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan])
array([nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", + " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan])
array([nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan,\n", - " nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan, nan])
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([0.])
array([-795137.125, -791137.125, -787137.125, ..., 780862.875, 784862.875,\n", - " 788862.875])
array([-805847.688, -801847.688, -797847.688, ..., 770152.312, 774152.312,\n", - " 778152.312])
[157609 values with dtype=float64]
[157609 values with dtype=float64]
array(b'', dtype='|S1')
array(['1996-12-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([0.])
array([-795137.125, -791137.125, -787137.125, ..., 780862.875, 784862.875,\n", + " 788862.875])
array([-805847.688, -801847.688, -797847.688, ..., 770152.312, 774152.312,\n", + " 778152.312])
[157609 values with dtype=float64]
[157609 values with dtype=float64]
array(b'', dtype='|S1')
43819 rows × 4 columns
\n", + "43818 rows × 4 columns
\n", "" ], "text/plain": [ " SO2-Darsena SO2-UM NO2-UM NO2-ZAL Prat\n", "Horario: UTC \n", - "2017-01-01 00:00:00.000000 0.18 NaN 63.38 46.01\n", - "2017-01-01 01:00:00.000003 0.20 NaN 64.64 49.08\n", - "2017-01-01 01:59:59.999997 0.18 NaN 68.16 53.00\n", - "2017-01-01 03:00:00.000000 0.04 NaN 68.29 46.75\n", - "2017-01-01 04:00:00.000003 0.00 NaN 62.91 39.65\n", + "2017-01-01 00:00:00.000000 0.20 NaN 64.64 49.08\n", + "2017-01-01 01:00:00.000003 0.18 NaN 68.16 53.00\n", + "2017-01-01 01:59:59.999997 0.04 NaN 68.29 46.75\n", + "2017-01-01 03:00:00.000000 0.00 NaN 62.91 39.65\n", + "2017-01-01 04:00:00.000003 0.00 NaN 47.13 28.86\n", "... ... ... ... ...\n", - "2021-12-31 13:59:59.999997 0.77 NaN 35.38 24.96\n", - "2021-12-31 15:00:00.000000 0.37 NaN 35.80 24.16\n", - "2021-12-31 16:00:00.000003 0.23 NaN 29.10 25.79\n", - "2021-12-31 16:59:59.999997 0.14 NaN 9.24 29.82\n", - "2021-12-31 18:00:00.000000 0.20 NaN 12.76 28.66\n", + "2021-12-31 13:00:00.000003 0.77 NaN 35.38 24.96\n", + "2021-12-31 13:59:59.999997 0.37 NaN 35.80 24.16\n", + "2021-12-31 15:00:00.000000 0.23 NaN 29.10 25.79\n", + "2021-12-31 16:00:00.000003 0.14 NaN 9.24 29.82\n", + "2021-12-31 16:59:59.999997 0.20 NaN 12.76 28.66\n", "\n", - "[43819 rows x 4 columns]" + "[43818 rows x 4 columns]" ] }, "execution_count": 2, @@ -185,12 +185,12 @@ " '2017-01-01 06:00:00', '2017-01-01 07:00:00',\n", " '2017-01-01 08:00:00', '2017-01-01 09:00:00',\n", " ...\n", - " '2021-12-31 09:00:00', '2021-12-31 10:00:00',\n", - " '2021-12-31 11:00:00', '2021-12-31 12:00:00',\n", - " '2021-12-31 13:00:00', '2021-12-31 14:00:00',\n", - " '2021-12-31 15:00:00', '2021-12-31 16:00:00',\n", - " '2021-12-31 17:00:00', '2021-12-31 18:00:00'],\n", - " dtype='datetime64[ns]', length=43819, freq=None)" + " '2021-12-31 08:00:00', '2021-12-31 09:00:00',\n", + " '2021-12-31 10:00:00', '2021-12-31 11:00:00',\n", + " '2021-12-31 12:00:00', '2021-12-31 13:00:00',\n", + " '2021-12-31 14:00:00', '2021-12-31 15:00:00',\n", + " '2021-12-31 16:00:00', '2021-12-31 17:00:00'],\n", + " dtype='datetime64[ns]', length=43818, freq=None)" ] }, "execution_count": 3, @@ -302,9 +302,9 @@ "array([datetime.datetime(2017, 1, 1, 0, 0),\n", " datetime.datetime(2017, 1, 1, 1, 0),\n", " datetime.datetime(2017, 1, 1, 2, 0), ...,\n", + " datetime.datetime(2021, 12, 31, 15, 0),\n", " datetime.datetime(2021, 12, 31, 16, 0),\n", - " datetime.datetime(2021, 12, 31, 17, 0),\n", - " datetime.datetime(2021, 12, 31, 18, 0)], dtype=object)" + " datetime.datetime(2021, 12, 31, 17, 0)], dtype=object)" ] }, "execution_count": 5, @@ -416,7 +416,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/points_nes.py:401: UserWarning: WARNING!!! Different data types for variable station_nameInput dtype=<xarray.Dataset>\n", - "Dimensions: (time: 43819, station: 2)\n", + "Dimensions: (time: 43818, station: 2)\n", "Coordinates:\n", - " * time (time) datetime64[ns] 2017-01-01 ... 2021-12-31T18:00:00\n", + " * time (time) datetime64[ns] 2017-01-01 ... 2021-12-31T17:00:00\n", " * station (station) float64 0.0 1.0\n", "Data variables:\n", " station_name (station) object 'NO2-UM' 'NO2-ZAL Prat'\n", - " sconcno2 (time, station) float64 63.38 46.01 64.64 ... 12.76 28.66\n", + " sconcno2 (time, station) float64 64.64 49.08 68.16 ... 12.76 28.66\n", " lat (station) float64 41.37 41.32\n", " lon (station) float64 2.185 2.135\n", "Attributes:\n", - " Conventions: CF-1.7
array(['2017-01-01T00:00:00.000000000', '2017-01-01T01:00:00.000000000',\n", - " '2017-01-01T02:00:00.000000000', ..., '2021-12-31T16:00:00.000000000',\n", - " '2021-12-31T17:00:00.000000000', '2021-12-31T18:00:00.000000000'],\n", - " dtype='datetime64[ns]')
array([0., 1.])
array(['NO2-UM', 'NO2-ZAL Prat'], dtype=object)
array([[63.38, 46.01],\n", - " [64.64, 49.08],\n", + " Conventions: CF-1.7
array(['2017-01-01T00:00:00.000000000', '2017-01-01T01:00:00.000000000',\n", + " '2017-01-01T02:00:00.000000000', ..., '2021-12-31T15:00:00.000000000',\n", + " '2021-12-31T16:00:00.000000000', '2021-12-31T17:00:00.000000000'],\n", + " dtype='datetime64[ns]')
array([0., 1.])
array(['NO2-UM', 'NO2-ZAL Prat'], dtype=object)
array([[64.64, 49.08],\n", " [68.16, 53. ],\n", + " [68.29, 46.75],\n", " ...,\n", " [29.1 , 25.79],\n", " [ 9.24, 29.82],\n", - " [12.76, 28.66]])
array([41.373777, 41.317277])
array([2.184514, 2.134501])
array([41.373777, 41.317277])
array([2.184514, 2.134501])
43819 rows × 6 columns
\n", + "43818 rows × 6 columns
\n", "" ], "text/plain": [ " SO2-Darsena SO2-UM NO2-UM NO2-ZAL Prat month year\n", - "2017-01-01 00:00:00 0.18 NaN 63.38 46.01 1 2017\n", - "2017-01-01 01:00:00 0.20 NaN 64.64 49.08 1 2017\n", - "2017-01-01 02:00:00 0.18 NaN 68.16 53.00 1 2017\n", - "2017-01-01 03:00:00 0.04 NaN 68.29 46.75 1 2017\n", - "2017-01-01 04:00:00 0.00 NaN 62.91 39.65 1 2017\n", + "2017-01-01 00:00:00 0.20 NaN 64.64 49.08 1 2017\n", + "2017-01-01 01:00:00 0.18 NaN 68.16 53.00 1 2017\n", + "2017-01-01 02:00:00 0.04 NaN 68.29 46.75 1 2017\n", + "2017-01-01 03:00:00 0.00 NaN 62.91 39.65 1 2017\n", + "2017-01-01 04:00:00 0.00 NaN 47.13 28.86 1 2017\n", "... ... ... ... ... ... ...\n", - "2021-12-31 14:00:00 0.77 NaN 35.38 24.96 12 2021\n", - "2021-12-31 15:00:00 0.37 NaN 35.80 24.16 12 2021\n", - "2021-12-31 16:00:00 0.23 NaN 29.10 25.79 12 2021\n", - "2021-12-31 17:00:00 0.14 NaN 9.24 29.82 12 2021\n", - "2021-12-31 18:00:00 0.20 NaN 12.76 28.66 12 2021\n", + "2021-12-31 13:00:00 0.77 NaN 35.38 24.96 12 2021\n", + "2021-12-31 14:00:00 0.37 NaN 35.80 24.16 12 2021\n", + "2021-12-31 15:00:00 0.23 NaN 29.10 25.79 12 2021\n", + "2021-12-31 16:00:00 0.14 NaN 9.24 29.82 12 2021\n", + "2021-12-31 17:00:00 0.20 NaN 12.76 28.66 12 2021\n", "\n", - "[43819 rows x 6 columns]" + "[43818 rows x 6 columns]" ] }, "execution_count": 14, @@ -1506,20 +1506,20 @@ "Data variables:\n", " station_name (station) |S75 b'NO2-UM' b'NO2-ZAL Prat'\n", " altitude (station) float64 nan nan\n", - " sconcno2 (time, station) float64 6.88 23.03 8.77 ... 26.42 0.07 27.11\n", + " sconcno2 (time, station) float64 8.77 23.49 1.76 ... 27.11 13.87 30.07\n", " lat (station) float64 41.37 41.32\n", " lon (station) float64 2.185 2.135\n", "Attributes:\n", - " Conventions: CF-1.7array(['2017-05-01T00:00:00.000000000', '2017-05-01T01:00:00.000000000',\n", + " Conventions: CF-1.7
array(['2017-05-01T00:00:00.000000000', '2017-05-01T01:00:00.000000000',\n", " '2017-05-01T02:00:00.000000000', ..., '2017-05-31T21:00:00.000000000',\n", " '2017-05-31T22:00:00.000000000', '2017-05-31T23:00:00.000000000'],\n", - " dtype='datetime64[ns]')
array([0., 1.])
array([b'NO2-UM', b'NO2-ZAL Prat'], dtype='|S75')
array([nan, nan])
array([[ 6.88, 23.03],\n", - " [ 8.77, 23.49],\n", + " dtype='datetime64[ns]')
array([0., 1.])
array([b'NO2-UM', b'NO2-ZAL Prat'], dtype='|S75')
array([nan, nan])
array([[ 8.77, 23.49],\n", " [ 1.76, 21.73],\n", + " [ 0.57, 11.44],\n", " ...,\n", - " [21.96, 41.2 ],\n", " [ 0. , 26.42],\n", - " [ 0.07, 27.11]])
array([41.373777, 41.317277])
array([2.184514, 2.134501])
array([41.373777, 41.317277])
array([2.184514, 2.134501])
<xarray.Dataset>\n", - "Dimensions: (time: 744, station: 2)\n", + "Dimensions: (time: 738, station: 2)\n", "Coordinates:\n", - " * time (time) datetime64[ns] 2019-12-01 ... 2019-12-31T23:00:00\n", + " * time (time) datetime64[ns] 2021-12-01 ... 2021-12-31T17:00:00\n", " * station (station) float64 0.0 1.0\n", "Data variables:\n", " station_name (station) |S75 b'NO2-UM' b'NO2-ZAL Prat'\n", " altitude (station) float64 nan nan\n", - " sconcno2 (time, station) float64 65.22 80.1 63.97 ... 20.7 nan 16.95\n", + " sconcno2 (time, station) float64 13.33 nan 13.42 ... 29.82 12.76 28.66\n", " lat (station) float64 41.37 41.32\n", " lon (station) float64 2.185 2.135\n", "Attributes:\n", - " Conventions: CF-1.7
array(['2019-12-01T00:00:00.000000000', '2019-12-01T01:00:00.000000000',\n", - " '2019-12-01T02:00:00.000000000', ..., '2019-12-31T21:00:00.000000000',\n", - " '2019-12-31T22:00:00.000000000', '2019-12-31T23:00:00.000000000'],\n", - " dtype='datetime64[ns]')
array([0., 1.])
array([b'NO2-UM', b'NO2-ZAL Prat'], dtype='|S75')
array([nan, nan])
array([[65.22, 80.1 ],\n", - " [63.97, 81.52],\n", - " [57.43, 65.34],\n", + " Conventions: CF-1.7
array(['2021-12-01T00:00:00.000000000', '2021-12-01T01:00:00.000000000',\n", + " '2021-12-01T02:00:00.000000000', ..., '2021-12-31T15:00:00.000000000',\n", + " '2021-12-31T16:00:00.000000000', '2021-12-31T17:00:00.000000000'],\n", + " dtype='datetime64[ns]')
array([0., 1.])
array([b'NO2-UM', b'NO2-ZAL Prat'], dtype='|S75')
array([nan, nan])
array([[13.33, nan],\n", + " [13.42, nan],\n", + " [ 0.74, nan],\n", " ...,\n", - " [ nan, 25.97],\n", - " [ nan, 20.7 ],\n", - " [ nan, 16.95]])
array([41.373777, 41.317277])
array([2.184514, 2.134501])
array([41.373777, 41.317277])
array([2.184514, 2.134501])